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), fFillEMCALRegionSSHistograms(0),
59 fFillUEBandSubtractHistograms(1), fFillCellHistograms(0),
60 fFillOverlapHistograms(0), fStudyFECCorrelation(0), fStudyTracksInCone(0),
61 fFillTaggedDecayHistograms(0), fNDecayBits(0),
62 fDecayBits(), fDecayTagsM02Cut(0),
63 fFillNLMHistograms(0),
64 fLeadingOnly(0), fCheckLeadingWithNeutralClusters(0),
65 fSelectPrimariesInCone(0), fMakePrimaryPi0DecayStudy(0),
66 fFillBackgroundBinHistograms(0), fNBkgBin(0),
67 fFillPtTrigBinHistograms(0), fNPtTrigBin(0),
68 fMinCellsAngleOverlap(0),
70 fNCones(0), fNPtThresFrac(0),
71 fConeSizes(), fPtThresholds(),
72 fPtFractions(), fSumPtThresholds(),
73 fStudyPtCutInCone(0), fNPtCutsInCone(0), fPtCutInCone(),
74 fStudyEtaCutInCone(0), fNEtaCutsInCone(0), fEtaCutInCone(),
75 fMomentum(), fMomIso(),
76 fMomDaugh1(), fMomDaugh2(),
77 fTrackVector(), fProdVertex(),
78 fCluster(0), fClustersArr(0),
80 fhEIso(0), fhPtIso(0),
81 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
83 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
85 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
87 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
88 fhPtTrackInConeOtherBCPileUpSPD(0), fhPtTrackInConeVtxBC0(0),
89 fhPtTrackInConeBC0PileUpSPD(0),
90 fhPtInConePileUp(), fhPtInConeCent(0),
91 fhPerpConeSumPt(0), fhPerpConeSumPtTOFBC0(0),
92 fhPtInPerpCone(0), fhPtInPerpConeTOFBC0(0),
93 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
94 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
95 fhEtaPhiInPerpCone(0), fhEtaPhiInPerpConeTOFBC0(0),
96 fhEtaBandCluster(0), fhPhiBandCluster(0),
97 fhEtaBandTrack(0), fhPhiBandTrack(0),
98 fhEtaBandCell(0), fhPhiBandCell(0),
99 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
100 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
101 fhConeSumPt(0), fhPtLambda0Eiso(0), fhConeSumPtCellTrack(0),
102 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
103 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
104 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
105 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
106 fhConeSumPtTrigEtaPhi(0),
107 fhConeSumPtCellTrackTrigEtaPhi(0),
108 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
109 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
110 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
111 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
112 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
113 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
114 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
115 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
116 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
117 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
118 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
119 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
120 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
121 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
122 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
123 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
124 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
125 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
126 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
127 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
128 fhConeSumPtCellvsTrack(0),
129 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
130 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
131 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
132 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
133 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
134 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
135 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
136 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
137 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
138 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
139 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
140 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
141 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
142 fhPtPrimMCPi0DecayPairOutOfCone(0),
143 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
144 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
145 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
146 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
147 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
148 fhPtPrimMCPi0DecayPairNoOverlap(0),
149 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
150 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
151 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
152 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
153 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
154 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
155 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
156 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
157 fhPtPrimMCEtaDecayPairOutOfCone(0),
158 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
159 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
160 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
161 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
162 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
163 fhPtPrimMCEtaDecayPairNoOverlap(0),
164 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
165 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
166 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
167 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
168 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
169 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
170 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
171 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
172 fhPtLeadConeBin(0), fhSumPtConeBin(0),
173 fhPtLeadConeBinMC(0), fhSumPtConeBinMC(0),
174 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
175 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
176 fhPtLeadConeBinLambda0MC(0), fhSumPtConeBinLambda0MC(0),
177 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
178 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
179 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
180 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
181 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
184 fhELambda0LocMax1(), fhELambda1LocMax1(),
185 fhELambda0LocMax2(), fhELambda1LocMax2(),
186 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
188 fhEIsoPileUp(), fhPtIsoPileUp(),
189 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
190 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
191 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
192 fhTimeNPileUpVertContributors(0),
193 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
194 fhConeSumPtClusterPerPtCut(0), fhConeSumPtClusterPerPtCutLargePtTrig(0),
195 fhConeSumPtTrackPerPtCut(0), fhConeSumPtTrackPerPtCutLargePtTrig(0),
196 fhConeSumPtTrackPerEtaCut(0), fhConeSumPtTrackPerEtaCutLargePtTrig(0),
197 fhConeSumPtClusterFECCorrPair(0), fhConeSumPtClusterFECCorrOdd(0),
198 fhConeSumPtClusterFECCorrPair2Max(0), fhConeSumPtClusterFECCorrOdd2Max(0),
199 fhConeSumPtClusterFECCorrPairHighCut(0), fhConeSumPtClusterFECCorrOddHighCut(0),
200 fhConeSumPtClusterFECCorrPair2MaxHighCut(0), fhConeSumPtClusterFECCorrOdd2MaxHighCut(0),
201 fhConeSumPtTrackTOFBC0(0), fhConeSumPtTrackTOFBCN(0), fhConeSumPtTrackTOFNo(0),
202 fhPtTrackInConeTOFBC0 (0), fhPtTrackInConeTOFBCN (0), fhPtTrackInConeTOFNo (0),
203 fhPhiTrackInCone(0), fhEtaTrackInCone(0), fhEtaPhiTrackInCone(0),
204 fhPhiTrackInConeTOFBC0(0), fhPhiTrackInConeTOFBCN(0), fhPhiTrackInConeTOFNo(0),
205 fhEtaTrackInConeTOFBC0(0), fhEtaTrackInConeTOFBCN(0), fhEtaTrackInConeTOFNo(0),
206 fhEtaPhiTrackInConeTOFBC0(0), fhEtaPhiTrackInConeTOFBCN(0), fhEtaPhiTrackInConeTOFNo(0),
207 fhTrackTOFInCone(0), fhTrackTOFInConeBC0(0),
208 fhConeSumPtTrackITSRefitOnSPDOn(0), fhConeSumPtTrackITSRefitOnSPDOff(0), fhConeSumPtTrackITSRefitOffSPDOff(0),
209 fhPtTrackInConeITSRefitOnSPDOn(0), fhPtTrackInConeITSRefitOnSPDOff(0) , fhPtTrackInConeITSRefitOffSPDOff(0),
210 fhPhiTrackInConeITSRefitOnSPDOn(0), fhPhiTrackInConeITSRefitOnSPDOff(0), fhPhiTrackInConeITSRefitOffSPDOff(0),
211 fhEtaTrackInConeITSRefitOnSPDOn(0), fhEtaTrackInConeITSRefitOnSPDOff(0), fhEtaTrackInConeITSRefitOffSPDOff(0),
212 fhEtaPhiTrackInConeITSRefitOnSPDOn(0),fhEtaPhiTrackInConeITSRefitOnSPDOff(0), fhEtaPhiTrackInConeITSRefitOffSPDOff(0),
213 fhConeSumPtTrackTOFBC0ITSRefitOnSPDOn(0), fhPtTrackInConeTOFBC0ITSRefitOnSPDOn(0),
214 fhPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0), fhEtaTrackInConeTOFBC0ITSRefitOnSPDOn(0),fhEtaPhiTrackInConeTOFBC0ITSRefitOnSPDOn(0),
215 fhPerpConeSumPtITSRefitOnSPDOn (0), fhPtInPerpConeITSRefitOnSPDOn(0), fhEtaPhiInPerpConeITSRefitOnSPDOn(0),
216 fhPerpConeSumPtTOFBC0ITSRefitOnSPDOn (0), fhPtInPerpConeTOFBC0ITSRefitOnSPDOn (0), fhEtaPhiInPerpConeTOFBC0ITSRefitOnSPDOn(0)
220 for(
Int_t i = 0; i < 5 ; i++)
227 for(
Int_t j = 0; j < 5 ; j++)
262 for(
Int_t iso =0; iso < 2; iso++)
272 for(
Int_t i = 0; i < 5 ; i++)
291 for(
Int_t i = 0; i < 2 ; i++)
307 for(
Int_t i = 0; i < 2 ; i++)
323 for(
Int_t ieta = 0; ieta < 4; ieta++)
325 for(
Int_t iphi = 0; iphi < 3; iphi++)
336 for(
Int_t ism =0; ism < 20; ism++)
357 for(
Int_t i = 0 ; i < 7 ; i++)
366 for(
Int_t icase = 0; icase < 4; icase++)
368 for(
Int_t icase = 0; icase < 6; icase++)
396 Float_t ptTrig = pCandidate->Pt() ;
397 Float_t phiTrig = pCandidate->Phi();
398 Float_t etaTrig = pCandidate->Eta();
400 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
402 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
406 AliWarning(
"Cluster not available?");
411 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
412 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
464 Float_t phiTrig = pCandidate->Phi();
465 if(phiTrig<0) phiTrig += TMath::TwoPi();
467 Float_t etaTrig = pCandidate->Eta();
469 if(pCandidate->GetDetectorTag()==
kEMCAL)
471 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
474 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
476 if(!eGeom->CheckAbsCellId(absId))
return ;
479 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
481 Int_t imEta=-1, imPhi=-1;
482 Int_t ieta =-1, iphi =-1;
484 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
486 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
488 Int_t colTrig = ieta;
489 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
491 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
493 Int_t sqrSize = int(conesize/0.0143);
497 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
498 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
502 Int_t irowmin = rowTrig-sqrSize;
503 if ( irowmin < 0 ) irowmin = 0 ;
505 Int_t irowmax = rowTrig+sqrSize;
506 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
508 for(
Int_t irow = irowmin; irow <irowmax; irow++)
510 for(
Int_t icol = 0; icol < nTotalCols; icol++)
512 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
513 if(inSector==5)
continue;
516 if(icol < AliEMCALGeoParams::fgkEMCALCols)
518 inSupMod = 2*inSector + 1;
521 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
523 inSupMod = 2*inSector;
524 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
527 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
530 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
534 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
535 if(!eGeom->CheckAbsCellId(iabsId))
continue;
537 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
545 Int_t icolmin = colTrig-sqrSize;
546 if ( icolmin < 0 ) icolmin = 0;
548 Int_t icolmax = colTrig+sqrSize;
549 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
552 for(
Int_t icol = icolmin; icol < icolmax; icol++)
554 for(
Int_t irow = 0; irow < nTotalRows; irow++)
556 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
557 if ( inSector == 5 ) continue ;
563 if(icol < AliEMCALGeoParams::fgkEMCALCols)
566 inSupMod = 2*inSector + 1;
569 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
572 inSupMod = 2*inSector;
573 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
576 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
579 if(TMath::Abs(icol-colTrig) < sqrSize) {
582 if(TMath::Abs(irow-rowTrig) < sqrSize) {
586 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
587 if(!eGeom->CheckAbsCellId(iabsId))
589 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
593 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
603 Float_t ptTrig = pCandidate->Pt();
625 Float_t ptTrig = pCandidate->Pt() ;
626 Float_t phiTrig = pCandidate->Phi();
627 Float_t etaTrig = pCandidate->Eta();
632 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
634 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
638 AliWarning(
"Track not available?");
645 if ( pCandidate->GetDetectorTag() ==
kCTS )
648 Bool_t contained = kFALSE;
650 for(
Int_t i = 0; i < 4; i++)
652 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
655 if ( contained ) continue ;
671 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
673 phiBandPtSum+=track->Pt();
678 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
680 etaBandPtSum+=track->Pt();
686 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
687 Double_t dEta = etaTrig - track->Eta();
688 Double_t arg = dPhi*dPhi + dEta*dEta;
690 if(TMath::Sqrt(arg) < conesize)
692 sumptPerp+=track->Pt();
700 ULong_t status = track->GetStatus();
701 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
702 Int_t trackBC = track->GetTOFBunchCrossing(bz);
705 if(okTOF && trackBC == 0)
710 sumptPerpBC0+=track->Pt();
713 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
720 sumptPerpITSSPD+=track->Pt();
723 if(okTOF && trackBC == 0 && !bConstrained)
728 sumptPerpBC0ITSSPD+=track->Pt();
734 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
735 arg = dPhi*dPhi + dEta*dEta;
737 if(TMath::Sqrt(arg) < conesize)
739 sumptPerp+=track->Pt();
747 ULong_t status = track->GetStatus();
748 Bool_t okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
749 Int_t trackBC = track->GetTOFBunchCrossing(bz);
752 if(okTOF && trackBC == 0)
757 sumptPerpBC0+=track->Pt();
760 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
767 sumptPerpITSSPD+=track->Pt();
770 if(okTOF && trackBC == 0 && !bConstrained)
775 sumptPerpBC0ITSSPD+=track->Pt();
805 Float_t etaUEptsumCluster = 0 ;
806 Float_t phiUEptsumCluster = 0 ;
815 Float_t coneA = conesize*conesize*TMath::Pi();
816 Float_t ptTrig = pCandidate->Pt() ;
817 Float_t phiTrig = pCandidate->Phi();
818 Float_t etaTrig = pCandidate->Eta();
824 Float_t phiUEptsumTrackNorm = 0 ;
825 Float_t etaUEptsumTrackNorm = 0 ;
826 Float_t coneptsumTrackSubPhi = 0 ;
827 Float_t coneptsumTrackSubEta = 0 ;
828 Float_t coneptsumTrackSubPhiNorm = 0 ;
829 Float_t coneptsumTrackSubEtaNorm = 0 ;
830 etaBandptsumTrackNorm = 0 ;
842 Float_t correctConeSumTrack = 1;
843 Float_t correctConeSumTrackPhi = 1;
846 phiUEptsumTrack,etaUEptsumTrack,
847 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
848 correctConeSumTrack,correctConeSumTrackPhi);
850 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
851 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
853 etaBandptsumTrackNorm = etaUEptsumTrackNorm;
863 if(coneptsumTrack > 0)
865 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
866 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
878 Float_t phiUEptsumClusterNorm = 0 ;
879 Float_t etaUEptsumClusterNorm = 0 ;
880 Float_t coneptsumClusterSubPhi = 0 ;
881 Float_t coneptsumClusterSubEta = 0 ;
882 Float_t coneptsumClusterSubPhiNorm = 0 ;
883 Float_t coneptsumClusterSubEtaNorm = 0 ;
884 Float_t phiUEptsumCellNorm = 0 ;
885 Float_t etaUEptsumCellNorm = 0 ;
886 Float_t coneptsumCellSubPhi = 0 ;
887 Float_t coneptsumCellSubEta = 0 ;
888 Float_t coneptsumCellSubPhiNorm = 0 ;
889 Float_t coneptsumCellSubEtaNorm = 0 ;
890 etaBandptsumClusterNorm = 0;
905 Float_t correctConeSumClusterEta = 1;
906 Float_t correctConeSumClusterPhi = 1;
909 phiUEptsumCluster,etaUEptsumCluster,
910 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
911 correctConeSumClusterEta,correctConeSumClusterPhi);
922 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
923 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
925 etaBandptsumClusterNorm = etaUEptsumClusterNorm;
937 if(coneptsumCluster!=0)
939 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
940 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
961 Float_t emcPhiSize = TMath::DegToRad()*100.;
963 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
964 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
968 Float_t correctConeSumCellEta = 1;
969 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
971 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
975 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
978 Float_t correctConeSumCellPhi = 1;
980 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
981 (phiTrig-conesize < 80*TMath::DegToRad()))
984 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
985 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
991 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
995 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
996 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
1007 if ( coneptsumCell > 0.01 )
1009 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
1010 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
1026 Double_t sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
1027 Double_t sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
1048 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
1049 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
1073 coneptLeadCluster = 0;
1074 coneptsumCluster = 0;
1078 Float_t ptTrig = aodParticle->Pt();
1099 for(
Int_t icase = 0; icase < 4; icase++)
1101 for(
Int_t icase = 0; icase < 6; icase++)
1106 if(coneptLeadCluster > 0 || coneptsumCluster > 0)
1107 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumCluster,coneptLeadCluster));
1118 Float_t coneptsumClusterPerCut[20];
1122 coneptsumClusterPerCut[icut] = 0;
1129 Float_t ptSumPairFECHighCut = 0;
1130 Float_t ptSumPairFEC2MaxHighCut = 0;
1131 Float_t ptSumOddFECHighCut = 0;
1132 Float_t ptSumOddFEC2MaxHighCut = 0;
1133 Float_t ptSumFECMaxCorr[] = {0,0,0,0};
1134 Float_t ptSumFECCorr[] = {0,0,0,0,0,0};
1136 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1138 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1159 coneptsumCluster+=ptcone;
1160 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1166 if ( ptcone >
fPtCutInCone[icut] ) coneptsumClusterPerCut[icut]+=ptcone;
1175 Int_t absIdMax = -1;
1176 Int_t absIdMax2 = -1;
1179 Int_t ncells = calo->GetNCells();
1180 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1182 Int_t absId = calo->GetCellsAbsId()[ipos];
1189 absIdMax2 = absIdMax;
1194 else if(eCell > emax2)
1204 Int_t absIdCorr[] = {-1,-1,-1,-1};
1210 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1212 Int_t absId = calo->GetCellsAbsId()[ipos];
1218 if( absId == absIdMax || weight < 0.01 )
continue;
1220 for(
Int_t id = 0;
id < 4;
id++)
1222 if ( absId == absIdCorr[
id] )
1225 Int_t diff = TMath::Abs(absId-absIdMax);
1226 if ( diff == 1 ) near = kTRUE;
1227 else if( diff >= 8 ) far = kTRUE;
1234 if ( nCorr == 0 ) index = 0;
1237 if ( near && !far ) index = 1;
1238 else if (!near && far ) index = 2;
1239 else if ( near && far ) index = 3;
1242 if ( index >= 0 ) ptSumFECMaxCorr[index]+=ptcone;
1247 Int_t absIdCorrSecondary [] = {-1,-1,-1,-1};
1249 for (
Int_t ipos = 0; ipos < ncells; ipos++)
1251 Int_t absId = calo->GetCellsAbsId()[ipos];
1256 if( absId == absIdCorr[0] || absId == absIdCorr[1] ||
1257 absId == absIdCorr[2] || absId == absIdCorr[3] || weight < 0.01 )
continue;
1261 if(!correlOK)
continue;
1263 for (
Int_t ipos2 = ipos+1; ipos2 < ncells; ipos2++)
1265 Int_t absId2 = calo->GetCellsAbsId()[ipos2];
1270 if( absId2 == absIdCorr[0] || absId2 == absIdCorr[1] || absId2 == absId ||
1271 absId2 == absIdCorr[2] || absId2 == absIdCorr[3] || weight < 0.01 )
continue;
1273 for(
Int_t id = 0;
id < 4;
id++)
1275 if ( absId2 == absIdCorrSecondary[
id] )
1282 if ( counter == 1 ) distCase = 1;
1283 else if ( counter >= 2 ) distCase = 2;
1285 if ( nCorr > 0 ) distCase+=3;
1287 ptSumFECCorr[distCase]+=ptcone;
1295 for(
Int_t id = 0;
id < 4;
id++)
1297 if ( absIdMax == absIdCorr[
id]) idMax = id;
1298 else if( absIdMax2 == absIdCorr[
id]) idMax2 = id;
1301 if ( idMax==0 || idMax == 3 ) ptSumPairFEC+=ptcone;
1302 else ptSumOddFEC +=ptcone;
1306 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2Max+=ptcone;
1307 else ptSumOddFEC2Max +=ptcone;
1312 if ( idMax==0 || idMax == 3 ) ptSumPairFECHighCut+=ptcone;
1313 else ptSumOddFECHighCut +=ptcone;
1317 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2MaxHighCut+=ptcone;
1318 else ptSumOddFEC2MaxHighCut +=ptcone;
1327 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1328 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1351 for(
Int_t icase = 0; icase < 4; icase++)
1353 for(
Int_t icase = 0; icase < 6; icase++)
1369 Float_t ptTrig = aodParticle->Pt();
1370 Float_t phiTrig = aodParticle->Phi();
1371 if(phiTrig<0) phiTrig += TMath::TwoPi();
1372 Float_t etaTrig = aodParticle->Eta();
1374 if(aodParticle->GetDetectorTag()==
kEMCAL)
1376 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1379 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1381 if(!eGeom->CheckAbsCellId(absId))
return ;
1384 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1386 Int_t imEta=-1, imPhi=-1;
1387 Int_t ieta =-1, iphi =-1;
1389 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1391 Int_t iEta=-1, iPhi=-1;
1392 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1394 Int_t colTrig = iEta;
1395 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1396 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1398 Int_t sqrSize = int(conesize/0.0143);
1403 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1405 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1407 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1408 if(inSector==5)
continue;
1410 Int_t inSupMod = -1;
1412 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1414 inSupMod = 2*inSector + 1;
1417 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1419 inSupMod = 2*inSector;
1420 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1423 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1425 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1426 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1429 coneptsumCell += cells->GetCellAmplitude(iabsId);
1445 coneptLeadTrack = 0;
1449 Float_t ptTrig = aodParticle->Pt();
1469 if(coneptLeadTrack > 0 || coneptsumTrack > 0)
1470 AliError(Form(
"No ref tracks!!! sum %f, lead %f",coneptsumTrack,coneptLeadTrack));
1479 Float_t coneptsumTrackPerCut[20];
1480 Float_t coneptsumTrackPerEtaCut[20];
1481 Float_t coneptsumTrackTOFBC0 = 0;
1482 Float_t coneptsumTrackTOFBCN = 0;
1483 Float_t coneptsumTrackTOFNo = 0;
1484 Float_t coneptsumTrackITSRefitOnSPDOn = 0;
1485 Float_t coneptsumTrackITSRefitOnSPDOff = 0;
1486 Float_t coneptsumTrackITSRefitOffSPDOff = 0;
1487 Float_t coneptsumTrackTOFBC0ITSRefitOnSPDOn = 0;
1492 coneptsumTrackPerCut[icut] = 0;
1495 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1497 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1499 pTtrack = track->Pt();
1506 coneptsumTrack+=pTtrack;
1507 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1513 if ( pTtrack >
fPtCutInCone[icut] ) coneptsumTrackPerCut[icut]+=pTtrack;
1521 if ( TMath::Abs(track->Eta()) <
fEtaCutInCone[icut] ) coneptsumTrackPerEtaCut[icut]+=pTtrack;
1529 phitrack = track->Phi();
1530 etatrack = track->Eta();
1537 ULong_t status = track->GetStatus();
1538 okTOF = (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ;
1539 trackBC = track->GetTOFBunchCrossing(bz);
1540 Double32_t tof = track->GetTOFsignal()*1e-3;
1556 coneptsumTrackTOFBC0 += pTtrack;
1564 coneptsumTrackTOFBCN += pTtrack;
1573 coneptsumTrackTOFNo += pTtrack;
1577 Bool_t bITSRefit = (status & AliVTrack::kITSrefit) == AliVTrack::kITSrefit;
1578 Bool_t bConstrained = (!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1));
1586 coneptsumTrackITSRefitOnSPDOff += pTtrack;
1594 coneptsumTrackITSRefitOffSPDOff += pTtrack;
1603 coneptsumTrackITSRefitOnSPDOn += pTtrack;
1610 if(okTOF && trackBC == 0 && !bConstrained)
1616 coneptsumTrackTOFBC0ITSRefitOnSPDOn += pTtrack;
1657 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
1658 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
1712 if(energy < 8)
return;
1718 Int_t nVerticesSPD = -1;
1719 Int_t nVerticesTracks = -1;
1723 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
1724 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
1729 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
1730 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
1742 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
1746 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
1747 ncont=pv->GetNContributors();
1748 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
1750 diamZ = esdEv->GetDiamondZ();
1754 AliAODVertex *pv=aodEv->GetVertex(iVert);
1755 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
1756 ncont=pv->GetNContributors();
1757 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
1759 diamZ = aodEv->GetDiamondZ();
1762 Double_t distZ = TMath::Abs(z2-z1);
1763 diamZ = TMath::Abs(z2-diamZ);
1781 Int_t clusterID = pCandidate->GetCaloLabel(0) ;
1782 Int_t nMaxima = pCandidate->GetNLM();
1783 Int_t mcTag = pCandidate->GetTag() ;
1784 Bool_t isolated = pCandidate->IsIsolated();
1786 if ( clusterID < 0 )
1788 AliWarning(Form(
"ID of cluster = %d, not possible!", clusterID));
1792 Float_t m02 = pCandidate->GetM02() ;
1794 Float_t pt = pCandidate->Pt();
1795 Float_t eta = pCandidate->Eta();
1796 Float_t phi = pCandidate->Phi();
1797 if(phi<0) phi+= TMath::TwoPi();
1816 Int_t decayTag = pCandidate->DecayTag();
1817 if(decayTag < 0) decayTag = 0;
1854 Int_t ptsumBin = -1;
1855 Int_t leadptBin = -1;
1857 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
1867 if ( leadptBin >= 0 )
1876 if ( leadptBin == 0 )
1877 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
1881 if ( ptsumBin >= 0 )
1893 Int_t decayTag = pCandidate->DecayTag();
1894 if(decayTag < 0) decayTag = 0;
1956 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1978 Int_t ptTrigBin = -1;
1987 if ( ptTrigBin >= 0 )
2004 Int_t decayTag = pCandidate->DecayTag();
2005 if(decayTag < 0) decayTag = 0;
2108 Int_t noverlaps = 0;
2112 Int_t overpdg[nlabels];
2113 Int_t overlab[nlabels];
2143 if ( noverlaps == 1 )
2170 else if (noverlaps > 1 )
2256 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
2262 Float_t dEdx = track->GetTPCsignal();
2302 const Int_t buffersize = 255;
2303 char onePar[buffersize] ;
2305 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
2309 snprintf(onePar, buffersize,
"Isolation Cand Detector: %s;",
fIsoDetectorString.Data()) ;
2311 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
2313 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
2315 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
2317 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
2322 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
2324 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
2329 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
2334 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
2339 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
2344 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
2355 return new TObjString(parList) ;
2365 outputContainer->SetName(
"IsolatedParticleHistos") ;
2414 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
2417 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
2422 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
2425 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
2429 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
2432 TString sParticle =
", x^{0,#pm}";
2436 TString parTitle = Form(
"#it{R} = %2.2f%s%s",
GetIsolationCut()->GetConeSize(), sThreshold.Data(),sParticle.Data());
2438 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
2441 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
2442 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
2443 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
2444 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
2446 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
2447 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
2448 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
2449 "OtherDecay",
"Electron",
"Hadron"} ;
2452 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
2453 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
2454 "#pi^{0}" ,
"#eta"} ;
2456 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
2457 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
2463 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2464 nptbins,ptmin,ptmax);
2465 fhENoIso->SetYTitle(
"#it{counts}");
2466 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
2470 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2471 nptbins,ptmin,ptmax);
2473 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2477 Form(
"Number of not isolated leading particles #eta vs #phi, %s",parTitle.Data()),
2478 netabins,etamin,etamax,nphibins,phimin,phimax);
2489 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
2490 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2491 nptbins,ptmin,ptmax);
2493 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2496 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
2497 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2498 nptbins,ptmin,ptmax);
2499 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
2500 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2504 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2505 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2507 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2511 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2512 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2514 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2522 TString isoName [] = {
"NoIso",
"Iso"};
2523 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2527 for(
Int_t iso = 0; iso < 2; iso++)
2531 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2532 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2533 nptbins,ptmin,ptmax);
2534 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2535 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2536 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2540 Form(
"Number of %s leading Pi0 decay particles #eta vs #phi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2541 netabins,etamin,etamax,nphibins,phimin,phimax);
2549 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2550 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2551 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2562 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2563 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2564 nptbins,ptmin,ptmax);
2565 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2566 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2567 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2576 TString isoName [] = {
"NoIso",
"Iso"};
2577 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2580 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2581 nptbins,ptmin,ptmax);
2582 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2583 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2584 outputContainer->Add(
fhEIso) ;
2587 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2588 nptbins,ptmin,ptmax);
2589 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
2590 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2591 outputContainer->Add(
fhPtIso) ;
2594 Form(
"Number of isolated particles vs #phi, %s",parTitle.Data()),
2595 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2597 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2601 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
2602 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2604 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2608 Form(
"Number of isolated particles #eta vs #phi, %s",parTitle.Data()),
2609 netabins,etamin,etamax,nphibins,phimin,phimax);
2617 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2618 nptbins,ptmin,ptmax, 100,0,100);
2624 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2625 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
2634 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2635 nptbins,ptmin,ptmax,10,0,10);
2640 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2641 nptbins,ptmin,ptmax,10,0,10);
2648 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2649 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2650 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
2651 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2655 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2656 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2658 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2662 (Form(
"hPtLambda0Eiso"),
2663 Form(
"ABCD 3D Matrix: #it{p}_{T} vs #lambda_{0}^{2} vs E_{T}^{iso}, %s",parTitle.Data()),70,0.,70.,200,0.,2.,110,-10.,100.);
2670 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2671 netabins,etamin,etamax,nphibins,phimin,phimax);
2678 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for #it{R} = %2.2f",r),
2679 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2680 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
2681 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2715 (Form(
"hPtLeadCone_Bin%d",ibin),
2716 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2723 (Form(
"hSumPtCone_Bin%d",ibin),
2724 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2737 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2738 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2745 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2746 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2760 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2761 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2768 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2769 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2780 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
2781 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2788 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
2789 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2801 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2802 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2809 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2810 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2850 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
2851 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2858 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
2859 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2872 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2873 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2880 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2881 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2895 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2896 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2903 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2904 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2915 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
2916 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2923 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
2924 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2936 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2937 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2944 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2945 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2959 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2960 100,0,100,nptinconebins,ptinconemin,ptinconemax);
2970 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2971 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2979 Form(
"Cluster #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f",r),
2988 Form(
"Cluster #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3000 {
"No correlation with main cell",
"Main cell corr. with near cell",
3001 "Main cell corr. with far cell",
"Main cell corr. with near and far"};
3004 {
"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",
3005 "Main cell corr. & no other inside cluster" ,
"Main cell corr. & 1 cell corr. inside cluster",
"Main cell corr. & >=2 cell corr. inside cluster"};
3007 for(
Int_t icase = 0; icase < 4; icase++)
3010 (Form(
"hConePtSumClusterFECCorrMax_Case%d",icase),
3011 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, %s",r, titleMaxFEC[icase].
Data()),
3012 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3018 for(
Int_t icase = 0; icase < 6; icase++)
3021 (Form(
"hConePtSumClusterFECCorr_Case%d",icase),
3022 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, %s",r, titleFEC[icase].
Data()),
3023 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3030 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC",r),
3031 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3037 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, (n,n+1) max cells",r),
3038 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3044 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, #it{E}_{T} > 0.7 GeV",r),
3045 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3051 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),
3052 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3058 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC",r),
3059 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3065 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, (n,n+1) max cells",r),
3066 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3072 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, #it{E}_{T} > 0.7 GeV",r),
3073 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3079 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),
3080 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3088 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
3089 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3098 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3099 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3108 "#Sigma cluster #it{p}_{T} in UE Eta Band",
3109 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3115 "#Sigma cluster #it{p}_{T} UE Phi Band",
3116 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3122 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} in UE Eta Band",
3123 netabins,etamin,etamax,nphibins,phimin,phimax);
3130 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} UE Phi Band",
3131 netabins,etamin,etamax,nphibins,phimin,phimax);
3140 "#Sigma cell #it{p}_{T} in UE Eta Band",
3141 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3147 "#Sigma cell #it{p}_{T} UE Phi Band",
3148 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3154 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} in UE Eta Band",
3155 netabins,etamin,etamax,nphibins,phimin,phimax);
3162 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} UE Phi Band",
3163 netabins,etamin,etamax,nphibins,phimin,phimax);
3171 Form(
"#eta vs #phi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
3172 netabins,-1,1,nphibins,0,TMath::TwoPi());
3178 Form(
"#eta vs #phi of clusters in #phi band isolation cone for #it{R} = %2.2f",r),
3179 netabins,-1,1,nphibins,0,TMath::TwoPi());
3185 Form(
"#eta vs #phi of clusters in cone for #it{R} = %2.2f",r),
3186 netabins,-1,1,nphibins,0,TMath::TwoPi());
3192 Form(
"#eta vs #phi of all clusters"),
3193 netabins,-1,1,nphibins,0,TMath::TwoPi());
3201 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
3202 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3210 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
3211 nptbins,ptmin,ptmax,1000,0,50);
3217 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
3224 Form(
"#col vs #row of cells in #phi band isolation cone for #it{R} = %2.2f",r),
3234 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3235 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3241 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3242 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3248 Form(
"Trigger #eta vs #phi, Clusters #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3249 netabins,etamin,etamax,nphibins,phimin,phimax);
3256 Form(
"Trigger #eta vs #phi, Clusters #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3257 netabins,etamin,etamax,nphibins,phimin,phimax);
3266 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3267 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3273 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3274 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3280 Form(
"Trigger #eta vs #phi, Cells #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3281 netabins,etamin,etamax,nphibins,phimin,phimax);
3288 Form(
"Trigger #eta vs #phi, Cells #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3289 netabins,etamin,etamax,nphibins,phimin,phimax);
3297 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
3298 nptbins,ptmin,ptmax,100,0,1);
3304 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#phi ",r),
3305 netabins,etamin,etamax,nphibins,phimin,phimax);
3312 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance",r),
3313 nptbins,ptmin,ptmax,100,0,1);
3319 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance, in trigger #eta-#phi ",r),
3320 netabins,etamin,etamax,nphibins,phimin,phimax);
3327 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3328 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3334 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3335 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3341 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),
3342 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3348 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),
3349 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3355 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for cluster (before normalization), R=%2.2f",r),
3356 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
3362 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for cluster (before normalization), R=%2.2f",r),
3363 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
3371 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
3372 nptbins,ptmin,ptmax,100,0,1);
3378 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#phi ",r),
3379 netabins,etamin,etamax,nphibins,phimin,phimax);
3386 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance",r),
3387 nptbins,ptmin,ptmax,100,0,1);
3393 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance, in trigger #eta-#phi ",r),
3394 netabins,etamin,etamax,nphibins,phimin,phimax);
3402 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3403 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3409 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3410 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3416 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),
3417 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3423 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),
3424 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3437 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
3438 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3444 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
3445 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3453 Form(
"Track #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f",r),
3462 Form(
"Track #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3474 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f",r),
3483 Form(
"Track #Sigma #it{p}_{T}, different #eta cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3493 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3494 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3501 Int_t ntofbins = 1000;
3502 Int_t mintof = -500;
3505 fhTrackTOFInCone =
new TH2F (
"hTrackTOFInCone",
"TOF signal vs track pT", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
3510 fhTrackTOFInConeBC0 =
new TH2F (
"hTrackTOFInConeBC0",
"TOF signal vs track pT, BC=0", nptbins,ptmin,ptmax,ntofbins,mintof,maxtof);
3516 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3517 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3523 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f",r),
3524 netabins,-1,1,nphibins,0,TMath::TwoPi());
3530 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f",r),
3531 nptbins,ptmin,ptmax,netabins,-1,1);
3537 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f",r),
3538 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3549 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0",r),
3550 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3556 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC!=0",r),
3557 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3563 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no TOF",r),
3564 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3570 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
3571 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3577 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
3578 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3584 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no TOF",r),
3585 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3592 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
3593 netabins,-1,1,nphibins,0,TMath::TwoPi());
3599 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
3600 netabins,-1,1,nphibins,0,TMath::TwoPi());
3606 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no TOF",r),
3607 netabins,-1,1,nphibins,0,TMath::TwoPi());
3613 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
3614 nptbins,ptmin,ptmax,netabins,-1,1);
3620 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
3621 nptbins,ptmin,ptmax,netabins,-1,1);
3627 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, no TOF",r),
3628 nptbins,ptmin,ptmax,netabins,-1,1);
3634 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC=0",r),
3635 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3641 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC!=0",r),
3642 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3648 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, no TOF",r),
3649 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3658 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD On",r),
3659 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3665 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track ITS Refit SPD Off",r),
3666 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3672 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track no ITS Refit SPD Off",r),
3673 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3679 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0",r),
3680 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3686 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
3687 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3693 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
3694 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3701 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
3702 netabins,-1,1,nphibins,0,TMath::TwoPi());
3708 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
3709 netabins,-1,1,nphibins,0,TMath::TwoPi());
3715 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
3716 netabins,-1,1,nphibins,0,TMath::TwoPi());
3722 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
3723 nptbins,ptmin,ptmax,netabins,-1,1);
3729 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
3730 nptbins,ptmin,ptmax,netabins,-1,1);
3736 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
3737 nptbins,ptmin,ptmax,netabins,-1,1);
3743 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD On",r),
3744 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3750 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, ITS Refit SPD Off",r),
3751 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3757 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, no ITS Refit SPD Off",r),
3758 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3767 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f, track TOF BC=0, track ITS Refit SPD On",r),
3768 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3774 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC=0, track ITS Refit SPD On",r),
3775 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3781 Form(
"#eta vs #varphi of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
3782 netabins,-1,1,nphibins,0,TMath::TwoPi());
3788 Form(
"#eta vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
3789 nptbins,ptmin,ptmax,netabins,-1,1);
3795 Form(
"#varphi vs pT of Tracks in cone for #it{R} = %2.2f, TOF BC=0, track ITS Refit SPD On",r),
3796 nptbins,ptmin,ptmax,nphibins,0,TMath::TwoPi());
3805 "#Sigma track #it{p}_{T} in UE Eta Band",
3806 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3812 "#Sigma track #it{p}_{T} in UE Phi Band",
3813 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax*8);
3820 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Eta Band",
3821 netabins,etamin,etamax,nphibins,phimin,phimax);
3828 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Phi Band",
3829 netabins,etamin,etamax,nphibins,phimin,phimax);
3836 Form(
"#eta vs #phi of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
3837 netabins,-1,1,nphibins,0,TMath::TwoPi());
3843 Form(
"#eta vs #phi of tracks in #phi band isolation cone for #it{R} = %2.2f",r),
3844 netabins,-1,1,nphibins,0,TMath::TwoPi());
3850 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3851 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3857 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3858 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3864 Form(
"Trigger #eta vs #phi, Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3865 netabins,etamin,etamax,nphibins,phimin,phimax);
3872 Form(
"Trigger #eta vs #phi, Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3873 netabins,etamin,etamax,nphibins,phimin,phimax);
3880 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance",r),
3881 nptbins,ptmin,ptmax,100,0,1);
3887 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance, in trigger #eta-#phi ",r),
3888 netabins,etamin,etamax,nphibins,phimin,phimax);
3895 Form(
"#Sigma #it{p}_{T} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3896 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3902 Form(
"#Sigma #it{p}_{T, norm} in cone after bkg sub from #phi band vs #Sigma #it{p}_{T} in cone before bkg sub, R=%2.2f",r),
3903 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3909 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),
3910 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3916 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),
3917 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3925 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
3926 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3932 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
3933 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3941 Form(
"#eta vs #varphi of all Tracks"),
3942 netabins,-1,1,nphibins,0,TMath::TwoPi());
3949 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, TOF BC=0",r),
3950 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3956 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, TOF BC=0",r),
3957 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3963 Form(
"#eta vs #varphi of all Tracks, TOF BC=0"),
3964 netabins,-1,1,nphibins,0,TMath::TwoPi());
3971 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, ITS Refit, SPD On",r),
3972 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3978 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, ITS Refit, SPD On",r),
3979 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3985 Form(
"#eta vs #varphi of all Tracks, ITS Refit, SPD On"),
3986 netabins,-1,1,nphibins,0,TMath::TwoPi());
3994 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, TOF BC=0, ITS refit, SPD on",r),
3995 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4001 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f, TOF BC=0, ITS refit, SPD on",r),
4002 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4008 Form(
"#eta vs #varphi of all Tracks, TOF BC=0, ITS refit, SPD on"),
4009 netabins,-1,1,nphibins,0,TMath::TwoPi());
4017 Form(
"#eta vs #phi of all Tracks"),
4018 netabins,-1,1,nphibins,0,TMath::TwoPi());
4024 Form(
"#eta vs #phi of Tracks in cone for #it{R} = %2.2f",r),
4025 netabins,-1,1,nphibins,0,TMath::TwoPi());
4031 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for tracks (before normalization), R=%2.2f",r),
4032 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
4038 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for tracks (before normalization), R=%2.2f",r),
4039 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
4049 Form(
"Track vs Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4050 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4056 Form(
"#Sigma #it{p}_{T}^{cluster}/#Sigma #it{p}_{T}^{track} in isolation cone for #it{R} = %2.2f",r),
4057 nptbins,ptmin,ptmax,200,0,5);
4064 Form(
"Track vs Cluster lead #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4065 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
4071 Form(
" #it{p}^{leading cluster}_{T}/#it{p}^{leading track}_{T} in isolation cone for #it{R} = %2.2f",r),
4072 nptbins,ptmin,ptmax,200,0,5);
4081 Form(
"Track vs cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4082 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4088 Form(
"Track and Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4089 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4095 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
4096 netabins,etamin,etamax,nphibins,phimin,phimax);
4106 Form(
"#Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
4107 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4113 Form(
"#Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
4114 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4120 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
4121 netabins,etamin,etamax,nphibins,phimin,phimax);
4128 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
4129 netabins,etamin,etamax,nphibins,phimin,phimax);
4136 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
4137 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4143 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
4144 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4150 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
4151 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4157 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
4158 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
4164 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
4165 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4171 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
4172 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4178 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
4179 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4185 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
4186 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4195 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
4196 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4202 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
4203 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
4209 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
4210 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4216 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
4217 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
4223 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
4224 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4230 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
4231 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
4237 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
4238 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4244 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
4245 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
4251 Form(
"Trigger #eta vs #phi, Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
4252 netabins,etamin,etamax,nphibins,phimin,phimax);
4259 Form(
"Trigger #eta vs #phi, Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
4260 netabins,etamin,etamax,nphibins,phimin,phimax);
4269 TString region[] = {
"ITS",
"TPC",
"TRD",
"TOF",
"Top EMCal",
"In EMCal"};
4271 for(
Int_t iso = 0; iso < 2; iso++)
4276 (Form(
"hTrackMatchedDEta%s",isoName[iso].
Data()),
4277 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
4278 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4283 (Form(
"hTrackMatchedDPhi%s",isoName[iso].
Data()),
4284 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
4285 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
4290 (Form(
"hTrackMatchedDEtaDPhi%s",isoName[iso].
Data()),
4291 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
4292 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
4305 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
4306 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
4311 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
4312 nptbins,ptmin,ptmax,nresetabins,resphimin,resphimax);
4317 (Form(
"hTrackMatchedDEtaDPhi%s_MC%s",isoName[iso].
Data(),mcPartName[imc].
Data()),
4318 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
4319 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
4331 (Form(
"hdEdx%s",isoName[iso].
Data()),
4332 Form(
"%s - Matched track <d#it{E}/d#it{x}> vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
4333 nptbins,ptmin,ptmax,ndedxbins, dedxmin, dedxmax);
4334 fhdEdx[iso]->SetXTitle(
"#it{E} (GeV)");
4335 fhdEdx[iso]->SetYTitle(
"<d#it{E}/d#it{x}>");
4336 outputContainer->Add(
fhdEdx[iso]);
4339 (Form(
"hEOverP%s",isoName[iso].
Data()),
4340 Form(
"%s - Matched track #it{E}/#it{p} vs cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
4341 nptbins,ptmin,ptmax,nPoverEbins,pOverEmin,pOverEmax);
4342 fhEOverP[iso]->SetXTitle(
"#it{E} (GeV)");
4343 fhEOverP[iso]->SetYTitle(
"#it{E}/#it{p}");
4344 outputContainer->Add(
fhEOverP[iso]);
4349 (Form(
"hTrackMatchedMCParticle%s",isoName[iso].
Data()),
4350 Form(
"%s - Origin of particle vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
4351 nptbins,ptmin,ptmax,8,0,8);
4371 (Form(
"hELambda0%s",isoName[iso].
Data()),
4372 Form(
"%s cluster : #it{E} vs #lambda_{0}, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4373 fhELambda0[iso]->SetYTitle(
"#lambda_{0}^{2}");
4385 (Form(
"hPtLambda0%s",isoName[iso].
Data()),
4386 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, %s",isoTitle[iso].
Data(), parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4388 fhPtLambda0[iso]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
4396 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4397 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4399 fhPtLambda0MC[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
4403 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4404 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4412 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4413 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4419 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4420 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4426 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4427 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4433 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4434 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4441 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4442 nptbins,ptmin,ptmax,10,0,10);
4444 fhPtNOverlap[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
4448 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
4449 nptbins,ptmin,ptmax,10,0,10);
4461 (Form(
"hPtLambda0TRD%s",isoName[iso].
Data()),
4462 Form(
"%s cluster: #it{p}_{T} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4468 (Form(
"hELambda0TRD%s",isoName[iso].
Data()),
4469 Form(
"%s cluster: #it{E} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4485 (Form(
"hNLocMax%s",isoName[iso].
Data()),
4486 Form(
"%s - Number of local maxima in cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
4487 nptbins,ptmin,ptmax,10,0,10);
4489 fhNLocMax[iso]->SetXTitle(
"#it{E} (GeV)");
4493 (Form(
"hELambda0LocMax1%s",isoName[iso].
Data()),
4494 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{0}, #it{NLM}=1, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4500 (Form(
"hELambda1LocMax1%s",isoName[iso].
Data()),
4501 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{1}, #it{NLM}=1, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4507 (Form(
"hELambda0LocMax2%s",isoName[iso].
Data()),
4508 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{0}, #it{NLM}=2, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4514 (Form(
"hELambda1LocMax2%s",isoName[iso].
Data()),
4515 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{1}, #it{NLM}=2, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4521 ( Form(
"hELambda0LocMaxN%s",isoName[iso].
Data()),
4522 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{0}, #it{NLM}>2, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4528 (Form(
"hELambda1LocMaxN%s",isoName[iso].
Data()),
4529 Form(
"%s cluster (#eta) pairs: #it{E} vs #lambda_{1}, #it{NLM}>2, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4539 for(
Int_t ieta = 0; ieta < 4; ieta++)
4541 for(
Int_t iphi = 0; iphi < 3; iphi++)
4567 new TH2F(Form(
"hLam0_%s_eta%d_phi%d_sm%d",isoName[iso].
Data(),ieta,iphi,ism),
4568 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, sm %d, region eta %d, phi %d",
4569 isoTitle[iso].
Data(),ism,ieta,iphi),
4570 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
4578 Float_t ptLimit[] = {2,3,4,5,6,8,10,12};
4579 for(
Int_t ipt = 0; ipt < 7; ipt++)
4582 (Form(
"hEtaPhiLam0BinPtBin%d%s",ipt,isoName[iso].
Data()),
4583 Form(
"%s, #eta vs #phi in #it{p}_{T}=[%2.1f,%2.1f] GeV/#it{c} and #lambda^{2}_{0}=[0.3,0.4]",
4584 isoTitle[iso].
Data(),ptLimit[ipt],ptLimit[ipt+1]),
4585 netabins,etamin,etamax,nphibins,phimin,phimax);
4595 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s",
4596 isoTitle[iso].
Data(),parTitle.Data()),
4597 nptbins,ptmin,ptmax,500,0,500);