17 #include <TClonesArray.h> 19 #include <TObjString.h> 24 #include "TDatabasePDG.h" 29 #include "AliMCEvent.h" 34 #include "AliVParticle.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 fFillPerSMHistograms(0), fFillPerTCardIndexHistograms(0), fTCardIndex(-1),
59 fFillEMCALRegionHistograms(0), fFillUEBandSubtractHistograms(1),
60 fFillCellHistograms(0),
61 fFillOverlapHistograms(0),
62 fStudyTracksInCone(0), fStudyMCConversionRadius(0),
63 fFillTaggedDecayHistograms(0), fNDecayBits(0),
64 fDecayBits(), fDecayTagsM02Cut(0),
65 fFillNLMHistograms(0),
66 fLeadingOnly(0), fCheckLeadingWithNeutralClusters(0),
67 fSelectPrimariesInCone(0), fMakePrimaryPi0DecayStudy(0),
68 fFillBackgroundBinHistograms(0), fNBkgBin(0),
69 fFillPtTrigBinHistograms(0), fNPtTrigBin(0),
70 fMinCellsAngleOverlap(0),
72 fNCones(0), fNPtThresFrac(0),
73 fConeSizes(), fPtThresholds(),
74 fPtFractions(), fSumPtThresholds(),
75 fStudyPtCutInCone(0), fNPtCutsInCone(0),
76 fMinPtCutInCone(), fMaxPtCutInCone(),
77 fStudyEtaCutInCone(0), fNEtaCutsInCone(0), fEtaCutInCone(),
78 fStudyRCutInCone(0), fNRCutsInCone(0), fRCutInCone(),
79 fStudyNCellsCut(0), fNNCellsInCandidate(0), fNCellsInCandidate(),
80 fNCellsWithWeight(0), fTrigSupMod(-1),
81 fStudyExoticTrigger(0), fNExoCutInCandidate(0), fExoCutInCandidate(),
82 fMomentum(), fMomIso(),
83 fMomDaugh1(), fMomDaugh2(),
84 fTrackVector(), fProdVertex(),
85 fCluster(0), fClustersArr(0), fCaloCells(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 fhSumPtConeAfterEtaBandUESubBin(0),
194 fhPtLeadConeBinMC(0),
195 fhSumPtConeAfterEtaBandUESubBinMC(0),
197 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
198 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
199 fhSumPtConeAfterEtaBandUESubBinLambda0(0), fhPtLeadConeBinLambda0MC(0),
200 fhSumPtConeAfterEtaBandUESubBinLambda0MC(0),fhSumPtConeBinLambda0MC(0),
201 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
203 fhPtTrigBinSumPtTrackCone(0), fhPtTrigBinSumPtClusterCone(0),
204 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
205 fhPtTrigBinSumPtTrackConeMC(0), fhPtTrigBinSumPtClusterConeMC(0),
206 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
207 fhPtTrigBinSumPtTrackConeDecay(0), fhPtTrigBinSumPtClusterConeDecay(0),
208 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
209 fhPtTrigBinLambda0vsSumPtTrackCone(0), fhPtTrigBinLambda0vsSumPtClusterCone(0),
210 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
211 fhPtTrigBinLambda0vsSumPtTrackConeMC(0), fhPtTrigBinLambda0vsSumPtClusterConeMC(0),
212 fhPtTrigBinLambda0vsSumPtConeMCNoOverlap(0),fhPtTrigBinLambda0vsSumPtTrackConeMCNoOverlap(0),
213 fhPtTrigBinLambda0vsSumPtClusterConeMCNoOverlap(0), fhPtTrigBinLambda0vsSumPtConeMC1Overlap(0),
214 fhPtTrigBinLambda0vsSumPtTrackConeMC1Overlap(0),fhPtTrigBinLambda0vsSumPtClusterConeMC1Overlap(0),
219 fhELambda0LocMax1(), fhELambda1LocMax1(),
220 fhELambda0LocMax2(), fhELambda1LocMax2(),
221 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
223 fhEIsoPileUp(), fhPtIsoPileUp(),
224 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
225 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
226 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
227 fhTimeNPileUpVertContributors(0),
228 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
230 fhPtClusterInConePerRCut(0), fhPtClusterInConePerRCutLargePtTrig(0),
231 fhPtTrackInConePerRCut(0), fhPtTrackInConePerRCutLargePtTrig(0),
232 fhConeSumPtClusterPerRCut(0), fhConeSumPtClusterPerRCutLargePtTrig(0),
233 fhConeSumPtTrackPerRCut(0), fhConeSumPtTrackPerRCutLargePtTrig(0),
234 fhConeNClusterPerMinPtCut(0), fhConeNClusterPerMinPtCutLargePtTrig(0),
235 fhConeNTrackPerMinPtCut(0), fhConeNTrackPerMinPtCutLargePtTrig(0),
236 fhPerpConeNTrackPerMinPtCut(0), fhPerpConeNTrackPerMinPtCutLargePtTrig(0),
237 fhConeSumPtClusterPerMinPtCut(0), fhConeSumPtClusterPerMinPtCutLargePtTrig(0),
238 fhConeSumPtTrackPerMinPtCut(0), fhConeSumPtTrackPerMinPtCutLargePtTrig(0),
239 fhPerpConeSumPtTrackPerMinPtCut(0), fhPerpConeSumPtTrackPerMinPtCutLargePtTrig(0),
240 fhConeSumPtClusterPerMaxPtCut(0), fhConeSumPtClusterPerMaxPtCutLargePtTrig(0),
241 fhConeSumPtTrackPerMaxPtCut(0), fhConeSumPtTrackPerMaxPtCutLargePtTrig(0),
242 fhConeSumPtTrackPerEtaCut(0), fhConeSumPtTrackPerEtaCutLargePtTrig(0),
244 fhPtClusterInConePerNCellCut(0), fhPtClusterInConePerNCellCutLargePtTrig(0),
245 fhPtTrackInConePerNCellCut(0), fhPtTrackInConePerNCellCutLargePtTrig(0),
246 fhConeSumPtClusterPerNCellCut(0), fhConeSumPtClusterPerNCellCutLargePtTrig(0),
247 fhConeSumPtTrackPerNCellCut(0), fhConeSumPtTrackPerNCellCutLargePtTrig(0),
248 fhPtClusterInConePerExoCut(0), fhPtClusterInConePerExoCutLargePtTrig(0),
249 fhPtTrackInConePerExoCut(0), fhPtTrackInConePerExoCutLargePtTrig(0),
250 fhConeSumPtClusterPerExoCut(0), fhConeSumPtClusterPerExoCutLargePtTrig(0),
251 fhConeSumPtTrackPerExoCut(0), fhConeSumPtTrackPerExoCutLargePtTrig(0),
253 fhConeSumPtTrackTOFBC0(0), fhConeSumPtTrackTOFBCN(0), fhConeSumPtTrackTOFNo(0),
254 fhPtTrackInConeTOFBC0 (0), fhPtTrackInConeTOFBCN (0), fhPtTrackInConeTOFNo (0),
255 fhPhiTrackInCone(0), fhEtaTrackInCone(0), fhEtaPhiTrackInCone(0),
256 fhPhiTrackInConeTOFBC0(0), fhPhiTrackInConeTOFBCN(0), fhPhiTrackInConeTOFNo(0),
257 fhEtaTrackInConeTOFBC0(0), fhEtaTrackInConeTOFBCN(0), fhEtaTrackInConeTOFNo(0),
258 fhEtaPhiTrackInConeTOFBC0(0), fhEtaPhiTrackInConeTOFBCN(0), fhEtaPhiTrackInConeTOFNo(0),
259 fhTrackTOFInCone(0), fhTrackTOFInConeBC0(0), fhTrackTOFInConeExoTrigger(0),
261 fhConeSumPtTrackITSRefitOnSPDOn(0), fhConeSumPtTrackITSRefitOnSPDOff(0), fhConeSumPtTrackITSRefitOffSPDOff(0),
262 fhPtTrackInConeITSRefitOnSPDOn(0), fhPtTrackInConeITSRefitOnSPDOff(0) , fhPtTrackInConeITSRefitOffSPDOff(0),
263 fhPhiTrackInConeITSRefitOnSPDOn(0), fhPhiTrackInConeITSRefitOnSPDOff(0), fhPhiTrackInConeITSRefitOffSPDOff(0),
264 fhEtaTrackInConeITSRefitOnSPDOn(0), fhEtaTrackInConeITSRefitOnSPDOff(0), fhEtaTrackInConeITSRefitOffSPDOff(0),
265 fhEtaPhiTrackInConeITSRefitOnSPDOn(0),fhEtaPhiTrackInConeITSRefitOnSPDOff(0), fhEtaPhiTrackInConeITSRefitOffSPDOff(0),
266 fhConeSumPtTrackTOFBC0ITSRefitOnSPDOn(0), fhPtTrackInConeTOFBC0ITSRefitOnSPDOn(0),
267 fhPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0), fhEtaTrackInConeTOFBC0ITSRefitOnSPDOn(0),fhEtaPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0),
268 fhPerpConeSumPtITSRefitOnSPDOn (0), fhPtInPerpConeITSRefitOnSPDOn(0), fhEtaPhiInPerpConeITSRefitOnSPDOn(0),
269 fhPerpConeSumPtTOFBC0ITSRefitOnSPDOn (0), fhPtInPerpConeTOFBC0ITSRefitOnSPDOn (0), fhEtaPhiInPerpConeTOFBC0ITSRefitOnSPDOn(0)
273 for(
Int_t i = 0; i < 5 ; i++)
280 for(
Int_t j = 0; j < 5 ; j++)
315 for(
Int_t iso =0; iso < 2; iso++)
325 for(
Int_t i = 0; i < 5 ; i++)
344 for(
Int_t i = 0; i < 2 ; i++)
360 for(
Int_t imc = 0; imc < 4; imc++)
368 for(
Int_t i = 0; i < 2 ; i++)
384 for(
Int_t ieta = 0; ieta < 4; ieta++)
386 for(
Int_t iphi = 0; iphi < 3; iphi++)
397 for(
Int_t ism =0; ism < 20; ism++)
420 for(
Int_t i = 0 ; i < 7 ; i++)
429 for(
Int_t i = 0 ; i < 3 ; i++)
435 for(
Int_t i = 0 ; i < 4 ; i++)
445 for(
Int_t ism =0; ism < 20; ism++)
461 for(
Int_t ism =0; ism < 20; ism++)
504 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
506 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
510 AliWarning(
"Cluster not available?");
516 cluster->GetID() == pCandidate->
GetCaloLabel(1) )
continue ;
580 if(phiTrig<0) phiTrig += TMath::TwoPi();
586 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
589 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
591 if(!eGeom->CheckAbsCellId(absId))
return ;
594 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
596 Int_t imEta=-1, imPhi=-1;
597 Int_t ieta =-1, iphi =-1;
599 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
601 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
603 Int_t colTrig = ieta;
604 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
606 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
608 Int_t sqrSize = int(conesize/0.0143);
612 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
613 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
617 Int_t irowmin = rowTrig-sqrSize;
618 if ( irowmin < 0 ) irowmin = 0 ;
620 Int_t irowmax = rowTrig+sqrSize;
621 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
623 for(
Int_t irow = irowmin; irow <irowmax; irow++)
625 for(
Int_t icol = 0; icol < nTotalCols; icol++)
627 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
628 if(inSector==5)
continue;
631 if(icol < AliEMCALGeoParams::fgkEMCALCols)
633 inSupMod = 2*inSector + 1;
636 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
638 inSupMod = 2*inSector;
639 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
642 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
645 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
649 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
650 if(!eGeom->CheckAbsCellId(iabsId))
continue;
652 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
660 Int_t icolmin = colTrig-sqrSize;
661 if ( icolmin < 0 ) icolmin = 0;
663 Int_t icolmax = colTrig+sqrSize;
664 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
667 for(
Int_t icol = icolmin; icol < icolmax; icol++)
669 for(
Int_t irow = 0; irow < nTotalRows; irow++)
671 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
672 if ( inSector == 5 ) continue ;
678 if(icol < AliEMCALGeoParams::fgkEMCALCols)
681 inSupMod = 2*inSector + 1;
684 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
687 inSupMod = 2*inSector;
688 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
691 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
694 if(TMath::Abs(icol-colTrig) < sqrSize) {
697 if(TMath::Abs(irow-rowTrig) < sqrSize) {
701 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
702 if(!eGeom->CheckAbsCellId(iabsId))
704 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
708 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
741 Float_t coneptsumPerpTrackPerMinCut[20];
742 Float_t coneNPerpTrackPerMinCut [20];
748 coneptsumPerpTrackPerMinCut[icut] = 0;
749 coneNPerpTrackPerMinCut [icut] = 0;
760 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
762 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
766 AliWarning(
"Track not available?");
776 Bool_t contained = kFALSE;
778 for(
Int_t i = 0; i < 4; i++)
780 if( trackID == pCandidate->
GetTrackLabel(i) ) contained = kTRUE;
783 if ( contained ) continue ;
801 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
803 phiBandPtSum+=track->Pt();
810 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
812 etaBandPtSum+=track->Pt();
820 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
821 Double_t dEta = etaTrig - track->Eta();
822 Double_t arg = dPhi*dPhi + dEta*dEta;
824 if(TMath::Sqrt(arg) < conesize)
826 sumptPerp+=track->Pt();
834 ULong_t status = track->GetStatus();
835 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
836 Int_t trackBC = track->GetTOFBunchCrossing(bz);
839 if(okTOF && trackBC == 0)
844 sumptPerpBC0+=track->Pt();
847 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
854 sumptPerpITSSPD+=track->Pt();
857 if(okTOF && trackBC == 0 && !bConstrained)
862 sumptPerpBC0ITSSPD+=track->Pt();
868 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
869 arg = dPhi*dPhi + dEta*dEta;
871 if(TMath::Sqrt(arg) < conesize)
873 sumptPerp+=track->Pt();
881 coneptsumPerpTrackPerMinCut[icut]+=track->Pt();
882 coneNPerpTrackPerMinCut [icut]++;
893 ULong_t status = track->GetStatus();
894 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
895 Int_t trackBC = track->GetTOFBunchCrossing(bz);
898 if(okTOF && trackBC == 0)
903 sumptPerpBC0+=track->Pt();
906 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
913 sumptPerpITSSPD+=track->Pt();
916 if(okTOF && trackBC == 0 && !bConstrained)
921 sumptPerpBC0ITSSPD+=track->Pt();
932 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
933 dcaCons = aodTrack->DCA();
936 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
997 Float_t etaUEptsumCluster = 0 ;
998 Float_t phiUEptsumCluster = 0 ;
1007 Float_t coneA = conesize*conesize*TMath::Pi();
1018 Float_t phiUEptsumTrackNorm = 0 ;
1019 Float_t etaUEptsumTrackNorm = 0 ;
1020 Float_t coneptsumTrackSubPhi = 0 ;
1021 Float_t coneptsumTrackSubEta = 0 ;
1022 Float_t coneptsumTrackSubPhiNorm = 0 ;
1023 Float_t coneptsumTrackSubEtaNorm = 0 ;
1037 Float_t correctConeSumTrack = 1;
1038 Float_t correctConeSumTrackPhi = 1;
1041 phiUEptsumTrack,etaUEptsumTrack,
1042 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
1043 correctConeSumTrack,correctConeSumTrackPhi);
1045 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
1046 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
1054 if(coneptsumTrack > 0)
1056 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
1057 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
1078 Float_t phiUEptsumClusterNorm = 0 ;
1079 Float_t etaUEptsumClusterNorm = 0 ;
1080 Float_t coneptsumClusterSubPhi = 0 ;
1081 Float_t coneptsumClusterSubEta = 0 ;
1082 Float_t coneptsumClusterSubPhiNorm = 0 ;
1083 Float_t coneptsumClusterSubEtaNorm = 0 ;
1084 Float_t phiUEptsumCellNorm = 0 ;
1085 Float_t etaUEptsumCellNorm = 0 ;
1086 Float_t coneptsumCellSubPhi = 0 ;
1087 Float_t coneptsumCellSubEta = 0 ;
1088 Float_t coneptsumCellSubPhiNorm = 0 ;
1089 Float_t coneptsumCellSubEtaNorm = 0 ;
1107 Float_t correctConeSumClusterEta = 1;
1108 Float_t correctConeSumClusterPhi = 1;
1111 phiUEptsumCluster,etaUEptsumCluster,
1112 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
1113 correctConeSumClusterEta,correctConeSumClusterPhi);
1124 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
1125 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
1132 if(coneptsumCluster!=0)
1134 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
1135 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
1167 Float_t emcPhiSize = TMath::DegToRad()*100.;
1169 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
1170 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
1174 Float_t correctConeSumCellEta = 1;
1175 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
1177 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
1181 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
1184 Float_t correctConeSumCellPhi = 1;
1186 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
1187 (phiTrig-conesize < 80*TMath::DegToRad()))
1190 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
1191 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
1197 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
1201 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
1202 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1213 if ( coneptsumCell > 0.01 )
1215 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1216 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1226 sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1227 sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1257 Int_t ptsumAfterEtaBandUESubBin = -1;
1259 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, n bins %d",pt,coneptsumTrack+coneptsumCluster,
fNBkgBin));
1268 if ( ptsumAfterEtaBandUESubBin >= 0 )
1270 AliDebug(1,Form(
"\t Sum bin %d [%2.2f,%2.2f]" , ptsumAfterEtaBandUESubBin ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin] ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin +1]));
1279 Int_t ptsumAfterEtaBandSubBinMC = ptsumAfterEtaBandUESubBin+mcIndex*
fNBkgBin;
1281 if( ptsumAfterEtaBandUESubBin >=0 )
1291 if( ptsumAfterEtaBandUESubBin >=0 )
1310 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1312 if( ptsumAfterEtaBandUESubBin >=0 )
1331 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1332 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1356 coneptLeadCluster = 0;
1357 coneptsumCluster = 0;
1378 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1379 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1390 Float_t coneNClusterPerMinCut [20];
1391 Float_t coneptsumClusterPerMinCut[20];
1392 Float_t coneptsumClusterPerMaxCut[20];
1393 Float_t coneptsumClusterPerRCut [10];
1395 Float_t coneptsumClusterPerNCellCut[20];
1396 Float_t coneptsumClusterPerExoCut [20];
1402 coneNClusterPerMinCut [icut] = 0;
1403 coneptsumClusterPerMinCut[icut] = 0;
1404 coneptsumClusterPerMaxCut[icut] = 0;
1412 coneptsumClusterPerRCut[icut] = 0;
1420 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1421 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1422 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1423 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1427 coneptsumClusterPerNCellCut[icut] = 0;
1435 coneptsumClusterPerExoCut[icut] = 0;
1439 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1441 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1480 coneptsumCluster+=ptcone;
1481 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1489 coneptsumClusterPerMinCut[icut]+=ptcone;
1490 coneNClusterPerMinCut [icut]++;
1493 if ( ptcone <
fMaxPtCutInCone[icut] ) coneptsumClusterPerMaxCut[icut]+=ptcone;
1504 coneptsumClusterPerRCut[icut]+=ptcone;
1513 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1520 coneptsumClusterPerNCellCut[icut]+=ptcone;
1533 coneptsumClusterPerExoCut[icut]+=ptcone;
1584 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1619 if(phiTrig<0) phiTrig += TMath::TwoPi();
1624 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1627 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1629 if(!eGeom->CheckAbsCellId(absId))
return ;
1632 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1634 Int_t imEta=-1, imPhi=-1;
1635 Int_t ieta =-1, iphi =-1;
1637 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1639 Int_t iEta=-1, iPhi=-1;
1640 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1642 Int_t colTrig = iEta;
1643 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1644 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1646 Int_t sqrSize = int(conesize/0.0143);
1651 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1653 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1655 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1656 if(inSector==5)
continue;
1658 Int_t inSupMod = -1;
1660 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1662 inSupMod = 2*inSector + 1;
1665 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1667 inSupMod = 2*inSector;
1668 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1671 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1673 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1674 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1677 coneptsumCell += cells->GetCellAmplitude(iabsId);
1693 coneptLeadTrack = 0;
1727 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1728 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1737 Float_t coneNTrackPerMinCut [20];
1738 Float_t coneptsumTrackPerMinCut[20];
1739 Float_t coneptsumTrackPerMaxCut[20];
1740 Float_t coneptsumTrackPerEtaCut[10];
1741 Float_t coneptsumTrackPerNCellCut[20];
1742 Float_t coneptsumTrackPerExoCut[20];
1743 Float_t coneptsumTrackPerRCut [10];
1744 Float_t coneptsumTrackTOFBC0 = 0;
1745 Float_t coneptsumTrackTOFBCN = 0;
1746 Float_t coneptsumTrackTOFNo = 0;
1747 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1748 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1749 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1750 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1756 coneptsumTrackPerMinCut[icut] = 0;
1757 coneptsumTrackPerMaxCut[icut] = 0;
1758 coneNTrackPerMinCut [icut] = 0;
1766 coneptsumTrackPerEtaCut[icut] = 0;
1774 coneptsumTrackPerRCut[icut] = 0;
1782 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1783 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1784 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1785 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1789 coneptsumTrackPerNCellCut[icut] = 0;
1797 coneptsumTrackPerExoCut[icut] = 0;
1801 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1803 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1805 pTtrack = track->Pt();
1824 Int_t trackLabel = TMath::Abs(track->GetLabel());
1826 AliVParticle * mcpart =
GetMC()->GetTrack(trackLabel);
1827 if( !mcpart )
continue;
1829 Int_t partInConeCharge = TMath::Abs(mcpart->Charge());
1830 Int_t partInConePDG = mcpart->PdgCode();
1831 Bool_t physPrimary = mcpart->IsPhysicalPrimary();
1833 if ( partInConeCharge > 0 && TMath::Abs(partInConePDG) != 11 )
1836 if ( TMath::Abs(partInConePDG) == 211 ) mcChTag = 0;
1837 else if ( TMath::Abs(partInConePDG) == 321 ) mcChTag = 1;
1838 else if ( TMath::Abs(partInConePDG) == 2212 ) mcChTag = 2;
1854 coneptsumTrack+=pTtrack;
1855 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1863 coneptsumTrackPerMinCut[icut]+=pTtrack;
1864 coneNTrackPerMinCut [icut]++;
1867 if ( pTtrack <
fMaxPtCutInCone[icut] ) coneptsumTrackPerMaxCut[icut]+=pTtrack;
1875 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1886 coneptsumTrackPerRCut[icut]+=pTtrack;
1897 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1902 coneptsumTrackPerNCellCut[icut]+=pTtrack;
1915 coneptsumTrackPerExoCut[icut]+=pTtrack;
1926 phitrack = track->Phi();
1927 etatrack = track->Eta();
1934 ULong_t status = track->GetStatus();
1935 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1936 trackBC = track->GetTOFBunchCrossing(bz);
1937 Double32_t tof = track->GetTOFsignal()*1e-3;
1956 coneptsumTrackTOFBC0 += pTtrack;
1964 coneptsumTrackTOFBCN += pTtrack;
1973 coneptsumTrackTOFNo += pTtrack;
1977 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1978 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1986 coneptsumTrackITSRefitOnSPDOff += pTtrack;
1994 coneptsumTrackITSRefitOffSPDOff += pTtrack;
2003 coneptsumTrackITSRefitOnSPDOn += pTtrack;
2010 if(okTOF && trackBC == 0 && !bConstrained)
2016 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
2029 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
2030 dcaCons = aodTrack->DCA();
2033 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
2133 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
2164 if(energy < 8)
return;
2170 Int_t nVerticesSPD = -1;
2171 Int_t nVerticesTracks = -1;
2175 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
2176 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
2181 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
2182 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
2194 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
2198 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
2199 ncont=pv->GetNContributors();
2200 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
2202 diamZ = esdEv->GetDiamondZ();
2206 AliAODVertex *pv=aodEv->GetVertex(iVert);
2207 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
2208 ncont=pv->GetNContributors();
2209 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
2211 diamZ = aodEv->GetDiamondZ();
2214 Double_t distZ = TMath::Abs(z2-z1);
2215 diamZ = TMath::Abs(z2-diamZ);
2250 if(decayTag < 0) decayTag = 0;
2287 Int_t ptsumBin = -1;
2288 Int_t leadptBin = -1;
2290 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
2300 if ( leadptBin >= 0 )
2309 if ( leadptBin == 0 )
2310 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
2314 if ( ptsumBin >= 0 )
2327 if(decayTag < 0) decayTag = 0;
2389 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
2409 Int_t ptTrigBin = -1;
2418 if ( ptTrigBin >= 0 )
2447 if(decayTag < 0) decayTag = 0;
2542 Int_t noverlaps = 0;
2546 Int_t overpdg[nlabels];
2547 Int_t overlab[nlabels];
2549 GetMC(), overpdg, overlab);
2577 if ( noverlaps == 1 )
2604 else if (noverlaps == 0 )
2643 if ( ptTrigBin >= 0 )
2692 if ( noverlaps == 0 )
2737 else if ( noverlaps == 1 )
2843 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2849 Float_t dEdx = track->GetTPCsignal();
2889 const Int_t buffersize = 255;
2890 char onePar[buffersize] ;
2892 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2896 snprintf(onePar, buffersize,
"Isolation Cand. Detector: %s;",
fIsoDetectorString.Data()) ;
2900 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2902 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2904 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2906 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2911 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2913 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2918 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2923 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2928 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2933 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2944 return new TObjString(parList) ;
2954 outputContainer->SetName(
"IsolatedParticleHistos") ;
2998 Int_t cellBins = 15;
3016 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
3019 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
3024 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
3027 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
3031 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
3034 TString sParticle =
", x^{0,#pm}";
3038 TString parTitle = Form(
"#it{R} = %2.2f%s%s" ,
GetIsolationCut()->GetConeSize(),sThreshold.Data(),sParticle.Data());
3045 parTitle = Form(
"%2.2f<#it{R}<%2.2f%s%s" ,
3047 parTitleR = Form(
"%2.2f<#it{R}<%2.2f%s" ,
3049 parTitleRCh = Form(
"%2.2f<#it{R}<%2.2f, x^{#pm}",
3051 parTitleRNe = Form(
"%2.2f<#it{R}<%2.2f, x^{0}" ,
3055 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
3058 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
3059 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
3060 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
3061 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
3063 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
3064 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
3065 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
3066 "OtherDecay",
"Electron",
"Hadron"} ;
3069 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
3070 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
3071 "#pi^{0}" ,
"#eta"} ;
3073 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
3074 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
3080 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
3081 nptbins,ptmin,ptmax);
3082 fhENoIso->SetYTitle(
"#it{counts}");
3083 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
3087 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
3088 nptbins,ptmin,ptmax);
3090 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3096 (
"hENoIsoExoTrigger",
3097 Form(
"Number of not isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3098 nptbins,ptmin,ptmax);
3104 (
"hPtNoIsoExoTrigger",
3105 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3106 nptbins,ptmin,ptmax);
3113 Form(
"Number of not isolated leading particles #eta vs #varphi, %s",parTitle.Data()),
3114 netabins,etamin,etamax,nphibins,phimin,phimax);
3125 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
3126 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3127 nptbins,ptmin,ptmax);
3129 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3132 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
3133 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3134 nptbins,ptmin,ptmax);
3135 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
3136 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3140 Form(
"#varphi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3141 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3143 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3147 Form(
"#eta vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3148 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3150 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3158 TString isoName [] = {
"NoIso",
"Iso"};
3159 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3163 for(
Int_t iso = 0; iso < 2; iso++)
3167 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3168 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3169 nptbins,ptmin,ptmax);
3170 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
3171 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3172 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
3176 Form(
"Number of %s leading Pi0 decay particles #eta vs #varphi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3177 netabins,etamin,etamax,nphibins,phimin,phimax);
3185 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3186 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
3187 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3198 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
3199 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
3200 nptbins,ptmin,ptmax);
3201 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
3202 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3203 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
3212 TString isoName [] = {
"NoIso",
"Iso"};
3213 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3216 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
3217 nptbins,ptmin,ptmax);
3218 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
3219 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
3220 outputContainer->Add(
fhEIso) ;
3223 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3224 nptbins,ptmin,ptmax);
3225 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
3226 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3227 outputContainer->Add(
fhPtIso) ;
3233 Form(
"Number of isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3234 nptbins,ptmin,ptmax);
3240 (
"hPtIsoExoTrigger",
3241 Form(
"Number of isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3242 nptbins,ptmin,ptmax);
3249 Form(
"Number of isolated particles vs #varphi, %s",parTitle.Data()),
3250 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3251 fhPhiIso->SetYTitle(
"#varphi (rad)");
3252 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3256 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
3257 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3259 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3263 Form(
"Number of isolated particles #eta vs #varphi, %s",parTitle.Data()),
3264 netabins,etamin,etamax,nphibins,phimin,phimax);
3272 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3273 nptbins,ptmin,ptmax, 100,0,100);
3279 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3280 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
3289 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3290 nptbins,ptmin,ptmax,10,0,10);
3295 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3296 nptbins,ptmin,ptmax,10,0,10);
3303 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3304 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3305 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
3306 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3311 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3312 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3313 fhConeSumPt->SetYTitle(
"#Sigma #it{p}_{T} (GeV/#it{c})");
3314 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3321 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3324 (Form(
"hConePtSum_SM%d",ism),
3325 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, SM %d",r,ism),
3326 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3335 for(
Int_t itc = 0; itc < 16; itc++)
3338 (Form(
"hConePtSum_TC%d",itc),
3339 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, TC index %d",r,itc),
3340 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3350 (
"hConePtSumExoTrigger",
3351 Form(
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3352 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3367 Form(
"Trigger #eta vs #varphi, #Sigma #it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3368 netabins,etamin,etamax,nphibins,phimin,phimax);
3375 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s",parTitleR.Data()),
3376 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3377 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
3378 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3385 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3388 (Form(
"hPtInCone_SM%d",ism),
3389 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s, SM %d", parTitleR.Data(),ism),
3390 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3399 for(
Int_t itc = 0; itc < 16; itc++)
3402 (Form(
"hPtInCone_TC%d",itc),
3403 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s, TC index %d", parTitleR.Data(),itc),
3404 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3414 Form(
"#it{p}_{T} in isolation cone for %s, exotic trigger",parTitleR.Data()),
3415 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3457 (Form(
"hPtLeadCone_Bin%d",ibin),
3458 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3465 (Form(
"hSumPtCone_Bin%d",ibin),
3466 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3475 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d",ibin),
3476 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3490 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3491 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3498 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3499 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3513 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3514 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3521 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3522 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3531 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3532 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3544 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
3545 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3552 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
3553 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3562 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d",ibin),
3563 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3576 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3577 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3584 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3585 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3594 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3595 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3674 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
3675 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3682 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
3683 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3692 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d",ibin),
3693 Form(
"#Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3700 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d",ibin),
3701 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3715 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3716 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3723 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3724 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3733 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3734 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3741 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3742 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3757 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3758 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3765 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3766 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3775 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3776 Form(
"#Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3783 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3784 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3796 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
3797 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3804 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
3805 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3814 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d",ibin),
3815 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3822 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d",ibin),
3823 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3836 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3837 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3839 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3845 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3846 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3848 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3856 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3857 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",
3859 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3865 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3866 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",
3868 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3877 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3878 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",
3880 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3886 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3887 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",
3889 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3897 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3898 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",
3900 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3906 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3907 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",
3909 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3915 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3916 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",
3918 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3924 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3925 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",
3927 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3943 Form(
"#it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3944 100,0,100,nptinconebins,ptinconemin,ptinconemax);
3954 (
"hConePtSumCluster",
3955 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3956 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3965 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3968 (Form(
"hConePtSumCluster_SM%d",ism),
3969 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, SM %d",r,ism),
3970 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3979 for(
Int_t itc = 0; itc < 16; itc++)
3982 (Form(
"hConePtSumCluster_TC%d",itc),
3983 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, TC index %d",r,itc),
3984 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3994 (
"hConePtSumClusterExoTrigger",
3995 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3996 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4005 (
"hConeNClusterPerMinPtCut",
4006 Form(
"N clusters, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4015 (
"hConeNClusterPerMinPtCutLargePtTrig",
4016 Form(
"N cluster, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4025 (
"hConePtSumClusterPerMinPtCut",
4026 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4035 (
"hConePtSumClusterPerMinPtCutLargePtTrig",
4036 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),
4046 (
"hConePtSumClusterPerMaxPtCut",
4047 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4056 (
"hConePtSumClusterPerMaxPtCutLargePtTrig",
4057 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),
4069 (
"hConePtSumClusterPerRCut",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts",
4078 (
"hConePtSumClusterPerRCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4087 (
"hPtClusterInConePerRCut",
"Cluster #it{p}_{T}, different #it{R} cuts",
4096 (
"hPtClusterInConePerRCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4108 (
"hConePtSumClusterPerNCellCut",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
4117 (
"hConePtSumClusterPerNCellCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4126 (
"hPtClusterInConePerNCellCut",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts",
4135 (
"hPtClusterInConePerNCellCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4143 for(
Int_t ishsh = 0; ishsh < 4; ishsh++)
4146 (Form(
"hPtClusterInConePerNCellPerSM_ShSh%d",ishsh),
4147 Form(
"Cluster #it{p}_{T} vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
4155 (Form(
"hConeSumPtClusterPerNCellPerSM_ShSh%d",ishsh),
4156 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),
4168 (
"hConePtSumClusterPerExoCut",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts",
4177 (
"hConePtSumClusterPerExoCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4186 (
"hPtClusterInConePerExoCut",
"Cluster #it{p}_{T}, different exoticity cuts",
4195 (
"hPtClusterInConePerExoCutLargePtTrig",
"Cluster #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4205 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
4206 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4215 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4216 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4225 "#Sigma cluster #it{p}_{T} in UE Eta Band",
4226 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4232 "#Sigma cluster #it{p}_{T} UE Phi Band",
4233 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4241 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} in UE Eta Band",
4242 netabins,etamin,etamax,nphibins,phimin,phimax);
4249 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} UE Phi Band",
4250 netabins,etamin,etamax,nphibins,phimin,phimax);
4257 Form(
"#eta vs #varphi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4258 netabins,-1,1,nphibins,0,TMath::TwoPi());
4264 Form(
"#eta vs #varphi of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4265 netabins,-1,1,nphibins,0,TMath::TwoPi());
4271 Form(
"#eta vs #varphi of clusters in cone for #it{R} = %2.2f",r),
4272 netabins,-1,1,nphibins,0,TMath::TwoPi());
4278 Form(
"#eta vs #varphi of all clusters"),
4279 netabins,-1,1,nphibins,0,TMath::TwoPi());
4288 "#Sigma cell #it{p}_{T} in UE Eta Band",
4289 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4295 "#Sigma cell #it{p}_{T} UE Phi Band",
4296 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4302 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} in UE Eta Band",
4303 netabins,etamin,etamax,nphibins,phimin,phimax);
4310 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} UE Phi Band",
4311 netabins,etamin,etamax,nphibins,phimin,phimax);
4319 Form(
"#it{p}_{T} of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4320 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4326 Form(
"#it{p}_{T} of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4327 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4334 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
4335 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4344 if ( ism < fFirstModule || ism >
fLastModule )
continue;
4347 (Form(
"hPtClusterInCone_SM%d",ism),
4348 Form(
"#it{p}_{T} of clusters in isolation cone for %s, SM %d", parTitleR.Data(),ism),
4349 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4358 for(
Int_t itc = 0; itc < 16; itc++)
4361 (Form(
"hPtClusterInCone_TC%d",itc),
4362 Form(
"#it{p}_{T} of clusters in isolation cone for %s, TC index %d", parTitleR.Data(),itc),
4363 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4373 (
"hPtClusterInConeExoTrigger",
4374 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4375 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4384 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
4385 nptbins,ptmin,ptmax,1000,0,50);
4391 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
4398 Form(
"#col vs #row of cells in #varphi band isolation cone for #it{R} = %2.2f",r),
4408 Form(
"Clusters #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
4409 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4415 Form(
"Clusters #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
4416 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4422 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4423 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4429 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4430 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4438 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),
4439 netabins,etamin,etamax,nphibins,phimin,phimax);
4446 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),
4447 netabins,etamin,etamax,nphibins,phimin,phimax);
4457 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4458 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4464 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4465 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4471 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),
4472 netabins,etamin,etamax,nphibins,phimin,phimax);
4479 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),
4480 netabins,etamin,etamax,nphibins,phimin,phimax);
4490 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
4491 nptbins,ptmin,ptmax,100,0,1);
4497 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#varphi ",r),
4498 netabins,etamin,etamax,nphibins,phimin,phimax);
4505 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance",r),
4506 nptbins,ptmin,ptmax,100,0,1);
4512 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance, in trigger #eta-#varphi ",r),
4513 netabins,etamin,etamax,nphibins,phimin,phimax);
4521 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),
4522 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4528 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),
4529 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4535 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),
4536 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4542 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),
4543 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4549 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #eta band for cluster (before normalization), R=%2.2f",r),
4550 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
4556 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #varphi band for cluster (before normalization), R=%2.2f",r),
4557 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
4566 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
4567 nptbins,ptmin,ptmax,100,0,1);
4573 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#varphi ",r),
4574 netabins,etamin,etamax,nphibins,phimin,phimax);
4581 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance",r),
4582 nptbins,ptmin,ptmax,100,0,1);
4588 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance, in trigger #eta-#varphi ",r),
4589 netabins,etamin,etamax,nphibins,phimin,phimax);
4597 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),
4598 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4604 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),
4605 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4611 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),
4612 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4618 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),
4619 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4631 (
"hConeLeadPtTrack",
4632 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
4633 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4640 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
4641 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4650 if ( ism < fFirstModule || ism >
fLastModule )
continue;
4653 (Form(
"hPtTrackInCone_SM%d",ism),
4654 Form(
"#it{p}_{T} of tracks in isolation cone for %s, SM %d", parTitleR.Data(),ism),
4655 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4664 for(
Int_t itc = 0; itc < 16; itc++)
4667 (Form(
"hPtTrackInCone_TC%d",itc),
4668 Form(
"#it{p}_{T} of tracks in isolation cone for %s, TC index %d", parTitleR.Data(),itc),
4669 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4678 TString mcChPartName[] = {
"Pion",
"Kaon",
"Proton",
"Other"};
4679 for(
Int_t imc = 0; imc < 4; imc++)
4682 (Form(
"hPtTrackInCone_Primary_%s",mcChPartName[imc].
Data()),
4683 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4684 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);