MakePurityCalculationAndComparisons.C File Reference

Calculate isolated photon purity. More...

#include <TFile.h>
#include <TDirectoryFile.h>
#include <TList.h>
#include <TString.h>
#include <TROOT.h>
#include <TStyle.h>
#include <TH1D.h>
#include <TH2F.h>
#include <TCanvas.h>
#include <TPad.h>
#include <TLegend.h>
#include <TObject.h>
#include <TAxis.h>
#include <TGaxis.h>
#include <TLine.h>
#include <TF1.h>
#include <TMath.h>
#include <TGraphErrors.h>
#include "PlotUtils.C"

Include dependency graph for MakePurityCalculationAndComparisons.C:

Go to the source code of this file.

Functions
void	MakeDoubleRatios (TString titleProd, TString gjProd, Int_t nCases, TString *mcCasesPath, TString *mcCases, TString *mcLeg, Bool_t addGJ=1, Float_t gjFactor=1, Bool_t usePi0=0, Int_t rebin=0, Float_t m02LowCutMin=0.10, Float_t m02LowCutMax=0.27, Float_t m02HigCutMin=0.4, Float_t m02HigCutMax=3.00, Int_t debug=kFALSE)
void	MakePurity (TString dataName, TString simuName, Int_t nCases, TString *mcCasesPath, TString *mcCases, TString *mcLeg, TString gj="", Bool_t usePi0=kTRUE, Float_t m02LowCutMin=0.10, Float_t m02LowCutMax=0.27, Float_t m02HigCutMin=0.40, Float_t m02HigCutMax=3.00, Bool_t debug=kFALSE)
void	MakePurityCalculationAndComparisons (Bool_t doDoRat=kTRUE, Bool_t doPurity=kTRUE, Bool_t debug=kFALSE)

Variables
Int_t	color [] = {1,4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2,4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2}
Int_t	lineStyle [] = {1,1,1,1,1,1, 1, 1, 1, 1, 1,2,2,2,2,2, 2, 2, 2, 2, 2,2,2,2,2,2,2,2,2,2,2,2}
Int_t	marker [] = {24,21,20,25,24,24,24,24,24,24,24}
TString	format = ".eps"
TString	opt = "100MeV_New"
TString	dataFileDir = "data3"

Detailed Description

Calculate isolated photon purity.

Example macro to calculate isolated photon purity for different data/MC productions and different shower shape selection regions with the ABCD method.

Author

Gustavo Conesa Balbastre Gusta.nosp@m.vo.C.nosp@m.onesa.nosp@m..Bal.nosp@m.bastr.nosp@m.e@ce.nosp@m.rn.ch, (LPSC-CNRS)

Definition in file MakePurityCalculationAndComparisons.C.

Function Documentation

void MakeDoubleRatios	(	TString	titleProd,
		TString	gjProd,
		Int_t	nCases,
		TString *	mcCasesPath,
		TString *	mcCases,
		TString *	mcLeg,
		Bool_t	addGJ = 1,
		Float_t	gjFactor = 1,
		Bool_t	usePi0 = 0,
		Int_t	rebin = 0,
		Float_t	m02LowCutMin = 0.10,
		Float_t	m02LowCutMax = 0.27,
		Float_t	m02HigCutMin = 0.4,
		Float_t	m02HigCutMax = 3.00,
		Int_t	debug = kFALSE
	)

Get the histograms defining the different purity regions ABDC, and do the double ratios needed for purity calculation in MakePurity() where data and MC double ratios are combined. This method is used the same for that and MC.

Parameters

titleProd	: Name of analized production
gjProd	: Name of gamma-jet analized production
nCases	: in case of MC, the different cases studied, in case of data it is 1
mcCasesPath	: path to MC/Data case input file
mcCases	: simplified name of MC/Data sample
mcLef	: name of MC/Data sample for legends in plots
addGJ	: add the GJ sample (P3) or not (P2)
gjFactor	: in case of addition, scale the GJ sample by this factor
usePi0	: Use as bkg region the pi0 shower shape band
rebin	: if 0 use the energy binning defined globally, if not, use simple binning number
m02LowCutMin	Signal region, minimum shower shape value
m02LowCutMax	Signal region, maximum shower shape value
m02HighCutMin	Bkg region, minimum shower shape value
m02HighCutMax	Bkg region, maximum shower shape value
debug	: Bool to activate printfs.

ABCD ///

ABCD Ratio of the different MC productions///

