17 #include <TClonesArray.h>
19 #include <TObjString.h>
24 #include "TParticle.h"
25 #include "TDatabasePDG.h"
35 #include "AliAODMCParticle.h"
36 #include "AliAODPWG4ParticleCorrelation.h"
38 #include "AliVTrack.h"
39 #include "AliVCluster.h"
40 #include "AliESDEvent.h"
41 #include "AliAODEvent.h"
44 #include "AliEMCALGeometry.h"
45 #include "AliPHOSGeoUtils.h"
56 fIsoDetector(-1), fIsoDetectorString(
""),
57 fReMakeIC(0), fMakeSeveralIC(0),
58 fFillTMHisto(0), fFillSSHisto(1),
59 fFillEMCALRegionHistograms(0), fFillUEBandSubtractHistograms(1),
60 fFillCellHistograms(0),
61 fFillOverlapHistograms(0), fStudyFECCorrelation(0),
62 fStudyTracksInCone(0), fStudyMCConversionRadius(0),
63 fStudyExoticTrigger(0),
64 fFillTaggedDecayHistograms(0), fNDecayBits(0),
65 fDecayBits(), fDecayTagsM02Cut(0),
66 fFillNLMHistograms(0),
67 fLeadingOnly(0), fCheckLeadingWithNeutralClusters(0),
68 fSelectPrimariesInCone(0), fMakePrimaryPi0DecayStudy(0),
69 fFillBackgroundBinHistograms(0), fNBkgBin(0),
70 fFillPtTrigBinHistograms(0), fNPtTrigBin(0),
71 fMinCellsAngleOverlap(0),
73 fNCones(0), fNPtThresFrac(0),
74 fConeSizes(), fPtThresholds(),
75 fPtFractions(), fSumPtThresholds(),
76 fStudyPtCutInCone(0), fNPtCutsInCone(0),
77 fMinPtCutInCone(), fMaxPtCutInCone(),
78 fStudyEtaCutInCone(0), fNEtaCutsInCone(0), fEtaCutInCone(),
79 fStudyRCutInCone(0), fNRCutsInCone(0), fRCutInCone(),
80 fNNCellsInCandidate(0), fNCellsInCandidate(),
81 fNExoCutInCandidate(0), fExoCutInCandidate(),
82 fMomentum(), fMomIso(),
83 fMomDaugh1(), fMomDaugh2(),
84 fTrackVector(), fProdVertex(),
85 fCluster(0), fClustersArr(0),
86 fIsExoticTrigger(0), fClusterExoticity(1),
88 fhEIso(0), fhPtIso(0),
89 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
91 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
93 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
94 fhEIsoExoTrigger(0), fhENoIsoExoTrigger(0),
95 fhPtIsoExoTrigger(0), fhPtNoIsoExoTrigger(0),
98 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
99 fhPtInConeExoTrigger(0), fhPtClusterInConeExoTrigger(0), fhPtTrackInConeExoTrigger(0),
100 fhPtTrackInConeOtherBCPileUpSPD(0), fhPtTrackInConeVtxBC0(0),
101 fhPtTrackInConeBC0PileUpSPD(0),
102 fhPtInConePileUp(), fhPtInConeCent(0),
104 fhPerpConeSumPt(0), fhPerpConeSumPtTOFBC0(0),
105 fhPtInPerpCone(0), fhPtInPerpConeTOFBC0(0),
106 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
107 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
108 fhEtaPhiInPerpCone(0), fhEtaPhiInPerpConeTOFBC0(0),
109 fhEtaBandClusterEtaPhi(0), fhPhiBandClusterEtaPhi(0),
110 fhEtaBandTrackEtaPhi(0), fhPhiBandTrackEtaPhi(0),
111 fhEtaBandClusterPt(0), fhPhiBandClusterPt(0),
112 fhEtaBandTrackPt(0), fhPhiBandTrackPt(0),
113 fhEtaBandCell(0), fhPhiBandCell(0),
114 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
115 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
117 fhConeSumPtCellTrack(0),
118 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
119 fhConeSumPtExoTrigger(0), fhConeSumPtClusterExoTrigger(0), fhConeSumPtTrackExoTrigger(0),
121 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
122 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
123 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
124 fhConeSumPtTrigEtaPhi(0),
125 fhConeSumPtCellTrackTrigEtaPhi(0),
126 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
127 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
128 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
129 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
130 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
131 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
132 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
133 fhConeSumPtEtaUENormCluster(0), fhConeSumPtPhiUENormCluster(0),
134 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
135 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
136 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
137 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
138 fhConeSumPtEtaUENormTrack(0), fhConeSumPtPhiUENormTrack(0),
139 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
140 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
141 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
142 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
143 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
144 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
145 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
146 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
147 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
148 fhConeSumPtCellvsTrack(0),
149 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
150 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
151 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
152 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
153 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
154 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
155 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
156 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
157 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
158 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
159 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
160 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
161 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
162 fhPtPrimMCPi0DecayPairOutOfCone(0),
163 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
164 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
165 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
166 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
167 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
168 fhPtPrimMCPi0DecayPairNoOverlap(0),
169 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
170 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
171 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
172 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
173 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
174 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
175 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
176 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
177 fhPtPrimMCEtaDecayPairOutOfCone(0),
178 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
179 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
180 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
181 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
182 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
183 fhPtPrimMCEtaDecayPairNoOverlap(0),
184 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
185 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
186 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
187 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
188 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
189 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
190 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
191 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
192 fhPtLeadConeBin(0), fhSumPtConeBin(0),
193 fhPtLeadConeBinMC(0), fhSumPtConeBinMC(0),
194 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
195 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
196 fhPtLeadConeBinLambda0MC(0), fhSumPtConeBinLambda0MC(0),
197 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
198 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
199 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
200 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
201 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
204 fhELambda0LocMax1(), fhELambda1LocMax1(),
205 fhELambda0LocMax2(), fhELambda1LocMax2(),
206 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
208 fhEIsoPileUp(), fhPtIsoPileUp(),
209 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
210 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
211 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
212 fhTimeNPileUpVertContributors(0),
213 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
215 fhPtClusterInConePerRCut(0), fhPtClusterInConePerRCutLargePtTrig(0),
216 fhPtTrackInConePerRCut(0), fhPtTrackInConePerRCutLargePtTrig(0),
217 fhConeSumPtClusterPerRCut(0), fhConeSumPtClusterPerRCutLargePtTrig(0),
218 fhConeSumPtTrackPerRCut(0), fhConeSumPtTrackPerRCutLargePtTrig(0),
219 fhConeSumPtClusterPerMinPtCut(0), fhConeSumPtClusterPerMinPtCutLargePtTrig(0),
220 fhConeSumPtTrackPerMinPtCut(0), fhConeSumPtTrackPerMinPtCutLargePtTrig(0),
221 fhConeSumPtClusterPerMaxPtCut(0), fhConeSumPtClusterPerMaxPtCutLargePtTrig(0),
222 fhConeSumPtTrackPerMaxPtCut(0), fhConeSumPtTrackPerMaxPtCutLargePtTrig(0),
223 fhConeSumPtTrackPerEtaCut(0), fhConeSumPtTrackPerEtaCutLargePtTrig(0),
225 fhPtClusterInConePerNCellCut(0), fhPtClusterInConePerNCellCutLargePtTrig(0),
226 fhPtTrackInConePerNCellCut(0), fhPtTrackInConePerNCellCutLargePtTrig(0),
227 fhConeSumPtClusterPerNCellCut(0), fhConeSumPtClusterPerNCellCutLargePtTrig(0),
228 fhConeSumPtTrackPerNCellCut(0), fhConeSumPtTrackPerNCellCutLargePtTrig(0),
229 fhPtClusterInConePerExoCut(0), fhPtClusterInConePerExoCutLargePtTrig(0),
230 fhPtTrackInConePerExoCut(0), fhPtTrackInConePerExoCutLargePtTrig(0),
231 fhConeSumPtClusterPerExoCut(0), fhConeSumPtClusterPerExoCutLargePtTrig(0),
232 fhConeSumPtTrackPerExoCut(0), fhConeSumPtTrackPerExoCutLargePtTrig(0),
234 fhConeSumPtClusterFECCorrPair(0), fhConeSumPtClusterFECCorrOdd(0),
235 fhConeSumPtClusterFECCorrPair2Max(0), fhConeSumPtClusterFECCorrOdd2Max(0),
236 fhConeSumPtClusterFECCorrPairHighCut(0), fhConeSumPtClusterFECCorrOddHighCut(0),
237 fhConeSumPtClusterFECCorrPair2MaxHighCut(0), fhConeSumPtClusterFECCorrOdd2MaxHighCut(0),
239 fhConeSumPtTrackTOFBC0(0), fhConeSumPtTrackTOFBCN(0), fhConeSumPtTrackTOFNo(0),
240 fhPtTrackInConeTOFBC0 (0), fhPtTrackInConeTOFBCN (0), fhPtTrackInConeTOFNo (0),
241 fhPhiTrackInCone(0), fhEtaTrackInCone(0), fhEtaPhiTrackInCone(0),
242 fhPhiTrackInConeTOFBC0(0), fhPhiTrackInConeTOFBCN(0), fhPhiTrackInConeTOFNo(0),
243 fhEtaTrackInConeTOFBC0(0), fhEtaTrackInConeTOFBCN(0), fhEtaTrackInConeTOFNo(0),
244 fhEtaPhiTrackInConeTOFBC0(0), fhEtaPhiTrackInConeTOFBCN(0), fhEtaPhiTrackInConeTOFNo(0),
245 fhTrackTOFInCone(0), fhTrackTOFInConeBC0(0), fhTrackTOFInConeExoTrigger(0),
247 fhConeSumPtTrackITSRefitOnSPDOn(0), fhConeSumPtTrackITSRefitOnSPDOff(0), fhConeSumPtTrackITSRefitOffSPDOff(0),
248 fhPtTrackInConeITSRefitOnSPDOn(0), fhPtTrackInConeITSRefitOnSPDOff(0) , fhPtTrackInConeITSRefitOffSPDOff(0),
249 fhPhiTrackInConeITSRefitOnSPDOn(0), fhPhiTrackInConeITSRefitOnSPDOff(0), fhPhiTrackInConeITSRefitOffSPDOff(0),
250 fhEtaTrackInConeITSRefitOnSPDOn(0), fhEtaTrackInConeITSRefitOnSPDOff(0), fhEtaTrackInConeITSRefitOffSPDOff(0),
251 fhEtaPhiTrackInConeITSRefitOnSPDOn(0),fhEtaPhiTrackInConeITSRefitOnSPDOff(0), fhEtaPhiTrackInConeITSRefitOffSPDOff(0),
252 fhConeSumPtTrackTOFBC0ITSRefitOnSPDOn(0), fhPtTrackInConeTOFBC0ITSRefitOnSPDOn(0),
253 fhPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0), fhEtaTrackInConeTOFBC0ITSRefitOnSPDOn(0),fhEtaPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0),
254 fhPerpConeSumPtITSRefitOnSPDOn (0), fhPtInPerpConeITSRefitOnSPDOn(0), fhEtaPhiInPerpConeITSRefitOnSPDOn(0),
255 fhPerpConeSumPtTOFBC0ITSRefitOnSPDOn (0), fhPtInPerpConeTOFBC0ITSRefitOnSPDOn (0), fhEtaPhiInPerpConeTOFBC0ITSRefitOnSPDOn(0)
259 for(
Int_t i = 0; i < 5 ; i++)
266 for(
Int_t j = 0; j < 5 ; j++)
301 for(
Int_t iso =0; iso < 2; iso++)
311 for(
Int_t i = 0; i < 5 ; i++)
330 for(
Int_t i = 0; i < 2 ; i++)
346 for(
Int_t imc = 0; imc < 4; imc++)
354 for(
Int_t i = 0; i < 2 ; i++)
370 for(
Int_t ieta = 0; ieta < 4; ieta++)
372 for(
Int_t iphi = 0; iphi < 3; iphi++)
383 for(
Int_t ism =0; ism < 20; ism++)
406 for(
Int_t i = 0 ; i < 7 ; i++)
415 for(
Int_t icase = 0; icase < 4; icase++)
417 for(
Int_t icase = 0; icase < 6; icase++)
445 Float_t ptTrig = pCandidate->Pt() ;
446 Float_t phiTrig = pCandidate->Phi();
447 Float_t etaTrig = pCandidate->Eta();
449 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
451 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
455 AliWarning(
"Cluster not available?");
460 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
461 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
514 Float_t phiTrig = pCandidate->Phi();
515 if(phiTrig<0) phiTrig += TMath::TwoPi();
517 Float_t etaTrig = pCandidate->Eta();
519 if(pCandidate->GetDetectorTag()==
kEMCAL)
521 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
524 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
526 if(!eGeom->CheckAbsCellId(absId))
return ;
529 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
531 Int_t imEta=-1, imPhi=-1;
532 Int_t ieta =-1, iphi =-1;
534 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
536 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
538 Int_t colTrig = ieta;
539 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
541 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
543 Int_t sqrSize = int(conesize/0.0143);
547 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
548 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
552 Int_t irowmin = rowTrig-sqrSize;
553 if ( irowmin < 0 ) irowmin = 0 ;
555 Int_t irowmax = rowTrig+sqrSize;
556 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
558 for(
Int_t irow = irowmin; irow <irowmax; irow++)
560 for(
Int_t icol = 0; icol < nTotalCols; icol++)
562 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
563 if(inSector==5)
continue;
566 if(icol < AliEMCALGeoParams::fgkEMCALCols)
568 inSupMod = 2*inSector + 1;
571 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
573 inSupMod = 2*inSector;
574 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
577 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
580 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
584 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
585 if(!eGeom->CheckAbsCellId(iabsId))
continue;
587 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
595 Int_t icolmin = colTrig-sqrSize;
596 if ( icolmin < 0 ) icolmin = 0;
598 Int_t icolmax = colTrig+sqrSize;
599 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
602 for(
Int_t icol = icolmin; icol < icolmax; icol++)
604 for(
Int_t irow = 0; irow < nTotalRows; irow++)
606 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
607 if ( inSector == 5 ) continue ;
613 if(icol < AliEMCALGeoParams::fgkEMCALCols)
616 inSupMod = 2*inSector + 1;
619 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
622 inSupMod = 2*inSector;
623 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
626 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
629 if(TMath::Abs(icol-colTrig) < sqrSize) {
632 if(TMath::Abs(irow-rowTrig) < sqrSize) {
636 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
637 if(!eGeom->CheckAbsCellId(iabsId))
639 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
643 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
653 Float_t ptTrig = pCandidate->Pt();
675 Float_t ptTrig = pCandidate->Pt() ;
676 Float_t phiTrig = pCandidate->Phi();
677 Float_t etaTrig = pCandidate->Eta();
682 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
684 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
688 AliWarning(
"Track not available?");
695 if ( pCandidate->GetDetectorTag() ==
kCTS )
698 Bool_t contained = kFALSE;
700 for(
Int_t i = 0; i < 4; i++)
702 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
705 if ( contained ) continue ;
721 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
723 phiBandPtSum+=track->Pt();
729 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
731 etaBandPtSum+=track->Pt();
738 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
739 Double_t dEta = etaTrig - track->Eta();
740 Double_t arg = dPhi*dPhi + dEta*dEta;
742 if(TMath::Sqrt(arg) < conesize)
744 sumptPerp+=track->Pt();
752 ULong_t status = track->GetStatus();
753 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
754 Int_t trackBC = track->GetTOFBunchCrossing(bz);
757 if(okTOF && trackBC == 0)
762 sumptPerpBC0+=track->Pt();
765 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
772 sumptPerpITSSPD+=track->Pt();
775 if(okTOF && trackBC == 0 && !bConstrained)
780 sumptPerpBC0ITSSPD+=track->Pt();
786 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
787 arg = dPhi*dPhi + dEta*dEta;
789 if(TMath::Sqrt(arg) < conesize)
791 sumptPerp+=track->Pt();
799 ULong_t status = track->GetStatus();
800 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
801 Int_t trackBC = track->GetTOFBunchCrossing(bz);
804 if(okTOF && trackBC == 0)
809 sumptPerpBC0+=track->Pt();
812 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
819 sumptPerpITSSPD+=track->Pt();
822 if(okTOF && trackBC == 0 && !bConstrained)
827 sumptPerpBC0ITSSPD+=track->Pt();
857 Float_t etaUEptsumCluster = 0 ;
858 Float_t phiUEptsumCluster = 0 ;
867 Float_t coneA = conesize*conesize*TMath::Pi();
868 Float_t ptTrig = pCandidate->Pt() ;
869 Float_t phiTrig = pCandidate->Phi();
870 Float_t etaTrig = pCandidate->Eta();
876 Float_t phiUEptsumTrackNorm = 0 ;
877 Float_t etaUEptsumTrackNorm = 0 ;
878 Float_t coneptsumTrackSubPhi = 0 ;
879 Float_t coneptsumTrackSubEta = 0 ;
880 Float_t coneptsumTrackSubPhiNorm = 0 ;
881 Float_t coneptsumTrackSubEtaNorm = 0 ;
893 Float_t correctConeSumTrack = 1;
894 Float_t correctConeSumTrackPhi = 1;
897 phiUEptsumTrack,etaUEptsumTrack,
898 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
899 correctConeSumTrack,correctConeSumTrackPhi);
901 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
902 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
914 if(coneptsumTrack > 0)
916 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
917 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
929 Float_t phiUEptsumClusterNorm = 0 ;
930 Float_t etaUEptsumClusterNorm = 0 ;
931 Float_t coneptsumClusterSubPhi = 0 ;
932 Float_t coneptsumClusterSubEta = 0 ;
933 Float_t coneptsumClusterSubPhiNorm = 0 ;
934 Float_t coneptsumClusterSubEtaNorm = 0 ;
935 Float_t phiUEptsumCellNorm = 0 ;
936 Float_t etaUEptsumCellNorm = 0 ;
937 Float_t coneptsumCellSubPhi = 0 ;
938 Float_t coneptsumCellSubEta = 0 ;
939 Float_t coneptsumCellSubPhiNorm = 0 ;
940 Float_t coneptsumCellSubEtaNorm = 0 ;
955 Float_t correctConeSumClusterEta = 1;
956 Float_t correctConeSumClusterPhi = 1;
959 phiUEptsumCluster,etaUEptsumCluster,
960 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
961 correctConeSumClusterEta,correctConeSumClusterPhi);
972 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
973 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
987 if(coneptsumCluster!=0)
989 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
990 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
1011 Float_t emcPhiSize = TMath::DegToRad()*100.;
1013 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
1014 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
1018 Float_t correctConeSumCellEta = 1;
1019 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
1021 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
1025 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
1028 Float_t correctConeSumCellPhi = 1;
1030 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
1031 (phiTrig-conesize < 80*TMath::DegToRad()))
1034 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
1035 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
1041 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
1045 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
1046 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1057 if ( coneptsumCell > 0.01 )
1059 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1060 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1070 sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1071 sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1098 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1099 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1123 coneptLeadCluster = 0;
1124 coneptsumCluster = 0;
1128 Float_t ptTrig = aodParticle->Pt();
1151 for(
Int_t icase = 0; icase < 4; icase++)
1153 for(
Int_t icase = 0; icase < 6; icase++)
1158 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1159 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1170 Float_t coneptsumClusterPerMinCut[20];
1171 Float_t coneptsumClusterPerMaxCut[20];
1172 Float_t coneptsumClusterPerRCut [10];
1174 Float_t coneptsumClusterPerNCellCut[20];
1175 Float_t coneptsumClusterPerExoCut [20];
1181 coneptsumClusterPerMinCut[icut] = 0;
1182 coneptsumClusterPerMaxCut[icut] = 0;
1190 coneptsumClusterPerRCut[icut] = 0;
1198 coneptsumClusterPerNCellCut[icut] = 0;
1203 coneptsumClusterPerExoCut[icut] = 0;
1211 Float_t ptSumPairFECHighCut = 0;
1212 Float_t ptSumPairFEC2MaxHighCut = 0;
1213 Float_t ptSumOddFECHighCut = 0;
1214 Float_t ptSumOddFEC2MaxHighCut = 0;
1215 Float_t ptSumFECMaxCorr[] = {0,0,0,0};
1216 Float_t ptSumFECCorr[] = {0,0,0,0,0,0};
1218 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1220 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1247 coneptsumCluster+=ptcone;
1248 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1254 if ( ptcone >
fMinPtCutInCone[icut] ) coneptsumClusterPerMinCut[icut]+=ptcone;
1255 if ( ptcone <
fMaxPtCutInCone[icut] ) coneptsumClusterPerMaxCut[icut]+=ptcone;
1266 coneptsumClusterPerRCut[icut]+=ptcone;
1279 coneptsumClusterPerNCellCut[icut]+=ptcone;
1289 coneptsumClusterPerExoCut[icut]+=ptcone;
1301 Int_t absIdMax = -1;
1302 Int_t absIdMax2 = -1;
1305 Int_t ncells = calo->GetNCells();
1306 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1308 Int_t absId = calo->GetCellsAbsId()[ipos];
1315 absIdMax2 = absIdMax;
1320 else if(eCell > emax2)
1330 Int_t absIdCorr[] = {-1,-1,-1,-1};
1336 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1338 Int_t absId = calo->GetCellsAbsId()[ipos];
1344 if( absId == absIdMax || weight < 0.01 )
continue;
1346 for(
Int_t id = 0;
id < 4;
id++)
1348 if ( absId == absIdCorr[
id] )
1351 Int_t diff = TMath::Abs(absId-absIdMax);
1352 if ( diff == 1 ) near = kTRUE;
1353 else if( diff >= 8 ) far = kTRUE;
1360 if ( nCorr == 0 ) index = 0;
1363 if ( near && !far ) index = 1;
1364 else if (!near && far ) index = 2;
1365 else if ( near && far ) index = 3;
1368 if ( index >= 0 ) ptSumFECMaxCorr[index]+=ptcone;
1373 Int_t absIdCorrSecondary [] = {-1,-1,-1,-1};
1375 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1377 Int_t absId = calo->GetCellsAbsId()[ipos];
1382 if( absId == absIdCorr[0] || absId == absIdCorr[1] ||
1383 absId == absIdCorr[2] || absId == absIdCorr[3] || weight < 0.01 )
continue;
1387 if(!correlOK)
continue;
1389 for (
Int_t ipos2 = ipos+1; ipos2 < ncells; ipos2++)
1391 Int_t absId2 = calo->GetCellsAbsId()[ipos2];
1396 if( absId2 == absIdCorr[0] || absId2 == absIdCorr[1] || absId2 == absId ||
1397 absId2 == absIdCorr[2] || absId2 == absIdCorr[3] || weight < 0.01 )
continue;
1399 for(
Int_t id = 0;
id < 4;
id++)
1401 if ( absId2 == absIdCorrSecondary[
id] )
1408 if ( counter == 1 ) distCase = 1;
1409 else if ( counter >= 2 ) distCase = 2;
1411 if ( nCorr > 0 ) distCase+=3;
1413 ptSumFECCorr[distCase]+=ptcone;
1421 for(
Int_t id = 0;
id < 4;
id++)
1423 if ( absIdMax == absIdCorr[
id]) idMax = id;
1424 else if( absIdMax2 == absIdCorr[
id]) idMax2 = id;
1427 if ( idMax==0 || idMax == 3 ) ptSumPairFEC+=ptcone;
1428 else ptSumOddFEC +=ptcone;
1432 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2Max+=ptcone;
1433 else ptSumOddFEC2Max +=ptcone;
1438 if ( idMax==0 || idMax == 3 ) ptSumPairFECHighCut+=ptcone;
1439 else ptSumOddFECHighCut +=ptcone;
1443 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2MaxHighCut+=ptcone;
1444 else ptSumOddFEC2MaxHighCut +=ptcone;
1455 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1456 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1506 for(
Int_t icase = 0; icase < 4; icase++)
1508 for(
Int_t icase = 0; icase < 6; icase++)
1526 Float_t ptTrig = aodParticle->Pt();
1527 Float_t phiTrig = aodParticle->Phi();
1528 if(phiTrig<0) phiTrig += TMath::TwoPi();
1529 Float_t etaTrig = aodParticle->Eta();
1531 if(aodParticle->GetDetectorTag()==
kEMCAL)
1533 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1536 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1538 if(!eGeom->CheckAbsCellId(absId))
return ;
1541 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1543 Int_t imEta=-1, imPhi=-1;
1544 Int_t ieta =-1, iphi =-1;
1546 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1548 Int_t iEta=-1, iPhi=-1;
1549 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1551 Int_t colTrig = iEta;
1552 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1553 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1555 Int_t sqrSize = int(conesize/0.0143);
1560 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1562 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1564 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1565 if(inSector==5)
continue;
1567 Int_t inSupMod = -1;
1569 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1571 inSupMod = 2*inSector + 1;
1574 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1576 inSupMod = 2*inSector;
1577 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1580 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1582 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1583 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1586 coneptsumCell += cells->GetCellAmplitude(iabsId);
1602 coneptLeadTrack = 0;
1606 Float_t ptTrig = aodParticle->Pt();
1628 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1629 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1638 Float_t coneptsumTrackPerMinCut[20];
1639 Float_t coneptsumTrackPerMaxCut[20];
1640 Float_t coneptsumTrackPerEtaCut[10];
1641 Float_t coneptsumTrackPerNCellCut[20];
1642 Float_t coneptsumTrackPerExoCut[20];
1643 Float_t coneptsumTrackPerRCut [10];
1644 Float_t coneptsumTrackTOFBC0 = 0;
1645 Float_t coneptsumTrackTOFBCN = 0;
1646 Float_t coneptsumTrackTOFNo = 0;
1647 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1648 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1649 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1650 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1656 coneptsumTrackPerMinCut[icut] = 0;
1657 coneptsumTrackPerMaxCut[icut] = 0;
1665 coneptsumTrackPerEtaCut[icut] = 0;
1673 coneptsumTrackPerRCut[icut] = 0;
1681 coneptsumTrackPerNCellCut[icut] = 0;
1686 coneptsumTrackPerExoCut[icut] = 0;
1690 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1692 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1694 pTtrack = track->Pt();
1701 Int_t partInConeCharge = 0, partInConePDG = 0;
1702 Bool_t physPrimary = kFALSE;
1703 Int_t trackLabel = TMath::Abs(track->GetLabel());
1707 TParticle * mcpartESD =
static_cast<TParticle*
>(
GetMCStack()->Particle(trackLabel));
1710 partInConePDG = mcpartESD->GetPdgCode();
1711 partInConeCharge = TMath::Abs((
Int_t) TDatabasePDG::Instance()->GetParticle(partInConePDG)->Charge());
1712 physPrimary = mcpartESD->IsPrimary();
1720 AliAODMCParticle * mcpartAOD = (AliAODMCParticle *) mcparticles->At(trackLabel);
1723 partInConeCharge = TMath::Abs(mcpartAOD->Charge());
1724 partInConePDG = mcpartAOD->GetPdgCode();
1725 physPrimary = mcpartAOD->IsPhysicalPrimary();
1730 if ( partInConeCharge > 0 && TMath::Abs(partInConePDG) != 11 )
1733 if ( TMath::Abs(partInConePDG) == 211 ) mcChTag = 0;
1734 else if ( TMath::Abs(partInConePDG) == 321 ) mcChTag = 1;
1735 else if ( TMath::Abs(partInConePDG) == 2212 ) mcChTag = 2;
1751 coneptsumTrack+=pTtrack;
1752 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1758 if ( pTtrack >
fMinPtCutInCone[icut] ) coneptsumTrackPerMinCut[icut]+=pTtrack;
1759 if ( pTtrack <
fMaxPtCutInCone[icut] ) coneptsumTrackPerMaxCut[icut]+=pTtrack;
1767 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1778 coneptsumTrackPerRCut[icut]+=pTtrack;
1791 coneptsumTrackPerNCellCut[icut]+=pTtrack;
1801 coneptsumTrackPerExoCut[icut]+=pTtrack;
1812 phitrack = track->Phi();
1813 etatrack = track->Eta();
1820 ULong_t status = track->GetStatus();
1821 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1822 trackBC = track->GetTOFBunchCrossing(bz);
1823 Double32_t tof = track->GetTOFsignal()*1e-3;
1842 coneptsumTrackTOFBC0 += pTtrack;
1850 coneptsumTrackTOFBCN += pTtrack;
1859 coneptsumTrackTOFNo += pTtrack;
1863 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1864 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1872 coneptsumTrackITSRefitOnSPDOff += pTtrack;
1880 coneptsumTrackITSRefitOffSPDOff += pTtrack;
1889 coneptsumTrackITSRefitOnSPDOn += pTtrack;
1896 if(okTOF && trackBC == 0 && !bConstrained)
1902 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
1945 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
1946 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
2005 if(energy < 8)
return;
2011 Int_t nVerticesSPD = -1;
2012 Int_t nVerticesTracks = -1;
2016 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
2017 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
2022 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
2023 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
2035 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
2039 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
2040 ncont=pv->GetNContributors();
2041 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
2043 diamZ = esdEv->GetDiamondZ();
2047 AliAODVertex *pv=aodEv->GetVertex(iVert);
2048 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
2049 ncont=pv->GetNContributors();
2050 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
2052 diamZ = aodEv->GetDiamondZ();
2055 Double_t distZ = TMath::Abs(z2-z1);
2056 diamZ = TMath::Abs(z2-diamZ);
2074 Int_t nMaxima = pCandidate->GetNLM();
2075 Int_t mcTag = pCandidate->GetTag() ;
2076 Bool_t isolated = pCandidate->IsIsolated();
2078 Float_t m02 = pCandidate->GetM02() ;
2080 Float_t pt = pCandidate->Pt();
2081 Float_t eta = pCandidate->Eta();
2082 Float_t phi = pCandidate->Phi();
2083 if(phi<0) phi+= TMath::TwoPi();
2089 Int_t decayTag = pCandidate->DecayTag();
2090 if(decayTag < 0) decayTag = 0;
2127 Int_t ptsumBin = -1;
2128 Int_t leadptBin = -1;
2130 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
2140 if ( leadptBin >= 0 )
2149 if ( leadptBin == 0 )
2150 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
2154 if ( ptsumBin >= 0 )
2166 Int_t decayTag = pCandidate->DecayTag();
2167 if(decayTag < 0) decayTag = 0;
2229 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
2251 Int_t ptTrigBin = -1;
2260 if ( ptTrigBin >= 0 )
2277 Int_t decayTag = pCandidate->DecayTag();
2278 if(decayTag < 0) decayTag = 0;
2383 Int_t noverlaps = 0;
2387 Int_t overpdg[nlabels];
2388 Int_t overlab[nlabels];
2418 if ( noverlaps == 1 )
2445 else if (noverlaps > 1 )
2531 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2537 Float_t dEdx = track->GetTPCsignal();
2577 const Int_t buffersize = 255;
2578 char onePar[buffersize] ;
2580 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2584 snprintf(onePar, buffersize,
"Isolation Cand. Detector: %s;",
fIsoDetectorString.Data()) ;
2588 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2590 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2592 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2594 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2599 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2601 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2606 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2611 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2616 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2621 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2632 return new TObjString(parList) ;
2642 outputContainer->SetName(
"IsolatedParticleHistos") ;
2691 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
2694 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
2699 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
2702 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
2706 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
2709 TString sParticle =
", x^{0,#pm}";
2713 TString parTitle = Form(
"#it{R} = %2.2f%s%s",
GetIsolationCut()->GetConeSize(), sThreshold.Data(),sParticle.Data());
2715 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
2718 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
2719 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
2720 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
2721 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
2723 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
2724 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
2725 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
2726 "OtherDecay",
"Electron",
"Hadron"} ;
2729 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
2730 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
2731 "#pi^{0}" ,
"#eta"} ;
2733 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
2734 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
2740 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2741 nptbins,ptmin,ptmax);
2742 fhENoIso->SetYTitle(
"#it{counts}");
2743 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
2747 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2748 nptbins,ptmin,ptmax);
2750 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2756 (
"hENoIsoExoTrigger",
2757 Form(
"Number of not isolated particles vs E, %s, exot>0.97",parTitle.Data()),
2758 nptbins,ptmin,ptmax);
2764 (
"hPtNoIsoExoTrigger",
2765 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
2766 nptbins,ptmin,ptmax);
2773 Form(
"Number of not isolated leading particles #eta vs #varphi, %s",parTitle.Data()),
2774 netabins,etamin,etamax,nphibins,phimin,phimax);
2785 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
2786 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2787 nptbins,ptmin,ptmax);
2789 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2792 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
2793 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2794 nptbins,ptmin,ptmax);
2795 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
2796 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2800 Form(
"#varphi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2801 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2803 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2807 Form(
"#eta vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2808 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2810 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2818 TString isoName [] = {
"NoIso",
"Iso"};
2819 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2823 for(
Int_t iso = 0; iso < 2; iso++)
2827 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2828 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2829 nptbins,ptmin,ptmax);
2830 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2831 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2832 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2836 Form(
"Number of %s leading Pi0 decay particles #eta vs #varphi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2837 netabins,etamin,etamax,nphibins,phimin,phimax);
2845 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2846 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2847 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2858 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2859 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2860 nptbins,ptmin,ptmax);
2861 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2862 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2863 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2872 TString isoName [] = {
"NoIso",
"Iso"};
2873 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2876 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2877 nptbins,ptmin,ptmax);
2878 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2879 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2880 outputContainer->Add(
fhEIso) ;
2883 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2884 nptbins,ptmin,ptmax);
2885 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
2886 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2887 outputContainer->Add(
fhPtIso) ;
2893 Form(
"Number of isolated particles vs E, %s, exot>0.97",parTitle.Data()),
2894 nptbins,ptmin,ptmax);
2900 (
"hPtIsoExoTrigger",
2901 Form(
"Number of isolated particles vs #it{p}_{T}, %s, exot>0.97",parTitle.Data()),
2902 nptbins,ptmin,ptmax);
2909 Form(
"Number of isolated particles vs #varphi, %s",parTitle.Data()),
2910 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2911 fhPhiIso->SetYTitle(
"#varphi (rad)");
2912 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2916 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
2917 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2919 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2923 Form(
"Number of isolated particles #eta vs #varphi, %s",parTitle.Data()),
2924 netabins,etamin,etamax,nphibins,phimin,phimax);
2932 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2933 nptbins,ptmin,ptmax, 100,0,100);
2939 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2940 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
2949 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2950 nptbins,ptmin,ptmax,10,0,10);
2955 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2956 nptbins,ptmin,ptmax,10,0,10);
2963 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2964 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2965 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
2966 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2971 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2972 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2973 fhConeSumPt->SetYTitle(
"#Sigma #it{p}_{T} (GeV/#it{c})");
2974 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2980 (
"hConePtSumExoTrigger",
2981 Form(
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
2982 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2997 Form(
"Trigger #eta vs #varphi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2998 netabins,etamin,etamax,nphibins,phimin,phimax);
3005 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for #it{R} = %2.2f",r),
3006 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3007 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
3008 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3014 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f, exotic trigger",r),
3015 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3052 (Form(
"hPtLeadCone_Bin%d",ibin),
3053 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3060 (Form(
"hSumPtCone_Bin%d",ibin),
3061 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3074 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3075 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3082 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3083 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3097 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3098 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3105 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3106 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3117 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
3118 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
3125 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
3126 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
3138 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3139 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
3146 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3147 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
3187 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
3188 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3195 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
3196 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3209 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3210 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3217 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
3218 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3232 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3233 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3240 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3241 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3252 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
3253 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
3260 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
3261 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
3273 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
3274 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
3281 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
3282 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
3296 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3297 100,0,100,nptinconebins,ptinconemin,ptinconemax);
3307 (
"hConePtSumCluster",
3308 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3309 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3317 (
"hConePtSumClusterExoTrigger",
3318 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
3319 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3328 (
"hConePtSumClusterPerMinPtCut",
3329 Form(
"Cluster #Sigma #it{p}_{T}, different min #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3338 (
"hConePtSumClusterPerMinPtCutLargePtTrig",
3339 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),
3348 (
"hConePtSumClusterPerMaxPtCut",
3349 Form(
"Cluster #Sigma #it{p}_{T}, different max #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
3358 (
"hConePtSumClusterPerMaxPtCutLargePtTrig",
3359 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),
3371 (
"hConePtSumClusterPerRCut",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts",
3380 (
"hConePtSumClusterPerRCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3389 (
"hPtClusterInConePerRCut",
"Cluster #it{p}_{T}, different #it{R} cuts",
3398 (
"hPtClusterInConePerRCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
3410 (
"hConePtSumClusterPerNCellCut",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
3419 (
"hConePtSumClusterPerNCellCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3428 (
"hPtClusterInConePerNCellCut",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts",
3437 (
"hPtClusterInConePerNCellCutLargePtTrig",
"Cluster #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
3447 (
"hConePtSumClusterPerExoCut",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts",
3456 (
"hConePtSumClusterPerExoCutLargePtTrig",
"Cluster #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
3465 (
"hPtClusterInConePerExoCut",
"Cluster #it{p}_{T}, different exoticity cuts",
3474 (
"hPtClusterInConePerExoCutLargePtTrig",
"Cluster #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
3487 {
"No correlation with main cell",
"Main cell corr. with near cell",
3488 "Main cell corr. with far cell",
"Main cell corr. with near and far"};
3491 {
"No main cell corr., no other inside cluster" ,
"No main cell corr., 1 cell correl in cluster" ,
"No main cell corr., >=2 cell correl in cluster",
3492 "Main cell corr. & no other inside cluster" ,
"Main cell corr. & 1 cell corr. inside cluster",
"Main cell corr. & >=2 cell corr. inside cluster"};
3494 for(
Int_t icase = 0; icase < 4; icase++)
3497 (Form(
"hConePtSumClusterFECCorrMax_Case%d",icase),
3498 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, %s",r, titleMaxFEC[icase].
Data()),
3499 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3505 for(
Int_t icase = 0; icase < 6; icase++)
3508 (Form(
"hConePtSumClusterFECCorr_Case%d",icase),
3509 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, %s",r, titleFEC[icase].
Data()),
3510 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3517 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC",r),
3518 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3524 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, (n,n+1) max cells",r),
3525 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3531 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, #it{E}_{T} > 0.7 GeV",r),
3532 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3538 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, (n,n+1) max cells, #it{E}_{T}> 0.7 GeV",r),
3539 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3545 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC",r),
3546 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3552 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, (n,n+1) max cells",r),
3553 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3559 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, #it{E}_{T} > 0.7 GeV",r),
3560 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3566 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, (n,n+1) max cells, #it{E}_{T}> 0.7 GeV",r),
3567 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3575 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
3576 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3585 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3586 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3595 "#Sigma cluster #it{p}_{T} in UE Eta Band",
3596 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3602 "#Sigma cluster #it{p}_{T} UE Phi Band",
3603 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3609 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} in UE Eta Band",
3610 netabins,etamin,etamax,nphibins,phimin,phimax);
3617 "Trigger #eta vs #varphi, #Sigma cluster #it{p}_{T} UE Phi Band",
3618 netabins,etamin,etamax,nphibins,phimin,phimax);
3627 "#Sigma cell #it{p}_{T} in UE Eta Band",
3628 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3634 "#Sigma cell #it{p}_{T} UE Phi Band",
3635 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3641 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} in UE Eta Band",
3642 netabins,etamin,etamax,nphibins,phimin,phimax);
3649 "Trigger #eta vs #varphi, #Sigma cell #it{p}_{T} UE Phi Band",
3650 netabins,etamin,etamax,nphibins,phimin,phimax);
3658 Form(
"#eta vs #varphi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
3659 netabins,-1,1,nphibins,0,TMath::TwoPi());
3665 Form(
"#eta vs #varphi of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
3666 netabins,-1,1,nphibins,0,TMath::TwoPi());
3673 Form(
"#it{p}_{T} of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
3674 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3680 Form(
"#it{p}_{T} of clusters in #varphi band isolation cone for #it{R} = %2.2f",r),
3681 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3687 Form(
"#eta vs #varphi of clusters in cone for #it{R} = %2.2f",r),
3688 netabins,-1,1,nphibins,0,TMath::TwoPi());
3694 Form(
"#eta vs #varphi of all clusters"),
3695 netabins,-1,1,nphibins,0,TMath::TwoPi());
3703 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
3704 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3712 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f, exotic trigger",r),
3713 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3722 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
3723 nptbins,ptmin,ptmax,1000,0,50);
3729 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
3736 Form(
"#col vs #row of cells in #varphi band isolation cone for #it{R} = %2.2f",r),
3746 Form(
"Clusters #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
3747 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3753 Form(
"Clusters #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
3754 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3760 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
3761 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3767 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
3768 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3774 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),
3775 netabins,etamin,etamax,nphibins,phimin,phimax);
3782 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),
3783 netabins,etamin,etamax,nphibins,phimin,phimax);
3792 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
3793 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3799 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
3800 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3806 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),
3807 netabins,etamin,etamax,nphibins,phimin,phimax);
3814 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),
3815 netabins,etamin,etamax,nphibins,phimin,phimax);
3823 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
3824 nptbins,ptmin,ptmax,100,0,1);
3830 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#varphi ",r),
3831 netabins,etamin,etamax,nphibins,phimin,phimax);
3838 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance",r),
3839 nptbins,ptmin,ptmax,100,0,1);
3845 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #varphi acceptance, in trigger #eta-#varphi ",r),
3846 netabins,etamin,etamax,nphibins,phimin,phimax);
3853 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),
3854 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3860 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),
3861 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3867 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),
3868 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3874 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),
3875 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3881 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #eta band for cluster (before normalization), R=%2.2f",r),
3882 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
3888 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in #varphi band for cluster (before normalization), R=%2.2f",r),
3889 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
3897 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
3898 nptbins,ptmin,ptmax,100,0,1);
3904 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#varphi ",r),
3905 netabins,etamin,etamax,nphibins,phimin,phimax);
3912 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance",r),
3913 nptbins,ptmin,ptmax,100,0,1);
3919 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #varphi acceptance, in trigger #eta-#varphi ",r),
3920 netabins,etamin,etamax,nphibins,phimin,phimax);
3928 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),
3929 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3935 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),
3936 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3942 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),
3943 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3949 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),
3950 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3962 (
"hConeLeadPtTrack",
3963 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
3964 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3971 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
3972 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3979 TString mcChPartName[] = {
"Pion",
"Kaon",
"Proton",
"Other"};
3980 for(
Int_t imc = 0; imc < 4; imc++)
3983 (Form(
"hPtTrackInCone_Primary_%s",mcChPartName[imc].
Data()),
3984 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
3985 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3991 (Form(
"hPtTrackInCone_Secondary_%s",mcChPartName[imc].
Data()),
3992 Form(
"reconstructed #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
3993 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3999 (Form(
"hPtTrackInCone_Gener_Primary_%s",mcChPartName[imc].
Data()),
4000 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4001 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4007 (Form(
"hPtTrackInCone_Gener_Secondary_%s",mcChPartName[imc].
Data()),
4008 Form(
"generated #it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, primary MC %s",r,mcChPartName[imc].
Data()),
4009 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4020 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, exotic trigger",r),
4021 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4030 (
"hConePtSumTrackPerMinPtCut",
4031 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4040 (
"hConePtSumTrackPerMinPtCutLargePtTrig",
4041 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4050 (
"hConePtSumTrackPerMaxPtCut",
4051 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f",r),
4060 (
"hConePtSumTrackPerMaxPtCutLargePtTrig",
4061 Form(
"Track #Sigma #it{p}_{T}, different #it{p}_{T} cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4073 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f",r),
4082 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
4094 (
"hConePtSumTrackPerRCut",
"Track #Sigma #it{p}_{T}, different #it{R} cuts",
4103 (
"hConePtSumTrackPerRCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4112 (
"hPtTrackInConePerRCut",
"Track #it{p}_{T}, different #it{R} cuts",
4121 (
"hPtTrackInConePerRCutLargePtTrig",
"Track #it{p}_{T}, different #it{R} cuts, #it{p}_{T}^{trig} > 10 GeV",
4133 (
"hConePtSumTrackPerNCellCut",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts",
4142 (
"hConePtSumTrackPerNCellCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4151 (
"hPtTrackInConePerNCellCut",
"Track #it{p}_{T}, different #it{N}_{cell} cuts",
4160 (
"hPtTrackInConePerNCellCutLargePtTrig",
"Track #it{p}_{T}, different #it{N}_{cell} cuts, #it{p}_{T}^{trig} > 10 GeV",
4169 (
"hConePtSumTrackPerExoCut",
"Track #Sigma #it{p}_{T}, different exoticity cuts",
4178 (
"hConePtSumTrackPerExoCutLargePtTrig",
"Track #Sigma #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4187 (
"hPtTrackInConePerExoCut",
"Track #it{p}_{T}, different exoticity cuts",
4196 (
"hPtTrackInConePerExoCutLargePtTrig",
"Track #it{p}_{T}, different exoticity cuts, #it{p}_{T}^{trig} > 10 GeV",
4208 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4209 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4217 (
"hConePtSumTrackExoTrigger",
4218 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, exo trigger",r),
4219 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4227 Int_t ntofbins = 1000;
4228 Int_t mintof = -500;
4232 (
"hTrackTOFInCone",
"TOF signal vs track #it{p}_{T}",
4233 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4241 (
"hTrackTOFInConeExoTrigger",
"TOF signal vs track #it{p}_{T}, exoticity > 0.97",
4242 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4249 (
"hTrackTOFInConeBC0",
"TOF signal vs track #it{p}_{T}, BC=0",
4250 nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
4256 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
4257 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4263 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f",r),
4264 netabins,-1,1,nphibins,0,TMath::TwoPi());
4270 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4271 nptbins,ptmin,ptmax,netabins,-1,1);
4277 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f",r),
4278 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4289 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0",r),
4290 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4296 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC!=0",r),
4297 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4303 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no TOF",r),
4304 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4310 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
4311 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4317 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
4318 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4324 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no TOF",r),
4325 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4332 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4333 netabins,-1,1,nphibins,0,TMath::TwoPi());
4339 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4340 netabins,-1,1,nphibins,0,TMath::TwoPi());
4346 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4347 netabins,-1,1,nphibins,0,TMath::TwoPi());
4353 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4354 nptbins,ptmin,ptmax,netabins,-1,1);
4360 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4361 nptbins,ptmin,ptmax,netabins,-1,1);
4367 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4368 nptbins,ptmin,ptmax,netabins,-1,1);
4374 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
4375 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4381 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
4382 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4388 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no TOF",r),
4389 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4398 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD On",r),
4399 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4405 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD Off",r),
4406 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4412 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no ITS Refit SPD Off",r),
4413 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4419 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
4420 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4426 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
4427 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4433 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4434 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4441 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4442 netabins,-1,1,nphibins,0,TMath::TwoPi());
4448 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4449 netabins,-1,1,nphibins,0,TMath::TwoPi());
4455 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4456 netabins,-1,1,nphibins,0,TMath::TwoPi());
4462 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4463 nptbins,ptmin,ptmax,netabins,-1,1);
4469 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4470 nptbins,ptmin,ptmax,netabins,-1,1);
4476 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4477 nptbins,ptmin,ptmax,netabins,-1,1);
4483 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
4484 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4490 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
4491 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4497 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
4498 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4507 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0, track ITS Refit SPD On",r),
4508 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4514 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0, track ITS Refit SPD On",r),
4515 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4521 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
4522 netabins,-1,1,nphibins,0,TMath::TwoPi());
4528 Form(
"#eta vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
4529 nptbins,ptmin,ptmax,netabins,-1,1);
4535 Form(
"#varphi vs #it{p}_{T} of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
4536 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
4545 "#Sigma track #it{p}_{T} in UE Eta Band",
4546 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4552 "#Sigma track #it{p}_{T} in UE Phi Band",
4553 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4560 "Trigger #eta vs #varphi, #Sigma track #it{p}_{T} in UE Eta Band",
4561 netabins,etamin,etamax,nphibins,phimin,phimax);
4568 "Trigger #eta vs #varphi, #Sigma track #it{p}_{T} in UE Phi Band",
4569 netabins,etamin,etamax,nphibins,phimin,phimax);
4576 Form(
"#eta vs #varphi of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
4577 netabins,-1,1,nphibins,0,TMath::TwoPi());
4583 Form(
"#eta vs #varphi of tracks in #varphi band isolation cone for #it{R} = %2.2f",r),
4584 netabins,-1,1,nphibins,0,TMath::TwoPi());
4590 Form(
"#it{p}_{T} of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
4591 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4597 Form(
"#eta vs #varphi of tracks in #varphi band isolation cone for #it{R} = %2.2f",r),
4598 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4604 Form(
"Tracks #Sigma #it{p}_{T} in normalized #eta band, #it{R} = %2.2f",r),
4605 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4611 Form(
"Tracks #Sigma #it{p}_{T} in normalized #varphi band, #it{R} = %2.2f",r),
4612 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4619 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4620 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4626 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from #varphi band in the isolation cone for #it{R} = %2.2f",r),
4627 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4633 Form(
"Trigger #eta vs #varphi, Tracks #Sigma #it{p}_{T} after bkg subtraction from #eta band in the isolation cone for #it{R} = %2.2f",r),
4634 netabins,etamin,etamax,nphibins,phimin,phimax);