17 #include <TClonesArray.h>
19 #include <TObjString.h>
24 #include "TDatabasePDG.h"
29 #include "AliMCEvent.h"
34 #include "AliVParticle.h"
35 #include "AliAODPWG4ParticleCorrelation.h"
37 #include "AliVTrack.h"
38 #include "AliVCluster.h"
39 #include "AliESDEvent.h"
40 #include "AliAODEvent.h"
43 #include "AliEMCALGeometry.h"
44 #include "AliPHOSGeoUtils.h"
55 fIsoDetector(-1), fIsoDetectorString(
""),
56 fReMakeIC(0), fMakeSeveralIC(0),
57 fFillTMHisto(0), fFillSSHisto(1),
58 fFillEMCALRegionHistograms(0), fFillUEBandSubtractHistograms(1),
59 fFillCellHistograms(0),
60 fFillOverlapHistograms(0),
61 fStudyTracksInCone(0), fStudyMCConversionRadius(0),
62 fFillTaggedDecayHistograms(0), fNDecayBits(0),
63 fDecayBits(), fDecayTagsM02Cut(0),
64 fFillNLMHistograms(0),
65 fLeadingOnly(0), fCheckLeadingWithNeutralClusters(0),
66 fSelectPrimariesInCone(0), fMakePrimaryPi0DecayStudy(0),
67 fFillBackgroundBinHistograms(0), fNBkgBin(0),
68 fFillPtTrigBinHistograms(0), fNPtTrigBin(0),
69 fMinCellsAngleOverlap(0),
71 fNCones(0), fNPtThresFrac(0),
72 fConeSizes(), fPtThresholds(),
73 fPtFractions(), fSumPtThresholds(),
74 fStudyPtCutInCone(0), fNPtCutsInCone(0),
75 fMinPtCutInCone(), fMaxPtCutInCone(),
76 fStudyEtaCutInCone(0), fNEtaCutsInCone(0), fEtaCutInCone(),
77 fStudyRCutInCone(0), fNRCutsInCone(0), fRCutInCone(),
78 fStudyNCellsCut(0), fNNCellsInCandidate(0), fNCellsInCandidate(),
79 fNCellsWithWeight(0), fTrigSupMod(-1),
80 fStudyExoticTrigger(0), fNExoCutInCandidate(0), fExoCutInCandidate(),
81 fMomentum(), fMomIso(),
82 fMomDaugh1(), fMomDaugh2(),
83 fTrackVector(), fProdVertex(),
84 fCluster(0), fClustersArr(0), fCaloCells(0),
85 fIsExoticTrigger(0), fClusterExoticity(1),
87 fhEIso(0), fhPtIso(0),
88 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
90 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
92 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
93 fhEIsoExoTrigger(0), fhENoIsoExoTrigger(0),
94 fhPtIsoExoTrigger(0), fhPtNoIsoExoTrigger(0),
97 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
98 fhPtInConeExoTrigger(0), fhPtClusterInConeExoTrigger(0), fhPtTrackInConeExoTrigger(0),
99 fhPtTrackInConeOtherBCPileUpSPD(0), fhPtTrackInConeVtxBC0(0),
100 fhPtTrackInConeBC0PileUpSPD(0),
101 fhPtInConePileUp(), fhPtInConeCent(0),
103 fhPerpConeSumPt(0), fhPerpConeSumPtTOFBC0(0),
104 fhPtInPerpCone(0), fhPtInPerpConeTOFBC0(0),
105 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
106 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
107 fhEtaPhiInPerpCone(0), fhEtaPhiInPerpConeTOFBC0(0),
108 fhEtaBandClusterEtaPhi(0), fhPhiBandClusterEtaPhi(0),
109 fhEtaBandTrackEtaPhi(0), fhPhiBandTrackEtaPhi(0),
110 fhEtaBandClusterPt(0), fhPhiBandClusterPt(0),
111 fhEtaBandTrackPt(0), fhPhiBandTrackPt(0),
112 fhEtaBandCell(0), fhPhiBandCell(0),
113 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
114 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
116 fhConeSumPtCellTrack(0),
117 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
118 fhConeSumPtExoTrigger(0), fhConeSumPtClusterExoTrigger(0), fhConeSumPtTrackExoTrigger(0),
120 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
121 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
122 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
123 fhConeSumPtTrigEtaPhi(0),
124 fhConeSumPtCellTrackTrigEtaPhi(0),
125 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
126 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
127 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
128 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
129 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
130 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
131 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
132 fhConeSumPtEtaUENormCluster(0), fhConeSumPtPhiUENormCluster(0),
133 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
134 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
135 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
136 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
137 fhConeSumPtEtaUENormTrack(0), fhConeSumPtPhiUENormTrack(0),
138 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
139 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
140 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
141 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
142 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
143 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
144 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
145 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
146 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
147 fhConeSumPtCellvsTrack(0),
148 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
149 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
150 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
151 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
152 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
153 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
154 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
155 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
156 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
157 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
158 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
159 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
160 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
161 fhPtPrimMCPi0DecayPairOutOfCone(0),
162 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
163 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
164 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
165 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
166 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
167 fhPtPrimMCPi0DecayPairNoOverlap(0),
168 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
169 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
170 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
171 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
172 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
173 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
174 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
175 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
176 fhPtPrimMCEtaDecayPairOutOfCone(0),
177 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
178 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
179 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
180 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
181 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
182 fhPtPrimMCEtaDecayPairNoOverlap(0),
183 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
184 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
185 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
186 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
187 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
188 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
189 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
190 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
191 fhPtLeadConeBin(0), fhSumPtConeBin(0),
192 fhSumPtConeAfterEtaBandUESubBin(0),
193 fhPtLeadConeBinMC(0),
194 fhSumPtConeAfterEtaBandUESubBinMC(0),
196 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
197 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
198 fhSumPtConeAfterEtaBandUESubBinLambda0(0), fhPtLeadConeBinLambda0MC(0),
199 fhSumPtConeAfterEtaBandUESubBinLambda0MC(0),fhSumPtConeBinLambda0MC(0),
200 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
202 fhPtTrigBinSumPtTrackCone(0), fhPtTrigBinSumPtClusterCone(0),
203 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
204 fhPtTrigBinSumPtTrackConeMC(0), fhPtTrigBinSumPtClusterConeMC(0),
205 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
206 fhPtTrigBinSumPtTrackConeDecay(0), fhPtTrigBinSumPtClusterConeDecay(0),
207 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
208 fhPtTrigBinLambda0vsSumPtTrackCone(0), fhPtTrigBinLambda0vsSumPtClusterCone(0),
209 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
210 fhPtTrigBinLambda0vsSumPtTrackConeMC(0), fhPtTrigBinLambda0vsSumPtClusterConeMC(0),
211 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap(0),fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap(0),
212 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap(0), fhPtTrigBinLambda0vsSumPtConeMC1Overlap(0),
213 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap(0),fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap(0),
218 fhELambda0LocMax1(), fhELambda1LocMax1(),
219 fhELambda0LocMax2(), fhELambda1LocMax2(),
220 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
222 fhEIsoPileUp(), fhPtIsoPileUp(),
223 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
224 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
225 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
226 fhTimeNPileUpVertContributors(0),
227 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
229 fhPtClusterInConePerRCut(0), fhPtClusterInConePerRCutLargePtTrig(0),
230 fhPtTrackInConePerRCut(0), fhPtTrackInConePerRCutLargePtTrig(0),
231 fhConeSumPtClusterPerRCut(0), fhConeSumPtClusterPerRCutLargePtTrig(0),
232 fhConeSumPtTrackPerRCut(0), fhConeSumPtTrackPerRCutLargePtTrig(0),
233 fhConeNClusterPerMinPtCut(0), fhConeNClusterPerMinPtCutLargePtTrig(0),
234 fhConeNTrackPerMinPtCut(0), fhConeNTrackPerMinPtCutLargePtTrig(0),
235 fhPerpConeNTrackPerMinPtCut(0), fhPerpConeNTrackPerMinPtCutLargePtTrig(0),
236 fhConeSumPtClusterPerMinPtCut(0), fhConeSumPtClusterPerMinPtCutLargePtTrig(0),
237 fhConeSumPtTrackPerMinPtCut(0), fhConeSumPtTrackPerMinPtCutLargePtTrig(0),
238 fhPerpConeSumPtTrackPerMinPtCut(0), fhPerpConeSumPtTrackPerMinPtCutLargePtTrig(0),
239 fhConeSumPtClusterPerMaxPtCut(0), fhConeSumPtClusterPerMaxPtCutLargePtTrig(0),
240 fhConeSumPtTrackPerMaxPtCut(0), fhConeSumPtTrackPerMaxPtCutLargePtTrig(0),
241 fhConeSumPtTrackPerEtaCut(0), fhConeSumPtTrackPerEtaCutLargePtTrig(0),
243 fhPtClusterInConePerNCellCut(0), fhPtClusterInConePerNCellCutLargePtTrig(0),
244 fhPtTrackInConePerNCellCut(0), fhPtTrackInConePerNCellCutLargePtTrig(0),
245 fhConeSumPtClusterPerNCellCut(0), fhConeSumPtClusterPerNCellCutLargePtTrig(0),
246 fhConeSumPtTrackPerNCellCut(0), fhConeSumPtTrackPerNCellCutLargePtTrig(0),
247 fhPtClusterInConePerExoCut(0), fhPtClusterInConePerExoCutLargePtTrig(0),
248 fhPtTrackInConePerExoCut(0), fhPtTrackInConePerExoCutLargePtTrig(0),
249 fhConeSumPtClusterPerExoCut(0), fhConeSumPtClusterPerExoCutLargePtTrig(0),
250 fhConeSumPtTrackPerExoCut(0), fhConeSumPtTrackPerExoCutLargePtTrig(0),
252 fhConeSumPtTrackTOFBC0(0), fhConeSumPtTrackTOFBCN(0), fhConeSumPtTrackTOFNo(0),
253 fhPtTrackInConeTOFBC0 (0), fhPtTrackInConeTOFBCN (0), fhPtTrackInConeTOFNo (0),
254 fhPhiTrackInCone(0), fhEtaTrackInCone(0), fhEtaPhiTrackInCone(0),
255 fhPhiTrackInConeTOFBC0(0), fhPhiTrackInConeTOFBCN(0), fhPhiTrackInConeTOFNo(0),
256 fhEtaTrackInConeTOFBC0(0), fhEtaTrackInConeTOFBCN(0), fhEtaTrackInConeTOFNo(0),
257 fhEtaPhiTrackInConeTOFBC0(0), fhEtaPhiTrackInConeTOFBCN(0), fhEtaPhiTrackInConeTOFNo(0),
258 fhTrackTOFInCone(0), fhTrackTOFInConeBC0(0), fhTrackTOFInConeExoTrigger(0),
260 fhConeSumPtTrackITSRefitOnSPDOn(0), fhConeSumPtTrackITSRefitOnSPDOff(0), fhConeSumPtTrackITSRefitOffSPDOff(0),
261 fhPtTrackInConeITSRefitOnSPDOn(0), fhPtTrackInConeITSRefitOnSPDOff(0) , fhPtTrackInConeITSRefitOffSPDOff(0),
262 fhPhiTrackInConeITSRefitOnSPDOn(0), fhPhiTrackInConeITSRefitOnSPDOff(0), fhPhiTrackInConeITSRefitOffSPDOff(0),
263 fhEtaTrackInConeITSRefitOnSPDOn(0), fhEtaTrackInConeITSRefitOnSPDOff(0), fhEtaTrackInConeITSRefitOffSPDOff(0),
264 fhEtaPhiTrackInConeITSRefitOnSPDOn(0),fhEtaPhiTrackInConeITSRefitOnSPDOff(0), fhEtaPhiTrackInConeITSRefitOffSPDOff(0),
265 fhConeSumPtTrackTOFBC0ITSRefitOnSPDOn(0), fhPtTrackInConeTOFBC0ITSRefitOnSPDOn(0),
266 fhPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0), fhEtaTrackInConeTOFBC0ITSRefitOnSPDOn(0),fhEtaPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0),
267 fhPerpConeSumPtITSRefitOnSPDOn (0), fhPtInPerpConeITSRefitOnSPDOn(0), fhEtaPhiInPerpConeITSRefitOnSPDOn(0),
268 fhPerpConeSumPtTOFBC0ITSRefitOnSPDOn (0), fhPtInPerpConeTOFBC0ITSRefitOnSPDOn (0), fhEtaPhiInPerpConeTOFBC0ITSRefitOnSPDOn(0)
272 for(
Int_t i = 0; i < 5 ; i++)
279 for(
Int_t j = 0; j < 5 ; j++)
314 for(
Int_t iso =0; iso < 2; iso++)
324 for(
Int_t i = 0; i < 5 ; i++)
343 for(
Int_t i = 0; i < 2 ; i++)
359 for(
Int_t imc = 0; imc < 4; imc++)
367 for(
Int_t i = 0; i < 2 ; i++)
383 for(
Int_t ieta = 0; ieta < 4; ieta++)
385 for(
Int_t iphi = 0; iphi < 3; iphi++)
396 for(
Int_t ism =0; ism < 20; ism++)
419 for(
Int_t i = 0 ; i < 7 ; i++)
428 for(
Int_t i = 0 ; i < 3 ; i++)
434 for(
Int_t i = 0 ; i < 4 ; i++)
467 Float_t ptTrig = pCandidate->Pt() ;
468 Float_t phiTrig = pCandidate->Phi();
469 Float_t etaTrig = pCandidate->Eta();
471 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
473 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
477 AliWarning(
"Cluster not available?");
482 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
483 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
546 Float_t phiTrig = pCandidate->Phi();
547 if(phiTrig<0) phiTrig += TMath::TwoPi();
549 Float_t etaTrig = pCandidate->Eta();
551 if(pCandidate->GetDetectorTag()==
kEMCAL)
553 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
556 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
558 if(!eGeom->CheckAbsCellId(absId))
return ;
561 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
563 Int_t imEta=-1, imPhi=-1;
564 Int_t ieta =-1, iphi =-1;
566 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
568 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
570 Int_t colTrig = ieta;
571 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
573 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
575 Int_t sqrSize = int(conesize/0.0143);
579 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
580 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
584 Int_t irowmin = rowTrig-sqrSize;
585 if ( irowmin < 0 ) irowmin = 0 ;
587 Int_t irowmax = rowTrig+sqrSize;
588 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
590 for(
Int_t irow = irowmin; irow <irowmax; irow++)
592 for(
Int_t icol = 0; icol < nTotalCols; icol++)
594 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
595 if(inSector==5)
continue;
598 if(icol < AliEMCALGeoParams::fgkEMCALCols)
600 inSupMod = 2*inSector + 1;
603 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
605 inSupMod = 2*inSector;
606 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
609 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
612 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
616 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
617 if(!eGeom->CheckAbsCellId(iabsId))
continue;
619 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
627 Int_t icolmin = colTrig-sqrSize;
628 if ( icolmin < 0 ) icolmin = 0;
630 Int_t icolmax = colTrig+sqrSize;
631 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
634 for(
Int_t icol = icolmin; icol < icolmax; icol++)
636 for(
Int_t irow = 0; irow < nTotalRows; irow++)
638 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
639 if ( inSector == 5 ) continue ;
645 if(icol < AliEMCALGeoParams::fgkEMCALCols)
648 inSupMod = 2*inSector + 1;
651 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
654 inSupMod = 2*inSector;
655 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
658 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
661 if(TMath::Abs(icol-colTrig) < sqrSize) {
664 if(TMath::Abs(irow-rowTrig) < sqrSize) {
668 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
669 if(!eGeom->CheckAbsCellId(iabsId))
671 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
675 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
685 Float_t ptTrig = pCandidate->Pt();
708 Float_t coneptsumPerpTrackPerMinCut[20];
709 Float_t coneNPerpTrackPerMinCut [20];
715 coneptsumPerpTrackPerMinCut[icut] = 0;
716 coneNPerpTrackPerMinCut [icut] = 0;
720 Float_t ptTrig = pCandidate->Pt() ;
721 Float_t phiTrig = pCandidate->Phi();
722 Float_t etaTrig = pCandidate->Eta();
727 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
729 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
733 AliWarning(
"Track not available?");
740 if ( pCandidate->GetDetectorTag() ==
kCTS )
743 Bool_t contained = kFALSE;
745 for(
Int_t i = 0; i < 4; i++)
747 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
750 if ( contained ) continue ;
768 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
770 phiBandPtSum+=track->Pt();
777 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
779 etaBandPtSum+=track->Pt();
787 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
788 Double_t dEta = etaTrig - track->Eta();
789 Double_t arg = dPhi*dPhi + dEta*dEta;
791 if(TMath::Sqrt(arg) < conesize)
793 sumptPerp+=track->Pt();
801 ULong_t status = track->GetStatus();
802 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
803 Int_t trackBC = track->GetTOFBunchCrossing(bz);
806 if(okTOF && trackBC == 0)
811 sumptPerpBC0+=track->Pt();
814 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
821 sumptPerpITSSPD+=track->Pt();
824 if(okTOF && trackBC == 0 && !bConstrained)
829 sumptPerpBC0ITSSPD+=track->Pt();
835 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
836 arg = dPhi*dPhi + dEta*dEta;
838 if(TMath::Sqrt(arg) < conesize)
840 sumptPerp+=track->Pt();
848 coneptsumPerpTrackPerMinCut[icut]+=track->Pt();
849 coneNPerpTrackPerMinCut [icut]++;
860 ULong_t status = track->GetStatus();
861 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
862 Int_t trackBC = track->GetTOFBunchCrossing(bz);
865 if(okTOF && trackBC == 0)
870 sumptPerpBC0+=track->Pt();
873 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
880 sumptPerpITSSPD+=track->Pt();
883 if(okTOF && trackBC == 0 && !bConstrained)
888 sumptPerpBC0ITSSPD+=track->Pt();
899 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
900 dcaCons = aodTrack->DCA();
903 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
964 Float_t etaUEptsumCluster = 0 ;
965 Float_t phiUEptsumCluster = 0 ;
974 Float_t coneA = conesize*conesize*TMath::Pi();
975 Float_t ptTrig = pCandidate->Pt() ;
976 Float_t phiTrig = pCandidate->Phi();
977 Float_t etaTrig = pCandidate->Eta();
979 Float_t m02 = pCandidate->GetM02();
980 Int_t mcTag = pCandidate->GetTag() ;
985 Float_t phiUEptsumTrackNorm = 0 ;
986 Float_t etaUEptsumTrackNorm = 0 ;
987 Float_t coneptsumTrackSubPhi = 0 ;
988 Float_t coneptsumTrackSubEta = 0 ;
989 Float_t coneptsumTrackSubPhiNorm = 0 ;
990 Float_t coneptsumTrackSubEtaNorm = 0 ;
1004 Float_t correctConeSumTrack = 1;
1005 Float_t correctConeSumTrackPhi = 1;
1008 phiUEptsumTrack,etaUEptsumTrack,
1009 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
1010 correctConeSumTrack,correctConeSumTrackPhi);
1012 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
1013 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
1021 if(coneptsumTrack > 0)
1023 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
1024 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
1045 Float_t phiUEptsumClusterNorm = 0 ;
1046 Float_t etaUEptsumClusterNorm = 0 ;
1047 Float_t coneptsumClusterSubPhi = 0 ;
1048 Float_t coneptsumClusterSubEta = 0 ;
1049 Float_t coneptsumClusterSubPhiNorm = 0 ;
1050 Float_t coneptsumClusterSubEtaNorm = 0 ;
1051 Float_t phiUEptsumCellNorm = 0 ;
1052 Float_t etaUEptsumCellNorm = 0 ;
1053 Float_t coneptsumCellSubPhi = 0 ;
1054 Float_t coneptsumCellSubEta = 0 ;
1055 Float_t coneptsumCellSubPhiNorm = 0 ;
1056 Float_t coneptsumCellSubEtaNorm = 0 ;
1074 Float_t correctConeSumClusterEta = 1;
1075 Float_t correctConeSumClusterPhi = 1;
1078 phiUEptsumCluster,etaUEptsumCluster,
1079 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
1080 correctConeSumClusterEta,correctConeSumClusterPhi);
1091 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
1092 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
1099 if(coneptsumCluster!=0)
1101 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
1102 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
1134 Float_t emcPhiSize = TMath::DegToRad()*100.;
1136 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
1137 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
1141 Float_t correctConeSumCellEta = 1;
1142 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
1144 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
1148 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
1151 Float_t correctConeSumCellPhi = 1;
1153 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
1154 (phiTrig-conesize < 80*TMath::DegToRad()))
1157 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
1158 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
1164 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
1168 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
1169 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1180 if ( coneptsumCell > 0.01 )
1182 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1183 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1193 sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1194 sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1224 Int_t ptsumAfterEtaBandUESubBin = -1;
1226 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, n bins %d",pt,coneptsumTrack+coneptsumCluster,
fNBkgBin));
1235 if ( ptsumAfterEtaBandUESubBin >= 0 )
1237 AliDebug(1,Form(
"\t Sum bin %d [%2.2f,%2.2f]" , ptsumAfterEtaBandUESubBin ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin] ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin +1]));
1246 Int_t ptsumAfterEtaBandSubBinMC = ptsumAfterEtaBandUESubBin+mcIndex*
fNBkgBin;
1248 if( ptsumAfterEtaBandUESubBin >=0 )
1258 if( ptsumAfterEtaBandUESubBin >=0 )
1277 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1279 if( ptsumAfterEtaBandUESubBin >=0 )
1298 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1299 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1323 coneptLeadCluster = 0;
1324 coneptsumCluster = 0;
1328 Float_t ptTrig = aodParticle->Pt();
1339 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1340 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1351 Float_t coneNClusterPerMinCut [20];
1352 Float_t coneptsumClusterPerMinCut[20];
1353 Float_t coneptsumClusterPerMaxCut[20];
1354 Float_t coneptsumClusterPerRCut [10];
1356 Float_t coneptsumClusterPerNCellCut[20];
1357 Float_t coneptsumClusterPerExoCut [20];
1363 coneNClusterPerMinCut [icut] = 0;
1364 coneptsumClusterPerMinCut[icut] = 0;
1365 coneptsumClusterPerMaxCut[icut] = 0;
1373 coneptsumClusterPerRCut[icut] = 0;
1380 Float_t m02 = aodParticle->GetM02();
1381 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1382 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1383 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1384 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1388 coneptsumClusterPerNCellCut[icut] = 0;
1396 coneptsumClusterPerExoCut[icut] = 0;
1400 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1402 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1429 coneptsumCluster+=ptcone;
1430 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1438 coneptsumClusterPerMinCut[icut]+=ptcone;
1439 coneNClusterPerMinCut [icut]++;
1442 if ( ptcone <
fMaxPtCutInCone[icut] ) coneptsumClusterPerMaxCut[icut]+=ptcone;
1453 coneptsumClusterPerRCut[icut]+=ptcone;
1462 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1469 coneptsumClusterPerNCellCut[icut]+=ptcone;
1482 coneptsumClusterPerExoCut[icut]+=ptcone;
1495 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1496 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1526 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1559 Float_t ptTrig = aodParticle->Pt();
1560 Float_t phiTrig = aodParticle->Phi();
1561 if(phiTrig<0) phiTrig += TMath::TwoPi();
1562 Float_t etaTrig = aodParticle->Eta();
1564 if(aodParticle->GetDetectorTag()==
kEMCAL)
1566 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1569 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1571 if(!eGeom->CheckAbsCellId(absId))
return ;
1574 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1576 Int_t imEta=-1, imPhi=-1;
1577 Int_t ieta =-1, iphi =-1;
1579 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1581 Int_t iEta=-1, iPhi=-1;
1582 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1584 Int_t colTrig = iEta;
1585 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1586 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1588 Int_t sqrSize = int(conesize/0.0143);
1593 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1595 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1597 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1598 if(inSector==5)
continue;
1600 Int_t inSupMod = -1;
1602 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1604 inSupMod = 2*inSector + 1;
1607 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1609 inSupMod = 2*inSector;
1610 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1613 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1615 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1616 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1619 coneptsumCell += cells->GetCellAmplitude(iabsId);
1635 coneptLeadTrack = 0;
1639 Float_t ptTrig = aodParticle->Pt();
1662 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1663 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1672 Float_t coneNTrackPerMinCut [20];
1673 Float_t coneptsumTrackPerMinCut[20];
1674 Float_t coneptsumTrackPerMaxCut[20];
1675 Float_t coneptsumTrackPerEtaCut[10];
1676 Float_t coneptsumTrackPerNCellCut[20];
1677 Float_t coneptsumTrackPerExoCut[20];
1678 Float_t coneptsumTrackPerRCut [10];
1679 Float_t coneptsumTrackTOFBC0 = 0;
1680 Float_t coneptsumTrackTOFBCN = 0;
1681 Float_t coneptsumTrackTOFNo = 0;
1682 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1683 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1684 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1685 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1691 coneptsumTrackPerMinCut[icut] = 0;
1692 coneptsumTrackPerMaxCut[icut] = 0;
1693 coneNTrackPerMinCut [icut] = 0;
1701 coneptsumTrackPerEtaCut[icut] = 0;
1709 coneptsumTrackPerRCut[icut] = 0;
1716 Float_t m02 = aodParticle->GetM02();
1717 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1718 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1719 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1720 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1724 coneptsumTrackPerNCellCut[icut] = 0;
1732 coneptsumTrackPerExoCut[icut] = 0;
1736 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1738 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1740 pTtrack = track->Pt();
1747 Int_t trackLabel = TMath::Abs(track->GetLabel());
1749 AliVParticle * mcpart =
GetMC()->GetTrack(trackLabel);
1750 if( !mcpart )
continue;
1752 Int_t partInConeCharge = TMath::Abs(mcpart->Charge());
1753 Int_t partInConePDG = mcpart->PdgCode();
1754 Bool_t physPrimary = mcpart->IsPhysicalPrimary();
1756 if ( partInConeCharge > 0 && TMath::Abs(partInConePDG) != 11 )
1759 if ( TMath::Abs(partInConePDG) == 211 ) mcChTag = 0;
1760 else if ( TMath::Abs(partInConePDG) == 321 ) mcChTag = 1;
1761 else if ( TMath::Abs(partInConePDG) == 2212 ) mcChTag = 2;
1777 coneptsumTrack+=pTtrack;
1778 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1786 coneptsumTrackPerMinCut[icut]+=pTtrack;
1787 coneNTrackPerMinCut [icut]++;
1790 if ( pTtrack <
fMaxPtCutInCone[icut] ) coneptsumTrackPerMaxCut[icut]+=pTtrack;
1798 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1809 coneptsumTrackPerRCut[icut]+=pTtrack;
1820 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1825 coneptsumTrackPerNCellCut[icut]+=pTtrack;
1838 coneptsumTrackPerExoCut[icut]+=pTtrack;
1849 phitrack = track->Phi();
1850 etatrack = track->Eta();
1857 ULong_t status = track->GetStatus();
1858 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1859 trackBC = track->GetTOFBunchCrossing(bz);
1860 Double32_t tof = track->GetTOFsignal()*1e-3;
1879 coneptsumTrackTOFBC0 += pTtrack;
1887 coneptsumTrackTOFBCN += pTtrack;
1896 coneptsumTrackTOFNo += pTtrack;
1900 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1901 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1909 coneptsumTrackITSRefitOnSPDOff += pTtrack;
1917 coneptsumTrackITSRefitOffSPDOff += pTtrack;
1926 coneptsumTrackITSRefitOnSPDOn += pTtrack;
1933 if(okTOF && trackBC == 0 && !bConstrained)
1939 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
1952 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
1953 dcaCons = aodTrack->DCA();
1956 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
2008 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
2009 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
2048 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
2079 if(energy < 8)
return;
2085 Int_t nVerticesSPD = -1;
2086 Int_t nVerticesTracks = -1;
2090 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
2091 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
2096 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
2097 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
2109 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
2113 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
2114 ncont=pv->GetNContributors();
2115 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
2117 diamZ = esdEv->GetDiamondZ();
2121 AliAODVertex *pv=aodEv->GetVertex(iVert);
2122 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
2123 ncont=pv->GetNContributors();
2124 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
2126 diamZ = aodEv->GetDiamondZ();
2129 Double_t distZ = TMath::Abs(z2-z1);
2130 diamZ = TMath::Abs(z2-diamZ);
2143 (AliAODPWG4ParticleCorrelation *pCandidate,
2147 if(!fFillTMHisto && !fFillSSHisto && !fFillBackgroundBinHistograms && !fFillTaggedDecayHistograms)
return;
2149 Int_t nMaxima = pCandidate->GetNLM();
2150 Int_t mcTag = pCandidate->GetTag() ;
2151 Bool_t isolated = pCandidate->IsIsolated();
2153 Float_t m02 = pCandidate->GetM02() ;
2155 Float_t pt = pCandidate->Pt();
2156 Float_t eta = pCandidate->Eta();
2157 Float_t phi = pCandidate->Phi();
2158 if(phi<0) phi+= TMath::TwoPi();
2162 if(fFillTaggedDecayHistograms)
2164 Int_t decayTag = pCandidate->DecayTag();
2165 if(decayTag < 0) decayTag = 0;
2167 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
2169 if(!GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
continue;
2171 if(fFillSSHisto) fhPtLambda0Decay[isolated][ibit]->Fill(pt, m02, GetEventWeight());
2175 if( m02 > fDecayTagsM02Cut )
continue;
2177 fhPtDecay [isolated][ibit]->Fill(pt, GetEventWeight());
2178 fhEtaPhiDecay[isolated][ibit]->Fill(eta, phi, GetEventWeight());
2182 fhPtDecayMC[isolated][ibit][mcIndex]->Fill(pt, GetEventWeight());
2185 fhPtDecayMC[isolated][ibit][kmcPhoton]->Fill(pt, GetEventWeight());
2189 if ( mcIndex == kmcPi0Decay ) fhPtDecayMC[isolated][ibit][kmcPi0DecayLostPair]->Fill(pt, GetEventWeight());
2190 else if( mcIndex == kmcEtaDecay ) fhPtDecayMC[isolated][ibit][kmcEtaDecayLostPair]->Fill(pt, GetEventWeight());
2199 if(fFillBackgroundBinHistograms)
2202 Int_t ptsumBin = -1;
2203 Int_t leadptBin = -1;
2205 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,fNBkgBin));
2207 for(
Int_t ibin = 0; ibin < fNBkgBin; ibin++)
2209 if( coneptsum >= fBkgBinLimit[ibin] && coneptsum < fBkgBinLimit[ibin+1]) ptsumBin = ibin;
2210 if( coneleadpt >= fBkgBinLimit[ibin] && coneleadpt < fBkgBinLimit[ibin+1]) leadptBin = ibin;
2215 if ( leadptBin >= 0 )
2217 AliDebug(1,Form(
"\t Lead bin %d [%2.2f,%2.2f]", leadptBin,fBkgBinLimit[leadptBin],fBkgBinLimit[leadptBin+1]));
2219 fhPtLeadConeBin[leadptBin]->Fill(pt, GetEventWeight());
2222 fhPtLeadConeBinLambda0[leadptBin]->Fill(pt, m02, GetEventWeight());
2224 if ( leadptBin == 0 )
2225 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
2226 pt, GetTrackMultiplicity(),GetEMCALClusters()->GetEntriesFast()));
2229 if ( ptsumBin >= 0 )
2231 AliDebug(1,Form(
"\t Sum bin %d [%2.2f,%2.2f]" , ptsumBin ,fBkgBinLimit[ptsumBin] ,fBkgBinLimit[ptsumBin +1]));
2233 fhSumPtConeBin[ptsumBin]->Fill(pt, GetEventWeight());
2235 if(fFillSSHisto) fhSumPtConeBinLambda0[ptsumBin]->Fill(pt, m02, GetEventWeight());
2239 if( fFillTaggedDecayHistograms && m02 < fDecayTagsM02Cut )
2241 Int_t decayTag = pCandidate->DecayTag();
2242 if(decayTag < 0) decayTag = 0;
2244 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
2246 if(GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
2248 Int_t leadptBinDecay = leadptBin+ibit*fNBkgBin;
2249 Int_t ptsumBinDecay = ptsumBin+ibit*fNBkgBin;
2250 if( leadptBin >=0 ) fhPtLeadConeBinDecay[leadptBinDecay]->Fill(pt, GetEventWeight());
2251 if( ptsumBin >=0 ) fhSumPtConeBinDecay [ ptsumBinDecay]->Fill(pt, GetEventWeight());
2258 Int_t leadptBinMC = leadptBin+mcIndex*fNBkgBin;
2259 Int_t ptsumBinMC = ptsumBin+mcIndex*fNBkgBin;
2263 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2264 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2269 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt, GetEventWeight());
2270 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2275 leadptBinMC = leadptBin+kmcPhoton*fNBkgBin;
2276 ptsumBinMC = ptsumBin+kmcPhoton*fNBkgBin;
2279 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2280 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2285 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt, GetEventWeight());
2286 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2293 if ( mcIndex == kmcPi0Decay )
2295 leadptBinMC = leadptBin+kmcPi0DecayLostPair*fNBkgBin;
2296 ptsumBinMC = ptsumBin+kmcPi0DecayLostPair*fNBkgBin;
2298 else if(mcIndex == kmcEtaDecay)
2300 leadptBinMC = leadptBin+kmcEtaDecayLostPair*fNBkgBin;
2301 ptsumBinMC = ptsumBin+kmcEtaDecayLostPair*fNBkgBin;
2304 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
2308 fhPtLeadConeBinMC[leadptBinMC]->Fill(pt, GetEventWeight());
2309 if(fFillSSHisto) fhPtLeadConeBinLambda0MC[leadptBinMC]->Fill(pt, m02, GetEventWeight());
2314 fhSumPtConeBinMC [ ptsumBinMC]->Fill(pt);
2315 if(fFillSSHisto) fhSumPtConeBinLambda0MC [ ptsumBinMC]->Fill(pt, m02, GetEventWeight());
2324 Int_t ptTrigBin = -1;
2325 if(fFillPtTrigBinHistograms)
2327 for(
Int_t ibin = 0; ibin < fNPtTrigBin; ibin++)
2329 if( pt >= fPtTrigBinLimit[ibin] && pt < fPtTrigBinLimit[ibin+1]) ptTrigBin = ibin;
2333 if ( ptTrigBin >= 0 )
2335 AliDebug(1,Form(
"Trigger pT %f, bin %d [%2.2f,%2.2f]",pt,ptTrigBin,fPtTrigBinLimit[ptTrigBin],fPtTrigBinLimit[ptTrigBin+1]));
2337 fhPtTrigBinPtLeadCone[ptTrigBin]->Fill(coneleadpt, GetEventWeight());
2338 fhPtTrigBinSumPtCone [ptTrigBin]->Fill(coneptsum , GetEventWeight());
2342 fhPtTrigBinSumPtTrackCone [ptTrigBin]->Fill(coneptsumTrack, GetEventWeight());
2343 fhPtTrigBinSumPtClusterCone[ptTrigBin]->Fill(coneptsumClust, GetEventWeight());
2348 fhPtTrigBinLambda0vsPtLeadCone[ptTrigBin]->Fill(coneleadpt, m02, GetEventWeight());
2349 fhPtTrigBinLambda0vsSumPtCone [ptTrigBin]->Fill(coneptsum , m02, GetEventWeight());
2353 fhPtTrigBinLambda0vsSumPtTrackCone [ptTrigBin]->Fill(coneptsumTrack, m02, GetEventWeight());
2354 fhPtTrigBinLambda0vsSumPtClusterCone[ptTrigBin]->Fill(coneptsumClust, m02, GetEventWeight());
2359 if( fFillTaggedDecayHistograms && m02 < fDecayTagsM02Cut )
2361 Int_t decayTag = pCandidate->DecayTag();
2362 if(decayTag < 0) decayTag = 0;
2364 for(
Int_t ibit = 0; ibit < fNDecayBits; ibit++)
2366 if(GetNeutralMesonSelection()->CheckDecayBit(decayTag,fDecayBits[ibit]))
2368 Int_t binDecay = ptTrigBin+ibit*fNPtTrigBin;
2371 fhPtTrigBinPtLeadConeDecay[binDecay]->Fill(coneleadpt, GetEventWeight());
2372 fhPtTrigBinSumPtConeDecay [binDecay]->Fill(coneptsum , GetEventWeight());
2376 fhPtTrigBinSumPtTrackConeDecay [binDecay]->Fill(coneptsumTrack, GetEventWeight());
2377 fhPtTrigBinSumPtClusterConeDecay[binDecay]->Fill(coneptsumClust, GetEventWeight());
2386 Int_t ptTrigBinMC = ptTrigBin+mcIndex*fNPtTrigBin;
2388 fhPtTrigBinPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt , GetEventWeight());
2389 fhPtTrigBinSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , GetEventWeight());
2393 fhPtTrigBinSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, GetEventWeight());
2394 fhPtTrigBinSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, GetEventWeight());
2408 fhELambda0 [isolated]->Fill(energy, m02, GetEventWeight());
2409 fhPtLambda0[isolated]->Fill(pt, m02, GetEventWeight());
2418 fhPtLambda0MC[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2422 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MC[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2423 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MC[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2426 fhPtLambda0MC[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2431 fhPtLambda0MCConv[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2435 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConv[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2436 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConv[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2439 fhPtLambda0MCConv[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2445 Int_t noverlaps = 0;
2446 if ( fFillOverlapHistograms && fCluster )
2448 const UInt_t nlabels = fCluster->GetNLabels();
2449 Int_t overpdg[nlabels];
2450 Int_t overlab[nlabels];
2451 noverlaps = GetMCAnalysisUtils()->GetNOverlaps(fCluster->GetLabels(), nlabels, mcTag, -1,
2452 GetMC(), overpdg, overlab);
2455 fhPtNOverlap[kmcPhoton][isolated]->Fill(pt, noverlaps, GetEventWeight());
2459 if ( mcIndex == kmcPi0Decay ) fhPtNOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2460 else if( mcIndex == kmcEtaDecay ) fhPtNOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2463 fhPtNOverlap[mcIndex][isolated]->Fill(pt, noverlaps, GetEventWeight());
2468 fhPtNOverlapConv[kmcPhoton][isolated]->Fill(pt, noverlaps, GetEventWeight());
2472 if ( mcIndex == kmcPi0Decay ) fhPtNOverlapConv[kmcPi0DecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2473 else if( mcIndex == kmcEtaDecay ) fhPtNOverlapConv[kmcEtaDecayLostPair][isolated]->Fill(pt, noverlaps, GetEventWeight());
2476 fhPtNOverlapConv[mcIndex][isolated]->Fill(pt, noverlaps, GetEventWeight());
2480 if ( noverlaps == 1 )
2483 fhPtLambda0MCWith1Overlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2487 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCWith1Overlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2488 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCWith1Overlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2491 fhPtLambda0MCWith1Overlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2496 fhPtLambda0MCConvWith1Overlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2500 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConvWith1Overlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2501 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConvWith1Overlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2504 fhPtLambda0MCConvWith1Overlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2507 else if (noverlaps == 0 )
2510 fhPtLambda0MCWithNoOverlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2514 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCWithNoOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2515 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCWithNoOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2518 fhPtLambda0MCWithNoOverlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2523 fhPtLambda0MCConvWithNoOverlap[kmcPhoton][isolated]->Fill(pt, m02, GetEventWeight());
2527 if ( mcIndex == kmcPi0Decay ) fhPtLambda0MCConvWithNoOverlap[kmcPi0DecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2528 else if( mcIndex == kmcEtaDecay ) fhPtLambda0MCConvWithNoOverlap[kmcEtaDecayLostPair][isolated]->Fill(pt, m02, GetEventWeight());
2531 fhPtLambda0MCConvWithNoOverlap[mcIndex][isolated]->Fill(pt, m02, GetEventWeight());
2539 if(fFillPtTrigBinHistograms && fFillSSHisto)
2541 Int_t ptTrigBinMC = ptTrigBin+mcIndex *fNPtTrigBin;
2542 Int_t ptTrigBinMCPhoton = ptTrigBin+kmcPhoton*fNPtTrigBin;
2543 Int_t ptTrigBinMCPi0Lost = ptTrigBin+kmcPi0DecayLostPair*fNPtTrigBin;
2544 Int_t ptTrigBinMCEtaLost = ptTrigBin+kmcEtaDecayLostPair*fNPtTrigBin;
2546 if ( ptTrigBin >= 0 )
2548 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMC]->Fill(coneptsum , m02, GetEventWeight());
2549 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMC]->Fill(coneleadpt, m02, GetEventWeight());
2553 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2554 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2559 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMCPhoton]->Fill(coneptsum , m02, GetEventWeight());
2560 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMCPhoton]->Fill(coneleadpt, m02, GetEventWeight());
2564 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2565 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2569 if( mcIndex == kmcPi0Decay )
2571 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMCPi0Lost]->Fill(coneptsum , m02, GetEventWeight());
2572 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMCPi0Lost]->Fill(coneleadpt, m02, GetEventWeight());
2576 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2577 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2581 if( mcIndex == kmcEtaDecay )
2583 fhPtTrigBinLambda0vsSumPtConeMC [ptTrigBinMCEtaLost]->Fill(coneptsum , m02, GetEventWeight());
2584 fhPtTrigBinLambda0vsPtLeadConeMC[ptTrigBinMCEtaLost]->Fill(coneleadpt, m02, GetEventWeight());
2588 fhPtTrigBinLambda0vsSumPtTrackConeMC [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2589 fhPtTrigBinLambda0vsSumPtClusterConeMC[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2593 if(fFillOverlapHistograms)
2595 if ( noverlaps == 0 )
2597 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMC]->Fill(coneptsum, m02, GetEventWeight());
2601 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2602 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2607 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCPhoton]->Fill(coneptsum, m02, GetEventWeight());
2611 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2612 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2616 if( mcIndex == kmcPi0Decay )
2618 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCPi0Lost]->Fill(coneptsum, m02, GetEventWeight());
2622 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2623 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2627 if( mcIndex == kmcEtaDecay )
2629 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap[ptTrigBinMCEtaLost]->Fill(coneptsum, m02, GetEventWeight());
2633 fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2634 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2640 else if ( noverlaps == 1 )
2642 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMC]->Fill(coneptsum, m02, GetEventWeight());
2645 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMC]->Fill(coneptsumTrack, m02, GetEventWeight());
2646 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMC]->Fill(coneptsumClust, m02, GetEventWeight());
2651 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCPhoton]->Fill(coneptsum, m02, GetEventWeight());
2655 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCPhoton]->Fill(coneptsumTrack, m02, GetEventWeight());
2656 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCPhoton]->Fill(coneptsumClust, m02, GetEventWeight());
2659 if( mcIndex == kmcPi0Decay )
2661 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCPi0Lost]->Fill(coneptsum, m02, GetEventWeight());
2665 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCPi0Lost]->Fill(coneptsumTrack, m02, GetEventWeight());
2666 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCPi0Lost]->Fill(coneptsumClust, m02, GetEventWeight());
2670 if( mcIndex == kmcEtaDecay )
2672 fhPtTrigBinLambda0vsSumPtConeMC1Overlap[ptTrigBinMCEtaLost]->Fill(coneptsum, m02, GetEventWeight());
2676 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap [ptTrigBinMCEtaLost]->Fill(coneptsumTrack, m02, GetEventWeight());
2677 fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap[ptTrigBinMCEtaLost]->Fill(coneptsumClust, m02, GetEventWeight());
2689 if(GetCalorimeter() ==
kEMCAL && GetFirstSMCoveredByTRD() >= 0 &&
2690 GetModuleNumber(pCandidate) >= GetFirstSMCoveredByTRD() )
2692 fhELambda0TRD [isolated]->Fill(energy, m02, GetEventWeight());
2693 fhPtLambda0TRD[isolated]->Fill(pt , m02, GetEventWeight());
2697 if(fFillNLMHistograms)
2699 fhNLocMax[isolated]->Fill(energy, nMaxima, GetEventWeight());
2702 fhELambda0LocMax1[isolated]->Fill(energy, m02, GetEventWeight());
2703 fhELambda1LocMax1[isolated]->Fill(energy, m02, GetEventWeight());
2707 fhELambda0LocMax2[isolated]->Fill(energy, m02, GetEventWeight());
2708 fhELambda1LocMax2[isolated]->Fill(energy, m02, GetEventWeight());
2712 fhELambda0LocMaxN[isolated]->Fill(energy, m02, GetEventWeight());
2713 fhELambda1LocMaxN[isolated]->Fill(energy, m02, GetEventWeight());
2719 if ( fFillTMHisto && fCluster )
2721 Float_t dZ = fCluster->GetTrackDz();
2722 Float_t dR = fCluster->GetTrackDx();
2731 if(fhTrackMatchedDEta[isolated] && TMath::Abs(dR) < 999)
2733 fhTrackMatchedDEta[isolated]->Fill(energy, dZ, GetEventWeight());
2734 fhTrackMatchedDPhi[isolated]->Fill(energy, dR, GetEventWeight());
2735 if(energy > 0.5) fhTrackMatchedDEtaDPhi[isolated]->Fill(dZ, dR, GetEventWeight());
2738 fhTrackMatchedDEtaMC[mcIndex][isolated]->Fill(energy, dZ, GetEventWeight());
2739 fhTrackMatchedDPhiMC[mcIndex][isolated]->Fill(energy, dR, GetEventWeight());
2740 if(energy > 0.5) fhTrackMatchedDEtaDPhiMC[mcIndex][isolated]->Fill(dZ, dR, GetEventWeight());
2746 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2748 AliVTrack *track = GetCaloUtils()->GetMatchedTrack(fCluster, GetReader()->GetInputEvent());
2752 Float_t dEdx = track->GetTPCsignal();
2753 fhdEdx[isolated]->Fill(fCluster->E(), dEdx, GetEventWeight());
2755 Float_t eOverp = fCluster->E()/track->P();
2756 fhEOverP[isolated]->Fill(fCluster->E(), eOverp, GetEventWeight());
2767 GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCEta) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 2.5, GetEventWeight());
2768 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCPhoton) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 0.5, GetEventWeight());
2769 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCElectron) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 1.5, GetEventWeight());
2770 else fhTrackMatchedMCParticle[isolated]->Fill(energy, 3.5, GetEventWeight());
2776 GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCEta) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 6.5, GetEventWeight());
2777 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCPhoton) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 4.5, GetEventWeight());
2778 else if ( GetMCAnalysisUtils()->CheckTagBit(mcTag,
AliMCAnalysisUtils::kMCElectron) ) fhTrackMatchedMCParticle[isolated]->Fill(energy, 5.5, GetEventWeight());
2779 else fhTrackMatchedMCParticle[isolated]->Fill(energy, 7.5, GetEventWeight());
2792 const Int_t buffersize = 255;
2793 char onePar[buffersize] ;
2795 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2799 snprintf(onePar, buffersize,
"Isolation Cand. Detector: %s;",
fIsoDetectorString.Data()) ;
2803 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2805 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2807 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2809 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2814 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2816 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2821 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2826 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2831 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2836 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2847 return new TObjString(parList) ;
2857 outputContainer->SetName(
"IsolatedParticleHistos") ;
2901 Int_t cellBins = 15;
2919 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
2922 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
2927 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
2930 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
2934 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
2937 TString sParticle =
", x^{0,#pm}";
2941 TString parTitle = Form(
"#it{R} = %2.2f%s%s" ,
GetIsolationCut()->GetConeSize(),sThreshold.Data(),sParticle.Data());
2946 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
2949 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
2950 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
2951 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
2952 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
2954 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
2955 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
2956 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
2957 "OtherDecay",
"Electron",
"Hadron"} ;
2960 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
2961 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
2962 "#pi^{0}" ,
"#eta"} ;
2964 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
2965 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
2971 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2972 nptbins,ptmin,ptmax);
2973 fhENoIso->SetYTitle(
"#it{counts}");
2974 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
2978 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2979 nptbins,ptmin,ptmax);
2981 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2987 (
"hENoIsoExoTrigger",
2988 Form(
"Number of not isolated particles vs E, %s, exot>0.97",parTitle.Data()),
2989 nptbins,ptmin,ptmax);
2995 (
"hPtNoIsoExoTrigger",
2996 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
2997 nptbins,ptmin,ptmax);
3004 Form(
"Number of not isolated leading particles #eta vs #varphi, %s",parTitle.Data()),
3005 netabins,etamin,etamax,nphibins,phimin,phimax);
3016 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
3017 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3018 nptbins,ptmin,ptmax);
3020 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3023 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
3024 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3025 nptbins,ptmin,ptmax);
3026 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
3027 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3031 Form(
"#varphi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3032 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3034 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3038 Form(
"#eta vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3039 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3041 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3049 TString isoName [] = {
"NoIso",
"Iso"};
3050 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3054 for(
Int_t iso = 0; iso < 2; iso++)
3058 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3059 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3060 nptbins,ptmin,ptmax);
3061 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
3062 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3063 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
3067 Form(
"Number of %s leading Pi0 decay particles #eta vs #varphi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3068 netabins,etamin,etamax,nphibins,phimin,phimax);
3076 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3077 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
3078 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3089 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
3090 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
3091 nptbins,ptmin,ptmax);
3092 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
3093 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3094 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
3103 TString isoName [] = {
"NoIso",
"Iso"};
3104 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3107 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
3108 nptbins,ptmin,ptmax);
3109 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
3110 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
3111 outputContainer->Add(
fhEIso) ;
3114 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3115 nptbins,ptmin,ptmax);
3116 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
3117 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3118 outputContainer->Add(
fhPtIso) ;
3124 Form(
"Number of isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3125 nptbins,ptmin,ptmax);
3131 (
"hPtIsoExoTrigger",
3132 Form(
"Number of isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3133 nptbins,ptmin,ptmax);
3140 Form(
"Number of isolated particles vs #varphi, %s",parTitle.Data()),
3141 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3142 fhPhiIso->SetYTitle(
"#varphi (rad)");
3143 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3147 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
3148 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3150 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3154 Form(
"Number of isolated particles #eta vs #varphi, %s",parTitle.Data()),
3155 netabins,etamin,etamax,nphibins,phimin,phimax);
3163 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3164 nptbins,ptmin,ptmax, 100,0,100);
3170 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3171 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
3180 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3181 nptbins,ptmin,ptmax,10,0,10);
3186 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3187 nptbins,ptmin,ptmax,10,0,10);
3194 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3195 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3196 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
3197 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3202 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3203 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3204 fhConeSumPt->SetYTitle(
"#Sigma #it{p}_{T} (GeV/#it{c})");
3205 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3211 (
"hConePtSumExoTrigger",
3212 Form(
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3213 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3228 Form(
"Trigger #eta vs #varphi, #Sigma #it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3229 netabins,etamin,etamax,nphibins,phimin,phimax);
3236 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s",parTitleR.Data()),
3237 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3238 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
3239 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3245 Form(
"#it{p}_{T} in isolation cone for %s, exotic trigger",parTitleR.Data()),
3246 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3288 (Form(
"hPtLeadCone_Bin%d",ibin),
3289 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3296 (Form(
"hSumPtCone_Bin%d",ibin),
3297 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3306 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d",ibin),
3307 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3321 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3322 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3329 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3330 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3344 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3345 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3352 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3353 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3362 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3363 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3375 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
3376 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3383 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
3384 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3393 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d",ibin),
3394 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3407 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3408 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3415 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3416 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3425 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3426 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3505 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
3506 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3513 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
3514 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3523 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d",ibin),
3524 Form(
"#Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3531 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d",ibin),
3532 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3546 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3547 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3554 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3555 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3564 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3565 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3572 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3573 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3588 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3589 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3596 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3597 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3606 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3607 Form(
"#Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3614 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3615 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3627 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
3628 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3635 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
3636 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3645 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d",ibin),
3646 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3653 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d",ibin),
3654 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3667 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3668 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3670 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3676 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3677 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3679 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3687 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3688 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",
3690 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3696 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3697 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",
3699 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3708 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3709 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",
3711 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3717 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3718 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",
3720 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3728 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3729 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",
3731 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3737 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3738 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",
3740 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3746 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3747 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",
3749 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3755 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3756 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",
3758 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3774 Form(
"#it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3775 100,0,100,nptinconebins,ptinconemin,ptinconemax);
3785 (
"hConePtSumCluster",
3786 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3787 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3795 (
"hConePtSumClusterExoTrigger",
3796 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3797 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3806 (
"hConeNClusterPerMinPtCut",
3807 Form(
"N clusters, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3816 (
"hConeNClusterPerMinPtCutLargePtTrig",
3817 Form(
"N cluster, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3826 (
"hConePtSumClusterPerMinPtCut",
3827 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3836 (
"hConePtSumClusterPerMinPtCutLargePtTrig",
3837 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),
3847 (
"hConePtSumClusterPerMaxPtCut",
3848 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3857 (
"hConePtSumClusterPerMaxPtCutLargePtTrig",
3858 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),
3870 (
"hConePtSumClusterPerRCut",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts",
3879 (
"hConePtSumClusterPerRCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3888 (
"hPtClusterInConePerRCut",
"Cluster #it{p}_{T}, different #it{R} cuts",
3897 (
"hPtClusterInConePerRCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3909 (
"hConePtSumClusterPerNCellCut",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
3918 (
"hConePtSumClusterPerNCellCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3927 (
"hPtClusterInConePerNCellCut",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts",
3936 (
"hPtClusterInConePerNCellCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3944 for(
Int_t ishsh = 0; ishsh < 4; ishsh++)
3947 (Form(
"hPtClusterInConePerNCellPerSM_ShSh%d",ishsh),
3948 Form(
"Cluster #it{p}_{T} vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
3956 (Form(
"hConeSumPtClusterPerNCellPerSM_ShSh%d",ishsh),
3957 Form(
"Cluster #Sigma #it{p}_{T} in cone vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
3969 (
"hConePtSumClusterPerExoCut",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts",
3978 (
"hConePtSumClusterPerExoCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
3987 (
"hPtClusterInConePerExoCut",
"Cluster #it{p}_{T}, different exoticity cuts",
3996 (
"hPtClusterInConePerExoCutLargePtTrig",
"Cluster #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4006 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
4007 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4016 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4017 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4026 "#Sigma cluster #it{p}_{T} in UE Eta Band",
4027 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4033 "#Sigma cluster #it{p}_{T} UE Phi Band",
4034 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4042 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} in UE Eta Band",
4043 netabins,etamin,etamax,nphibins,phimin,phimax);
4050 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} UE Phi Band",
4051 netabins,etamin,etamax,nphibins,phimin,phimax);
4058 Form(
"#eta vs #varphi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4059 netabins,-1,1,nphibins,0,TMath::TwoPi());
4065 Form(
"#eta vs #varphi of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4066 netabins,-1,1,nphibins,0,TMath::TwoPi());
4072 Form(
"#eta vs #varphi of clusters in cone for #it{R} = %2.2f",r),
4073 netabins,-1,1,nphibins,0,TMath::TwoPi());
4079 Form(
"#eta vs #varphi of all clusters"),
4080 netabins,-1,1,nphibins,0,TMath::TwoPi());
4089 "#Sigma cell #it{p}_{T} in UE Eta Band",
4090 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4096 "#Sigma cell #it{p}_{T} UE Phi Band",
4097 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4103 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} in UE Eta Band",
4104 netabins,etamin,etamax,nphibins,phimin,phimax);
4111 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} UE Phi Band",
4112 netabins,etamin,etamax,nphibins,phimin,phimax);
4120 Form(
"#it{p}_{T} of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4121 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4127 Form(
"#it{p}_{T} of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4128 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4135 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
4136 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4144 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4145 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4154 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
4155 nptbins,ptmin,ptmax,1000,0,50);
4161 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
4168 Form(
"#col vs #row of cells in #varphi band isolation cone for #it{R} = %2.2f",r),
4178 Form(
"Clusters #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
4179 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4185 Form(
"Clusters #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
4186 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4192 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4193 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4199 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4200 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4208 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),
4209 netabins,etamin,etamax,nphibins,phimin,phimax);
4216 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),
4217 netabins,etamin,etamax,nphibins,phimin,phimax);
4227 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4228 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4234 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4235 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4241 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),
4242 netabins,etamin,etamax,nphibins,phimin,phimax);
4249 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),
4250 netabins,etamin,etamax,nphibins,phimin,phimax);
4260 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
4261 nptbins,ptmin,ptmax,100,0,1);
4267 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#varphi ",r),
4268 netabins,etamin,etamax,nphibins,phimin,phimax);
4275 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance",r),
4276 nptbins,ptmin,ptmax,100,0,1);
4282 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance, in trigger #eta-#varphi ",r),
4283 netabins,etamin,etamax,nphibins,phimin,phimax);
4291 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),
4292 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4298 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),
4299 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4305 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),
4306 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4312 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),
4313 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4319 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #eta band for cluster (before normalization), R=%2.2f",r),
4320 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
4326 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #varphi band for cluster (before normalization), R=%2.2f",r),
4327 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
4336 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
4337 nptbins,ptmin,ptmax,100,0,1);
4343 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#varphi ",r),
4344 netabins,etamin,etamax,nphibins,phimin,phimax);
4351 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance",r),
4352 nptbins,ptmin,ptmax,100,0,1);
4358 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance, in trigger #eta-#varphi ",r),
4359 netabins,etamin,etamax,nphibins,phimin,phimax);
4367 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),
4368 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4374 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),
4375 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4381 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),
4382 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4388 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),
4389 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4401 (
"hConeLeadPtTrack",
4402 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
4403 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4410 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
4411 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4418 TString mcChPartName[] = {
"Pion",
"Kaon",
"Proton",
"Other"};
4419 for(
Int_t imc = 0; imc < 4; imc++)
4422 (Form(
"hPtTrackInCone_Primary_%s",mcChPartName[imc].
Data()),
4423 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4424 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4430 (Form(
"hPtTrackInCone_Secondary_%s",mcChPartName[imc].
Data()),
4431 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, secondary MC %s",r,mcChPartName[imc].
Data()),
4432 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4438 (Form(
"hPtTrackInCone_Gener_Primary_%s",mcChPartName[imc].
Data()),
4439 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4440 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4446 (Form(
"hPtTrackInCone_Gener_Secondary_%s",mcChPartName[imc].
Data()),
4447 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, secondary MC %s",r,mcChPartName[imc].
Data()),
4448 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4459 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4460 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4469 (
"hConeNTrackPerMinPtCut",
4470 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4479 (
"hConeNTrackPerMinPtCutLargePtTrig",
4480 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4489 (
"hConePtSumTrackPerMinPtCut",
4490 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4499 (
"hConePtSumTrackPerMinPtCutLargePtTrig",
4500 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),
4511 (
"hPerpConeNTrackPerMinPtCut",
4512 Form(
"N tracks, different #it{p}_{T} cuts in perpendicular cone for #it{R} = %2.2f",r),
4521 (
"hPerpConeNTrackPerMinPtCutLargePtTrig",
4522 Form(
"N tracks, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4531 (
"hPerpConePtSumTrackPerMinPtCut",
4532 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in perpendicular cone for #it{R} = %2.2f",r),
4541 (
"hPerpConePtSumTrackPerMinPtCutLargePtTrig",
4542 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),
4552 (
"hConePtSumTrackPerMaxPtCut",
4553 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4562 (
"hConePtSumTrackPerMaxPtCutLargePtTrig",
4563 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),
4575 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f",r),
4584 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4596 (
"hConePtSumTrackPerRCut",
"Track #Sigma #it{p}_{T}, different #it{R} cuts",
4605 (
"hConePtSumTrackPerRCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4614 (
"hPtTrackInConePerRCut",
"Track #it{p}_{T}, different #it{R} cuts",
4623 (
"hPtTrackInConePerRCutLargePtTrig",
"Track #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4635 (
"hConePtSumTrackPerNCellCut",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
4644 (
"hConePtSumTrackPerNCellCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4653 (
"hPtTrackInConePerNCellCut",
"Track #it{p}_{T}, different #it{N}_{cell} cuts",
4662 (
"hPtTrackInConePerNCellCutLargePtTrig",
"Track #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4670 for(
Int_t ishsh = 0; ishsh < 4; ishsh++)
4673 (Form(
"hPtTrackInConePerNCellPerSM_ShSh%d",ishsh),
4674 Form(
"Track #it{p}_{T} vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
4682 (Form(
"hConeSumPtTrackPerNCellPerSM_ShSh%d",ishsh),
4683 Form(
"Track #Sigma #it{p}_{T} in cone vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
4695 (
"hConePtSumTrackPerExoCut",
"Track #Sigma #it{p}_{T}, different exoticity cuts",
4704 (
"hConePtSumTrackPerExoCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4713 (
"hPtTrackInConePerExoCut",
"Track #it{p}_{T}, different exoticity cuts",
4722 (
"hPtTrackInConePerExoCutLargePtTrig",
"Track #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4734 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4735 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4743 (
"hConePtSumTrackExoTrigger",
4744 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
4745 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4753 Int_t ntofbins = 1000;
4754 Int_t mintof = -500;
4758 (
"hTrackTOFInCone",
"TOF signal vs track #it{p}_{T}",
4759 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4767 (
"hTrackTOFInConeExoTrigger",
"TOF signal vs track #it{p}_{T}, exoticity > 0.97",
4768 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4775 (
"hTrackTOFInConeBC0",
"TOF signal vs track #it{p}_{T}, BC=0",
4776 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4782 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
4783 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4789 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f",r),
4790 netabins,-1,1,nphibins,0,TMath::TwoPi());
4796 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4797 nptbins,ptmin,ptmax,netabins,-1,1);
4803 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4804 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4815 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0",r),
4816 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4822 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC!=0",r),
4823 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4829 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no TOF",r),
4830 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);