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),
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 fMomentum(), fMomIso(),
75 fMomDaugh1(), fMomDaugh2(),
76 fTrackVector(), fProdVertex(),
77 fCluster(0), fClustersArr(0),
79 fhEIso(0), fhPtIso(0),
80 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
82 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
84 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
86 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
87 fhPtTrackInConeOtherBC(0), fhPtTrackInConeOtherBCPileUpSPD(0),
88 fhPtTrackInConeBC0(0), fhPtTrackInConeVtxBC0(0),
89 fhPtTrackInConeBC0PileUpSPD(0),
90 fhPtInConePileUp(), fhPtInConeCent(0),
91 fhPerpConeSumPt(0), fhPtInPerpCone(0),
92 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
93 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
94 fhEtaBandCluster(0), fhPhiBandCluster(0),
95 fhEtaBandTrack(0), fhPhiBandTrack(0),
96 fhEtaBandCell(0), fhPhiBandCell(0),
97 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
98 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
99 fhConeSumPt(0), fhPtLambda0Eiso(0), fhConeSumPtCellTrack(0),
100 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
101 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
102 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
103 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
104 fhConeSumPtTrigEtaPhi(0),
105 fhConeSumPtCellTrackTrigEtaPhi(0),
106 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
107 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
108 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
109 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
110 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
111 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
112 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
113 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
114 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
115 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
116 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
117 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
118 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
119 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
120 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
121 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
122 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
123 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
124 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
125 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
126 fhConeSumPtCellvsTrack(0),
127 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
128 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
129 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
130 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
131 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
132 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
133 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
134 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
135 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
136 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
137 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
138 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
139 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
140 fhPtPrimMCPi0DecayPairOutOfCone(0),
141 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
142 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
143 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
144 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
145 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
146 fhPtPrimMCPi0DecayPairNoOverlap(0),
147 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
148 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
149 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
150 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
151 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
152 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
153 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
154 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
155 fhPtPrimMCEtaDecayPairOutOfCone(0),
156 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
157 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
158 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
159 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
160 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
161 fhPtPrimMCEtaDecayPairNoOverlap(0),
162 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
163 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
164 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
165 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
166 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
167 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
168 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
169 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
170 fhPtLeadConeBin(0), fhSumPtConeBin(0),
171 fhPtLeadConeBinMC(0), fhSumPtConeBinMC(0),
172 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
173 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
174 fhPtLeadConeBinLambda0MC(0), fhSumPtConeBinLambda0MC(0),
175 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
176 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
177 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
178 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
179 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
182 fhELambda0LocMax1(), fhELambda1LocMax1(),
183 fhELambda0LocMax2(), fhELambda1LocMax2(),
184 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
186 fhEIsoPileUp(), fhPtIsoPileUp(),
187 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
188 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
189 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
190 fhTimeNPileUpVertContributors(0),
191 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0),
192 fhConeSumPtClusterPerPtCut(0), fhConeSumPtClusterPerPtCutLargePtTrig(0),
193 fhConeSumPtTrackPerPtCut(0), fhConeSumPtTrackPerPtCutLargePtTrig(0),
194 fhConeSumPtClusterFECCorrPair(0), fhConeSumPtClusterFECCorrOdd(0),
195 fhConeSumPtClusterFECCorrPair2Max(0), fhConeSumPtClusterFECCorrOdd2Max(0),
196 fhConeSumPtClusterFECCorrPairHighCut(0), fhConeSumPtClusterFECCorrOddHighCut(0),
197 fhConeSumPtClusterFECCorrPair2MaxHighCut(0), fhConeSumPtClusterFECCorrOdd2MaxHighCut(0)
201 for(
Int_t i = 0; i < 5 ; i++)
208 for(
Int_t j = 0; j < 5 ; j++)
243 for(
Int_t iso =0; iso < 2; iso++)
253 for(
Int_t i = 0; i < 5 ; i++)
272 for(
Int_t i = 0; i < 2 ; i++)
288 for(
Int_t i = 0; i < 2 ; i++)
304 for(
Int_t ieta = 0; ieta < 4; ieta++)
306 for(
Int_t iphi = 0; iphi < 3; iphi++)
317 for(
Int_t ism =0; ism < 20; ism++)
338 for(
Int_t i = 0 ; i < 7 ; i++)
372 Float_t ptTrig = pCandidate->Pt() ;
373 Float_t phiTrig = pCandidate->Phi();
374 Float_t etaTrig = pCandidate->Eta();
376 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
378 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
382 AliWarning(
"Cluster not available?");
387 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
388 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
440 Float_t phiTrig = pCandidate->Phi();
441 if(phiTrig<0) phiTrig += TMath::TwoPi();
443 Float_t etaTrig = pCandidate->Eta();
445 if(pCandidate->GetDetectorTag()==
kEMCAL)
447 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
450 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
452 if(!eGeom->CheckAbsCellId(absId))
return ;
455 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
457 Int_t imEta=-1, imPhi=-1;
458 Int_t ieta =-1, iphi =-1;
460 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
462 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
464 Int_t colTrig = ieta;
465 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
467 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
469 Int_t sqrSize = int(conesize/0.0143);
473 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
474 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
478 Int_t irowmin = rowTrig-sqrSize;
479 if ( irowmin < 0 ) irowmin = 0 ;
481 Int_t irowmax = rowTrig+sqrSize;
482 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
484 for(
Int_t irow = irowmin; irow <irowmax; irow++)
486 for(
Int_t icol = 0; icol < nTotalCols; icol++)
488 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
489 if(inSector==5)
continue;
492 if(icol < AliEMCALGeoParams::fgkEMCALCols)
494 inSupMod = 2*inSector + 1;
497 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
499 inSupMod = 2*inSector;
500 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
503 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
506 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
510 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
511 if(!eGeom->CheckAbsCellId(iabsId))
continue;
513 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
521 Int_t icolmin = colTrig-sqrSize;
522 if ( icolmin < 0 ) icolmin = 0;
524 Int_t icolmax = colTrig+sqrSize;
525 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
528 for(
Int_t icol = icolmin; icol < icolmax; icol++)
530 for(
Int_t irow = 0; irow < nTotalRows; irow++)
532 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
533 if ( inSector == 5 ) continue ;
539 if(icol < AliEMCALGeoParams::fgkEMCALCols)
542 inSupMod = 2*inSector + 1;
545 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
548 inSupMod = 2*inSector;
549 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
552 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
555 if(TMath::Abs(icol-colTrig) < sqrSize) {
558 if(TMath::Abs(irow-rowTrig) < sqrSize) {
562 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
563 if(!eGeom->CheckAbsCellId(iabsId))
565 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
569 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
579 Float_t ptTrig = pCandidate->Pt();
598 Float_t ptTrig = pCandidate->Pt() ;
599 Float_t phiTrig = pCandidate->Phi();
600 Float_t etaTrig = pCandidate->Eta();
603 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
605 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
609 AliWarning(
"Track not available?");
616 if ( pCandidate->GetDetectorTag() ==
kCTS )
619 Bool_t contained = kFALSE;
621 for(
Int_t i = 0; i < 4; i++)
623 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
626 if ( contained ) continue ;
642 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
644 phiBandPtSum+=track->Pt();
649 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
651 etaBandPtSum+=track->Pt();
656 Double_t pTrack = TMath::Sqrt(track->Px()*track->Px()+track->Py()*track->Py());
657 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
658 Double_t dEta = etaTrig - track->Eta();
659 Double_t arg = dPhi*dPhi + dEta*dEta;
661 if(TMath::Sqrt(arg) < conesize)
664 sumptPerp+=track->Pt();
667 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
668 arg = dPhi*dPhi + dEta*dEta;
670 if(TMath::Sqrt(arg) < conesize)
673 sumptPerp+=track->Pt();
693 Float_t etaUEptsumCluster = 0 ;
694 Float_t phiUEptsumCluster = 0 ;
703 Float_t coneA = conesize*conesize*TMath::Pi();
704 Float_t ptTrig = pCandidate->Pt() ;
705 Float_t phiTrig = pCandidate->Phi();
706 Float_t etaTrig = pCandidate->Eta();
712 Float_t phiUEptsumTrackNorm = 0 ;
713 Float_t etaUEptsumTrackNorm = 0 ;
714 Float_t coneptsumTrackSubPhi = 0 ;
715 Float_t coneptsumTrackSubEta = 0 ;
716 Float_t coneptsumTrackSubPhiNorm = 0 ;
717 Float_t coneptsumTrackSubEtaNorm = 0 ;
718 etaBandptsumTrackNorm = 0 ;
730 Float_t correctConeSumTrack = 1;
731 Float_t correctConeSumTrackPhi = 1;
734 phiUEptsumTrack,etaUEptsumTrack,
735 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
736 correctConeSumTrack,correctConeSumTrackPhi);
738 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
739 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
741 etaBandptsumTrackNorm = etaUEptsumTrackNorm;
751 if(coneptsumTrack > 0)
753 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
754 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
766 Float_t phiUEptsumClusterNorm = 0 ;
767 Float_t etaUEptsumClusterNorm = 0 ;
768 Float_t coneptsumClusterSubPhi = 0 ;
769 Float_t coneptsumClusterSubEta = 0 ;
770 Float_t coneptsumClusterSubPhiNorm = 0 ;
771 Float_t coneptsumClusterSubEtaNorm = 0 ;
772 Float_t phiUEptsumCellNorm = 0 ;
773 Float_t etaUEptsumCellNorm = 0 ;
774 Float_t coneptsumCellSubPhi = 0 ;
775 Float_t coneptsumCellSubEta = 0 ;
776 Float_t coneptsumCellSubPhiNorm = 0 ;
777 Float_t coneptsumCellSubEtaNorm = 0 ;
778 etaBandptsumClusterNorm = 0;
793 Float_t correctConeSumClusterEta = 1;
794 Float_t correctConeSumClusterPhi = 1;
797 phiUEptsumCluster,etaUEptsumCluster,
798 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
799 correctConeSumClusterEta,correctConeSumClusterPhi);
810 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
811 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
813 etaBandptsumClusterNorm = etaUEptsumClusterNorm;
825 if(coneptsumCluster!=0)
827 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
828 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
849 Float_t emcPhiSize = TMath::DegToRad()*100.;
851 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
852 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
856 Float_t correctConeSumCellEta = 1;
857 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
859 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
863 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
866 Float_t correctConeSumCellPhi = 1;
868 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
869 (phiTrig-conesize < 80*TMath::DegToRad()))
872 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
873 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
879 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
883 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
884 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
895 if ( coneptsumCell > 0.01 )
897 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
898 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
914 Double_t sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
915 Double_t sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
936 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
937 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
961 coneptLeadCluster = 0;
962 coneptsumCluster = 0;
966 Float_t ptTrig = aodParticle->Pt();
984 Float_t coneptsumClusterPerCut[20];
988 coneptsumClusterPerCut[icut] = 0;
995 Float_t ptSumPairFECHighCut = 0;
996 Float_t ptSumPairFEC2MaxHighCut = 0;
997 Float_t ptSumOddFECHighCut = 0;
998 Float_t ptSumOddFEC2MaxHighCut = 0;
1000 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
1002 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
1023 coneptsumCluster+=ptcone;
1024 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1030 if ( ptcone >
fPtCutInCone[icut] ) coneptsumClusterPerCut[icut]+=ptcone;
1037 Int_t absIdMax = -1;
1038 Int_t absIdMax2 = -1;
1041 for (
Int_t ipos = 0; ipos < calo->GetNCells(); ipos++)
1043 Int_t absId = calo->GetCellsAbsId()[ipos];
1050 absIdMax2 = absIdMax;
1055 else if(eCell > emax2)
1062 Int_t absIdCorr[] = {-1,-1,-1,-1};
1068 for(
Int_t id = 0;
id < 4;
id++)
1070 if ( absIdMax == absIdCorr[
id]) idMax = id;
1071 else if( absIdMax2 == absIdCorr[
id]) idMax2 = id;
1074 if ( idMax==0 || idMax == 3 ) ptSumPairFEC+=ptcone;
1075 else ptSumOddFEC +=ptcone;
1079 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2Max+=ptcone;
1080 else ptSumOddFEC2Max +=ptcone;
1085 if ( idMax==0 || idMax == 3 ) ptSumPairFECHighCut+=ptcone;
1086 else ptSumOddFECHighCut +=ptcone;
1090 if ( idMax==0 || idMax == 3 ) ptSumPairFEC2MaxHighCut+=ptcone;
1091 else ptSumOddFEC2MaxHighCut +=ptcone;
1100 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1101 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1137 Float_t ptTrig = aodParticle->Pt();
1138 Float_t phiTrig = aodParticle->Phi();
1139 if(phiTrig<0) phiTrig += TMath::TwoPi();
1140 Float_t etaTrig = aodParticle->Eta();
1142 if(aodParticle->GetDetectorTag()==
kEMCAL)
1144 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1147 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1149 if(!eGeom->CheckAbsCellId(absId))
return ;
1152 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1154 Int_t imEta=-1, imPhi=-1;
1155 Int_t ieta =-1, iphi =-1;
1157 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1159 Int_t iEta=-1, iPhi=-1;
1160 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1162 Int_t colTrig = iEta;
1163 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1164 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1166 Int_t sqrSize = int(conesize/0.0143);
1171 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1173 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1175 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1176 if(inSector==5)
continue;
1178 Int_t inSupMod = -1;
1180 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1182 inSupMod = 2*inSector + 1;
1185 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1187 inSupMod = 2*inSector;
1188 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1191 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1193 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1194 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1197 coneptsumCell += cells->GetCellAmplitude(iabsId);
1213 coneptLeadTrack = 0;
1217 Float_t ptTrig = aodParticle->Pt();
1231 Float_t coneptsumTrackPerCut[20];
1235 coneptsumTrackPerCut[icut] = 0;
1238 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1240 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1242 pTtrack = track->Pt();
1249 ULong_t status = track->GetStatus();
1250 Bool_t okTOF = ( (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ) ;
1252 Int_t trackBC = track->GetTOFBunchCrossing(bz);
1273 coneptsumTrack+=pTtrack;
1274 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1280 if ( pTtrack >
fPtCutInCone[icut] ) coneptsumTrackPerCut[icut]+=pTtrack;
1288 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
1289 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
1334 if(energy < 8)
return;
1340 Int_t nVerticesSPD = -1;
1341 Int_t nVerticesTracks = -1;
1345 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
1346 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
1351 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
1352 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
1364 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
1368 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
1369 ncont=pv->GetNContributors();
1370 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
1372 diamZ = esdEv->GetDiamondZ();
1376 AliAODVertex *pv=aodEv->GetVertex(iVert);
1377 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
1378 ncont=pv->GetNContributors();
1379 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
1381 diamZ = aodEv->GetDiamondZ();
1384 Double_t distZ = TMath::Abs(z2-z1);
1385 diamZ = TMath::Abs(z2-diamZ);
1403 Int_t clusterID = pCandidate->GetCaloLabel(0) ;
1404 Int_t nMaxima = pCandidate->GetNLM();
1405 Int_t mcTag = pCandidate->GetTag() ;
1406 Bool_t isolated = pCandidate->IsIsolated();
1408 if ( clusterID < 0 )
1410 AliWarning(Form(
"ID of cluster = %d, not possible!", clusterID));
1414 Float_t m02 = pCandidate->GetM02() ;
1416 Float_t pt = pCandidate->Pt();
1417 Float_t eta = pCandidate->Eta();
1418 Float_t phi = pCandidate->Phi();
1419 if(phi<0) phi+= TMath::TwoPi();
1438 Int_t decayTag = pCandidate->DecayTag();
1439 if(decayTag < 0) decayTag = 0;
1476 Int_t ptsumBin = -1;
1477 Int_t leadptBin = -1;
1479 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
1489 if ( leadptBin >= 0 )
1498 if ( leadptBin == 0 )
1499 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
1503 if ( ptsumBin >= 0 )
1515 Int_t decayTag = pCandidate->DecayTag();
1516 if(decayTag < 0) decayTag = 0;
1578 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1600 Int_t ptTrigBin = -1;
1609 if ( ptTrigBin >= 0 )
1626 Int_t decayTag = pCandidate->DecayTag();
1627 if(decayTag < 0) decayTag = 0;
1730 Int_t noverlaps = 0;
1734 Int_t overpdg[nlabels];
1735 Int_t overlab[nlabels];
1765 if ( noverlaps == 1 )
1792 else if (noverlaps > 1 )
1878 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
1884 Float_t dEdx = track->GetTPCsignal();
1924 const Int_t buffersize = 255;
1925 char onePar[buffersize] ;
1927 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
1931 snprintf(onePar, buffersize,
"Isolation Cand Detector: %s;",
fIsoDetectorString.Data()) ;
1933 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
1935 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
1937 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
1939 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
1944 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
1946 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
1951 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
1956 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
1961 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
1966 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
1977 return new TObjString(parList) ;
1987 outputContainer->SetName(
"IsolatedParticleHistos") ;
2036 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
2039 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
2044 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
2047 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
2051 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
2054 TString sParticle =
", x^{0,#pm}";
2058 TString parTitle = Form(
"#it{R} = %2.2f%s%s",
GetIsolationCut()->GetConeSize(), sThreshold.Data(),sParticle.Data());
2060 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
2063 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
2064 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
2065 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
2066 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
2068 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
2069 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
2070 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
2071 "OtherDecay",
"Electron",
"Hadron"} ;
2074 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
2075 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
2076 "#pi^{0}" ,
"#eta"} ;
2078 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
2079 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
2085 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2086 nptbins,ptmin,ptmax);
2087 fhENoIso->SetYTitle(
"#it{counts}");
2088 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
2092 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
2093 nptbins,ptmin,ptmax);
2095 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2099 Form(
"Number of not isolated leading particles #eta vs #phi, %s",parTitle.Data()),
2100 netabins,etamin,etamax,nphibins,phimin,phimax);
2111 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
2112 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2113 nptbins,ptmin,ptmax);
2115 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2118 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
2119 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2120 nptbins,ptmin,ptmax);
2121 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
2122 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2126 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2127 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2129 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2133 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
2134 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2136 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2144 TString isoName [] = {
"NoIso",
"Iso"};
2145 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2149 for(
Int_t iso = 0; iso < 2; iso++)
2153 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2154 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2155 nptbins,ptmin,ptmax);
2156 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2157 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2158 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2162 Form(
"Number of %s leading Pi0 decay particles #eta vs #phi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2163 netabins,etamin,etamax,nphibins,phimin,phimax);
2171 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2172 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2173 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2184 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2185 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2186 nptbins,ptmin,ptmax);
2187 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2188 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2189 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2198 TString isoName [] = {
"NoIso",
"Iso"};
2199 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2202 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2203 nptbins,ptmin,ptmax);
2204 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2205 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2206 outputContainer->Add(
fhEIso) ;
2209 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2210 nptbins,ptmin,ptmax);
2211 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
2212 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2213 outputContainer->Add(
fhPtIso) ;
2216 Form(
"Number of isolated particles vs #phi, %s",parTitle.Data()),
2217 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2219 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2223 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
2224 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2226 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2230 Form(
"Number of isolated particles #eta vs #phi, %s",parTitle.Data()),
2231 netabins,etamin,etamax,nphibins,phimin,phimax);
2239 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2240 nptbins,ptmin,ptmax, 100,0,100);
2246 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2247 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
2256 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2257 nptbins,ptmin,ptmax,10,0,10);
2262 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2263 nptbins,ptmin,ptmax,10,0,10);
2270 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2271 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2272 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
2273 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2277 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2278 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2280 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2284 (Form(
"hPtLambda0Eiso"),
2285 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.);
2292 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2293 netabins,etamin,etamax,nphibins,phimin,phimax);
2300 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for #it{R} = %2.2f",r),
2301 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2302 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
2303 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2337 (Form(
"hPtLeadCone_Bin%d",ibin),
2338 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2345 (Form(
"hSumPtCone_Bin%d",ibin),
2346 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2359 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2360 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2367 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2368 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2382 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2383 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2390 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2391 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2402 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
2403 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2410 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
2411 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2423 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2424 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2431 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2432 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2472 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
2473 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2480 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
2481 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2494 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2495 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2502 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2503 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2517 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2518 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2525 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2526 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2537 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
2538 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2545 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
2546 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2558 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2559 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2566 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2567 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2581 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2582 100,0,100,nptinconebins,ptinconemin,ptinconemax);
2592 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2593 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2601 Form(
"Cluster #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f",r),
2610 Form(
"Cluster #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
2622 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC",r),
2623 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2629 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, (n,n+1) max cells",r),
2630 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2636 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, odd column in FEC, #it{E}_{T} > 0.7 GeV",r),
2637 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2643 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),
2644 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2650 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC",r),
2651 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2657 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, (n,n+1) max cells",r),
2658 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2664 Form(
"Cluster #Sigma #it{p}_{T} for #it{R} = %2.2f, pair column in FEC, #it{E}_{T} > 0.7 GeV",r),
2665 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2671 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),
2672 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2680 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
2681 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2690 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2691 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2700 "#Sigma cluster #it{p}_{T} in UE Eta Band",
2701 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2707 "#Sigma cluster #it{p}_{T} UE Phi Band",
2708 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2714 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} in UE Eta Band",
2715 netabins,etamin,etamax,nphibins,phimin,phimax);
2722 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} UE Phi Band",
2723 netabins,etamin,etamax,nphibins,phimin,phimax);
2732 "#Sigma cell #it{p}_{T} in UE Eta Band",
2733 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2739 "#Sigma cell #it{p}_{T} UE Phi Band",
2740 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2746 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} in UE Eta Band",
2747 netabins,etamin,etamax,nphibins,phimin,phimax);
2754 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} UE Phi Band",
2755 netabins,etamin,etamax,nphibins,phimin,phimax);
2763 Form(
"#eta vs #phi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
2764 netabins,-1,1,nphibins,0,TMath::TwoPi());
2770 Form(
"#eta vs #phi of clusters in #phi band isolation cone for #it{R} = %2.2f",r),
2771 netabins,-1,1,nphibins,0,TMath::TwoPi());
2777 Form(
"#eta vs #phi of clusters in cone for #it{R} = %2.2f",r),
2778 netabins,-1,1,nphibins,0,TMath::TwoPi());
2784 Form(
"#eta vs #phi of all clusters"),
2785 netabins,-1,1,nphibins,0,TMath::TwoPi());
2793 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
2794 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2802 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
2803 nptbins,ptmin,ptmax,1000,0,50);
2809 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
2816 Form(
"#col vs #row of cells in #phi band isolation cone for #it{R} = %2.2f",r),
2826 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2827 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2833 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2834 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2840 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),
2841 netabins,etamin,etamax,nphibins,phimin,phimax);
2848 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),
2849 netabins,etamin,etamax,nphibins,phimin,phimax);
2858 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2859 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2865 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2866 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2872 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),
2873 netabins,etamin,etamax,nphibins,phimin,phimax);
2880 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),
2881 netabins,etamin,etamax,nphibins,phimin,phimax);
2889 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
2890 nptbins,ptmin,ptmax,100,0,1);
2896 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#phi ",r),
2897 netabins,etamin,etamax,nphibins,phimin,phimax);
2904 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance",r),
2905 nptbins,ptmin,ptmax,100,0,1);
2911 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance, in trigger #eta-#phi ",r),
2912 netabins,etamin,etamax,nphibins,phimin,phimax);
2919 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),
2920 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2926 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),
2927 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2933 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),
2934 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2940 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),
2941 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2947 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for cluster (before normalization), R=%2.2f",r),
2948 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
2954 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for cluster (before normalization), R=%2.2f",r),
2955 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
2963 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
2964 nptbins,ptmin,ptmax,100,0,1);
2970 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#phi ",r),
2971 netabins,etamin,etamax,nphibins,phimin,phimax);
2978 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance",r),
2979 nptbins,ptmin,ptmax,100,0,1);
2985 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance, in trigger #eta-#phi ",r),
2986 netabins,etamin,etamax,nphibins,phimin,phimax);
2994 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),
2995 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3001 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),
3002 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3008 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),
3009 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3015 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),
3016 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3028 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3029 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3037 Form(
"Track #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f",r),
3046 Form(
"Track #Sigma #it{p}_{T}, different pT cuts in isolation cone for #it{R} = %2.2f, #it{p}_{T}^{trig} > 10 GeV",r),
3056 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
3057 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3063 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
3064 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3073 "#Sigma track #it{p}_{T} in UE Eta Band",
3074 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3080 "#Sigma track #it{p}_{T} in UE Phi Band",
3081 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax*8);
3088 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Eta Band",
3089 netabins,etamin,etamax,nphibins,phimin,phimax);
3096 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Phi Band",
3097 netabins,etamin,etamax,nphibins,phimin,phimax);
3104 Form(
"#eta vs #phi of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
3105 netabins,-1,1,nphibins,0,TMath::TwoPi());
3111 Form(
"#eta vs #phi of tracks in #phi band isolation cone for #it{R} = %2.2f",r),
3112 netabins,-1,1,nphibins,0,TMath::TwoPi());
3118 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3119 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3125 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3126 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3132 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),
3133 netabins,etamin,etamax,nphibins,phimin,phimax);
3140 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),
3141 netabins,etamin,etamax,nphibins,phimin,phimax);
3148 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance",r),
3149 nptbins,ptmin,ptmax,100,0,1);
3155 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance, in trigger #eta-#phi ",r),
3156 netabins,etamin,etamax,nphibins,phimin,phimax);
3163 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),
3164 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3170 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),
3171 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3177 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),
3178 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3184 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),
3185 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3193 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
3194 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3200 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
3201 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3207 Form(
"#eta vs #phi of all Tracks"),
3208 netabins,-1,1,nphibins,0,TMath::TwoPi());
3214 Form(
"#eta vs #phi of Tracks in cone for #it{R} = %2.2f",r),
3215 netabins,-1,1,nphibins,0,TMath::TwoPi());
3221 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for tracks (before normalization), R=%2.2f",r),
3222 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
3228 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for tracks (before normalization), R=%2.2f",r),
3229 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
3239 Form(
"Track vs Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3240 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3246 Form(
"#Sigma #it{p}_{T}^{cluster}/#Sigma #it{p}_{T}^{track} in isolation cone for #it{R} = %2.2f",r),
3247 nptbins,ptmin,ptmax,200,0,5);
3254 Form(
"Track vs Cluster lead #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3255 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3261 Form(
" #it{p}^{leading cluster}_{T}/#it{p}^{leading track}_{T} in isolation cone for #it{R} = %2.2f",r),
3262 nptbins,ptmin,ptmax,200,0,5);
3271 Form(
"Track vs cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3272 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3278 Form(
"Track and Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3279 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3285 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3286 netabins,etamin,etamax,nphibins,phimin,phimax);
3296 Form(
"#Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3297 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3303 Form(
"#Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3304 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3310 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),
3311 netabins,etamin,etamax,nphibins,phimin,phimax);
3318 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),
3319 netabins,etamin,etamax,nphibins,phimin,phimax);
3326 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3327 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3333 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3334 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3340 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3341 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3347 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3348 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3354 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3355 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3361 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3362 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3368 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3369 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3375 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3376 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3385 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3386 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3392 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3393 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3399 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3400 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3406 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3407 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3413 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3414 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3420 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3421 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3427 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3428 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3434 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3435 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3441 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),
3442 netabins,etamin,etamax,nphibins,phimin,phimax);
3449 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),
3450 netabins,etamin,etamax,nphibins,phimin,phimax);
3459 TString region[] = {
"ITS",
"TPC",
"TRD",
"TOF",
"Top EMCal",
"In EMCal"};
3461 for(
Int_t iso = 0; iso < 2; iso++)
3466 (Form(
"hTrackMatchedDEta%s",isoName[iso].
Data()),
3467 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3468 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3473 (Form(
"hTrackMatchedDPhi%s",isoName[iso].
Data()),
3474 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3475 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3480 (Form(
"hTrackMatchedDEtaDPhi%s",isoName[iso].
Data()),
3481 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3482 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3495 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3496 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3501 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3502 nptbins,ptmin,ptmax,nresetabins,resphimin,resphimax);
3507 (Form(
"hTrackMatchedDEtaDPhi%s_MC%s",isoName[iso].
Data(),mcPartName[imc].
Data()),
3508 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3509 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3521 (Form(
"hdEdx%s",isoName[iso].
Data()),
3522 Form(
"%s - Matched track <d#it{E}/d#it{x}> vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3523 nptbins,ptmin,ptmax,ndedxbins, dedxmin, dedxmax);
3524 fhdEdx[iso]->SetXTitle(
"#it{E} (GeV)");
3525 fhdEdx[iso]->SetYTitle(
"<d#it{E}/d#it{x}>");
3526 outputContainer->Add(
fhdEdx[iso]);
3529 (Form(
"hEOverP%s",isoName[iso].
Data()),
3530 Form(
"%s - Matched track #it{E}/#it{p} vs cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3531 nptbins,ptmin,ptmax,nPoverEbins,pOverEmin,pOverEmax);
3532 fhEOverP[iso]->SetXTitle(
"#it{E} (GeV)");
3533 fhEOverP[iso]->SetYTitle(
"#it{E}/#it{p}");
3534 outputContainer->Add(
fhEOverP[iso]);
3539 (Form(
"hTrackMatchedMCParticle%s",isoName[iso].
Data()),
3540 Form(
"%s - Origin of particle vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3541 nptbins,ptmin,ptmax,8,0,8);
3561 (Form(
"hELambda0%s",isoName[iso].
Data()),
3562 Form(
"%s cluster : #it{E} vs #lambda_{0}, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3563 fhELambda0[iso]->SetYTitle(
"#lambda_{0}^{2}");
3575 (Form(
"hPtLambda0%s",isoName[iso].
Data()),
3576 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, %s",isoTitle[iso].
Data(), parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3578 fhPtLambda0[iso]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3586 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3587 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3589 fhPtLambda0MC[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3593 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3594 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3602 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3603 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3609 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3610 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3616 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3617 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3623 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3624 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3631 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3632 nptbins,ptmin,ptmax,10,0,10);
3634 fhPtNOverlap[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3638 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3639 nptbins,ptmin,ptmax,10,0,10);
3651 (Form(
"hPtLambda0TRD%s",isoName[iso].
Data()),
3652 Form(
"%s cluster: #it{p}_{T} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3658 (Form(
"hELambda0TRD%s",isoName[iso].
Data()),
3659 Form(
"%s cluster: #it{E} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3675 (Form(
"hNLocMax%s",isoName[iso].
Data()),
3676 Form(
"%s - Number of local maxima in cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3677 nptbins,ptmin,ptmax,10,0,10);
3679 fhNLocMax[iso]->SetXTitle(
"#it{E} (GeV)");
3683 (Form(
"hELambda0LocMax1%s",isoName[iso].
Data()),
3684 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);
3690 (Form(
"hELambda1LocMax1%s",isoName[iso].
Data()),
3691 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);
3697 (Form(
"hELambda0LocMax2%s",isoName[iso].
Data()),
3698 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);
3704 (Form(
"hELambda1LocMax2%s",isoName[iso].
Data()),
3705 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);
3711 ( Form(
"hELambda0LocMaxN%s",isoName[iso].
Data()),
3712 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);
3718 (Form(
"hELambda1LocMaxN%s",isoName[iso].
Data()),
3719 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);
3729 for(
Int_t ieta = 0; ieta < 4; ieta++)
3731 for(
Int_t iphi = 0; iphi < 3; iphi++)
3757 new TH2F(Form(
"hLam0_%s_eta%d_phi%d_sm%d",isoName[iso].
Data(),ieta,iphi,ism),
3758 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, sm %d, region eta %d, phi %d",
3759 isoTitle[iso].
Data(),ism,ieta,iphi),
3760 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3768 Float_t ptLimit[] = {2,3,4,5,6,8,10,12};
3769 for(
Int_t ipt = 0; ipt < 7; ipt++)
3772 (Form(
"hEtaPhiLam0BinPtBin%d%s",ipt,isoName[iso].
Data()),
3773 Form(
"%s, #eta vs #phi in #it{p}_{T}=[%2.1f,%2.1f] GeV/#it{c} and #lambda^{2}_{0}=[0.3,0.4]",
3774 isoTitle[iso].
Data(),ptLimit[ipt],ptLimit[ipt+1]),
3775 netabins,etamin,etamax,nphibins,phimin,phimax);
3785 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s",
3786 isoTitle[iso].
Data(),parTitle.Data()),
3787 nptbins,ptmin,ptmax,500,0,500);
3795 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s, SM covered by TRD",
3796 isoTitle[iso].
Data(),parTitle.Data()),
3797 nptbins,ptmin,ptmax,500,0,500);
3805 for(
Int_t iR = 0; iR < 6; iR++)
3808 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s",
3809 isoTitle[iso].
Data(),region[iR].
Data()),
3810 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3818 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s, SM covered by TRD",
3819 isoTitle[iso].
Data(),region[iR].
Data()),
3820 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3870 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
3871 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3877 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0, pile-up from SPD",r),
3878 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3884 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3885 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3891 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3892 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3899 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0, pile-up from SPD",r),
3900 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3906 for (
Int_t i = 0; i < 7 ; i++)
3909 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f, from pile-up (%s)",r,pileUpName[i].
Data()),
3910 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3923 fhEPrimMC[i] =
new TH1F(Form(
"hEPrim_MC%s",ppname[i].
Data()),
3924 Form(
"primary photon %s : #it{E}, %s",pptype[i].
Data(),parTitle.Data()),
3925 nptbins,ptmin,ptmax);
3926 fhEPrimMC[i]->SetXTitle(
"#it{E} (GeV)");
3930 Form(
"primary photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3931 nptbins,ptmin,ptmax);
3932 fhPtPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3936 Form(
"primary isolated photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3937 nptbins,ptmin,ptmax);
3942 Form(
"primary photon %s : #eta vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3943 nptbins,ptmin,ptmax,200,-2,2);
3945 fhEtaPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3949 Form(
"primary photon %s : #phi vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3950 nptbins,ptmin,ptmax,200,0.,TMath::TwoPi());
3952 fhPhiPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3959 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3960 nptbins,ptmin,ptmax);
3965 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
3967 nptbins,ptmin,ptmax);
3972 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
3974 nptbins,ptmin,ptmax);
3979 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
3981 nptbins,ptmin,ptmax);
3986 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
3988 nptbins,ptmin,ptmax);
3993 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
3995 nptbins,ptmin,ptmax);
4001 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
4003 nptbins,ptmin,ptmax);
4008 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
4010 nptbins,ptmin,ptmax);
4015 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
4016 nptbins,ptmin,ptmax);
4021 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
4023 nptbins,ptmin,ptmax);
4028 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
4029 nptbins,ptmin,ptmax);
4034 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
4035 nptbins,ptmin,ptmax);
4040 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
4042 nptbins,ptmin,ptmax);
4047 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
4049 nptbins,ptmin,ptmax);
4054 Form(
"primary %s, overlap: #it{p}_{T}, %s",
4056 nptbins,ptmin,ptmax);
4061 Form(
"primary %s, overlap: #it{p}_{T}, %s",
4063 nptbins,ptmin,ptmax);
4075 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
4076 nptbins,ptmin,ptmax);
4081 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
4083 nptbins,ptmin,ptmax);
4088 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
4090 nptbins,ptmin,ptmax);
4095 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
4097 nptbins,ptmin,ptmax);
4102 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
4104 nptbins,ptmin,ptmax);
4109 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
4111 nptbins,ptmin,ptmax);
4117 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
4119 nptbins,ptmin,ptmax);
4124 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
4126 nptbins,ptmin,ptmax);
4131 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
4132 nptbins,ptmin,ptmax);
4137 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
4139 nptbins,ptmin,ptmax);
4144 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
4145 nptbins,ptmin,ptmax);
4150 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
4151 nptbins,ptmin,ptmax);
4156 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
4158 nptbins,ptmin,ptmax);
4163 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
4165 nptbins,ptmin,ptmax);
4170 Form(
"primary %s, overlap: #it{p}_{T}, %s",
4172 nptbins,ptmin,ptmax);
4177 Form(
"primary %s, overlap: #it{p}_{T}, %s",
4179 nptbins,ptmin,ptmax);
4191 const Int_t buffersize = 255;
4192 char name[buffersize];
4193 char title[buffersize];
4197 snprintf(name, buffersize,
"hSumPtLeadingPt_Cone_%d",icone);
4198 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
4199 fhSumPtLeadingPt[icone] =
new TH2F(name, title,nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
4205 snprintf(name, buffersize,
"hPtLeadingPt_Cone_%d",icone);
4206 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
4207 fhPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4208 fhPtLeadingPt[icone] ->SetYTitle(
"#it{p}_{T}^{cone} (GeV/#it{c})");
4209 fhPtLeadingPt[icone] ->SetXTitle(
"#it{p}_{T}^{leading} (GeV/#it{c})");
4213 snprintf(name, buffersize,
"hPerpSumPtLeadingPt_Cone_%d",icone);
4214 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
4221 snprintf(name, buffersize,
"hPerpPtLeadingPt_Cone_%d",icone);
4222 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
4223 fhPerpPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
4232 snprintf(name , buffersize,
"hSumPtLeadingPt_MC%s_Cone_%d",mcPartName[imc].
Data(),icone);
4233 snprintf(title, buffersize,
"Candidate %s #it{p}_{T} vs cone #Sigma #it{p}_{T} for #it{R}=%2.2f",mcPartType[imc].
Data(),
fConeSizes[icone]);
4243 snprintf(name, buffersize,
"hPtThres_Cone_%d_Pt%d",icone,ipt);
4244 snprintf(title, buffersize,
"Isolated candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtThresholds[ipt]);
4249 snprintf(name, buffersize,
"hPtFrac_Cone_%d_Pt%d",icone,ipt);
4250 snprintf(title, buffersize,
"Isolated candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4255 snprintf(name, buffersize,
"hSumPt_Cone_%d_Pt%d",icone,ipt);
4256 snprintf(title, buffersize,
"Isolated candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4257 fhSumPtIsolated[icone][ipt] =
new TH1F(name, title,nptbins,ptmin,ptmax);
4262 snprintf(name, buffersize,
"hPtSumDensity_Cone_%d_Pt%d",icone,ipt);
4263 snprintf(title, buffersize,
"Isolated candidate #it{p}_{T} distribution for density in #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4269 snprintf(name, buffersize,
"hPtFracPtSum_Cone_%d_Pt%d",icone,ipt);
4270 snprintf(title, buffersize,
"Isolated candidate #it{p}_{T} distribution for PtFracPtSum in #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4277 snprintf(name, buffersize,
"hEtaPhiPtThres_Cone_%d_Pt%d",icone,ipt);
4278 snprintf(title, buffersize,
"Isolated candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtThresholds[ipt]);
4284 snprintf(name, buffersize,
"hEtaPhiPtFrac_Cone_%d_Pt%d",icone,ipt);
4285 snprintf(title, buffersize,
"Isolated candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4286 fhEtaPhiPtFracIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4291 snprintf(name, buffersize,
"hEtaPhiPtSum_Cone_%d_Pt%d",icone,ipt);
4292 snprintf(title, buffersize,
"Isolated candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4293 fhEtaPhiPtSumIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4298 snprintf(name, buffersize,
"hEtaPhiSumDensity_Cone_%d_Pt%d",icone,ipt);
4299 snprintf(title, buffersize,
"Isolated candidate #eta:#phi distribution for density #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4305 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Cone_%d_Pt%d",icone,ipt);
4306 snprintf(title, buffersize,
"Isolated candidate #eta:#phi distribution for FracPtSum #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4315 snprintf(name, buffersize,
"hPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4316 snprintf(title, buffersize,
"Isolated decay candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtThresholds[ipt]);
4321 snprintf(name, buffersize,
"hPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4322 snprintf(title, buffersize,
"Isolated decay candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4327 snprintf(name, buffersize,
"hPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4328 snprintf(title, buffersize,
"Isolated decay candidate #it{p}_{T} distribution for #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4334 snprintf(name, buffersize,
"hPtSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4335 snprintf(title, buffersize,
"Isolated decay candidate #it{p}_{T} distribution for density in #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4341 snprintf(name, buffersize,
"hPtFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4342 snprintf(title, buffersize,
"Isolated decay candidate #it{p}_{T} distribution for PtFracPtSum in #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4349 snprintf(name, buffersize,
"hEtaPhiPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4350 snprintf(title, buffersize,
"Isolated decay candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtThresholds[ipt]);
4356 snprintf(name, buffersize,
"hEtaPhiPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4357 snprintf(title, buffersize,
"Isolated decay candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4364 snprintf(name, buffersize,
"hEtaPhiPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4365 snprintf(title, buffersize,
"Isolated decay candidate #eta:#phi distribution for #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4371 snprintf(name, buffersize,
"hEtaPhiSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4372 snprintf(title, buffersize,
"Isolated decay candidate #eta:#phi distribution for density #it{R} = %2.2f and #it{p}_{T}^{sum} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fSumPtThresholds[ipt]);
4378 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4379 snprintf(title, buffersize,
"Isolated decay candidate #eta:#phi distribution for FracPtSum #it{R} = %2.2f and #it{p}_{T}^{fr} = %2.2f GeV/#it{c}",
fConeSizes[icone],
fPtFractions[ipt]);
4391 snprintf(name , buffersize,
"hPtThreshMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4392 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #it{p}_{T}^{th}=%2.2f",
4400 snprintf(name , buffersize,
"hPtFracMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4401 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig}=%2.2f",
4408 snprintf(name , buffersize,
"hSumPtMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4409 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #Sigma #it{p}_{T}^{in cone}=%2.2f",
4423 for (
Int_t i = 0; i < 7 ; i++)
4426 Form(
"Number of isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4427 nptbins,ptmin,ptmax);
4433 Form(
"Number of isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4434 nptbins,ptmin,ptmax);
4440 Form(
"Number of not isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4441 nptbins,ptmin,ptmax);
4447 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4448 nptbins,ptmin,ptmax);
4454 fhTimeENoCut =
new TH2F (
"hTimeE_NoCut",
"time of cluster vs E of clusters, no cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4459 fhTimeESPD =
new TH2F (
"hTimeE_SPD",
"time of cluster vs E of clusters, SPD cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4464 fhTimeESPDMulti =
new TH2F (
"hTimeE_SPDMulti",
"time of cluster vs E of clusters, SPD multi cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4469 fhTimeNPileUpVertSPD =
new TH2F (
"hTime_NPileUpVertSPD",
"time of cluster vs N pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4474 fhTimeNPileUpVertTrack =
new TH2F (
"hTime_NPileUpVertTracks",
"time of cluster vs N pile-up Tracks vertex", ntimebins,timemin,timemax, 50,0,50 );
4479 fhTimeNPileUpVertContributors =
new TH2F (
"hTime_NPileUpVertContributors",
"time of cluster vs N constributors to pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4484 fhTimePileUpMainVertexZDistance =
new TH2F (
"hTime_PileUpMainVertexZDistance",
"time of cluster vs distance in Z pile-up SPD vertex - main SPD vertex",ntimebins,timemin,timemax,100,0,50);
4489 fhTimePileUpMainVertexZDiamond =
new TH2F (
"hTime_PileUpMainVertexZDiamond",
"time of cluster vs distance in Z pile-up SPD vertex - z diamond",ntimebins,timemin,timemax,100,0,50);
4495 return outputContainer ;