17 #include <TClonesArray.h>
19 #include <TObjString.h>
24 #include "TParticle.h"
25 #include "TDatabasePDG.h"
35 #include "AliAODMCParticle.h"
36 #include "AliAODPWG4ParticleCorrelation.h"
38 #include "AliVTrack.h"
39 #include "AliVCluster.h"
40 #include "AliESDEvent.h"
41 #include "AliAODEvent.h"
44 #include "AliEMCALGeometry.h"
45 #include "AliPHOSGeoUtils.h"
56 fIsoDetector(-1), fIsoDetectorString(
""),
57 fReMakeIC(0), fMakeSeveralIC(0),
58 fFillTMHisto(0), fFillSSHisto(1),
59 fFillEMCALRegionHistograms(0), fFillUEBandSubtractHistograms(1),
60 fFillCellHistograms(0),
61 fFillOverlapHistograms(0),
62 fStudyTracksInCone(0), fStudyMCConversionRadius(0),
63 fStudyExoticTrigger(0),
64 fFillTaggedDecayHistograms(0), fNDecayBits(0),
65 fDecayBits(), fDecayTagsM02Cut(0),
66 fFillNLMHistograms(0),
67 fLeadingOnly(0), fCheckLeadingWithNeutralClusters(0),
68 fSelectPrimariesInCone(0), fMakePrimaryPi0DecayStudy(0),
69 fFillBackgroundBinHistograms(0), fNBkgBin(0),
70 fFillPtTrigBinHistograms(0), fNPtTrigBin(0),
71 fMinCellsAngleOverlap(0),
73 fNCones(0), fNPtThresFrac(0),
74 fConeSizes(), fPtThresholds(),
75 fPtFractions(), fSumPtThresholds(),
76 fStudyPtCutInCone(0), fNPtCutsInCone(0),
77 fMinPtCutInCone(), fMaxPtCutInCone(),
78 fStudyEtaCutInCone(0), fNEtaCutsInCone(0), fEtaCutInCone(),
79 fStudyRCutInCone(0), fNRCutsInCone(0), fRCutInCone(),
80 fNNCellsInCandidate(0), fNCellsInCandidate(),
81 fNExoCutInCandidate(0), fExoCutInCandidate(),
82 fMomentum(), fMomIso(),
83 fMomDaugh1(), fMomDaugh2(),
84 fTrackVector(), fProdVertex(),
85 fCluster(0), fClustersArr(0),
86 fIsExoticTrigger(0), fClusterExoticity(1),
88 fhEIso(0), fhPtIso(0),
89 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
91 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
93 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
94 fhEIsoExoTrigger(0), fhENoIsoExoTrigger(0),
95 fhPtIsoExoTrigger(0), fhPtNoIsoExoTrigger(0),
98 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
99 fhPtInConeExoTrigger(0), fhPtClusterInConeExoTrigger(0), fhPtTrackInConeExoTrigger(0),
100 fhPtTrackInConeOtherBCPileUpSPD(0), fhPtTrackInConeVtxBC0(0),
101 fhPtTrackInConeBC0PileUpSPD(0),
102 fhPtInConePileUp(), fhPtInConeCent(0),
104 fhPerpConeSumPt(0), fhPerpConeSumPtTOFBC0(0),
105 fhPtInPerpCone(0), fhPtInPerpConeTOFBC0(0),
106 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
107 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
108 fhEtaPhiInPerpCone(0), fhEtaPhiInPerpConeTOFBC0(0),
109 fhEtaBandClusterEtaPhi(0), fhPhiBandClusterEtaPhi(0),
110 fhEtaBandTrackEtaPhi(0), fhPhiBandTrackEtaPhi(0),
111 fhEtaBandClusterPt(0), fhPhiBandClusterPt(0),
112 fhEtaBandTrackPt(0), fhPhiBandTrackPt(0),
113 fhEtaBandCell(0), fhPhiBandCell(0),
114 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
115 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
117 fhConeSumPtCellTrack(0),
118 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
119 fhConeSumPtExoTrigger(0), fhConeSumPtClusterExoTrigger(0), fhConeSumPtTrackExoTrigger(0),
121 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
122 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
123 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
124 fhConeSumPtTrigEtaPhi(0),
125 fhConeSumPtCellTrackTrigEtaPhi(0),
126 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
127 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
128 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
129 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
130 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
131 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
132 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
133 fhConeSumPtEtaUENormCluster(0), fhConeSumPtPhiUENormCluster(0),
134 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
135 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
136 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
137 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
138 fhConeSumPtEtaUENormTrack(0), fhConeSumPtPhiUENormTrack(0),
139 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
140 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
141 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
142 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
143 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
144 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
145 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
146 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
147 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
148 fhConeSumPtCellvsTrack(0),
149 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
150 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
151 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
152 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
153 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
154 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
155 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
156 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
157 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
158 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
159 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
160 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
161 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
162 fhPtPrimMCPi0DecayPairOutOfCone(0),
163 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
164 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
165 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
166 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
167 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
168 fhPtPrimMCPi0DecayPairNoOverlap(0),
169 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
170 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
171 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
172 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
173 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
174 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
175 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
176 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
177 fhPtPrimMCEtaDecayPairOutOfCone(0),
178 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
179 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
180 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
181 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
182 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
183 fhPtPrimMCEtaDecayPairNoOverlap(0),
184 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
185 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
186 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
187 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
188 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
189 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
190 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
191 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
192 fhPtLeadConeBin(0), fhSumPtConeBin(0),
193 fhPtLeadConeBinMC(0), fhSumPtConeBinMC(0),
194 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
195 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
196 fhPtLeadConeBinLambda0MC(0), fhSumPtConeBinLambda0MC(0),
197 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
198 fhPtTrigBinSumPtTrackCone(0), fhPtTrigBinSumPtClusterCone(0),
199 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
200 fhPtTrigBinSumPtTrackConeMC(0), fhPtTrigBinSumPtClusterConeMC(0),
201 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
202 fhPtTrigBinSumPtTrackConeDecay(0), fhPtTrigBinSumPtClusterConeDecay(0),
203 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
204 fhPtTrigBinLambda0vsSumPtTrackCone(0), fhPtTrigBinLambda0vsSumPtClusterCone(0),
205 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
206 fhPtTrigBinLambda0vsSumPtTrackConeMC(0), fhPtTrigBinLambda0vsSumPtClusterConeMC(0),
207 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap(0),fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap(0),
208 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap(0), fhPtTrigBinLambda0vsSumPtConeMC1Overlap(0),
209 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap(0),fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap(0),
212 fhELambda0LocMax1(), fhELambda1LocMax1(),
213 fhELambda0LocMax2(), fhELambda1LocMax2(),
214 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
216 fhEIsoPileUp(), fhPtIsoPileUp(),
217 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
218 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
219 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
220 fhTimeNPileUpVertContributors(0),
221 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
223 fhPtClusterInConePerRCut(0), fhPtClusterInConePerRCutLargePtTrig(0),
224 fhPtTrackInConePerRCut(0), fhPtTrackInConePerRCutLargePtTrig(0),
225 fhConeSumPtClusterPerRCut(0), fhConeSumPtClusterPerRCutLargePtTrig(0),
226 fhConeSumPtTrackPerRCut(0), fhConeSumPtTrackPerRCutLargePtTrig(0),
227 fhConeNClusterPerMinPtCut(0), fhConeNClusterPerMinPtCutLargePtTrig(0),
228 fhConeNTrackPerMinPtCut(0), fhConeNTrackPerMinPtCutLargePtTrig(0),
229 fhPerpConeNTrackPerMinPtCut(0), fhPerpConeNTrackPerMinPtCutLargePtTrig(0),
230 fhConeSumPtClusterPerMinPtCut(0), fhConeSumPtClusterPerMinPtCutLargePtTrig(0),
231 fhConeSumPtTrackPerMinPtCut(0), fhConeSumPtTrackPerMinPtCutLargePtTrig(0),
232 fhPerpConeSumPtTrackPerMinPtCut(0), fhPerpConeSumPtTrackPerMinPtCutLargePtTrig(0),
233 fhConeSumPtClusterPerMaxPtCut(0), fhConeSumPtClusterPerMaxPtCutLargePtTrig(0),
234 fhConeSumPtTrackPerMaxPtCut(0), fhConeSumPtTrackPerMaxPtCutLargePtTrig(0),
235 fhConeSumPtTrackPerEtaCut(0), fhConeSumPtTrackPerEtaCutLargePtTrig(0),
237 fhPtClusterInConePerNCellCut(0), fhPtClusterInConePerNCellCutLargePtTrig(0),
238 fhPtTrackInConePerNCellCut(0), fhPtTrackInConePerNCellCutLargePtTrig(0),
239 fhConeSumPtClusterPerNCellCut(0), fhConeSumPtClusterPerNCellCutLargePtTrig(0),
240 fhConeSumPtTrackPerNCellCut(0), fhConeSumPtTrackPerNCellCutLargePtTrig(0),
241 fhPtClusterInConePerExoCut(0), fhPtClusterInConePerExoCutLargePtTrig(0),
242 fhPtTrackInConePerExoCut(0), fhPtTrackInConePerExoCutLargePtTrig(0),
243 fhConeSumPtClusterPerExoCut(0), fhConeSumPtClusterPerExoCutLargePtTrig(0),
244 fhConeSumPtTrackPerExoCut(0), fhConeSumPtTrackPerExoCutLargePtTrig(0),
246 fhConeSumPtTrackTOFBC0(0), fhConeSumPtTrackTOFBCN(0), fhConeSumPtTrackTOFNo(0),
247 fhPtTrackInConeTOFBC0 (0), fhPtTrackInConeTOFBCN (0), fhPtTrackInConeTOFNo (0),
248 fhPhiTrackInCone(0), fhEtaTrackInCone(0), fhEtaPhiTrackInCone(0),
249 fhPhiTrackInConeTOFBC0(0), fhPhiTrackInConeTOFBCN(0), fhPhiTrackInConeTOFNo(0),
250 fhEtaTrackInConeTOFBC0(0), fhEtaTrackInConeTOFBCN(0), fhEtaTrackInConeTOFNo(0),
251 fhEtaPhiTrackInConeTOFBC0(0), fhEtaPhiTrackInConeTOFBCN(0), fhEtaPhiTrackInConeTOFNo(0),
252 fhTrackTOFInCone(0), fhTrackTOFInConeBC0(0), fhTrackTOFInConeExoTrigger(0),
254 fhConeSumPtTrackITSRefitOnSPDOn(0), fhConeSumPtTrackITSRefitOnSPDOff(0), fhConeSumPtTrackITSRefitOffSPDOff(0),
255 fhPtTrackInConeITSRefitOnSPDOn(0), fhPtTrackInConeITSRefitOnSPDOff(0) , fhPtTrackInConeITSRefitOffSPDOff(0),
256 fhPhiTrackInConeITSRefitOnSPDOn(0), fhPhiTrackInConeITSRefitOnSPDOff(0), fhPhiTrackInConeITSRefitOffSPDOff(0),
257 fhEtaTrackInConeITSRefitOnSPDOn(0), fhEtaTrackInConeITSRefitOnSPDOff(0), fhEtaTrackInConeITSRefitOffSPDOff(0),
258 fhEtaPhiTrackInConeITSRefitOnSPDOn(0),fhEtaPhiTrackInConeITSRefitOnSPDOff(0), fhEtaPhiTrackInConeITSRefitOffSPDOff(0),
259 fhConeSumPtTrackTOFBC0ITSRefitOnSPDOn(0), fhPtTrackInConeTOFBC0ITSRefitOnSPDOn(0),
260 fhPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0), fhEtaTrackInConeTOFBC0ITSRefitOnSPDOn(0),fhEtaPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0),
261 fhPerpConeSumPtITSRefitOnSPDOn (0), fhPtInPerpConeITSRefitOnSPDOn(0), fhEtaPhiInPerpConeITSRefitOnSPDOn(0),
262 fhPerpConeSumPtTOFBC0ITSRefitOnSPDOn (0), fhPtInPerpConeTOFBC0ITSRefitOnSPDOn (0), fhEtaPhiInPerpConeTOFBC0ITSRefitOnSPDOn(0)
266 for(
Int_t i = 0; i < 5 ; i++)
273 for(
Int_t j = 0; j < 5 ; j++)
308 for(
Int_t iso =0; iso < 2; iso++)
318 for(
Int_t i = 0; i < 5 ; i++)
337 for(
Int_t i = 0; i < 2 ; i++)
353 for(
Int_t imc = 0; imc < 4; imc++)
361 for(
Int_t i = 0; i < 2 ; i++)
377 for(
Int_t ieta = 0; ieta < 4; ieta++)
379 for(
Int_t iphi = 0; iphi < 3; iphi++)
390 for(
Int_t ism =0; ism < 20; ism++)
413 for(
Int_t i = 0 ; i < 7 ; i++)
448 Float_t ptTrig = pCandidate->Pt() ;
449 Float_t phiTrig = pCandidate->Phi();
450 Float_t etaTrig = pCandidate->Eta();
452 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
454 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
458 AliWarning(
"Cluster not available?");
463 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
464 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
527 Float_t phiTrig = pCandidate->Phi();
528 if(phiTrig<0) phiTrig += TMath::TwoPi();
530 Float_t etaTrig = pCandidate->Eta();
532 if(pCandidate->GetDetectorTag()==
kEMCAL)
534 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
537 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
539 if(!eGeom->CheckAbsCellId(absId))
return ;
542 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
544 Int_t imEta=-1, imPhi=-1;
545 Int_t ieta =-1, iphi =-1;
547 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
549 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
551 Int_t colTrig = ieta;
552 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
554 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
556 Int_t sqrSize = int(conesize/0.0143);
560 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
561 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
565 Int_t irowmin = rowTrig-sqrSize;
566 if ( irowmin < 0 ) irowmin = 0 ;
568 Int_t irowmax = rowTrig+sqrSize;
569 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
571 for(
Int_t irow = irowmin; irow <irowmax; irow++)
573 for(
Int_t icol = 0; icol < nTotalCols; icol++)
575 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
576 if(inSector==5)
continue;
579 if(icol < AliEMCALGeoParams::fgkEMCALCols)
581 inSupMod = 2*inSector + 1;
584 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
586 inSupMod = 2*inSector;
587 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
590 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
593 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
597 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
598 if(!eGeom->CheckAbsCellId(iabsId))
continue;
600 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
608 Int_t icolmin = colTrig-sqrSize;
609 if ( icolmin < 0 ) icolmin = 0;
611 Int_t icolmax = colTrig+sqrSize;
612 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
615 for(
Int_t icol = icolmin; icol < icolmax; icol++)
617 for(
Int_t irow = 0; irow < nTotalRows; irow++)
619 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
620 if ( inSector == 5 ) continue ;
626 if(icol < AliEMCALGeoParams::fgkEMCALCols)
629 inSupMod = 2*inSector + 1;
632 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
635 inSupMod = 2*inSector;
636 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
639 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
642 if(TMath::Abs(icol-colTrig) < sqrSize) {
645 if(TMath::Abs(irow-rowTrig) < sqrSize) {
649 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
650 if(!eGeom->CheckAbsCellId(iabsId))
652 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
656 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
666 Float_t ptTrig = pCandidate->Pt();
689 Float_t coneptsumPerpTrackPerMinCut[20];
690 Float_t coneNPerpTrackPerMinCut [20];
696 coneptsumPerpTrackPerMinCut[icut] = 0;
697 coneNPerpTrackPerMinCut [icut] = 0;
701 Float_t ptTrig = pCandidate->Pt() ;
702 Float_t phiTrig = pCandidate->Phi();
703 Float_t etaTrig = pCandidate->Eta();
708 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
710 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
714 AliWarning(
"Track not available?");
721 if ( pCandidate->GetDetectorTag() ==
kCTS )
724 Bool_t contained = kFALSE;
726 for(
Int_t i = 0; i < 4; i++)
728 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
731 if ( contained ) continue ;
749 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
751 phiBandPtSum+=track->Pt();
758 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
760 etaBandPtSum+=track->Pt();
768 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
769 Double_t dEta = etaTrig - track->Eta();
770 Double_t arg = dPhi*dPhi + dEta*dEta;
772 if(TMath::Sqrt(arg) < conesize)
774 sumptPerp+=track->Pt();
782 ULong_t status = track->GetStatus();
783 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
784 Int_t trackBC = track->GetTOFBunchCrossing(bz);
787 if(okTOF && trackBC == 0)
792 sumptPerpBC0+=track->Pt();
795 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
802 sumptPerpITSSPD+=track->Pt();
805 if(okTOF && trackBC == 0 && !bConstrained)
810 sumptPerpBC0ITSSPD+=track->Pt();
816 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
817 arg = dPhi*dPhi + dEta*dEta;
819 if(TMath::Sqrt(arg) < conesize)
821 sumptPerp+=track->Pt();
829 coneptsumPerpTrackPerMinCut[icut]+=track->Pt();
830 coneNPerpTrackPerMinCut [icut]++;
841 ULong_t status = track->GetStatus();
842 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
843 Int_t trackBC = track->GetTOFBunchCrossing(bz);
846 if(okTOF && trackBC == 0)
851 sumptPerpBC0+=track->Pt();
854 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
861 sumptPerpITSSPD+=track->Pt();
864 if(okTOF && trackBC == 0 && !bConstrained)
869 sumptPerpBC0ITSSPD+=track->Pt();
919 Float_t etaUEptsumCluster = 0 ;
920 Float_t phiUEptsumCluster = 0 ;
929 Float_t coneA = conesize*conesize*TMath::Pi();
930 Float_t ptTrig = pCandidate->Pt() ;
931 Float_t phiTrig = pCandidate->Phi();
932 Float_t etaTrig = pCandidate->Eta();
938 Float_t phiUEptsumTrackNorm = 0 ;
939 Float_t etaUEptsumTrackNorm = 0 ;
940 Float_t coneptsumTrackSubPhi = 0 ;
941 Float_t coneptsumTrackSubEta = 0 ;
942 Float_t coneptsumTrackSubPhiNorm = 0 ;
943 Float_t coneptsumTrackSubEtaNorm = 0 ;
957 Float_t correctConeSumTrack = 1;
958 Float_t correctConeSumTrackPhi = 1;
961 phiUEptsumTrack,etaUEptsumTrack,
962 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
963 correctConeSumTrack,correctConeSumTrackPhi);
965 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
966 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
974 if(coneptsumTrack > 0)
976 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
977 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
998 Float_t phiUEptsumClusterNorm = 0 ;
999 Float_t etaUEptsumClusterNorm = 0 ;
1000 Float_t coneptsumClusterSubPhi = 0 ;
1001 Float_t coneptsumClusterSubEta = 0 ;
1002 Float_t coneptsumClusterSubPhiNorm = 0 ;
1003 Float_t coneptsumClusterSubEtaNorm = 0 ;
1004 Float_t phiUEptsumCellNorm = 0 ;
1005 Float_t etaUEptsumCellNorm = 0 ;
1006 Float_t coneptsumCellSubPhi = 0 ;
1007 Float_t coneptsumCellSubEta = 0 ;
1008 Float_t coneptsumCellSubPhiNorm = 0 ;
1009 Float_t coneptsumCellSubEtaNorm = 0 ;
1027 Float_t correctConeSumClusterEta = 1;
1028 Float_t correctConeSumClusterPhi = 1;
1031 phiUEptsumCluster,etaUEptsumCluster,
1032 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
1033 correctConeSumClusterEta,correctConeSumClusterPhi);
1044 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
1045 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
1052 if(coneptsumCluster!=0)
1054 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
1055 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
1087 Float_t emcPhiSize = TMath::DegToRad()*100.;
1089 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
1090 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
1094 Float_t correctConeSumCellEta = 1;
1095 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
1097 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
1101 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
1104 Float_t correctConeSumCellPhi = 1;
1106 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
1107 (phiTrig-conesize < 80*TMath::DegToRad()))
1110 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
1111 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
1117 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
1121 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
1122 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1133 if ( coneptsumCell > 0.01 )
1135 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1136 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1146 sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1147 sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1177 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1178 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1202 coneptLeadCluster = 0;
1203 coneptsumCluster = 0;
1207 Float_t ptTrig = aodParticle->Pt();
1218 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1219 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1230 Float_t coneNClusterPerMinCut [20];
1231 Float_t coneptsumClusterPerMinCut[20];
1232 Float_t coneptsumClusterPerMaxCut[20];
1233 Float_t coneptsumClusterPerRCut [10];
1235 Float_t coneptsumClusterPerNCellCut[20];
1236 Float_t coneptsumClusterPerExoCut [20];
1242 coneNClusterPerMinCut [icut] = 0;
1243 coneptsumClusterPerMinCut[icut] = 0;
1244 coneptsumClusterPerMaxCut[icut] = 0;
1252 coneptsumClusterPerRCut[icut] = 0;
1260 coneptsumClusterPerNCellCut[icut] = 0;
1265 coneptsumClusterPerExoCut[icut] = 0;
1270 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1272 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1299 coneptsumCluster+=ptcone;
1300 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1308 coneptsumClusterPerMinCut[icut]+=ptcone;
1309 coneNClusterPerMinCut [icut]++;
1312 if ( ptcone <
fMaxPtCutInCone[icut] ) coneptsumClusterPerMaxCut[icut]+=ptcone;
1323 coneptsumClusterPerRCut[icut]+=ptcone;
1336 coneptsumClusterPerNCellCut[icut]+=ptcone;
1346 coneptsumClusterPerExoCut[icut]+=ptcone;
1359 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1360 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1417 Float_t ptTrig = aodParticle->Pt();
1418 Float_t phiTrig = aodParticle->Phi();
1419 if(phiTrig<0) phiTrig += TMath::TwoPi();
1420 Float_t etaTrig = aodParticle->Eta();
1422 if(aodParticle->GetDetectorTag()==
kEMCAL)
1424 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1427 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1429 if(!eGeom->CheckAbsCellId(absId))
return ;
1432 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1434 Int_t imEta=-1, imPhi=-1;
1435 Int_t ieta =-1, iphi =-1;
1437 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1439 Int_t iEta=-1, iPhi=-1;
1440 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1442 Int_t colTrig = iEta;
1443 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1444 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1446 Int_t sqrSize = int(conesize/0.0143);
1451 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1453 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1455 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1456 if(inSector==5)
continue;
1458 Int_t inSupMod = -1;
1460 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1462 inSupMod = 2*inSector + 1;
1465 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1467 inSupMod = 2*inSector;
1468 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1471 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1473 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1474 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1477 coneptsumCell += cells->GetCellAmplitude(iabsId);
1493 coneptLeadTrack = 0;
1497 Float_t ptTrig = aodParticle->Pt();
1519 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1520 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1529 Float_t coneNTrackPerMinCut [20];
1530 Float_t coneptsumTrackPerMinCut[20];
1531 Float_t coneptsumTrackPerMaxCut[20];
1532 Float_t coneptsumTrackPerEtaCut[10];
1533 Float_t coneptsumTrackPerNCellCut[20];
1534 Float_t coneptsumTrackPerExoCut[20];
1535 Float_t coneptsumTrackPerRCut [10];
1536 Float_t coneptsumTrackTOFBC0 = 0;
1537 Float_t coneptsumTrackTOFBCN = 0;
1538 Float_t coneptsumTrackTOFNo = 0;
1539 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1540 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1541 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1542 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1548 coneptsumTrackPerMinCut[icut] = 0;
1549 coneptsumTrackPerMaxCut[icut] = 0;
1550 coneNTrackPerMinCut [icut] = 0;
1558 coneptsumTrackPerEtaCut[icut] = 0;
1566 coneptsumTrackPerRCut[icut] = 0;
1574 coneptsumTrackPerNCellCut[icut] = 0;
1579 coneptsumTrackPerExoCut[icut] = 0;
1583 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1585 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1587 pTtrack = track->Pt();
1594 Int_t partInConeCharge = 0, partInConePDG = 0;
1595 Bool_t physPrimary = kFALSE;
1596 Int_t trackLabel = TMath::Abs(track->GetLabel());
1600 TParticle * mcpartESD =
static_cast<TParticle*
>(
GetMCStack()->Particle(trackLabel));
1603 partInConePDG = mcpartESD->GetPdgCode();
1604 partInConeCharge = TMath::Abs((
Int_t) TDatabasePDG::Instance()->GetParticle(partInConePDG)->Charge());
1605 physPrimary = mcpartESD->IsPrimary();
1613 AliAODMCParticle * mcpartAOD = (AliAODMCParticle *) mcparticles->At(trackLabel);
1616 partInConeCharge = TMath::Abs(mcpartAOD->Charge());
1617 partInConePDG = mcpartAOD->GetPdgCode();
1618 physPrimary = mcpartAOD->IsPhysicalPrimary();
1623 if ( partInConeCharge > 0 && TMath::Abs(partInConePDG) != 11 )
1626 if ( TMath::Abs(partInConePDG) == 211 ) mcChTag = 0;
1627 else if ( TMath::Abs(partInConePDG) == 321 ) mcChTag = 1;
1628 else if ( TMath::Abs(partInConePDG) == 2212 ) mcChTag = 2;
1644 coneptsumTrack+=pTtrack;
1645 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1653 coneptsumTrackPerMinCut[icut]+=pTtrack;
1654 coneNTrackPerMinCut [icut]++;
1657 if ( pTtrack <
fMaxPtCutInCone[icut] ) coneptsumTrackPerMaxCut[icut]+=pTtrack;
1665 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1676 coneptsumTrackPerRCut[icut]+=pTtrack;
1689 coneptsumTrackPerNCellCut[icut]+=pTtrack;
1699 coneptsumTrackPerExoCut[icut]+=pTtrack;
1710 phitrack = track->Phi();
1711 etatrack = track->Eta();
1718 ULong_t status = track->GetStatus();
1719 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1720 trackBC = track->GetTOFBunchCrossing(bz);
1721 Double32_t tof = track->GetTOFsignal()*1e-3;
1740 coneptsumTrackTOFBC0 += pTtrack;
1748 coneptsumTrackTOFBCN += pTtrack;
1757 coneptsumTrackTOFNo += pTtrack;
1761 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1762 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1770 coneptsumTrackITSRefitOnSPDOff += pTtrack;
1778 coneptsumTrackITSRefitOffSPDOff += pTtrack;
1787 coneptsumTrackITSRefitOnSPDOn += pTtrack;
1794 if(okTOF && trackBC == 0 && !bConstrained)
1800 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
1843 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
1844 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
1908 if(energy < 8)
return;
1914 Int_t nVerticesSPD = -1;
1915 Int_t nVerticesTracks = -1;
1919 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
1920 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
1925 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
1926 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
1938 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
1942 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
1943 ncont=pv->GetNContributors();
1944 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
1946 diamZ = esdEv->GetDiamondZ();
1950 AliAODVertex *pv=aodEv->GetVertex(iVert);
1951 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
1952 ncont=pv->GetNContributors();
1953 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
1955 diamZ = aodEv->GetDiamondZ();
1958 Double_t distZ = TMath::Abs(z2-z1);
1959 diamZ = TMath::Abs(z2-diamZ);
1972 (AliAODPWG4ParticleCorrelation *pCandidate,
1976 if(!fFillTMHisto && !fFillSSHisto && !fFillBackgroundBinHistograms && !fFillTaggedDecayHistograms)
return;
1978 Int_t nMaxima = pCandidate->GetNLM();
1979 Int_t mcTag = pCandidate->GetTag() ;
1980 Bool_t isolated = pCandidate->IsIsolated();
1982 Float_t m02 = pCandidate->GetM02() ;
1984 Float_t pt = pCandidate->Pt();
1985 Float_t eta = pCandidate->Eta();
1986 Float_t phi = pCandidate->Phi();
1987 if(phi<0) phi+= TMath::TwoPi();
1991 if(fFillTaggedDecayHistograms)
1993 Int_t decayTag = pCandidate->DecayTag();
1994 if(decayTag < 0) decayTag = 0;
1996 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
1998 if(!GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
continue;
2000 if(fFillSSHisto) fhPtLambda0Decay[isolated][ibit]->Fill(pt, m02, GetEventWeight());
2004 if( m02 > fDecayTagsM02Cut )
continue;
2006 fhPtDecay [isolated][ibit]->Fill(pt, GetEventWeight());
2007 fhEtaPhiDecay[isolated][ibit]->Fill(eta, phi, GetEventWeight());
2011 fhPtDecayMC[isolated][ibit][mcIndex]->Fill(pt, GetEventWeight());
2014 fhPtDecayMC[isolated][ibit][kmcPhoton]->Fill(pt, GetEventWeight());
2018 if ( mcIndex == kmcPi0Decay ) fhPtDecayMC[isolated][ibit][kmcPi0DecayLostPair]->Fill(pt, GetEventWeight());
2019 else if( mcIndex == kmcEtaDecay ) fhPtDecayMC[isolated][ibit][kmcEtaDecayLostPair]->Fill(pt, GetEventWeight());
2028 if(fFillBackgroundBinHistograms)
2031 Int_t ptsumBin = -1;
2032 Int_t leadptBin = -1;
2034 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,fNBkgBin));
2036 for(
Int_t ibin = 0; ibin < fNBkgBin; ibin++)
2038 if( coneptsum >= fBkgBinLimit[ibin] && coneptsum < fBkgBinLimit[ibin+1]) ptsumBin = ibin;
2039 if( coneleadpt >= fBkgBinLimit[ibin] && coneleadpt < fBkgBinLimit[ibin+1]) leadptBin = ibin;
2044 if ( leadptBin >= 0 )
2046 AliDebug(1,Form(
"\t Lead bin %d [%2.2f,%2.2f]", leadptBin,fBkgBinLimit[leadptBin],fBkgBinLimit[leadptBin+1]));
2048 fhPtLeadConeBin[leadptBin]->Fill(pt, GetEventWeight());
2051 fhPtLeadConeBinLambda0[leadptBin]->Fill(pt, m02, GetEventWeight());
2053 if ( leadptBin == 0 )
2054 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
2055 pt, GetTrackMultiplicity(),GetEMCALClusters()->GetEntriesFast()));
2058 if ( ptsumBin >= 0 )
2060 AliDebug(1,Form(
"\t Sum bin %d [%2.2f,%2.2f]" , ptsumBin ,fBkgBinLimit[ptsumBin] ,fBkgBinLimit[ptsumBin +1]));
2062 fhSumPtConeBin[ptsumBin]->Fill(pt, GetEventWeight());
2064 if(fFillSSHisto) fhSumPtConeBinLambda0[ptsumBin]->Fill(pt, m02, GetEventWeight());
2068 if( fFillTaggedDecayHistograms && m02 < fDecayTagsM02Cut )
2070 Int_t decayTag = pCandidate->DecayTag();
2071 if(decayTag < 0) decayTag = 0;
2073 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
2075 if(GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
2077 Int_t leadptBinDecay = leadptBin+ibit*fNBkgBin;
2078 Int_t ptsumBinDecay = ptsumBin+ibit*fNBkgBin;
2079 if( leadptBin >=0 ) fhPtLeadConeBinDecay[leadptBinDecay]->Fill(pt, GetEventWeight());
2080 if( ptsumBin >=0 ) fhSumPtConeBinDecay [ ptsumBinDecay]->Fill(pt, GetEventWeight());
2087 Int_t leadptBinMC = leadptBin+mcIndex*fNBkgBin;
2088 Int_t ptsumBinMC = ptsumBin+mcIndex*fNBkgBin;
2092 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2093 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2098 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt, GetEventWeight());
2099 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2104 leadptBinMC = leadptBin+kmcPhoton*fNBkgBin;
2105 ptsumBinMC = ptsumBin+kmcPhoton*fNBkgBin;
2108 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2109 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2114 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt, GetEventWeight());
2115 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2122 if ( mcIndex == kmcPi0Decay )
2124 leadptBinMC = leadptBin+kmcPi0DecayLostPair*fNBkgBin;
2125 ptsumBinMC = ptsumBin+kmcPi0DecayLostPair*fNBkgBin;
2127 else if(mcIndex == kmcEtaDecay)
2129 leadptBinMC = leadptBin+kmcEtaDecayLostPair*fNBkgBin;
2130 ptsumBinMC = ptsumBin+kmcEtaDecayLostPair*fNBkgBin;
2133 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
2137 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2138 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2143 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt);
2144 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2153 Int_t ptTrigBin = -1;
2154 if(fFillPtTrigBinHistograms)
2156 for(
Int_t ibin = 0; ibin < fNPtTrigBin; ibin++)
2158 if( pt >= fPtTrigBinLimit[ibin] && pt < fPtTrigBinLimit[ibin+1]) ptTrigBin = ibin;
2162 if ( ptTrigBin >= 0 )
2164 AliDebug(1,Form(
"Trigger pT %f, bin %d [%2.2f,%2.2f]",pt,ptTrigBin,fPtTrigBinLimit[ptTrigBin],fPtTrigBinLimit[ptTrigBin+1]));
2166 fhPtTrigBinPtLeadCone[ptTrigBin]->Fill(coneleadpt, GetEventWeight());
2167 fhPtTrigBinSumPtCone [ptTrigBin]->Fill(coneptsum , GetEventWeight());
2171 fhPtTrigBinSumPtTrackCone [ptTrigBin]->Fill(coneptsumTrack, GetEventWeight());
2172 fhPtTrigBinSumPtClusterCone[ptTrigBin]->Fill(coneptsumClust, GetEventWeight());
2177 fhPtTrigBinLambda0vsPtLeadCone[ptTrigBin]->Fill(coneleadpt, m02, GetEventWeight());
2178 fhPtTrigBinLambda0vsSumPtCone [ptTrigBin]->Fill(coneptsum , m02, GetEventWeight());
2182 fhPtTrigBinLambda0vsSumPtTrackCone [ptTrigBin]->Fill(coneptsumTrack, m02, GetEventWeight());
2183 fhPtTrigBinLambda0vsSumPtClusterCone[ptTrigBin]->Fill(coneptsumClust, m02, GetEventWeight());
2188 if( fFillTaggedDecayHistograms && m02 < fDecayTagsM02Cut )
2190 Int_t decayTag = pCandidate->DecayTag();
2191 if(decayTag < 0) decayTag = 0;
2193 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
2195 if(GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
2197 Int_t binDecay = ptTrigBin+ibit*fNPtTrigBin;
2200 fhPtTrigBinPtLeadConeDecay[binDecay]->Fill(coneleadpt, GetEventWeight());
2201 fhPtTrigBinSumPtConeDecay [binDecay]->Fill(coneptsum , GetEventWeight());
2205 fhPtTrigBinSumPtTrackConeDecay [binDecay]->Fill(coneptsumTrack, GetEventWeight());
2206 fhPtTrigBinSumPtClusterConeDecay[binDecay]->Fill(coneptsumClust, GetEventWeight());
2215 Int_t ptTrigBinMC = ptTrigBin+mcIndex*fNPtTrigBin;
2217 fhPtTrigBinPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt , GetEventWeight());
2218 fhPtTrigBinSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , GetEventWeight());
2222 fhPtTrigBinSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, GetEventWeight());
2223 fhPtTrigBinSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, GetEventWeight());
2228 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt , m02, GetEventWeight());
2229 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , m02, GetEventWeight());
2233 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2234 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2240 ptTrigBinMC = ptTrigBin+kmcPhoton*fNPtTrigBin;
2242 fhPtTrigBinPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt, GetEventWeight());
2243 fhPtTrigBinSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , GetEventWeight());
2247 fhPtTrigBinSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, GetEventWeight());
2248 fhPtTrigBinSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, GetEventWeight());
2253 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt, m02, GetEventWeight());
2254 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , m02, GetEventWeight());
2258 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2259 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2267 if ( mcIndex == kmcPi0Decay ) ptTrigBinMC = ptTrigBin+kmcPi0DecayLostPair*fNPtTrigBin;
2268 else if( mcIndex == kmcEtaDecay ) ptTrigBinMC = ptTrigBin+kmcEtaDecayLostPair*fNPtTrigBin;
2270 fhPtTrigBinPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt , GetEventWeight());
2271 fhPtTrigBinSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , GetEventWeight());
2275 fhPtTrigBinSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, GetEventWeight());
2276 fhPtTrigBinSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, GetEventWeight());
2281 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt , m02, GetEventWeight());
2282 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , m02, GetEventWeight());
2286 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2287 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2304 fhELambda0 [isolated]->Fill(energy, m02, GetEventWeight());
2305 fhPtLambda0[isolated]->Fill(pt, m02, GetEventWeight());
2314 fhPtLambda0MC[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2318 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MC[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2319 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MC[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2322 fhPtLambda0MC[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2327 fhPtLambda0MCConv[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2331 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConv[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2332 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConv[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2335 fhPtLambda0MCConv[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2341 Int_t noverlaps = 0;
2342 if ( fFillOverlapHistograms && fCluster )
2344 const UInt_t nlabels = fCluster->GetNLabels();
2345 Int_t overpdg[nlabels];
2346 Int_t overlab[nlabels];
2347 noverlaps = GetMCAnalysisUtils()->GetNOverlaps(fCluster->GetLabels(), nlabels, mcTag, -1,
2348 GetReader(), overpdg, overlab);
2351 fhPtNOverlap[kmcPhoton][isolated]->Fill(pt, noverlaps, GetEventWeight());
2355 if ( mcIndex == kmcPi0Decay ) fhPtNOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2356 else if( mcIndex == kmcEtaDecay ) fhPtNOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2359 fhPtNOverlap[mcIndex][isolated]->Fill(pt, noverlaps, GetEventWeight());
2364 fhPtNOverlapConv[kmcPhoton][isolated]->Fill(pt, noverlaps, GetEventWeight());
2368 if ( mcIndex == kmcPi0Decay ) fhPtNOverlapConv[kmcPi0DecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2369 else if( mcIndex == kmcEtaDecay ) fhPtNOverlapConv[kmcEtaDecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2372 fhPtNOverlapConv[mcIndex][isolated]->Fill(pt, noverlaps, GetEventWeight());
2376 if ( noverlaps == 1 )
2379 fhPtLambda0MCWith1Overlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2383 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCWith1Overlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2384 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCWith1Overlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2387 fhPtLambda0MCWith1Overlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2392 fhPtLambda0MCConvWith1Overlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2396 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConvWith1Overlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2397 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConvWith1Overlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2400 fhPtLambda0MCConvWith1Overlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2403 else if (noverlaps == 0 )
2406 fhPtLambda0MCWithNoOverlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2410 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCWithNoOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2411 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCWithNoOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2414 fhPtLambda0MCWithNoOverlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2419 fhPtLambda0MCConvWithNoOverlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2423 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConvWithNoOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2424 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConvWithNoOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2427 fhPtLambda0MCConvWithNoOverlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2435 if(fFillPtTrigBinHistograms && fFillSSHisto)
2437 Int_t ptTrigBinMC = ptTrigBin+mcIndex *fNPtTrigBin;
2438 Int_t ptTrigBinMCPhoton = ptTrigBin+kmcPhoton*fNPtTrigBin;
2439 Int_t ptTrigBinMCPi0Lost = ptTrigBin+kmcPi0DecayLostPair*fNPtTrigBin;
2440 Int_t ptTrigBinMCEtaLost = ptTrigBin+kmcEtaDecayLostPair*fNPtTrigBin;
2442 if ( ptTrigBin >= 0 )
2444 fhPtTrigBinLambda0vsSumPtConeMC[ptTrigBinMC]->Fill(coneptsum, m02, GetEventWeight());
2448 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2449 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2454 fhPtTrigBinLambda0vsSumPtConeMC[ptTrigBinMCPhoton]->Fill(coneptsum, m02, GetEventWeight());
2458 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2459 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2463 if( mcIndex == kmcPi0Decay )
2465 fhPtTrigBinLambda0vsSumPtConeMC[ptTrigBinMCPi0Lost]->Fill(coneptsum, m02, GetEventWeight());
2469 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2470 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2474 if( mcIndex == kmcEtaDecay )
2476 fhPtTrigBinLambda0vsSumPtConeMC[ptTrigBinMCEtaLost]->Fill(coneptsum, m02, GetEventWeight());
2480 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2481 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2485 if(fFillOverlapHistograms)
2487 if ( noverlaps == 0 )
2489 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMC]->Fill(coneptsum, m02, GetEventWeight());
2493 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2494 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2499 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCPhoton]->Fill(coneptsum, m02, GetEventWeight());
2503 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2504 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2508 if( mcIndex == kmcPi0Decay )
2510 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCPi0Lost]->Fill(coneptsum, m02, GetEventWeight());
2514 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2515 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2519 if( mcIndex == kmcEtaDecay )
2521 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCEtaLost]->Fill(coneptsum, m02, GetEventWeight());
2525 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2526 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2532 else if ( noverlaps == 1 )
2534 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMC]->Fill(coneptsum, m02, GetEventWeight());
2537 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2538 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2543 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCPhoton]->Fill(coneptsum, m02, GetEventWeight());
2547 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2548 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2551 if( mcIndex == kmcPi0Decay )
2553 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCPi0Lost]->Fill(coneptsum, m02, GetEventWeight());
2557 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2558 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2562 if( mcIndex == kmcEtaDecay )
2564 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCEtaLost]->Fill(coneptsum, m02, GetEventWeight());
2568 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2569 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2581 if(GetCalorimeter() == kEMCAL && GetFirstSMCoveredByTRD() >= 0 &&
2582 GetModuleNumber(pCandidate) >= GetFirstSMCoveredByTRD() )
2584 fhELambda0TRD [isolated]->Fill(energy, m02, GetEventWeight());
2585 fhPtLambda0TRD[isolated]->Fill(pt , m02, GetEventWeight());
2589 if(fFillNLMHistograms)
2591 fhNLocMax[isolated]->Fill(energy, nMaxima, GetEventWeight());
2594 fhELambda0LocMax1[isolated]->Fill(energy, m02, GetEventWeight());
2595 fhELambda1LocMax1[isolated]->Fill(energy, m02, GetEventWeight());
2599 fhELambda0LocMax2[isolated]->Fill(energy, m02, GetEventWeight());
2600 fhELambda1LocMax2[isolated]->Fill(energy, m02, GetEventWeight());
2604 fhELambda0LocMaxN[isolated]->Fill(energy, m02, GetEventWeight());
2605 fhELambda1LocMaxN[isolated]->Fill(energy, m02, GetEventWeight());
2611 if ( fFillTMHisto && fCluster )
2613 Float_t dZ = fCluster->GetTrackDz();
2614 Float_t dR = fCluster->GetTrackDx();
2623 if(fhTrackMatchedDEta[isolated] && TMath::Abs(dR) < 999)
2625 fhTrackMatchedDEta[isolated]->Fill(energy, dZ, GetEventWeight());
2626 fhTrackMatchedDPhi[isolated]->Fill(energy, dR, GetEventWeight());
2627 if(energy > 0.5) fhTrackMatchedDEtaDPhi[isolated]->Fill(dZ, dR, GetEventWeight());
2630 fhTrackMatchedDEtaMC[mcIndex][isolated]->Fill(energy, dZ, GetEventWeight());
2631 fhTrackMatchedDPhiMC[mcIndex][isolated]->Fill(energy, dR, GetEventWeight());
2632 if(energy > 0.5) fhTrackMatchedDEtaDPhiMC[mcIndex][isolated]->Fill(dZ, dR, GetEventWeight());
2638 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2640 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(fCluster, GetReader()->GetInputEvent());
2644 Float_t dEdx = track->GetTPCsignal();
2645 fhdEdx[isolated]->Fill(fCluster->E(), dEdx, GetEventWeight());
2647 Float_t eOverp = fCluster->E()/track->P();
2648 fhEOverP[isolated]->Fill(fCluster->E(), eOverp, GetEventWeight());
2659 GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCEta) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 2.5, GetEventWeight());
2660 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCPhoton) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 0.5, GetEventWeight());
2661 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCElectron) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 1.5, GetEventWeight());
2662 else fhTrackMatchedMCParticle[isolated]->Fill(energy, 3.5, GetEventWeight());
2668 GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCEta) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 6.5, GetEventWeight());
2669 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCPhoton) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 4.5, GetEventWeight());
2670 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCElectron) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 5.5, GetEventWeight());
2671 else fhTrackMatchedMCParticle[isolated]->Fill(energy, 7.5, GetEventWeight());
2684 const Int_t buffersize = 255;
2685 char onePar[buffersize] ;
2687 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2691 snprintf(onePar, buffersize,
"Isolation Cand. Detector: %s;",
fIsoDetectorString.Data()) ;
2695 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2697 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2699 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2701 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2706 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2708 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2713 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2718 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2723 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2728 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2739 return new TObjString(parList) ;
2749 outputContainer->SetName(
"IsolatedParticleHistos") ;
2802 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
2805 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
2810 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
2813 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
2817 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
2820 TString sParticle =
", x^{0,#pm}";
2824 TString parTitle = Form(
"#it{R} = %2.2f%s%s" ,
GetIsolationCut()->GetConeSize(),sThreshold.Data(),sParticle.Data());
2829 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
2832 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
2833 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
2834 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
2835 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
2837 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
2838 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
2839 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
2840 "OtherDecay",
"Electron",
"Hadron"} ;
2843 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
2844 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
2845 "#pi^{0}" ,
"#eta"} ;
2847 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
2848 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
2854 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2855 nptbins,ptmin,ptmax);
2856 fhENoIso->SetYTitle(
"#it{counts}");
2857 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
2861 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2862 nptbins,ptmin,ptmax);
2864 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2870 (
"hENoIsoExoTrigger",
2871 Form(
"Number of not isolated particles vs E, %s, exot>0.97",parTitle.Data()),
2872 nptbins,ptmin,ptmax);
2878 (
"hPtNoIsoExoTrigger",
2879 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
2880 nptbins,ptmin,ptmax);
2887 Form(
"Number of not isolated leading particles #eta vs #varphi, %s",parTitle.Data()),
2888 netabins,etamin,etamax,nphibins,phimin,phimax);
2899 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
2900 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2901 nptbins,ptmin,ptmax);
2903 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2906 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
2907 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2908 nptbins,ptmin,ptmax);
2909 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
2910 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2914 Form(
"#varphi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2915 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2917 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2921 Form(
"#eta vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2922 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2924 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2932 TString isoName [] = {
"NoIso",
"Iso"};
2933 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2937 for(
Int_t iso = 0; iso < 2; iso++)
2941 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2942 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2943 nptbins,ptmin,ptmax);
2944 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2945 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2946 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2950 Form(
"Number of %s leading Pi0 decay particles #eta vs #varphi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2951 netabins,etamin,etamax,nphibins,phimin,phimax);
2959 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2960 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2961 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2972 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2973 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2974 nptbins,ptmin,ptmax);
2975 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2976 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2977 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2986 TString isoName [] = {
"NoIso",
"Iso"};
2987 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2990 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2991 nptbins,ptmin,ptmax);
2992 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2993 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2994 outputContainer->Add(
fhEIso) ;
2997 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2998 nptbins,ptmin,ptmax);
2999 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
3000 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3001 outputContainer->Add(
fhPtIso) ;
3007 Form(
"Number of isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3008 nptbins,ptmin,ptmax);
3014 (
"hPtIsoExoTrigger",
3015 Form(
"Number of isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3016 nptbins,ptmin,ptmax);
3023 Form(
"Number of isolated particles vs #varphi, %s",parTitle.Data()),
3024 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3025 fhPhiIso->SetYTitle(
"#varphi (rad)");
3026 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3030 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
3031 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3033 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3037 Form(
"Number of isolated particles #eta vs #varphi, %s",parTitle.Data()),
3038 netabins,etamin,etamax,nphibins,phimin,phimax);
3046 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3047 nptbins,ptmin,ptmax, 100,0,100);
3053 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3054 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
3063 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3064 nptbins,ptmin,ptmax,10,0,10);
3069 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3070 nptbins,ptmin,ptmax,10,0,10);
3077 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3078 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3079 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
3080 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3085 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3086 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3087 fhConeSumPt->SetYTitle(
"#Sigma #it{p}_{T} (GeV/#it{c})");
3088 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3094 (
"hConePtSumExoTrigger",
3095 Form(
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3096 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3111 Form(
"Trigger #eta vs #varphi, #Sigma #it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3112 netabins,etamin,etamax,nphibins,phimin,phimax);
3119 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s",parTitleR.Data()),
3120 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3121 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
3122 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3128 Form(
"#it{p}_{T} in isolation cone for %s, exotic trigger",parTitleR.Data()),
3129 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3166 (Form(
"hPtLeadCone_Bin%d",ibin),
3167 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3174 (Form(
"hSumPtCone_Bin%d",ibin),
3175 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3188 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3189 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3196 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3197 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3211 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3212 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3219 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3220 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3231 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
3232 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3239 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
3240 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3252 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3253 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3260 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3261 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3338 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
3339 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3346 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
3347 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3356 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d",ibin),
3357 Form(
"#Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3364 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d",ibin),
3365 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3379 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3380 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3387 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3388 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3397 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3398 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3405 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3406 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3421 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3422 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3429 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3430 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3439 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3440 Form(
"#Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3447 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3448 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3460 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
3461 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3468 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
3469 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3478 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d",ibin),
3479 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3486 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d",ibin),
3487 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3500 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3501 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3503 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3509 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3510 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3512 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3520 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3521 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3523 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3529 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3530 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3532 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3541 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3542 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, No Overlaps %s, %s",
3544 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3550 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3551 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, 1 Overlap %s, %s",
3553 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3561 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3562 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, No Overlaps %s, %s",
3564 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3570 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3571 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, No Overlaps %s, %s",
3573 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3579 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3580 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, 1 Overlap %s, %s",
3582 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3588 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3589 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC, 1 Overlap %s, %s",
3591 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3607 Form(
"#it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3608 100,0,100,nptinconebins,ptinconemin,ptinconemax);
3618 (
"hConePtSumCluster",
3619 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3620 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3628 (
"hConePtSumClusterExoTrigger",
3629 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3630 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3639 (
"hConeNClusterPerMinPtCut",
3640 Form(
"N clusters, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3649 (
"hConeNClusterPerMinPtCutLargePtTrig",
3650 Form(
"N cluster, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3659 (
"hConePtSumClusterPerMinPtCut",
3660 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3669 (
"hConePtSumClusterPerMinPtCutLargePtTrig",
3670 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3680 (
"hConePtSumClusterPerMaxPtCut",
3681 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3690 (
"hConePtSumClusterPerMaxPtCutLargePtTrig",
3691 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3703 (
"hConePtSumClusterPerRCut",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts",
3712 (
"hConePtSumClusterPerRCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3721 (
"hPtClusterInConePerRCut",
"Cluster #it{p}_{T}, different #it{R} cuts",
3730 (
"hPtClusterInConePerRCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3742 (
"hConePtSumClusterPerNCellCut",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
3751 (
"hConePtSumClusterPerNCellCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3760 (
"hPtClusterInConePerNCellCut",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts",
3769 (
"hPtClusterInConePerNCellCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3779 (
"hConePtSumClusterPerExoCut",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts",
3788 (
"hConePtSumClusterPerExoCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
3797 (
"hPtClusterInConePerExoCut",
"Cluster #it{p}_{T}, different exoticity cuts",
3806 (
"hPtClusterInConePerExoCutLargePtTrig",
"Cluster #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
3816 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
3817 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3826 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3827 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3836 "#Sigma cluster #it{p}_{T} in UE Eta Band",
3837 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3843 "#Sigma cluster #it{p}_{T} UE Phi Band",
3844 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3852 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} in UE Eta Band",
3853 netabins,etamin,etamax,nphibins,phimin,phimax);
3860 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} UE Phi Band",
3861 netabins,etamin,etamax,nphibins,phimin,phimax);
3868 Form(
"#eta vs #varphi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
3869 netabins,-1,1,nphibins,0,TMath::TwoPi());
3875 Form(
"#eta vs #varphi of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
3876 netabins,-1,1,nphibins,0,TMath::TwoPi());
3882 Form(
"#eta vs #varphi of clusters in cone for #it{R} = %2.2f",r),
3883 netabins,-1,1,nphibins,0,TMath::TwoPi());
3889 Form(
"#eta vs #varphi of all clusters"),
3890 netabins,-1,1,nphibins,0,TMath::TwoPi());
3899 "#Sigma cell #it{p}_{T} in UE Eta Band",
3900 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3906 "#Sigma cell #it{p}_{T} UE Phi Band",
3907 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3913 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} in UE Eta Band",
3914 netabins,etamin,etamax,nphibins,phimin,phimax);
3921 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} UE Phi Band",
3922 netabins,etamin,etamax,nphibins,phimin,phimax);
3930 Form(
"#it{p}_{T} of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
3931 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3937 Form(
"#it{p}_{T} of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
3938 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3945 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
3946 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3954 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f, exotic trigger",r),
3955 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3964 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
3965 nptbins,ptmin,ptmax,1000,0,50);
3971 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
3978 Form(
"#col vs #row of cells in #varphi band isolation cone for #it{R} = %2.2f",r),
3988 Form(
"Clusters #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
3989 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
3995 Form(
"Clusters #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
3996 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4002 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4003 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4009 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4010 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4018 Form(
"Trigger #eta vs #varphi, Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4019 netabins,etamin,etamax,nphibins,phimin,phimax);
4026 Form(
"Trigger #eta vs #varphi, Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4027 netabins,etamin,etamax,nphibins,phimin,phimax);
4037 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4038 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4044 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4045 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4051 Form(
"Trigger #eta vs #varphi, Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4052 netabins,etamin,etamax,nphibins,phimin,phimax);
4059 Form(
"Trigger #eta vs #varphi, Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4060 netabins,etamin,etamax,nphibins,phimin,phimax);
4070 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
4071 nptbins,ptmin,ptmax,100,0,1);
4077 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#varphi ",r),
4078 netabins,etamin,etamax,nphibins,phimin,phimax);
4085 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance",r),
4086 nptbins,ptmin,ptmax,100,0,1);
4092 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance, in trigger #eta-#varphi ",r),
4093 netabins,etamin,etamax,nphibins,phimin,phimax);
4101 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #varphi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4102 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4108 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #varphi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4109 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4115 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #eta band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4116 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4122 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #eta band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4123 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4129 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #eta band for cluster (before normalization), R=%2.2f",r),
4130 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
4136 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #varphi band for cluster (before normalization), R=%2.2f",r),
4137 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
4146 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
4147 nptbins,ptmin,ptmax,100,0,1);
4153 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#varphi ",r),
4154 netabins,etamin,etamax,nphibins,phimin,phimax);
4161 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance",r),
4162 nptbins,ptmin,ptmax,100,0,1);
4168 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance, in trigger #eta-#varphi ",r),
4169 netabins,etamin,etamax,nphibins,phimin,phimax);
4177 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #varphi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4178 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4184 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #varphi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4185 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4191 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #eta band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4192 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4198 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #eta band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
4199 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4211 (
"hConeLeadPtTrack",
4212 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
4213 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4220 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
4221 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4228 TString mcChPartName[] = {
"Pion",
"Kaon",
"Proton",
"Other"};
4229 for(
Int_t imc = 0; imc < 4; imc++)
4232 (Form(
"hPtTrackInCone_Primary_%s",mcChPartName[imc].
Data()),
4233 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4234 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4240 (Form(
"hPtTrackInCone_Secondary_%s",mcChPartName[imc].
Data()),
4241 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4242 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4248 (Form(
"hPtTrackInCone_Gener_Primary_%s",mcChPartName[imc].
Data()),
4249 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4250 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4256 (Form(
"hPtTrackInCone_Gener_Secondary_%s",mcChPartName[imc].
Data()),
4257 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4258 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4269 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4270 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4279 (
"hConeNTrackPerMinPtCut",
4280 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4289 (
"hConeNTrackPerMinPtCutLargePtTrig",
4290 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4299 (
"hConePtSumTrackPerMinPtCut",
4300 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4309 (
"hConePtSumTrackPerMinPtCutLargePtTrig",
4310 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4321 (
"hPerpConeNTrackPerMinPtCut",
4322 Form(
"N tracks, different #it{p}_{T} cuts in perpendicular cone for #it{R} = %2.2f",r),
4331 (
"hPerpConeNTrackPerMinPtCutLargePtTrig",
4332 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4341 (
"hPerpConePtSumTrackPerMinPtCut",
4342 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in perpendicular cone for #it{R} = %2.2f",r),
4351 (
"hPerpConePtSumTrackPerMinPtCutLargePtTrig",
4352 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in perpendicular cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4362 (
"hConePtSumTrackPerMaxPtCut",
4363 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4372 (
"hConePtSumTrackPerMaxPtCutLargePtTrig",
4373 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4385 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f",r),
4394 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4406 (
"hConePtSumTrackPerRCut",
"Track #Sigma #it{p}_{T}, different #it{R} cuts",
4415 (
"hConePtSumTrackPerRCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4424 (
"hPtTrackInConePerRCut",
"Track #it{p}_{T}, different #it{R} cuts",
4433 (
"hPtTrackInConePerRCutLargePtTrig",
"Track #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4445 (
"hConePtSumTrackPerNCellCut",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
4454 (
"hConePtSumTrackPerNCellCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4463 (
"hPtTrackInConePerNCellCut",
"Track #it{p}_{T}, different #it{N}_{cell} cuts",
4472 (
"hPtTrackInConePerNCellCutLargePtTrig",
"Track #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4481 (
"hConePtSumTrackPerExoCut",
"Track #Sigma #it{p}_{T}, different exoticity cuts",
4490 (
"hConePtSumTrackPerExoCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4499 (
"hPtTrackInConePerExoCut",
"Track #it{p}_{T}, different exoticity cuts",
4508 (
"hPtTrackInConePerExoCutLargePtTrig",
"Track #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4520 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4521 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4529 (
"hConePtSumTrackExoTrigger",
4530 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
4531 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4539 Int_t ntofbins = 1000;
4540 Int_t mintof = -500;
4544 (
"hTrackTOFInCone",
"TOF signal vs track #it{p}_{T}",
4545 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4553 (
"hTrackTOFInConeExoTrigger",
"TOF signal vs track #it{p}_{T}, exoticity > 0.97",
4554 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4561 (
"hTrackTOFInConeBC0",
"TOF signal vs track #it{p}_{T}, BC=0",
4562 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4568 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
4569 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4575 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f",r),
4576 netabins,-1,1,nphibins,0,TMath::TwoPi());
4582 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4583 nptbins,ptmin,ptmax,netabins,-1,1);
4589 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4590 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4601 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0",r),
4602 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4608 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC!=0",r),
4609 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4615 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no TOF",r),
4616 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4622 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
4623 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4629 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
4630 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4636 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no TOF",r),
4637 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4644 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4645 netabins,-1,1,nphibins,0,TMath::TwoPi());
4651 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4652 netabins,-1,1,nphibins,0,TMath::TwoPi());
4658 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4659 netabins,-1,1,nphibins,0,TMath::TwoPi());
4665 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4666 nptbins,ptmin,ptmax,netabins,-1,1);
4672 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4673 nptbins,ptmin,ptmax,netabins,-1,1);
4679 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4680 nptbins,ptmin,ptmax,netabins,-1,1);
4686 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4687 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4693 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4694 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4700 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4701 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4710 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD On",r),
4711 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4717 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD Off",r),
4718 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4724 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no ITS Refit SPD Off",r),
4725 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4731 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
4732 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4738 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
4739 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4745 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4746 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4753 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4754 netabins,-1,1,nphibins,0,TMath::TwoPi());
4760 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4761 netabins,-1,1,nphibins,0,TMath::TwoPi());
4767 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4768 netabins,-1,1,nphibins,0,TMath::TwoPi());
4774 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4775 nptbins,ptmin,ptmax,netabins,-1,1);
4781 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4782 nptbins,ptmin,ptmax,netabins,-1,1);
4788 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4789 nptbins,ptmin,ptmax,netabins,-1,1);
4795 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4796 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4802 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4803 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4809 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4810 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4819 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0, track ITS Refit SPD On",r),
4820 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);