ABCD Double Ratio ///

ABCD Double Ratio Cases ///

Definition at line 67 of file MakePurityCalculationAndComparisons.C.

Referenced by MakePurityCalculationAndComparisons().

void MakePurity	(	TString	dataName,
		TString	simuName,
		Int_t	nCases,
		TString *	mcCasesPath,
		TString *	mcCases,
		TString *	mcLeg,
		TString	gj = "",
		Bool_t	usePi0 = kTRUE,
		Float_t	m02LowCutMin = 0.10,
		Float_t	m02LowCutMax = 0.27,
		Float_t	m02HigCutMin = 0.40,
		Float_t	m02HigCutMax = 3.00,
		Bool_t	debug = kFALSE
	)

Get the histograms defining the different purity regions ABDC, and do the double ratios needed for purity calculation in MakePurity() where data and MC double ratios are combined. This method is used the same for that and MC.

Parameters

dataName	: Name analyzed data sample
simuName	: Name analyzed MC sample
nCases	: in case of MC, the different cases studied, in case of data it is 1
mcCasesPath	: path to MC/Data case input file
mcCases	: simplified name of MC/Data sample
mcLef	: name of MC/Data sample for legends in plots
gj	: string to select the addition and scaling of the GJ sample in MC - "","_GJx0.50","_GJx1.00","_GJx2.00"
usePi0	: Use as bkg region the pi0 shower shape band
m02LowCutMin	Signal region, minimum shower shape value
m02LowCutMax	Signal region, maximum shower shape value
m02HighCutMin	Bkg region, minimum shower shape value
m02HighCutMax	Bkg region, maximum shower shape value
debug	: Bool to activate printfs.

Purity ///

Purity Ratio Cases ///

Definition at line 821 of file MakePurityCalculationAndComparisons.C.

Referenced by MakeDoubleRatios(), and MakePurityCalculationAndComparisons().

void MakePurityCalculationAndComparisons	(	Bool_t	doDoRat = kTRUE,
		Bool_t	doPurity = kTRUE,
		Bool_t	debug = kFALSE
	)

Main steering method. Execute here the different selection options

Parameters

doDoRat	: Execute method MakeDoubleRatios()
doPurity	: Execute method MakePurity()
debug	: Bool to activate printfs.

Definition at line 1055 of file MakePurityCalculationAndComparisons.C.

Referenced by MakePurity().

Variable Documentation

Int_t color[] = {1,4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2,4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2}

Definition at line 35 of file MakePurityCalculationAndComparisons.C.

Referenced by MakeDoubleRatios(), and MakePurity().

TString dataFileDir = "data3"

Definition at line 42 of file MakePurityCalculationAndComparisons.C.

TString format = ".eps"

Definition at line 38 of file MakePurityCalculationAndComparisons.C.

Referenced by MakeDoubleRatios(), and MakePurity().

Int_t lineStyle[] = {1,1,1,1,1,1, 1, 1, 1, 1, 1,2,2,2,2,2, 2, 2, 2, 2, 2,2,2,2,2,2,2,2,2,2,2,2}

Definition at line 36 of file MakePurityCalculationAndComparisons.C.

Referenced by SetHistogramRangesAndTitles().

Int_t marker[] = {24,21,20,25,24,24,24,24,24,24,24}

Definition at line 37 of file MakePurityCalculationAndComparisons.C.

Referenced by Drawer::GetDataOther(), Drawer::GetStack(), MakeDoubleRatios(), MakePurity(), AliForwarddNdetaTask::CentralityBin::operator=(), AliMCTruthdNdetaTask::CentralityBin::operator=(), AliBasedNdetaTask::PreEvent(), QAPlotter::Ring::SetAtt(), SetAttr(), SetAttributes(), and SetHistogramRangesAndTitles().

TString opt = "100MeV_New"

Definition at line 41 of file MakePurityCalculationAndComparisons.C.

