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++)
447 for(
Int_t iso =0; iso < 2; iso++)
467 for(
Int_t ism =0; ism < 20; ism++)
510 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
512 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
516 AliWarning(
"Cluster not available?");
522 cluster->GetID() == pCandidate->
GetCaloLabel(1) )
continue ;
586 if(phiTrig<0) phiTrig += TMath::TwoPi();
592 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
595 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
597 if(!eGeom->CheckAbsCellId(absId))
return ;
600 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
602 Int_t imEta=-1, imPhi=-1;
603 Int_t ieta =-1, iphi =-1;
605 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
607 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
609 Int_t colTrig = ieta;
610 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
612 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
614 Int_t sqrSize = int(conesize/0.0143);
618 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
619 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
623 Int_t irowmin = rowTrig-sqrSize;
624 if ( irowmin < 0 ) irowmin = 0 ;
626 Int_t irowmax = rowTrig+sqrSize;
627 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
629 for(
Int_t irow = irowmin; irow <irowmax; irow++)
631 for(
Int_t icol = 0; icol < nTotalCols; icol++)
633 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
634 if(inSector==5)
continue;
637 if(icol < AliEMCALGeoParams::fgkEMCALCols)
639 inSupMod = 2*inSector + 1;
642 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
644 inSupMod = 2*inSector;
645 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
648 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
651 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
655 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
656 if(!eGeom->CheckAbsCellId(iabsId))
continue;
658 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
666 Int_t icolmin = colTrig-sqrSize;
667 if ( icolmin < 0 ) icolmin = 0;
669 Int_t icolmax = colTrig+sqrSize;
670 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
673 for(
Int_t icol = icolmin; icol < icolmax; icol++)
675 for(
Int_t irow = 0; irow < nTotalRows; irow++)
677 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
678 if ( inSector == 5 ) continue ;
684 if(icol < AliEMCALGeoParams::fgkEMCALCols)
687 inSupMod = 2*inSector + 1;
690 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
693 inSupMod = 2*inSector;
694 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
697 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
700 if(TMath::Abs(icol-colTrig) < sqrSize) {
703 if(TMath::Abs(irow-rowTrig) < sqrSize) {
707 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
708 if(!eGeom->CheckAbsCellId(iabsId))
710 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
714 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
747 Float_t coneptsumPerpTrackPerMinCut[20];
748 Float_t coneNPerpTrackPerMinCut [20];
754 coneptsumPerpTrackPerMinCut[icut] = 0;
755 coneNPerpTrackPerMinCut [icut] = 0;
766 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
768 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
772 AliWarning(
"Track not available?");
782 Bool_t contained = kFALSE;
784 for(
Int_t i = 0; i < 4; i++)
786 if( trackID == pCandidate->
GetTrackLabel(i) ) contained = kTRUE;
789 if ( contained ) continue ;
807 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
809 phiBandPtSum+=track->Pt();
816 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
818 etaBandPtSum+=track->Pt();
826 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
827 Double_t dEta = etaTrig - track->Eta();
828 Double_t arg = dPhi*dPhi + dEta*dEta;
830 if(TMath::Sqrt(arg) < conesize)
832 sumptPerp+=track->Pt();
840 ULong_t status = track->GetStatus();
841 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
842 Int_t trackBC = track->GetTOFBunchCrossing(bz);
845 if(okTOF && trackBC == 0)
850 sumptPerpBC0+=track->Pt();
853 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
860 sumptPerpITSSPD+=track->Pt();
863 if(okTOF && trackBC == 0 && !bConstrained)
868 sumptPerpBC0ITSSPD+=track->Pt();
874 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
875 arg = dPhi*dPhi + dEta*dEta;
877 if(TMath::Sqrt(arg) < conesize)
879 sumptPerp+=track->Pt();
887 coneptsumPerpTrackPerMinCut[icut]+=track->Pt();
888 coneNPerpTrackPerMinCut [icut]++;
899 ULong_t status = track->GetStatus();
900 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
901 Int_t trackBC = track->GetTOFBunchCrossing(bz);
904 if(okTOF && trackBC == 0)
909 sumptPerpBC0+=track->Pt();
912 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
919 sumptPerpITSSPD+=track->Pt();
922 if(okTOF && trackBC == 0 && !bConstrained)
927 sumptPerpBC0ITSSPD+=track->Pt();
938 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
939 dcaCons = aodTrack->DCA();
942 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
1003 Float_t etaUEptsumCluster = 0 ;
1004 Float_t phiUEptsumCluster = 0 ;
1013 Float_t coneA = conesize*conesize*TMath::Pi();
1024 Float_t phiUEptsumTrackNorm = 0 ;
1025 Float_t etaUEptsumTrackNorm = 0 ;
1026 Float_t coneptsumTrackSubPhi = 0 ;
1027 Float_t coneptsumTrackSubEta = 0 ;
1028 Float_t coneptsumTrackSubPhiNorm = 0 ;
1029 Float_t coneptsumTrackSubEtaNorm = 0 ;
1043 Float_t correctConeSumTrack = 1;
1044 Float_t correctConeSumTrackPhi = 1;
1047 phiUEptsumTrack,etaUEptsumTrack,
1048 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
1049 correctConeSumTrack,correctConeSumTrackPhi);
1051 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
1052 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
1060 if(coneptsumTrack > 0)
1062 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
1063 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
1084 Float_t phiUEptsumClusterNorm = 0 ;
1085 Float_t etaUEptsumClusterNorm = 0 ;
1086 Float_t coneptsumClusterSubPhi = 0 ;
1087 Float_t coneptsumClusterSubEta = 0 ;
1088 Float_t coneptsumClusterSubPhiNorm = 0 ;
1089 Float_t coneptsumClusterSubEtaNorm = 0 ;
1090 Float_t phiUEptsumCellNorm = 0 ;
1091 Float_t etaUEptsumCellNorm = 0 ;
1092 Float_t coneptsumCellSubPhi = 0 ;
1093 Float_t coneptsumCellSubEta = 0 ;
1094 Float_t coneptsumCellSubPhiNorm = 0 ;
1095 Float_t coneptsumCellSubEtaNorm = 0 ;
1113 Float_t correctConeSumClusterEta = 1;
1114 Float_t correctConeSumClusterPhi = 1;
1117 phiUEptsumCluster,etaUEptsumCluster,
1118 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
1119 correctConeSumClusterEta,correctConeSumClusterPhi);
1130 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
1131 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
1138 if(coneptsumCluster!=0)
1140 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
1141 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
1173 Float_t emcPhiSize = TMath::DegToRad()*100.;
1175 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
1176 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
1180 Float_t correctConeSumCellEta = 1;
1181 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
1183 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
1187 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
1190 Float_t correctConeSumCellPhi = 1;
1192 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
1193 (phiTrig-conesize < 80*TMath::DegToRad()))
1196 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
1197 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
1203 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
1207 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
1208 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1219 if ( coneptsumCell > 0.01 )
1221 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1222 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1232 sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1233 sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1263 Int_t ptsumAfterEtaBandUESubBin = -1;
1265 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, n bins %d",pt,coneptsumTrack+coneptsumCluster,
fNBkgBin));
1274 if ( ptsumAfterEtaBandUESubBin >= 0 )
1276 AliDebug(1,Form(
"\t Sum bin %d [%2.2f,%2.2f]" , ptsumAfterEtaBandUESubBin ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin] ,
fBkgBinLimit[ptsumAfterEtaBandUESubBin +1]));
1285 Int_t ptsumAfterEtaBandSubBinMC = ptsumAfterEtaBandUESubBin+mcIndex*
fNBkgBin;
1287 if( ptsumAfterEtaBandUESubBin >=0 )
1297 if( ptsumAfterEtaBandUESubBin >=0 )
1316 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1318 if( ptsumAfterEtaBandUESubBin >=0 )
1337 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1338 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1362 coneptLeadCluster = 0;
1363 coneptsumCluster = 0;
1384 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1385 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1396 Float_t coneNClusterPerMinCut [20];
1397 Float_t coneptsumClusterPerMinCut[20];
1398 Float_t coneptsumClusterPerMaxCut[20];
1399 Float_t coneptsumClusterPerRCut [10];
1401 Float_t coneptsumClusterPerNCellCut[20];
1402 Float_t coneptsumClusterPerExoCut [20];
1408 coneNClusterPerMinCut [icut] = 0;
1409 coneptsumClusterPerMinCut[icut] = 0;
1410 coneptsumClusterPerMaxCut[icut] = 0;
1418 coneptsumClusterPerRCut[icut] = 0;
1426 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1427 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1428 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1429 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1433 coneptsumClusterPerNCellCut[icut] = 0;
1441 coneptsumClusterPerExoCut[icut] = 0;
1445 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1447 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1486 coneptsumCluster+=ptcone;
1487 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1495 coneptsumClusterPerMinCut[icut]+=ptcone;
1496 coneNClusterPerMinCut [icut]++;
1499 if ( ptcone <
fMaxPtCutInCone[icut] ) coneptsumClusterPerMaxCut[icut]+=ptcone;
1510 coneptsumClusterPerRCut[icut]+=ptcone;
1519 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1526 coneptsumClusterPerNCellCut[icut]+=ptcone;
1539 coneptsumClusterPerExoCut[icut]+=ptcone;
1590 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1625 if(phiTrig<0) phiTrig += TMath::TwoPi();
1630 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1633 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1635 if(!eGeom->CheckAbsCellId(absId))
return ;
1638 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1640 Int_t imEta=-1, imPhi=-1;
1641 Int_t ieta =-1, iphi =-1;
1643 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1645 Int_t iEta=-1, iPhi=-1;
1646 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1648 Int_t colTrig = iEta;
1649 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1650 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1652 Int_t sqrSize = int(conesize/0.0143);
1657 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1659 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1661 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1662 if(inSector==5)
continue;
1664 Int_t inSupMod = -1;
1666 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1668 inSupMod = 2*inSector + 1;
1671 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1673 inSupMod = 2*inSector;
1674 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1677 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1679 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1680 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1683 coneptsumCell += cells->GetCellAmplitude(iabsId);
1699 coneptLeadTrack = 0;
1733 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1734 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1743 Float_t coneNTrackPerMinCut [20];
1744 Float_t coneptsumTrackPerMinCut[20];
1745 Float_t coneptsumTrackPerMaxCut[20];
1746 Float_t coneptsumTrackPerEtaCut[10];
1747 Float_t coneptsumTrackPerNCellCut[20];
1748 Float_t coneptsumTrackPerExoCut[20];
1749 Float_t coneptsumTrackPerRCut [10];
1750 Float_t coneptsumTrackTOFBC0 = 0;
1751 Float_t coneptsumTrackTOFBCN = 0;
1752 Float_t coneptsumTrackTOFNo = 0;
1753 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1754 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1755 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1756 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1762 coneptsumTrackPerMinCut[icut] = 0;
1763 coneptsumTrackPerMaxCut[icut] = 0;
1764 coneNTrackPerMinCut [icut] = 0;
1772 coneptsumTrackPerEtaCut[icut] = 0;
1780 coneptsumTrackPerRCut[icut] = 0;
1788 if ( m02 > 0.1 && m02 <= 0.3 ) ishsh = 0;
1789 else if ( m02 > 0.3 && m02 <= 0.4 ) ishsh = 1;
1790 else if ( m02 > 0.4 && m02 <= 1.0 ) ishsh = 2;
1791 else if ( m02 > 1.0 && m02 <= 3.0 ) ishsh = 3;
1795 coneptsumTrackPerNCellCut[icut] = 0;
1803 coneptsumTrackPerExoCut[icut] = 0;
1807 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1809 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1811 pTtrack = track->Pt();
1830 Int_t trackLabel = TMath::Abs(track->GetLabel());
1832 AliVParticle * mcpart =
GetMC()->GetTrack(trackLabel);
1833 if( !mcpart )
continue;
1835 Int_t partInConeCharge = TMath::Abs(mcpart->Charge());
1836 Int_t partInConePDG = mcpart->PdgCode();
1837 Bool_t physPrimary = mcpart->IsPhysicalPrimary();
1839 if ( partInConeCharge > 0 && TMath::Abs(partInConePDG) != 11 )
1842 if ( TMath::Abs(partInConePDG) == 211 ) mcChTag = 0;
1843 else if ( TMath::Abs(partInConePDG) == 321 ) mcChTag = 1;
1844 else if ( TMath::Abs(partInConePDG) == 2212 ) mcChTag = 2;
1860 coneptsumTrack+=pTtrack;
1861 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1869 coneptsumTrackPerMinCut[icut]+=pTtrack;
1870 coneNTrackPerMinCut [icut]++;
1873 if ( pTtrack <
fMaxPtCutInCone[icut] ) coneptsumTrackPerMaxCut[icut]+=pTtrack;
1881 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1892 coneptsumTrackPerRCut[icut]+=pTtrack;
1903 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
1908 coneptsumTrackPerNCellCut[icut]+=pTtrack;
1921 coneptsumTrackPerExoCut[icut]+=pTtrack;
1932 phitrack = track->Phi();
1933 etatrack = track->Eta();
1940 ULong_t status = track->GetStatus();
1941 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1942 trackBC = track->GetTOFBunchCrossing(bz);
1943 Double32_t tof = track->GetTOFsignal()*1e-3;
1962 coneptsumTrackTOFBC0 += pTtrack;
1970 coneptsumTrackTOFBCN += pTtrack;
1979 coneptsumTrackTOFNo += pTtrack;
1983 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1984 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1992 coneptsumTrackITSRefitOnSPDOff += pTtrack;
2000 coneptsumTrackITSRefitOffSPDOff += pTtrack;
2009 coneptsumTrackITSRefitOnSPDOn += pTtrack;
2016 if(okTOF && trackBC == 0 && !bConstrained)
2022 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
2035 AliAODTrack * aodTrack =
dynamic_cast<AliAODTrack*
>(track);
2036 dcaCons = aodTrack->DCA();
2039 track->PropagateToDCA(
GetReader()->GetInputEvent()->GetPrimaryVertex(),bz,100.,dca,covar);
2139 if ( ptTrig > 8 && ptTrig < 12 && ishsh >=0 )
2170 if(energy < 8)
return;
2176 Int_t nVerticesSPD = -1;
2177 Int_t nVerticesTracks = -1;
2181 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
2182 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
2187 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
2188 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
2200 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
2204 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
2205 ncont=pv->GetNContributors();
2206 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
2208 diamZ = esdEv->GetDiamondZ();
2212 AliAODVertex *pv=aodEv->GetVertex(iVert);
2213 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
2214 ncont=pv->GetNContributors();
2215 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
2217 diamZ = aodEv->GetDiamondZ();
2220 Double_t distZ = TMath::Abs(z2-z1);
2221 diamZ = TMath::Abs(z2-diamZ);
2256 if(decayTag < 0) decayTag = 0;
2293 Int_t ptsumBin = -1;
2294 Int_t leadptBin = -1;
2296 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
2306 if ( leadptBin >= 0 )
2315 if ( leadptBin == 0 )
2316 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
2320 if ( ptsumBin >= 0 )
2333 if(decayTag < 0) decayTag = 0;
2395 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
2415 Int_t ptTrigBin = -1;
2424 if ( ptTrigBin >= 0 )
2453 if(decayTag < 0) decayTag = 0;
2570 Int_t noverlaps = 0;
2574 Int_t overpdg[nlabels];
2575 Int_t overlab[nlabels];
2577 GetMC(), overpdg, overlab);
2605 if ( noverlaps == 1 )
2632 else if (noverlaps == 0 )
2671 if ( ptTrigBin >= 0 )
2720 if ( noverlaps == 0 )
2765 else if ( noverlaps == 1 )
2871 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2877 Float_t dEdx = track->GetTPCsignal();
2917 const Int_t buffersize = 255;
2918 char onePar[buffersize] ;
2920 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2924 snprintf(onePar, buffersize,
"Isolation Cand. Detector: %s;",
fIsoDetectorString.Data()) ;
2928 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2930 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2932 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2934 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2939 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2941 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2946 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2951 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2956 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2961 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2972 return new TObjString(parList) ;
2982 outputContainer->SetName(
"IsolatedParticleHistos") ;
3026 Int_t cellBins = 15;
3044 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
3047 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
3052 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
3055 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
3059 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
3062 TString sParticle =
", x^{0,#pm}";
3066 TString parTitle = Form(
"#it{R} = %2.2f%s%s" ,
GetIsolationCut()->GetConeSize(),sThreshold.Data(),sParticle.Data());
3073 parTitle = Form(
"%2.2f<#it{R}<%2.2f%s%s" ,
3075 parTitleR = Form(
"%2.2f<#it{R}<%2.2f%s" ,
3077 parTitleRCh = Form(
"%2.2f<#it{R}<%2.2f, x^{#pm}",
3079 parTitleRNe = Form(
"%2.2f<#it{R}<%2.2f, x^{0}" ,
3083 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
3086 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
3087 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
3088 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
3089 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
3091 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
3092 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
3093 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
3094 "OtherDecay",
"Electron",
"Hadron"} ;
3097 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
3098 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
3099 "#pi^{0}" ,
"#eta"} ;
3101 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
3102 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
3108 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
3109 nptbins,ptmin,ptmax);
3110 fhENoIso->SetYTitle(
"#it{counts}");
3111 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
3115 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
3116 nptbins,ptmin,ptmax);
3118 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3124 (
"hENoIsoExoTrigger",
3125 Form(
"Number of not isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3126 nptbins,ptmin,ptmax);
3132 (
"hPtNoIsoExoTrigger",
3133 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3134 nptbins,ptmin,ptmax);
3141 Form(
"Number of not isolated leading particles #eta vs #varphi, %s",parTitle.Data()),
3142 netabins,etamin,etamax,nphibins,phimin,phimax);
3153 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
3154 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3155 nptbins,ptmin,ptmax);
3157 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3160 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
3161 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3162 nptbins,ptmin,ptmax);
3163 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
3164 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3168 Form(
"#varphi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3169 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3171 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3175 Form(
"#eta vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
3176 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3178 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3186 TString isoName [] = {
"NoIso",
"Iso"};
3187 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3191 for(
Int_t iso = 0; iso < 2; iso++)
3195 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3196 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3197 nptbins,ptmin,ptmax);
3198 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
3199 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3200 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
3204 Form(
"Number of %s leading Pi0 decay particles #eta vs #varphi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
3205 netabins,etamin,etamax,nphibins,phimin,phimax);
3213 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
3214 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
3215 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3226 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
3227 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
3228 nptbins,ptmin,ptmax);
3229 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
3230 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3231 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
3240 TString isoName [] = {
"NoIso",
"Iso"};
3241 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
3244 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
3245 nptbins,ptmin,ptmax);
3246 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
3247 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
3248 outputContainer->Add(
fhEIso) ;
3251 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3252 nptbins,ptmin,ptmax);
3253 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
3254 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3255 outputContainer->Add(
fhPtIso) ;
3261 Form(
"Number of isolated particles vs E, %s, exot>0.97",parTitle.Data()),
3262 nptbins,ptmin,ptmax);
3268 (
"hPtIsoExoTrigger",
3269 Form(
"Number of isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
3270 nptbins,ptmin,ptmax);
3277 Form(
"Number of isolated particles vs #varphi, %s",parTitle.Data()),
3278 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
3279 fhPhiIso->SetYTitle(
"#varphi (rad)");
3280 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3284 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
3285 nptbins,ptmin,ptmax,netabins,etamin,etamax);
3287 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3291 Form(
"Number of isolated particles #eta vs #varphi, %s",parTitle.Data()),
3292 netabins,etamin,etamax,nphibins,phimin,phimax);
3300 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3301 nptbins,ptmin,ptmax, 100,0,100);
3307 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
3308 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
3317 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3318 nptbins,ptmin,ptmax,10,0,10);
3323 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
3324 nptbins,ptmin,ptmax,10,0,10);
3331 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3332 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3333 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
3334 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3339 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3340 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3341 fhConeSumPt->SetYTitle(
"#Sigma #it{p}_{T} (GeV/#it{c})");
3342 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
3349 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3352 (Form(
"hConePtSum_SM%d",ism),
3353 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, SM %d",r,ism),
3354 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3363 for(
Int_t itc = 0; itc < 16; itc++)
3366 (Form(
"hConePtSum_TC%d",itc),
3367 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, TC index %d",r,itc),
3368 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3378 (
"hConePtSumExoTrigger",
3379 Form(
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3380 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3395 Form(
"Trigger #eta vs #varphi, #Sigma #it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3396 netabins,etamin,etamax,nphibins,phimin,phimax);
3403 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s",parTitleR.Data()),
3404 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3405 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
3406 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3413 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3416 (Form(
"hPtInCone_SM%d",ism),
3417 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s, SM %d", parTitleR.Data(),ism),
3418 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3427 for(
Int_t itc = 0; itc < 16; itc++)
3430 (Form(
"hPtInCone_TC%d",itc),
3431 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for %s, TC index %d", parTitleR.Data(),itc),
3432 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3442 Form(
"#it{p}_{T} in isolation cone for %s, exotic trigger",parTitleR.Data()),
3443 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3485 (Form(
"hPtLeadCone_Bin%d",ibin),
3486 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3493 (Form(
"hSumPtCone_Bin%d",ibin),
3494 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3503 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d",ibin),
3504 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3518 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3519 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3526 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3527 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3541 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3542 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3549 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3550 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3559 (Form(
"hSumPtConeAfterEtaBandUESub_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3560 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3572 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
3573 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3580 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
3581 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3590 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d",ibin),
3591 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3604 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3605 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3612 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3613 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3622 (Form(
"hSumPtConeAfterEtaBandUESubLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3623 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3702 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
3703 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3710 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
3711 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3720 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d",ibin),
3721 Form(
"#Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3728 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d",ibin),
3729 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3743 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3744 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3751 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3752 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3761 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3762 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3769 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3770 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3785 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3786 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3793 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3794 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3803 (Form(
"hPtTrigBin_SumPtTrackCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3804 Form(
"#Sigma #it{p}_{T}^{in cone}_{track}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3811 (Form(
"hPtTrigBin_SumPtClusterCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3812 Form(
"#Sigma #it{p}_{T}^{in cone}_{cluster}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3824 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
3825 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3832 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
3833 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3842 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d",ibin),
3843 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{track} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3850 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d",ibin),
3851 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}_{cluster} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3864 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3865 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3867 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3873 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3874 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3876 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3884 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3885 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",
3887 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3893 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3894 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",
3896 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3905 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3906 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",
3908 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3914 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3915 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",
3917 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3925 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3926 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",
3928 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3934 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_NoOverlap%s",ibin,mcPartName[imc].
Data()),
3935 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",
3937 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3943 (Form(
"hPtTrigBin_SumPtTrackConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3944 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",
3946 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3952 (Form(
"hPtTrigBin_SumPtClusterConeVSLambda0_Bin%d_MC_1Overlap%s",ibin,mcPartName[imc].
Data()),
3953 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",
3955 nptsumbins,ptsummin,ptsummax,ssbins,ssmin,ssmax);
3971 Form(
"#it{p}_{T} in isolation cone for %s",parTitleR.Data()),
3972 100,0,100,nptinconebins,ptinconemin,ptinconemax);
3982 (
"hConePtSumCluster",
3983 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3984 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3993 if ( ism < fFirstModule || ism >
fLastModule )
continue;
3996 (Form(
"hConePtSumCluster_SM%d",ism),
3997 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, SM %d",r,ism),
3998 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4007 for(
Int_t itc = 0; itc < 16; itc++)
4010 (Form(
"hConePtSumCluster_TC%d",itc),
4011 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, TC index %d",r,itc),
4012 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4022 (
"hConePtSumClusterExoTrigger",
4023 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
4024 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4033 (
"hConeNClusterPerMinPtCut",
4034 Form(
"N clusters, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4043 (
"hConeNClusterPerMinPtCutLargePtTrig",
4044 Form(
"N cluster, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4053 (
"hConePtSumClusterPerMinPtCut",
4054 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4063 (
"hConePtSumClusterPerMinPtCutLargePtTrig",
4064 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),
4074 (
"hConePtSumClusterPerMaxPtCut",
4075 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4084 (
"hConePtSumClusterPerMaxPtCutLargePtTrig",
4085 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),
4097 (
"hConePtSumClusterPerRCut",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts",
4106 (
"hConePtSumClusterPerRCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4115 (
"hPtClusterInConePerRCut",
"Cluster #it{p}_{T}, different #it{R} cuts",
4124 (
"hPtClusterInConePerRCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4136 (
"hConePtSumClusterPerNCellCut",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
4145 (
"hConePtSumClusterPerNCellCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4154 (
"hPtClusterInConePerNCellCut",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts",
4163 (
"hPtClusterInConePerNCellCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4171 for(
Int_t ishsh = 0; ishsh < 4; ishsh++)
4174 (Form(
"hPtClusterInConePerNCellPerSM_ShSh%d",ishsh),
4175 Form(
"Cluster #it{p}_{T} vs #it{n}_{cell} vs SM , 8 < #it{p}_{T}^{trig} < 12 GeV, sh. shape bin %d",ishsh),
4183 (Form(
"hConeSumPtClusterPerNCellPerSM_ShSh%d",ishsh),
4184 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),
4196 (
"hConePtSumClusterPerExoCut",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts",
4205 (
"hConePtSumClusterPerExoCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4214 (
"hPtClusterInConePerExoCut",
"Cluster #it{p}_{T}, different exoticity cuts",
4223 (
"hPtClusterInConePerExoCutLargePtTrig",
"Cluster #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4233 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
4234 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4243 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4244 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4253 "#Sigma cluster #it{p}_{T} in UE Eta Band",
4254 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4260 "#Sigma cluster #it{p}_{T} UE Phi Band",
4261 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4269 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} in UE Eta Band",
4270 netabins,etamin,etamax,nphibins,phimin,phimax);
4277 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} UE Phi Band",
4278 netabins,etamin,etamax,nphibins,phimin,phimax);
4285 Form(
"#eta vs #varphi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4286 netabins,-1,1,nphibins,0,TMath::TwoPi());
4292 Form(
"#eta vs #varphi of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4293 netabins,-1,1,nphibins,0,TMath::TwoPi());
4299 Form(
"#eta vs #varphi of clusters in cone for #it{R} = %2.2f",r),
4300 netabins,-1,1,nphibins,0,TMath::TwoPi());
4306 Form(
"#eta vs #varphi of all clusters"),
4307 netabins,-1,1,nphibins,0,TMath::TwoPi());
4316 "#Sigma cell #it{p}_{T} in UE Eta Band",
4317 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4323 "#Sigma cell #it{p}_{T} UE Phi Band",
4324 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4330 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} in UE Eta Band",
4331 netabins,etamin,etamax,nphibins,phimin,phimax);
4338 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} UE Phi Band",
4339 netabins,etamin,etamax,nphibins,phimin,phimax);
4347 Form(
"#it{p}_{T} of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
4348 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4354 Form(
"#it{p}_{T} of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
4355 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4362 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
4363 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4372 if ( ism < fFirstModule || ism >
fLastModule )
continue;
4375 (Form(
"hPtClusterInCone_SM%d",ism),
4376 Form(
"#it{p}_{T} of clusters in isolation cone for %s, SM %d", parTitleR.Data(),ism),
4377 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4386 for(
Int_t itc = 0; itc < 16; itc++)
4389 (Form(
"hPtClusterInCone_TC%d",itc),
4390 Form(
"#it{p}_{T} of clusters in isolation cone for %s, TC index %d", parTitleR.Data(),itc),
4391 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4401 (
"hPtClusterInConeExoTrigger",
4402 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4403 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4412 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
4413 nptbins,ptmin,ptmax,1000,0,50);
4419 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
4426 Form(
"#col vs #row of cells in #varphi band isolation cone for #it{R} = %2.2f",r),
4436 Form(
"Clusters #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
4437 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4443 Form(
"Clusters #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
4444 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4450 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4451 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4457 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4458 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4466 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),
4467 netabins,etamin,etamax,nphibins,phimin,phimax);
4474 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),
4475 netabins,etamin,etamax,nphibins,phimin,phimax);
4485 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4486 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4492 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4493 nptbins,ptmin,ptmax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4499 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),
4500 netabins,etamin,etamax,nphibins,phimin,phimax);
4507 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),
4508 netabins,etamin,etamax,nphibins,phimin,phimax);
4518 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
4519 nptbins,ptmin,ptmax,100,0,1);
4525 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#varphi ",r),
4526 netabins,etamin,etamax,nphibins,phimin,phimax);
4533 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance",r),
4534 nptbins,ptmin,ptmax,100,0,1);
4540 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance, in trigger #eta-#varphi ",r),
4541 netabins,etamin,etamax,nphibins,phimin,phimax);
4549 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),
4550 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4556 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),
4557 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4563 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),
4564 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4570 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),
4571 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4577 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #eta band for cluster (before normalization), R=%2.2f",r),
4578 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
4584 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #varphi band for cluster (before normalization), R=%2.2f",r),
4585 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
4594 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
4595 nptbins,ptmin,ptmax,100,0,1);
4601 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#varphi ",r),
4602 netabins,etamin,etamax,nphibins,phimin,phimax);
4609 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance",r),
4610 nptbins,ptmin,ptmax,100,0,1);
4616 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance, in trigger #eta-#varphi ",r),
4617 netabins,etamin,etamax,nphibins,phimin,phimax);
4625 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),
4626 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4632 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),
4633 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4639 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),
4640 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4646 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),
4647 nptsumbins,ptsummin,ptsummax,1.2*nptsumbins,-ptsummax*0.2,ptsummax);
4659 (
"hConeLeadPtTrack",
4660 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
4661 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4668 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
4669 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);