From 81f2660fa11e1ab5eaf6b6c5819c656ecc6f3775 Mon Sep 17 00:00:00 2001 From: mcosenti Date: Wed, 19 Dec 2012 10:22:35 +0000 Subject: [PATCH] including histos to tag the isolation distribution of true gammas --- PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx | 14 ++++++++++++++ PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h | 2 ++ 2 files changed, 16 insertions(+) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx index 72af76b3941..7610eb3f502 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx @@ -78,6 +78,8 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton() : fMcPtInConeSBG(0), fMcPtInConeBGnoUE(0), fMcPtInConeSBGnoUE(0), + fAllIsoEtMcGamma(0), + fAllIsoNoUeEtMcGamma(0), fMCDirPhotonPtEtaPhiNoClus(0), fHnOutput(0) { @@ -130,6 +132,8 @@ AliAnalysisTaskEMCALIsoPhoton::AliAnalysisTaskEMCALIsoPhoton(const char *name) : fMcPtInConeSBG(0), fMcPtInConeBGnoUE(0), fMcPtInConeSBGnoUE(0), + fAllIsoEtMcGamma(0), + fAllIsoNoUeEtMcGamma(0), fMCDirPhotonPtEtaPhiNoClus(0), fHnOutput(0) { @@ -210,6 +214,12 @@ void AliAnalysisTaskEMCALIsoPhoton::UserCreateOutputObjects() fMcPtInConeSBGnoUE = new TH2F("hMcPtInConeSBGnoUE","#sum_{in-cone}p_{T}^{mc-primaries} vs. ISO^{TRK+EMC} (SBG range);ISO^{TRK+EMC} (GeV);#sum_{in-cone}p_{T}^{mc-primaries}",600,-10,50,1000,0,100); fOutputList->Add(fMcPtInConeSBGnoUE); + fAllIsoEtMcGamma = new TH2F("hAllIsoEtMcGammaE","ISO^{TRK+EMC} vs. E_{T}^{clus} for clusters comming from a MC prompt #gamma; E_{T}^{clus} (GeV);ISO^{TRK+EMC} (GeV);",1000,0,100,600,-10,50); + fOutputList->Add(fAllIsoEtMcGamma); + + fAllIsoNoUeEtMcGamma = new TH2F("hAllIsoNoUeEtMcGammaE","ISO^{TRK+EMC}_{noue} vs. E_{T}^{clus} for clusters comming from a MC prompt #gamma; E_{T}^{clus} (GeV);ISO^{TRK+EMC}_{noue} (GeV);",1000,0,100,600,-10,50); + fOutputList->Add(fAllIsoNoUeEtMcGamma); + fMCDirPhotonPtEtaPhiNoClus = new TH3F("hMCDirPhotonPhiEtaNoClus","p_{T}, #eta and #phi of prompt photons with no reco clusters;p_{T};#eta;#phi",1000,0,100,154,-0.77,0.77,130,1.38,3.20); fOutputList->Add(fMCDirPhotonPtEtaPhiNoClus); @@ -416,6 +426,10 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists() if(c->GetM02()>0.1 && c->GetM02()<0.3){ fMcPtInConeSBG->Fill(alliso-Et-allisoue, mcptsum); fMcPtInConeSBGnoUE->Fill(alliso-Et, mcptsum); + if(fMcIdFamily.Contains((Form("%d",c->GetLabel())))){ + fAllIsoEtMcGamma->Fill(Et, alliso-Et/*cecore*/-allisoue); + fAllIsoNoUeEtMcGamma->Fill(Et, alliso-Et/*cecore*/); + } } const Int_t ndims = fNDimensions; Double_t outputValues[ndims]; diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h index a357ca29268..cbaaeb8b0ca 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h @@ -102,6 +102,8 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE { TH2F *fMcPtInConeSBG; //!sum of mc-pt of "primary" particles inside de cone, as a function of NET-ISO in SBG range TH2F *fMcPtInConeBGnoUE; //!sum of mc-pt of "primary" particles inside de cone, as a function of NET-ISO in BG template no UE sub TH2F *fMcPtInConeSBGnoUE; //!sum of mc-pt of "primary" particles inside de cone, as a function of NET-ISO in SBG range no UE sub + TH2F *fAllIsoEtMcGamma; //!all iso distribution vs. Et clus for clusters comming from a MC prompt photon + TH2F *fAllIsoNoUeEtMcGamma; //!all iso distribution (without UE subtraction) vs. Et clus for clusters comming from a MC prompt photon TH3F *fMCDirPhotonPtEtaPhiNoClus; //!pt x eta x phi for prompt photons that didn't produce clusters THnSparse *fHnOutput; //!Output matrix with 7 dimensions -- 2.43.0