AliPhysics  ef3b16e (ef3b16e)

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
m02LowCutMinSignal region, minimum shower shape value
m02LowCutMaxSignal region, maximum shower shape value
m02HighCutMinBkg region, minimum shower shape value
m02HighCutMaxBkg 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
m02LowCutMinSignal region, minimum shower shape value
m02LowCutMaxSignal region, maximum shower shape value
m02HighCutMinBkg region, minimum shower shape value
m02HighCutMaxBkg 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().

TString opt = "100MeV_New"

Definition at line 41 of file MakePurityCalculationAndComparisons.C.

Referenced by AddChi2JJDecLowHighGraphs(), AddTaskCaloFilter(), AliEmcalPhysicsSelectionTask::AliEmcalPhysicsSelectionTask(), EGCfg::Ampt(), AnalyseTuple(), anaTree(), 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().