Referenced by AddChi2JJDecLowHighGraphs(), AddTaskCaloFilter(), AliEmcalPhysicsSelectionTask::AliEmcalPhysicsSelectionTask(), EGCfg::Ampt(), AnalyseTuple(), anaTree(), anaTreeV3(), AliRDHFCutsDstoKKpi::ApplyCutOnResonances(), AutoZoom(), BuildScript(), AliMTRChEffAnalysis::CompareMergedEfficiencies(), CompareTH3DataAndMCProd(), ConfigurePi0EbEAnalysis(), ConfigurePi0Flow(), EMCalTriggerPtAnalysis::AliEMCalTriggerAnaTriggerDecision::ConfigureTriggerDecision(), AliJetContainer::ConnectClusterContainer(), GraphSysErr::CopyKeys(), LInfo::CreateHistograms(), TrainSetup::CreatePhysicsSelection(), THistManager::CreateTProfile(), AliPWGHistoTools::DoIntegral(), AliFMDCorrELossFit::ELossFit::Draw(), AliTrackletPtPidStrWeights::Draw(), Ring::Draw(), AliFMDCorrELossFit::Draw(), GraphSysErr::Draw(), AliTrackletDeltaWeights::Draw(), SummaryMCTrackDrawer::DrawBetaGammadEdx(), DrawCorrCentralSecMap2(), SummarydNdetaDrawer::DrawResTitle(), GraphSysErr::Export(), FastLoad(), AliAnaCaloTrackCorrMaker::GetCaloUtils(), AliClusterContainer::GetClusPtCut(), AliFiducialCut::GetDCALFiducialCut(), AliAODCentralMult::GetHistogram(), AliMCParticleContainer::GetLeadingMCParticle(), AliTrackContainer::GetLeadingTrack(), AliMCAnalysisUtils::GetMCGeneratorString(), GraphSysErr::GetMinMax(), AliIntSpotEstimator::GetTrackMaxP(), AliAnaChargedParticles::GetTrackSector(), AliEmcalCorrectionTask::GetYAMLConfiguration(), OptionList::Help(), EGCfg::Hijing(), AliFMDMCHitHandler::Init(), AliMCAuxHandler::Init(), TDraw::IsGood(), ForwardOADBGUI::LoadCorrDraw(), LoadOne(), AliAnalysisTaskSEPicoV0Filter::LocalInit(), AliAnalysisTaskEmcalTmpSparseMaker::LocalInit(), AliAnalysisTaskEmcalJetV0Filter::LocalInit(), AliAnalysisTaskEmcalJetSparseMaker::LocalInit(), AliAnalysisTaskEmcalJetV0CF::LocalInit(), AliAnalysisTaskSEPicoV0MakerMC::LocalInit(), AliAnalysisTaskSEPicoV0Maker::LocalInit(), AliAnaScale::LocalInit(), AliFMDCorrAcceptance::ls(), MakeDataMCComparisonPerSMClusterEbin(), GraphSysErr::MakeMulti(), MakeNDataNMCComparisonPerSMClusterEbin(), AliForwardUtil::Histos::operator=(), GraphSysErr::Holder::operator=(), Railway::OutputLocation(), OptionList::Parse(), AliBasedNdetaTask::ParseNormalizationScheme(), PlotGraphs(), TupleSelector::PlotOne(), AliFMDEncodedEdx::Spec::Print(), AliCorrectionManagerBase::Print(), AliFMDDensityCalculator::Print(), Fast1DCentEstimator::Print(), AliBasedNdetaTask::Print(), AliCorrectionManagerBase::Correction::Print(), GraphSysErr::Print(), FastNchCentEstimator::Print(), V0CentEstimator::Print(), AliFMDCorrELossFit::Print(), ZNCentEstimator::Print(), AliBasedNdetaTask::CentralityBin::Print(), GraphSysErr::Holder::Print(), GraphSysErr::HolderP2P::Print(), GraphSysErr::HolderCommon::Print(), ProcessCandidates(), ProcessFast(), TupleSelector::Proof(), FastSim::ProofRun(), QuickDraw(), THistManager::ReleaseOwner(), FastAnalysis::Run(), FastSim::Run(), EPosSim::Run(), RunEPos(), RunFast(), SaveCanvas(), GraphSysErr::SavePrimitive(), AliAnalysisTaskDmesonMCPerform::SetAODMismatchProtection(), AliAnalysisTaskTrackingSysPropagation::SetAODMismatchProtection(), AliAnalysisTaskCheckEvSel::SetAODMismatchProtection(), AliAnalysisTaskSESignificance::SetAODMismatchProtection(), AliAnalysisTaskSEDStarSpectra::SetAODMismatchProtection(), AliAnalysisTaskSEDs::SetAODMismatchProtection(), AliAnalysisTaskSEHFQA::SetAODMismatchProtection(), AliAnalysisTaskSEDplus::SetAODMismatchProtection(), AliAnalysisTaskSEDstoK0sK::SetAODMismatchProtection(), AliAnalysisTaskSED0MassNonPromptFraction::SetAODMismatchProtection(), AliAnalysisTaskSED0Mass::SetAODMismatchProtection(), AliAnalysisTaskSECharmFraction::SetAODMismatchProtection(), AliAnalysisTaskSELc2V0bachelor::SetAODMismatchProtection(), AliAnalysisTaskEmcalEmbeddingHelper::SetAutoConfigureIdentifier(), AliCaloTrackReader::SetCentralityOpt(), GraphSysErr::SetCommonSumOption(), PWGJE::EMCALJetTasks::AliAnalysisTaskEmcalJetHPerformance::SetConfigurationPath(), AliCopyHeaderTask::SetCopyOptions(), AliAnalysisTaskCheckEvSel::SetCutOnzVertexSPD(), AliRDHFCuts::SetCutOnzVertexSPD(), GraphSysErr::SetDataOption(), AliAnalysisTaskSEDs::SetDoCutVarHistos(), GraphSysErr::Holder::SetDOption(), AliHFMultiTrials::SetDrawIndividualFits(), AliHFInvMassMultiTrialFit::SetDrawIndividualFits(), AliAnalysisTaskDmesonMCPerform::SetFillExtraHistos(), AliCFVertexingHF3Prong::SetGeneratedDsOption(), AliAnalysisTaskCombinHF::SetGoUpToQuark(), AliHFVnVsMassFitter::SetInitialFrac2Gaus(), AliHFVnVsMassFitter::SetInitialGaussianMean(), AliHFVnVsMassFitter::SetInitialGaussianSigma(), AliHFVnVsMassFitter::SetInitialGaussianSigma2Gaus(), AliFMDMCTrackELoss::SetMaxConsequtiveStrips(), AliAnalysisTaskSEImproveITS::SetMimicData(), AliAnalysisTaskLocalRho::SetModulationFitOptions(), AliAnalysisTaskJetV2::SetModulationFitOptions(), AliAnalysisTaskJetV3::SetModulationFitOptions(), AliHFInvMassFitter::SetNSigma4SideBands(), AliRDHFCuts::SetOptPileup(), AliParticleContainer::SetParticlePhiLimits(), AliRDHFCutsDstoKKpi::SetPidOption(), AliRDHFCutsDstoK0sK::SetPidOption(), AliAODPidHF::SetppLowEn2011(), AliJetRandomizerTask::SetRandomizeEta(), AliAnalysisTaskCheckHFMCProd::SetReadMC(), AliHFVnVsMassFitter::SetReflVnOption(), AliCFVertexingHF::SetRejectCandidateIfNotFromQuark(), AliCFTaskVertexingHFCutVarFDSub::SetRejectCandidateIfNotFromQuark(), AliCFTaskVertexingHF::SetRejectCandidateIfNotFromQuark(), AliRDHFCuts::SetRemoveTrackletOutliers(), AliHFMultiTrials::SetSaveBkgValue(), AliHFInvMassMultiTrialFit::SetSaveBkgValue(), AliRDHFCutsDplustoKpipi::SetScaleNormDLxyBypOverPt(), AliAnalysisTaskCheckHFMCProd::SetSearchUpToQuark(), AliHFInvMassFitter::SetSmoothReflectionTemplate(), AliAnalysisTaskSEDvsEventShapes::SetSubtractTrackletsFromDaughters(), AliAnalysisTaskSEDvsMultiplicity::SetSubtractTrackletsFromDaughters(), GraphSysErr::SetSumOption(), AliHFVnVsMassFitter::SetTemplateReflections(), AliAnalysisTaskSEImproveITS::SetUpdatePulls(), AliAnalysisTaskSEImproveITS::SetUpdateSTCovMatrix(), EPosSim::SetupProof(), TrackletdNdetaTrain::SetupReweighting(), MakeTrackletTrain::SetupReweighting(), AliHFInvMassMultiTrialFit::SetUse2GausSigmaRatioSignal(), AliHFInvMassMultiTrialFit::SetUse2GausSignal(), AliRDHFCuts::SetUseCentralityFlatteningInMC(), AliRDHFCuts::SetUseCutGeoNcrNcl(), AliHFMultiTrials::SetUseExpoBackground(), AliHFInvMassMultiTrialFit::SetUseExpoBackground(), AliHFMultiTrials::SetUseFixedMeanFreeS(), AliHFInvMassMultiTrialFit::SetUseFixedMeanFreeS(), AliHFMultiTrials::SetUseFixSigDownFreeMean(), AliHFInvMassMultiTrialFit::SetUseFixSigDownFreeMean(), AliHFMultiTrials::SetUseFixSigFixMean(), AliHFInvMassMultiTrialFit::SetUseFixSigFixMean(), AliHFMultiTrials::SetUseFixSigFreeMean(), AliHFInvMassMultiTrialFit::SetUseFixSigFreeMean(), AliHFMultiTrials::SetUseFixSigUpFreeMean(), AliHFInvMassMultiTrialFit::SetUseFixSigUpFreeMean(), AliHFMultiTrials::SetUseFreeS(), AliHFInvMassMultiTrialFit::SetUseFreeS(), AliAnalysisVertexingHF::SetUseKaonPIDfor3Prong(), AliAnalysisVertexingHF::SetUseKaonPIDforDs(), AliHFMultiTrials::SetUseLinBackground(), AliHFInvMassMultiTrialFit::SetUseLinBackground(), AliCFTaskVertexingHFCutVarFDSub::SetUseMCVertex(), AliCFTaskVertexingHF::SetUseMCVertex(), AliRDHFCuts::SetUseMultDepPileupCut(), AliAnalysisTaskSELambdacTMVA::SetUseNchWeight(), AliAnalysisVertexingHF::SetUsePidTag(), AliHFMultiTrials::SetUsePol2Background(), AliHFInvMassMultiTrialFit::SetUsePol2Background(), AliHFMultiTrials::SetUsePol3Background(), AliHFInvMassMultiTrialFit::SetUsePol3Background(), AliHFMultiTrials::SetUsePol4Background(), AliHFInvMassMultiTrialFit::SetUsePol4Background(), AliHFMultiTrials::SetUsePol5Background(), AliHFInvMassMultiTrialFit::SetUsePol5Background(), AliHFMultiTrials::SetUsePowerLawBackground(), AliHFInvMassMultiTrialFit::SetUsePowerLawBackground(), AliHFMultiTrials::SetUsePowerLawTimesExpoBackground(), AliHFInvMassMultiTrialFit::SetUsePowerLawTimesExpoBackground(), AliAnalysisTaskSEDvsEventShapes::SetUseQuarkLevel(), AliAnalysisTaskSEB0toDminuspi::SetUseQuarkLevelTag(), AliAnalysisTaskSEDplus::SetUseQuarkLevelTag(), AliAnalysisTaskSED0MassNonPromptFraction::SetUseQuarkLevelTag(), AliAnalysisTaskSED0Mass::SetUseQuarkLevelTag(), AliAnalysisTaskSEDs::SetUseSelectionBit(), AliAnalysisVertexingHF::SetUseTOFPID(), AliAnalysisVertexingHF::SetUseTPCPID(), AliAnalysisVertexingHF::SetUseTPCPIDOnlyIfNoTOF(), AliAnalysisTaskSEDs::SetWriteOnlySignalInNtuple(), OptionList::Show(), AliHFVnVsMassFitter::SimultaneusFit(), OptionList::Store(), AliCaloTrackReader::SwitchOnLEDEventsRemoval(), AliAnalysisTaskEMCALClusterize::SwitchOnUseClusterMCLabelForCell(), OptionList::Widest(), WithSysError(), AliAnaCalorimeterQA::~AliAnaCalorimeterQA(), AliAnaClusterPileUp::~AliAnaClusterPileUp(), AliAnaClusterShapeCorrelStudies::~AliAnaClusterShapeCorrelStudies(), AliAnaElectron::~AliAnaElectron(), AliAnaEMCALTriggerClusters::~AliAnaEMCALTriggerClusters(), AliAnaInsideClusterInvariantMass::~AliAnaInsideClusterInvariantMass(), AliAnaParticleIsolation::~AliAnaParticleIsolation(), AliAnaParticleJetLeadingConeCorrelation::~AliAnaParticleJetLeadingConeCorrelation(), AliAnaParticlePartonCorrelation::~AliAnaParticlePartonCorrelation(), AliAnaPhoton::~AliAnaPhoton(), AliAnaPhotonConvInCalo::~AliAnaPhotonConvInCalo(), AliAnaPi0EbE::~AliAnaPi0EbE(), AliAnaRandomTrigger::~AliAnaRandomTrigger(), AliNeutralMesonSelection::~AliNeutralMesonSelection(), and LInfo::~LInfo().