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),
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 fMomentum(), fMomIso(),
74 fMomDaugh1(), fMomDaugh2(),
75 fTrackVector(), fProdVertex(),
76 fCluster(0), fClustersArr(0),
78 fhEIso(0), fhPtIso(0),
79 fhPtCentralityIso(0), fhPtEventPlaneIso(0),
81 fhPhiIso(0), fhEtaIso(0), fhEtaPhiIso(0),
83 fhENoIso(0), fhPtNoIso(0), fhPtNLocMaxNoIso(0),
85 fhPtClusterInCone(0), fhPtCellInCone(0), fhPtTrackInCone(0),
86 fhPtTrackInConeOtherBC(0), fhPtTrackInConeOtherBCPileUpSPD(0),
87 fhPtTrackInConeBC0(0), fhPtTrackInConeVtxBC0(0),
88 fhPtTrackInConeBC0PileUpSPD(0),
89 fhPtInConePileUp(), fhPtInConeCent(0),
90 fhPerpConeSumPt(0), fhPtInPerpCone(0),
91 fhEtaPhiInConeCluster(0), fhEtaPhiCluster(0),
92 fhEtaPhiInConeTrack(0), fhEtaPhiTrack(0),
93 fhEtaBandCluster(0), fhPhiBandCluster(0),
94 fhEtaBandTrack(0), fhPhiBandTrack(0),
95 fhEtaBandCell(0), fhPhiBandCell(0),
96 fhConePtLead(0), fhConePtLeadCluster(0), fhConePtLeadTrack(0),
97 fhConePtLeadClustervsTrack(0), fhConePtLeadClusterTrackFrac(0),
98 fhConeSumPt(0), fhPtLambda0Eiso(0), fhConeSumPtCellTrack(0),
99 fhConeSumPtCell(0), fhConeSumPtCluster(0), fhConeSumPtTrack(0),
100 fhConeSumPtEtaBandUECluster(0), fhConeSumPtPhiBandUECluster(0),
101 fhConeSumPtEtaBandUETrack(0), fhConeSumPtPhiBandUETrack(0),
102 fhConeSumPtEtaBandUECell(0), fhConeSumPtPhiBandUECell(0),
103 fhConeSumPtTrigEtaPhi(0),
104 fhConeSumPtCellTrackTrigEtaPhi(0),
105 fhConeSumPtEtaBandUEClusterTrigEtaPhi(0), fhConeSumPtPhiBandUEClusterTrigEtaPhi(0),
106 fhConeSumPtEtaBandUETrackTrigEtaPhi(0), fhConeSumPtPhiBandUETrackTrigEtaPhi(0),
107 fhConeSumPtEtaBandUECellTrigEtaPhi(0), fhConeSumPtPhiBandUECellTrigEtaPhi(0),
108 fhConeSumPtEtaUESub(0), fhConeSumPtPhiUESub(0),
109 fhConeSumPtEtaUESubTrigEtaPhi(0), fhConeSumPtPhiUESubTrigEtaPhi(0),
110 fhConeSumPtEtaUESubTrackCell(0), fhConeSumPtPhiUESubTrackCell(0),
111 fhConeSumPtEtaUESubTrackCellTrigEtaPhi(0), fhConeSumPtPhiUESubTrackCellTrigEtaPhi(0),
112 fhConeSumPtEtaUESubCluster(0), fhConeSumPtPhiUESubCluster(0),
113 fhConeSumPtEtaUESubClusterTrigEtaPhi(0), fhConeSumPtPhiUESubClusterTrigEtaPhi(0),
114 fhConeSumPtEtaUESubCell(0), fhConeSumPtPhiUESubCell(0),
115 fhConeSumPtEtaUESubCellTrigEtaPhi(0), fhConeSumPtPhiUESubCellTrigEtaPhi(0),
116 fhConeSumPtEtaUESubTrack(0), fhConeSumPtPhiUESubTrack(0),
117 fhConeSumPtEtaUESubTrackTrigEtaPhi(0), fhConeSumPtPhiUESubTrackTrigEtaPhi(0),
118 fhFractionTrackOutConeEta(0), fhFractionTrackOutConeEtaTrigEtaPhi(0),
119 fhFractionClusterOutConeEta(0), fhFractionClusterOutConeEtaTrigEtaPhi(0),
120 fhFractionClusterOutConePhi(0), fhFractionClusterOutConePhiTrigEtaPhi(0),
121 fhFractionCellOutConeEta(0), fhFractionCellOutConeEtaTrigEtaPhi(0),
122 fhFractionCellOutConePhi(0), fhFractionCellOutConePhiTrigEtaPhi(0),
123 fhConeSumPtClustervsTrack(0), fhConeSumPtClusterTrackFrac(0),
124 fhConeSumPtEtaUESubClustervsTrack(0), fhConeSumPtPhiUESubClustervsTrack(0),
125 fhConeSumPtCellvsTrack(0),
126 fhConeSumPtEtaUESubCellvsTrack(0), fhConeSumPtPhiUESubCellvsTrack(0),
127 fhEtaBandClustervsTrack(0), fhPhiBandClustervsTrack(0),
128 fhEtaBandNormClustervsTrack(0), fhPhiBandNormClustervsTrack(0),
129 fhEtaBandCellvsTrack(0), fhPhiBandCellvsTrack(0),
130 fhEtaBandNormCellvsTrack(0), fhPhiBandNormCellvsTrack(0),
131 fhConeSumPtSubvsConeSumPtTotPhiTrack(0), fhConeSumPtSubNormvsConeSumPtTotPhiTrack(0),
132 fhConeSumPtSubvsConeSumPtTotEtaTrack(0), fhConeSumPtSubNormvsConeSumPtTotEtaTrack(0),
133 fhConeSumPtSubvsConeSumPtTotPhiCluster(0), fhConeSumPtSubNormvsConeSumPtTotPhiCluster(0),
134 fhConeSumPtSubvsConeSumPtTotEtaCluster(0), fhConeSumPtSubNormvsConeSumPtTotEtaCluster(0),
135 fhConeSumPtSubvsConeSumPtTotPhiCell(0), fhConeSumPtSubNormvsConeSumPtTotPhiCell(0),
136 fhConeSumPtSubvsConeSumPtTotEtaCell(0), fhConeSumPtSubNormvsConeSumPtTotEtaCell(0),
137 fhConeSumPtVSUETracksEtaBand(0), fhConeSumPtVSUETracksPhiBand(0),
138 fhConeSumPtVSUEClusterEtaBand(0), fhConeSumPtVSUEClusterPhiBand(0),
139 fhPtPrimMCPi0DecayPairOutOfCone(0),
140 fhPtPrimMCPi0DecayPairOutOfAcceptance(0),
141 fhPtPrimMCPi0DecayPairOutOfAcceptanceNoOverlap(0),
142 fhPtPrimMCPi0DecayPairAcceptInConeLowPt(0),
143 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlap(0),
144 fhPtPrimMCPi0DecayPairAcceptInConeLowPtNoOverlapCaloE(0),
145 fhPtPrimMCPi0DecayPairNoOverlap(0),
146 fhPtPrimMCPi0DecayIsoPairOutOfCone(0),
147 fhPtPrimMCPi0DecayIsoPairOutOfAcceptance(0),
148 fhPtPrimMCPi0DecayIsoPairOutOfAcceptanceNoOverlap(0),
149 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPt(0),
150 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlap(0),
151 fhPtPrimMCPi0DecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
152 fhPtPrimMCPi0DecayIsoPairNoOverlap(0),
153 fhPtPrimMCPi0Overlap(0), fhPtPrimMCPi0IsoOverlap(0),
154 fhPtPrimMCEtaDecayPairOutOfCone(0),
155 fhPtPrimMCEtaDecayPairOutOfAcceptance(0),
156 fhPtPrimMCEtaDecayPairOutOfAcceptanceNoOverlap(0),
157 fhPtPrimMCEtaDecayPairAcceptInConeLowPt(0),
158 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlap(0),
159 fhPtPrimMCEtaDecayPairAcceptInConeLowPtNoOverlapCaloE(0),
160 fhPtPrimMCEtaDecayPairNoOverlap(0),
161 fhPtPrimMCEtaDecayIsoPairOutOfCone(0),
162 fhPtPrimMCEtaDecayIsoPairOutOfAcceptance(0),
163 fhPtPrimMCEtaDecayIsoPairOutOfAcceptanceNoOverlap(0),
164 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPt(0),
165 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlap(0),
166 fhPtPrimMCEtaDecayIsoPairAcceptInConeLowPtNoOverlapCaloE(0),
167 fhPtPrimMCEtaDecayIsoPairNoOverlap(0),
168 fhPtPrimMCEtaOverlap(0), fhPtPrimMCEtaIsoOverlap(0),
169 fhPtLeadConeBin(0), fhSumPtConeBin(0),
170 fhPtLeadConeBinMC(0), fhSumPtConeBinMC(0),
171 fhPtLeadConeBinDecay(0), fhSumPtConeBinDecay(0),
172 fhPtLeadConeBinLambda0(0), fhSumPtConeBinLambda0(0),
173 fhPtLeadConeBinLambda0MC(0), fhSumPtConeBinLambda0MC(0),
174 fhPtTrigBinPtLeadCone(0), fhPtTrigBinSumPtCone(0),
175 fhPtTrigBinPtLeadConeMC(0), fhPtTrigBinSumPtConeMC(0),
176 fhPtTrigBinPtLeadConeDecay(0), fhPtTrigBinSumPtConeDecay(0),
177 fhPtTrigBinLambda0vsPtLeadCone(0), fhPtTrigBinLambda0vsSumPtCone(0),
178 fhPtTrigBinLambda0vsPtLeadConeMC(0), fhPtTrigBinLambda0vsSumPtConeMC(0),
181 fhELambda0LocMax1(), fhELambda1LocMax1(),
182 fhELambda0LocMax2(), fhELambda1LocMax2(),
183 fhELambda0LocMaxN(), fhELambda1LocMaxN(),
185 fhEIsoPileUp(), fhPtIsoPileUp(),
186 fhENoIsoPileUp(), fhPtNoIsoPileUp(),
187 fhTimeENoCut(0), fhTimeESPD(0), fhTimeESPDMulti(0),
188 fhTimeNPileUpVertSPD(0), fhTimeNPileUpVertTrack(0),
189 fhTimeNPileUpVertContributors(0),
190 fhTimePileUpMainVertexZDistance(0), fhTimePileUpMainVertexZDiamond(0)
194 for(
Int_t i = 0; i < 5 ; i++)
201 for(
Int_t j = 0; j < 5 ; j++)
236 for(
Int_t iso =0; iso < 2; iso++)
246 for(
Int_t i = 0; i < 5 ; i++)
265 for(
Int_t i = 0; i < 2 ; i++)
281 for(
Int_t i = 0; i < 2 ; i++)
297 for(
Int_t ieta = 0; ieta < 4; ieta++)
299 for(
Int_t iphi = 0; iphi < 3; iphi++)
310 for(
Int_t ism =0; ism < 20; ism++)
331 for(
Int_t i = 0 ; i < 7 ; i++)
365 Float_t ptTrig = pCandidate->Pt() ;
366 Float_t phiTrig = pCandidate->Phi();
367 Float_t etaTrig = pCandidate->Eta();
369 for(
Int_t icluster=0; icluster < pl->GetEntriesFast(); icluster++)
371 AliVCluster* cluster = (AliVCluster *) pl->At(icluster);
375 AliWarning(
"Cluster not available?");
380 if(cluster->GetID() == pCandidate->GetCaloLabel(0) ||
381 cluster->GetID() == pCandidate->GetCaloLabel(1) )
continue ;
433 Float_t phiTrig = pCandidate->Phi();
434 if(phiTrig<0) phiTrig += TMath::TwoPi();
436 Float_t etaTrig = pCandidate->Eta();
438 if(pCandidate->GetDetectorTag()==
kEMCAL)
440 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
443 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
445 if(!eGeom->CheckAbsCellId(absId))
return ;
448 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
450 Int_t imEta=-1, imPhi=-1;
451 Int_t ieta =-1, iphi =-1;
453 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
455 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
457 Int_t colTrig = ieta;
458 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + ieta ;
460 Int_t rowTrig = iphi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
462 Int_t sqrSize = int(conesize/0.0143);
466 Int_t nTotalRows = AliEMCALGeoParams::fgkEMCALRows*16/3 ;
467 Int_t nTotalCols = 2*AliEMCALGeoParams::fgkEMCALCols;
471 Int_t irowmin = rowTrig-sqrSize;
472 if ( irowmin < 0 ) irowmin = 0 ;
474 Int_t irowmax = rowTrig+sqrSize;
475 if ( irowmax > AliEMCALGeoParams::fgkEMCALRows ) irowmax = AliEMCALGeoParams::fgkEMCALRows;
477 for(
Int_t irow = irowmin; irow <irowmax; irow++)
479 for(
Int_t icol = 0; icol < nTotalCols; icol++)
481 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
482 if(inSector==5)
continue;
485 if(icol < AliEMCALGeoParams::fgkEMCALCols)
487 inSupMod = 2*inSector + 1;
490 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
492 inSupMod = 2*inSector;
493 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
496 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
499 if(TMath::Abs(icol-colTrig) < sqrSize || TMath::Abs(irow-rowTrig) < sqrSize){
503 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
504 if(!eGeom->CheckAbsCellId(iabsId))
continue;
506 etaBandPtSumCells += cells->GetCellAmplitude(iabsId);
514 Int_t icolmin = colTrig-sqrSize;
515 if ( icolmin < 0 ) icolmin = 0;
517 Int_t icolmax = colTrig+sqrSize;
518 if ( icolmax > AliEMCALGeoParams::fgkEMCALCols ) icolmax = AliEMCALGeoParams::fgkEMCALCols;
521 for(
Int_t icol = icolmin; icol < icolmax; icol++)
523 for(
Int_t irow = 0; irow < nTotalRows; irow++)
525 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
526 if ( inSector == 5 ) continue ;
532 if(icol < AliEMCALGeoParams::fgkEMCALCols)
535 inSupMod = 2*inSector + 1;
538 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
541 inSupMod = 2*inSector;
542 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
545 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
548 if(TMath::Abs(icol-colTrig) < sqrSize) {
551 if(TMath::Abs(irow-rowTrig) < sqrSize) {
555 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
556 if(!eGeom->CheckAbsCellId(iabsId))
558 AliWarning(Form(
"!eGeom->CheckAbsCellId(iabsId=%i) inSupMod %i irowLoc %i icolLoc %i",iabsId,inSupMod, irowLoc, icolLoc));
562 phiBandPtSumCells += cells->GetCellAmplitude(iabsId);
572 Float_t ptTrig = pCandidate->Pt();
591 Float_t ptTrig = pCandidate->Pt() ;
592 Float_t phiTrig = pCandidate->Phi();
593 Float_t etaTrig = pCandidate->Eta();
596 for(
Int_t itrack=0; itrack < trackList->GetEntriesFast(); itrack++)
598 AliVTrack* track = (AliVTrack *) trackList->At(itrack);
602 AliWarning(
"Track not available?");
609 if ( pCandidate->GetDetectorTag() ==
kCTS )
612 Bool_t contained = kFALSE;
614 for(
Int_t i = 0; i < 4; i++)
616 if( trackID == pCandidate->GetTrackLabel(i) ) contained = kTRUE;
619 if ( contained ) continue ;
635 if(track->Eta() > (etaTrig-conesize) && track->Eta() < (etaTrig+conesize))
637 phiBandPtSum+=track->Pt();
642 if(track->Phi() > (phiTrig-conesize) && track->Phi() < (phiTrig+conesize))
644 etaBandPtSum+=track->Pt();
649 Double_t pTrack = TMath::Sqrt(track->Px()*track->Px()+track->Py()*track->Py());
650 Double_t dPhi = phiTrig - track->Phi() + TMath::PiOver2();
651 Double_t dEta = etaTrig - track->Eta();
652 Double_t arg = dPhi*dPhi + dEta*dEta;
654 if(TMath::Sqrt(arg) < conesize)
657 sumptPerp+=track->Pt();
660 dPhi = phiTrig - track->Phi() - TMath::PiOver2();
661 arg = dPhi*dPhi + dEta*dEta;
663 if(TMath::Sqrt(arg) < conesize)
666 sumptPerp+=track->Pt();
686 Float_t etaUEptsumCluster = 0 ;
687 Float_t phiUEptsumCluster = 0 ;
696 Float_t coneA = conesize*conesize*TMath::Pi();
697 Float_t ptTrig = pCandidate->Pt() ;
698 Float_t phiTrig = pCandidate->Phi();
699 Float_t etaTrig = pCandidate->Eta();
705 Float_t phiUEptsumTrackNorm = 0 ;
706 Float_t etaUEptsumTrackNorm = 0 ;
707 Float_t coneptsumTrackSubPhi = 0 ;
708 Float_t coneptsumTrackSubEta = 0 ;
709 Float_t coneptsumTrackSubPhiNorm = 0 ;
710 Float_t coneptsumTrackSubEtaNorm = 0 ;
711 etaBandptsumTrackNorm = 0 ;
723 Float_t correctConeSumTrack = 1;
724 Float_t correctConeSumTrackPhi = 1;
727 phiUEptsumTrack,etaUEptsumTrack,
728 phiUEptsumTrackNorm,etaUEptsumTrackNorm,
729 correctConeSumTrack,correctConeSumTrackPhi);
731 coneptsumTrackSubPhi = coneptsumTrack - phiUEptsumTrackNorm;
732 coneptsumTrackSubEta = coneptsumTrack - etaUEptsumTrackNorm;
734 etaBandptsumTrackNorm = etaUEptsumTrackNorm;
744 if(coneptsumTrack > 0)
746 coneptsumTrackSubPhiNorm = coneptsumTrackSubPhi/coneptsumTrack;
747 coneptsumTrackSubEtaNorm = coneptsumTrackSubEta/coneptsumTrack;
759 Float_t phiUEptsumClusterNorm = 0 ;
760 Float_t etaUEptsumClusterNorm = 0 ;
761 Float_t coneptsumClusterSubPhi = 0 ;
762 Float_t coneptsumClusterSubEta = 0 ;
763 Float_t coneptsumClusterSubPhiNorm = 0 ;
764 Float_t coneptsumClusterSubEtaNorm = 0 ;
765 Float_t phiUEptsumCellNorm = 0 ;
766 Float_t etaUEptsumCellNorm = 0 ;
767 Float_t coneptsumCellSubPhi = 0 ;
768 Float_t coneptsumCellSubEta = 0 ;
769 Float_t coneptsumCellSubPhiNorm = 0 ;
770 Float_t coneptsumCellSubEtaNorm = 0 ;
771 etaBandptsumClusterNorm = 0;
786 Float_t correctConeSumClusterEta = 1;
787 Float_t correctConeSumClusterPhi = 1;
790 phiUEptsumCluster,etaUEptsumCluster,
791 phiUEptsumClusterNorm,etaUEptsumClusterNorm,
792 correctConeSumClusterEta,correctConeSumClusterPhi);
803 coneptsumClusterSubPhi = coneptsumCluster - phiUEptsumClusterNorm;
804 coneptsumClusterSubEta = coneptsumCluster - etaUEptsumClusterNorm;
806 etaBandptsumClusterNorm = etaUEptsumClusterNorm;
818 if(coneptsumCluster!=0)
820 coneptsumClusterSubPhiNorm = coneptsumClusterSubPhi/coneptsumCluster;
821 coneptsumClusterSubEtaNorm = coneptsumClusterSubEta/coneptsumCluster;
842 Float_t emcPhiSize = TMath::DegToRad()*100.;
844 if(((2*conesize*emcPhiSize)-coneA)!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*conesize*emcPhiSize)-coneA));
845 if(((2*conesize*emcEtaSize)-coneA)!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*conesize*emcEtaSize)-coneA));
849 Float_t correctConeSumCellEta = 1;
850 if(TMath::Abs(etaTrig)+conesize > emcEtaSize/2.)
852 Float_t excess = TMath::Abs(etaTrig) + conesize - emcEtaSize/2.;
856 if(((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta))!=0)phiUEptsumCellNorm = phiUEptsumCell*(coneA / ((2*(conesize-excess)*emcPhiSize)-(coneA-correctConeSumCellEta)));
859 Float_t correctConeSumCellPhi = 1;
861 if((phiTrig+conesize > 180*TMath::DegToRad()) ||
862 (phiTrig-conesize < 80*TMath::DegToRad()))
865 if( phiTrig+conesize > 180*TMath::DegToRad() ) excess = conesize + phiTrig - 180*TMath::DegToRad() ;
866 else excess = conesize - phiTrig + 80*TMath::DegToRad() ;
872 if(((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi))!=0)etaUEptsumCellNorm = etaUEptsumCell*(coneA / ((2*(conesize-excess)*emcEtaSize)-(coneA-correctConeSumCellPhi)));
876 coneptsumCellSubPhi = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - phiUEptsumCellNorm;
877 coneptsumCellSubEta = coneptsumCell*correctConeSumCellEta*correctConeSumCellPhi - etaUEptsumCellNorm;
888 if ( coneptsumCell > 0.01 )
890 coneptsumCellSubPhiNorm = coneptsumCellSubPhi/coneptsumCell;
891 coneptsumCellSubEtaNorm = coneptsumCellSubEta/coneptsumCell;
907 Double_t sumPhiUESub = coneptsumClusterSubPhi + coneptsumTrackSubPhi;
908 Double_t sumEtaUESub = coneptsumClusterSubEta + coneptsumTrackSubEta;
929 Double_t sumPhiUESubTrackCell = coneptsumCellSubPhi + coneptsumTrackSubPhi;
930 Double_t sumEtaUESubTrackCell = coneptsumCellSubEta + coneptsumTrackSubEta;
954 coneptLeadCluster = 0;
955 coneptsumCluster = 0;
959 Float_t ptTrig = aodParticle->Pt();
976 for(
Int_t icalo=0; icalo < refclusters->GetEntriesFast(); icalo++)
978 AliVCluster* calo = (AliVCluster *) refclusters->At(icalo);
999 coneptsumCluster+=ptcone;
1000 if(ptcone > coneptLeadCluster) coneptLeadCluster = ptcone;
1006 aodParticle->SetNeutralLeadPtInCone(coneptLeadCluster);
1007 aodParticle->SetNeutralPtSumInCone(coneptsumCluster);
1021 Float_t ptTrig = aodParticle->Pt();
1022 Float_t phiTrig = aodParticle->Phi();
1023 if(phiTrig<0) phiTrig += TMath::TwoPi();
1024 Float_t etaTrig = aodParticle->Eta();
1026 if(aodParticle->GetDetectorTag()==
kEMCAL)
1028 AliEMCALGeometry* eGeom = AliEMCALGeometry::GetInstance();
1031 if (eGeom->GetAbsCellIdFromEtaPhi(etaTrig,phiTrig,absId))
1033 if(!eGeom->CheckAbsCellId(absId))
return ;
1036 Int_t nSupMod = eGeom->GetSuperModuleNumber(absId);
1038 Int_t imEta=-1, imPhi=-1;
1039 Int_t ieta =-1, iphi =-1;
1041 if (eGeom->GetCellIndex(absId,nSupMod,nModule,imPhi,imEta))
1043 Int_t iEta=-1, iPhi=-1;
1044 eGeom->GetCellPhiEtaIndexInSModule(nSupMod,nModule,imPhi,imEta,iphi,ieta);
1046 Int_t colTrig = iEta;
1047 if (nSupMod % 2) colTrig = AliEMCALGeoParams::fgkEMCALCols + iEta ;
1048 Int_t rowTrig = iPhi + AliEMCALGeoParams::fgkEMCALRows*int(nSupMod/2);
1050 Int_t sqrSize = int(conesize/0.0143);
1055 for(
Int_t irow = rowTrig-sqrSize; irow < rowTrig+sqrSize; irow++)
1057 for(
Int_t icol = colTrig-sqrSize; icol < colTrig+sqrSize; icol++)
1059 Int_t inSector = int(irow/AliEMCALGeoParams::fgkEMCALRows);
1060 if(inSector==5)
continue;
1062 Int_t inSupMod = -1;
1064 if(icol < AliEMCALGeoParams::fgkEMCALCols)
1066 inSupMod = 2*inSector + 1;
1069 else if(icol > AliEMCALGeoParams::fgkEMCALCols - 1)
1071 inSupMod = 2*inSector;
1072 icolLoc = icol-AliEMCALGeoParams::fgkEMCALCols;
1075 Int_t irowLoc = irow - AliEMCALGeoParams::fgkEMCALRows*inSector ;
1077 Int_t iabsId = eGeom->GetAbsCellIdFromCellIndexes(inSupMod,irowLoc,icolLoc);
1078 if(!eGeom->CheckAbsCellId(iabsId))
continue;
1081 coneptsumCell += cells->GetCellAmplitude(iabsId);
1097 coneptLeadTrack = 0;
1101 Float_t ptTrig = aodParticle->Pt();
1114 for(
Int_t itrack=0; itrack < reftracks->GetEntriesFast(); itrack++)
1116 AliVTrack* track = (AliVTrack *) reftracks->At(itrack);
1118 pTtrack = track->Pt();
1125 ULong_t status = track->GetStatus();
1126 Bool_t okTOF = ( (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ) ;
1128 Int_t trackBC = track->GetTOFBunchCrossing(bz);
1149 coneptsumTrack+=pTtrack;
1150 if(pTtrack > coneptLeadTrack) coneptLeadTrack = pTtrack;
1156 aodParticle->SetChargedLeadPtInCone(coneptLeadTrack);
1157 aodParticle->SetChargedPtSumInCone(coneptsumTrack);
1193 if(energy < 8)
return;
1199 Int_t nVerticesSPD = -1;
1200 Int_t nVerticesTracks = -1;
1204 nVerticesSPD = esdEv->GetNumberOfPileupVerticesSPD();
1205 nVerticesTracks = esdEv->GetNumberOfPileupVerticesTracks();
1210 nVerticesSPD = aodEv->GetNumberOfPileupVerticesSPD();
1211 nVerticesTracks = aodEv->GetNumberOfPileupVerticesTracks();
1223 for(
Int_t iVert=0; iVert<nVerticesSPD;iVert++)
1227 const AliESDVertex* pv=esdEv->GetPileupVertexSPD(iVert);
1228 ncont=pv->GetNContributors();
1229 z1 = esdEv->GetPrimaryVertexSPD()->GetZ();
1231 diamZ = esdEv->GetDiamondZ();
1235 AliAODVertex *pv=aodEv->GetVertex(iVert);
1236 if(pv->GetType()!=AliAODVertex::kPileupSPD)
continue;
1237 ncont=pv->GetNContributors();
1238 z1=aodEv->GetPrimaryVertexSPD()->GetZ();
1240 diamZ = aodEv->GetDiamondZ();
1243 Double_t distZ = TMath::Abs(z2-z1);
1244 diamZ = TMath::Abs(z2-diamZ);
1262 Int_t clusterID = pCandidate->GetCaloLabel(0) ;
1263 Int_t nMaxima = pCandidate->GetNLM();
1264 Int_t mcTag = pCandidate->GetTag() ;
1265 Bool_t isolated = pCandidate->IsIsolated();
1267 if ( clusterID < 0 )
1269 AliWarning(Form(
"ID of cluster = %d, not possible!", clusterID));
1273 Float_t m02 = pCandidate->GetM02() ;
1275 Float_t pt = pCandidate->Pt();
1276 Float_t eta = pCandidate->Eta();
1277 Float_t phi = pCandidate->Phi();
1278 if(phi<0) phi+= TMath::TwoPi();
1297 Int_t decayTag = pCandidate->DecayTag();
1298 if(decayTag < 0) decayTag = 0;
1335 Int_t ptsumBin = -1;
1336 Int_t leadptBin = -1;
1338 AliDebug(1,Form(
"pT cand: %2.2f, In cone pT: Sum %2.2f, Lead %2.2f, n bins %d",pt,coneptsum,coneleadpt,
fNBkgBin));
1348 if ( leadptBin >= 0 )
1357 if ( leadptBin == 0 )
1358 AliDebug(1,Form(
"No track/clusters in isolation cone: cand pt %2.2f GeV/c, track multiplicity %d, N clusters %d",
1362 if ( ptsumBin >= 0 )
1374 Int_t decayTag = pCandidate->DecayTag();
1375 if(decayTag < 0) decayTag = 0;
1437 AliFatal(Form(
"Lost decay Bit assigned to bad case, mcIndex %d",mcIndex));
1459 Int_t ptTrigBin = -1;
1468 if ( ptTrigBin >= 0 )
1485 Int_t decayTag = pCandidate->DecayTag();
1486 if(decayTag < 0) decayTag = 0;
1589 Int_t noverlaps = 0;
1593 Int_t overpdg[nlabels];
1622 if ( noverlaps == 1 )
1649 else if (noverlaps > 1 )
1735 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
1741 Float_t dEdx = track->GetTPCsignal();
1781 const Int_t buffersize = 255;
1782 char onePar[buffersize] ;
1784 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
1788 snprintf(onePar, buffersize,
"Isolation Cand Detector: %s;",
fIsoDetectorString.Data()) ;
1790 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
1792 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
1794 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
1796 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
1801 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
1803 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
1808 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
1813 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
1818 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
1823 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
1834 return new TObjString(parList) ;
1844 outputContainer->SetName(
"IsolatedParticleHistos") ;
1893 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
1896 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
1901 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
1904 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
1908 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
1911 TString sParticle =
", x^{0,#pm}";
1915 TString parTitle = Form(
"#it{R} = %2.2f%s%s",
GetIsolationCut()->GetConeSize(), sThreshold.Data(),sParticle.Data());
1917 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
1920 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
1921 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
1922 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
1923 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
1925 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
1926 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
1927 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
1928 "OtherDecay",
"Electron",
"Hadron"} ;
1931 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
1932 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
1933 "#pi^{0}" ,
"#eta"} ;
1935 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
1936 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
1942 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
1943 nptbins,ptmin,ptmax);
1944 fhENoIso->SetYTitle(
"#it{counts}");
1945 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
1949 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
1950 nptbins,ptmin,ptmax);
1952 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
1956 Form(
"Number of not isolated leading particles #eta vs #phi, %s",parTitle.Data()),
1957 netabins,etamin,etamax,nphibins,phimin,phimax);
1968 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
1969 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1970 nptbins,ptmin,ptmax);
1972 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1975 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
1976 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1977 nptbins,ptmin,ptmax);
1978 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
1979 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1983 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1984 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
1986 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1990 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1991 nptbins,ptmin,ptmax,netabins,etamin,etamax);
1993 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2001 TString isoName [] = {
"NoIso",
"Iso"};
2002 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2006 for(
Int_t iso = 0; iso < 2; iso++)
2010 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2011 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2012 nptbins,ptmin,ptmax);
2013 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2014 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2015 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2019 Form(
"Number of %s leading Pi0 decay particles #eta vs #phi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2020 netabins,etamin,etamax,nphibins,phimin,phimax);
2028 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2029 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2030 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2041 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2042 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2043 nptbins,ptmin,ptmax);
2044 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2045 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2046 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2055 TString isoName [] = {
"NoIso",
"Iso"};
2056 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2059 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2060 nptbins,ptmin,ptmax);
2061 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2062 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2063 outputContainer->Add(
fhEIso) ;
2066 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2067 nptbins,ptmin,ptmax);
2068 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
2069 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2070 outputContainer->Add(
fhPtIso) ;
2073 Form(
"Number of isolated particles vs #phi, %s",parTitle.Data()),
2074 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2076 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2080 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
2081 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2083 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2087 Form(
"Number of isolated particles #eta vs #phi, %s",parTitle.Data()),
2088 netabins,etamin,etamax,nphibins,phimin,phimax);
2096 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2097 nptbins,ptmin,ptmax, 100,0,100);
2103 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2104 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
2113 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2114 nptbins,ptmin,ptmax,10,0,10);
2119 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2120 nptbins,ptmin,ptmax,10,0,10);
2127 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2128 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2129 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
2130 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2134 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2135 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2137 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2141 (Form(
"hPtLambda0Eiso"),
2142 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.);
2149 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2150 netabins,etamin,etamax,nphibins,phimin,phimax);
2157 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for #it{R} = %2.2f",r),
2158 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2159 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
2160 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2194 (Form(
"hPtLeadCone_Bin%d",ibin),
2195 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2202 (Form(
"hSumPtCone_Bin%d",ibin),
2203 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2216 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2217 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2224 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2225 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2239 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2240 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2247 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2248 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2259 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
2260 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2267 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
2268 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2280 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2281 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2288 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2289 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2329 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
2330 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2337 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
2338 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2351 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2352 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2359 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2360 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2374 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2375 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2382 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2383 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2394 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
2395 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2402 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
2403 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2415 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2416 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2423 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2424 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2438 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2439 100,0,100,nptinconebins,ptinconemin,ptinconemax);
2449 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2450 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2456 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
2457 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2466 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2467 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2476 "#Sigma cluster #it{p}_{T} in UE Eta Band",
2477 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2483 "#Sigma cluster #it{p}_{T} UE Phi Band",
2484 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2490 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} in UE Eta Band",
2491 netabins,etamin,etamax,nphibins,phimin,phimax);
2498 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} UE Phi Band",
2499 netabins,etamin,etamax,nphibins,phimin,phimax);
2508 "#Sigma cell #it{p}_{T} in UE Eta Band",
2509 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2515 "#Sigma cell #it{p}_{T} UE Phi Band",
2516 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2522 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} in UE Eta Band",
2523 netabins,etamin,etamax,nphibins,phimin,phimax);
2530 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} UE Phi Band",
2531 netabins,etamin,etamax,nphibins,phimin,phimax);
2539 Form(
"#eta vs #phi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
2540 netabins,-1,1,nphibins,0,TMath::TwoPi());
2546 Form(
"#eta vs #phi of clusters in #phi band isolation cone for #it{R} = %2.2f",r),
2547 netabins,-1,1,nphibins,0,TMath::TwoPi());
2553 Form(
"#eta vs #phi of clusters in cone for #it{R} = %2.2f",r),
2554 netabins,-1,1,nphibins,0,TMath::TwoPi());
2560 Form(
"#eta vs #phi of all clusters"),
2561 netabins,-1,1,nphibins,0,TMath::TwoPi());
2569 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
2570 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2578 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
2579 nptbins,ptmin,ptmax,1000,0,50);
2585 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
2592 Form(
"#col vs #row of cells in #phi band isolation cone for #it{R} = %2.2f",r),
2602 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2603 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2609 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2610 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2616 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),
2617 netabins,etamin,etamax,nphibins,phimin,phimax);
2624 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),
2625 netabins,etamin,etamax,nphibins,phimin,phimax);
2634 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2635 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2641 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2642 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2648 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),
2649 netabins,etamin,etamax,nphibins,phimin,phimax);
2656 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),
2657 netabins,etamin,etamax,nphibins,phimin,phimax);
2665 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
2666 nptbins,ptmin,ptmax,100,0,1);
2672 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#phi ",r),
2673 netabins,etamin,etamax,nphibins,phimin,phimax);
2680 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance",r),
2681 nptbins,ptmin,ptmax,100,0,1);
2687 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance, in trigger #eta-#phi ",r),
2688 netabins,etamin,etamax,nphibins,phimin,phimax);
2695 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),
2696 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2702 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),
2703 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2709 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),
2710 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2716 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),
2717 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2723 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for cluster (before normalization), R=%2.2f",r),
2724 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
2730 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for cluster (before normalization), R=%2.2f",r),
2731 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
2739 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
2740 nptbins,ptmin,ptmax,100,0,1);
2746 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#phi ",r),
2747 netabins,etamin,etamax,nphibins,phimin,phimax);
2754 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance",r),
2755 nptbins,ptmin,ptmax,100,0,1);
2761 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance, in trigger #eta-#phi ",r),
2762 netabins,etamin,etamax,nphibins,phimin,phimax);
2770 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),
2771 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2777 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),
2778 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2784 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),
2785 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2791 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),
2792 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2804 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2805 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2811 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
2812 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2818 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
2819 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2828 "#Sigma track #it{p}_{T} in UE Eta Band",
2829 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2835 "#Sigma track #it{p}_{T} in UE Phi Band",
2836 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax*8);
2843 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Eta Band",
2844 netabins,etamin,etamax,nphibins,phimin,phimax);
2851 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Phi Band",
2852 netabins,etamin,etamax,nphibins,phimin,phimax);
2859 Form(
"#eta vs #phi of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
2860 netabins,-1,1,nphibins,0,TMath::TwoPi());
2866 Form(
"#eta vs #phi of tracks in #phi band isolation cone for #it{R} = %2.2f",r),
2867 netabins,-1,1,nphibins,0,TMath::TwoPi());
2873 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2874 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2880 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2881 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2887 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),
2888 netabins,etamin,etamax,nphibins,phimin,phimax);
2895 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),
2896 netabins,etamin,etamax,nphibins,phimin,phimax);
2903 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance",r),
2904 nptbins,ptmin,ptmax,100,0,1);
2910 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance, in trigger #eta-#phi ",r),
2911 netabins,etamin,etamax,nphibins,phimin,phimax);
2918 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),
2919 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2925 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),
2926 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2932 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),
2933 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2939 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),
2940 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2948 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
2949 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2955 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
2956 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2962 Form(
"#eta vs #phi of all Tracks"),
2963 netabins,-1,1,nphibins,0,TMath::TwoPi());
2969 Form(
"#eta vs #phi of Tracks in cone for #it{R} = %2.2f",r),
2970 netabins,-1,1,nphibins,0,TMath::TwoPi());
2976 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for tracks (before normalization), R=%2.2f",r),
2977 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
2983 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for tracks (before normalization), R=%2.2f",r),
2984 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
2994 Form(
"Track vs Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2995 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3001 Form(
"#Sigma #it{p}_{T}^{cluster}/#Sigma #it{p}_{T}^{track} in isolation cone for #it{R} = %2.2f",r),
3002 nptbins,ptmin,ptmax,200,0,5);
3009 Form(
"Track vs Cluster lead #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3010 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3016 Form(
" #it{p}^{leading cluster}_{T}/#it{p}^{leading track}_{T} in isolation cone for #it{R} = %2.2f",r),
3017 nptbins,ptmin,ptmax,200,0,5);
3026 Form(
"Track vs cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3027 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3033 Form(
"Track and Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3034 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3040 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3041 netabins,etamin,etamax,nphibins,phimin,phimax);
3051 Form(
"#Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3052 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3058 Form(
"#Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3059 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3065 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),
3066 netabins,etamin,etamax,nphibins,phimin,phimax);
3073 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),
3074 netabins,etamin,etamax,nphibins,phimin,phimax);
3081 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3082 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3088 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3089 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3095 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3096 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3102 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3103 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3109 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3110 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3116 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3117 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3123 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3124 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3130 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3131 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3140 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3141 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3147 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3148 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3154 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3155 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3161 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3162 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3168 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3169 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3175 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3176 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3182 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3183 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3189 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3190 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3196 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),
3197 netabins,etamin,etamax,nphibins,phimin,phimax);
3204 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),
3205 netabins,etamin,etamax,nphibins,phimin,phimax);
3214 TString region[] = {
"ITS",
"TPC",
"TRD",
"TOF",
"Top EMCal",
"In EMCal"};
3216 for(
Int_t iso = 0; iso < 2; iso++)
3221 (Form(
"hTrackMatchedDEta%s",isoName[iso].
Data()),
3222 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3223 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3228 (Form(
"hTrackMatchedDPhi%s",isoName[iso].
Data()),
3229 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3230 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3235 (Form(
"hTrackMatchedDEtaDPhi%s",isoName[iso].
Data()),
3236 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3237 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3250 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3251 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3256 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3257 nptbins,ptmin,ptmax,nresetabins,resphimin,resphimax);
3262 (Form(
"hTrackMatchedDEtaDPhi%s_MC%s",isoName[iso].
Data(),mcPartName[imc].
Data()),
3263 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3264 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3276 (Form(
"hdEdx%s",isoName[iso].
Data()),
3277 Form(
"%s - Matched track <d#it{E}/d#it{x}> vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3278 nptbins,ptmin,ptmax,ndedxbins, dedxmin, dedxmax);
3279 fhdEdx[iso]->SetXTitle(
"#it{E} (GeV)");
3280 fhdEdx[iso]->SetYTitle(
"<d#it{E}/d#it{x}>");
3281 outputContainer->Add(
fhdEdx[iso]);
3284 (Form(
"hEOverP%s",isoName[iso].
Data()),
3285 Form(
"%s - Matched track #it{E}/#it{p} vs cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3286 nptbins,ptmin,ptmax,nPoverEbins,pOverEmin,pOverEmax);
3287 fhEOverP[iso]->SetXTitle(
"#it{E} (GeV)");
3288 fhEOverP[iso]->SetYTitle(
"#it{E}/#it{p}");
3289 outputContainer->Add(
fhEOverP[iso]);
3294 (Form(
"hTrackMatchedMCParticle%s",isoName[iso].
Data()),
3295 Form(
"%s - Origin of particle vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3296 nptbins,ptmin,ptmax,8,0,8);
3316 (Form(
"hELambda0%s",isoName[iso].
Data()),
3317 Form(
"%s cluster : #it{E} vs #lambda_{0}, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3318 fhELambda0[iso]->SetYTitle(
"#lambda_{0}^{2}");
3330 (Form(
"hPtLambda0%s",isoName[iso].
Data()),
3331 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, %s",isoTitle[iso].
Data(), parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3333 fhPtLambda0[iso]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3341 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3342 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3344 fhPtLambda0MC[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3348 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3349 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3357 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3358 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3364 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3365 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3371 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3372 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3378 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3379 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3386 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3387 nptbins,ptmin,ptmax,10,0,10);
3389 fhPtNOverlap[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3393 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3394 nptbins,ptmin,ptmax,10,0,10);
3406 (Form(
"hPtLambda0TRD%s",isoName[iso].
Data()),
3407 Form(
"%s cluster: #it{p}_{T} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3413 (Form(
"hELambda0TRD%s",isoName[iso].
Data()),
3414 Form(
"%s cluster: #it{E} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3430 (Form(
"hNLocMax%s",isoName[iso].
Data()),
3431 Form(
"%s - Number of local maxima in cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3432 nptbins,ptmin,ptmax,10,0,10);
3434 fhNLocMax[iso]->SetXTitle(
"#it{E} (GeV)");
3438 (Form(
"hELambda0LocMax1%s",isoName[iso].
Data()),
3439 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);
3445 (Form(
"hELambda1LocMax1%s",isoName[iso].
Data()),
3446 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);
3452 (Form(
"hELambda0LocMax2%s",isoName[iso].
Data()),
3453 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);
3459 (Form(
"hELambda1LocMax2%s",isoName[iso].
Data()),
3460 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);
3466 ( Form(
"hELambda0LocMaxN%s",isoName[iso].
Data()),
3467 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);
3473 (Form(
"hELambda1LocMaxN%s",isoName[iso].
Data()),
3474 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);
3484 for(
Int_t ieta = 0; ieta < 4; ieta++)
3486 for(
Int_t iphi = 0; iphi < 3; iphi++)
3512 new TH2F(Form(
"hLam0_%s_eta%d_phi%d_sm%d",isoName[iso].
Data(),ieta,iphi,ism),
3513 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, sm %d, region eta %d, phi %d",
3514 isoTitle[iso].
Data(),ism,ieta,iphi),
3515 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3523 Float_t ptLimit[] = {2,3,4,5,6,8,10,12};
3524 for(
Int_t ipt = 0; ipt < 7; ipt++)
3527 (Form(
"hEtaPhiLam0BinPtBin%d%s",ipt,isoName[iso].
Data()),
3528 Form(
"%s, #eta vs #phi in #it{p}_{T}=[%2.1f,%2.1f] GeV/#it{c} and #lambda^{2}_{0}=[0.3,0.4]",
3529 isoTitle[iso].
Data(),ptLimit[ipt],ptLimit[ipt+1]),
3530 netabins,etamin,etamax,nphibins,phimin,phimax);
3540 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s",
3541 isoTitle[iso].
Data(),parTitle.Data()),
3542 nptbins,ptmin,ptmax,500,0,500);
3550 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s, SM covered by TRD",
3551 isoTitle[iso].
Data(),parTitle.Data()),
3552 nptbins,ptmin,ptmax,500,0,500);
3560 for(
Int_t iR = 0; iR < 6; iR++)
3563 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s",
3564 isoTitle[iso].
Data(),region[iR].
Data()),
3565 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3573 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s, SM covered by TRD",
3574 isoTitle[iso].
Data(),region[iR].
Data()),
3575 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3625 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
3626 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3632 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0, pile-up from SPD",r),
3633 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3639 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3640 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3646 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3647 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3654 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0, pile-up from SPD",r),
3655 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3661 for (
Int_t i = 0; i < 7 ; i++)
3664 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f, from pile-up (%s)",r,pileUpName[i].
Data()),
3665 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3678 fhEPrimMC[i] =
new TH1F(Form(
"hEPrim_MC%s",ppname[i].
Data()),
3679 Form(
"primary photon %s : #it{E}, %s",pptype[i].
Data(),parTitle.Data()),
3680 nptbins,ptmin,ptmax);
3681 fhEPrimMC[i]->SetXTitle(
"#it{E} (GeV)");
3685 Form(
"primary photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3686 nptbins,ptmin,ptmax);
3687 fhPtPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3691 Form(
"primary isolated photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3692 nptbins,ptmin,ptmax);
3697 Form(
"primary photon %s : #eta vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3698 nptbins,ptmin,ptmax,200,-2,2);
3700 fhEtaPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3704 Form(
"primary photon %s : #phi vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3705 nptbins,ptmin,ptmax,200,0.,TMath::TwoPi());
3707 fhPhiPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3714 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3715 nptbins,ptmin,ptmax);
3720 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
3722 nptbins,ptmin,ptmax);
3727 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
3729 nptbins,ptmin,ptmax);
3734 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
3736 nptbins,ptmin,ptmax);
3741 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
3743 nptbins,ptmin,ptmax);
3748 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
3750 nptbins,ptmin,ptmax);
3756 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
3758 nptbins,ptmin,ptmax);
3763 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
3765 nptbins,ptmin,ptmax);
3770 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3771 nptbins,ptmin,ptmax);
3776 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
3778 nptbins,ptmin,ptmax);
3783 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3784 nptbins,ptmin,ptmax);
3789 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3790 nptbins,ptmin,ptmax);
3795 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
3797 nptbins,ptmin,ptmax);
3802 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
3804 nptbins,ptmin,ptmax);
3809 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3811 nptbins,ptmin,ptmax);
3816 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3818 nptbins,ptmin,ptmax);
3830 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3831 nptbins,ptmin,ptmax);
3836 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
3838 nptbins,ptmin,ptmax);
3843 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
3845 nptbins,ptmin,ptmax);
3850 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
3852 nptbins,ptmin,ptmax);
3857 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
3859 nptbins,ptmin,ptmax);
3864 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
3866 nptbins,ptmin,ptmax);
3872 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
3874 nptbins,ptmin,ptmax);
3879 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
3881 nptbins,ptmin,ptmax);
3886 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3887 nptbins,ptmin,ptmax);
3892 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
3894 nptbins,ptmin,ptmax);
3899 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3900 nptbins,ptmin,ptmax);
3905 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3906 nptbins,ptmin,ptmax);
3911 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
3913 nptbins,ptmin,ptmax);
3918 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
3920 nptbins,ptmin,ptmax);
3925 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3927 nptbins,ptmin,ptmax);
3932 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3934 nptbins,ptmin,ptmax);
3946 const Int_t buffersize = 255;
3947 char name[buffersize];
3948 char title[buffersize];
3952 snprintf(name, buffersize,
"hSumPtLeadingPt_Cone_%d",icone);
3953 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3954 fhSumPtLeadingPt[icone] =
new TH2F(name, title,nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3960 snprintf(name, buffersize,
"hPtLeadingPt_Cone_%d",icone);
3961 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3962 fhPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3963 fhPtLeadingPt[icone] ->SetYTitle(
"#it{p}_{T}^{cone} (GeV/#it{c})");
3964 fhPtLeadingPt[icone] ->SetXTitle(
"#it{p}_{T}^{leading} (GeV/#it{c})");
3968 snprintf(name, buffersize,
"hPerpSumPtLeadingPt_Cone_%d",icone);
3969 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3976 snprintf(name, buffersize,
"hPerpPtLeadingPt_Cone_%d",icone);
3977 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3978 fhPerpPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3987 snprintf(name , buffersize,
"hSumPtLeadingPt_MC%s_Cone_%d",mcPartName[imc].
Data(),icone);
3988 snprintf(title, buffersize,
"Candidate %s #it{p}_{T} vs cone #Sigma #it{p}_{T} for #it{R}=%2.2f",mcPartType[imc].
Data(),
fConeSizes[icone]);
3998 snprintf(name, buffersize,
"hPtThres_Cone_%d_Pt%d",icone,ipt);
3999 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]);
4004 snprintf(name, buffersize,
"hPtFrac_Cone_%d_Pt%d",icone,ipt);
4005 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]);
4010 snprintf(name, buffersize,
"hSumPt_Cone_%d_Pt%d",icone,ipt);
4011 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]);
4012 fhSumPtIsolated[icone][ipt] =
new TH1F(name, title,nptbins,ptmin,ptmax);
4017 snprintf(name, buffersize,
"hPtSumDensity_Cone_%d_Pt%d",icone,ipt);
4018 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]);
4024 snprintf(name, buffersize,
"hPtFracPtSum_Cone_%d_Pt%d",icone,ipt);
4025 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]);
4032 snprintf(name, buffersize,
"hEtaPhiPtThres_Cone_%d_Pt%d",icone,ipt);
4033 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]);
4039 snprintf(name, buffersize,
"hEtaPhiPtFrac_Cone_%d_Pt%d",icone,ipt);
4040 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]);
4041 fhEtaPhiPtFracIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4046 snprintf(name, buffersize,
"hEtaPhiPtSum_Cone_%d_Pt%d",icone,ipt);
4047 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]);
4048 fhEtaPhiPtSumIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4053 snprintf(name, buffersize,
"hEtaPhiSumDensity_Cone_%d_Pt%d",icone,ipt);
4054 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]);
4060 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Cone_%d_Pt%d",icone,ipt);
4061 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]);
4070 snprintf(name, buffersize,
"hPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4071 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]);
4076 snprintf(name, buffersize,
"hPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4077 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]);
4082 snprintf(name, buffersize,
"hPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4083 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]);
4089 snprintf(name, buffersize,
"hPtSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4090 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]);
4096 snprintf(name, buffersize,
"hPtFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4097 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]);
4104 snprintf(name, buffersize,
"hEtaPhiPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4105 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]);
4111 snprintf(name, buffersize,
"hEtaPhiPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4112 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]);
4119 snprintf(name, buffersize,
"hEtaPhiPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4120 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]);
4126 snprintf(name, buffersize,
"hEtaPhiSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4127 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]);
4133 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4134 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]);
4146 snprintf(name , buffersize,
"hPtThreshMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4147 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #it{p}_{T}^{th}=%2.2f",
4155 snprintf(name , buffersize,
"hPtFracMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4156 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",
4163 snprintf(name , buffersize,
"hSumPtMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4164 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #Sigma #it{p}_{T}^{in cone}=%2.2f",
4178 for (
Int_t i = 0; i < 7 ; i++)
4181 Form(
"Number of isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4182 nptbins,ptmin,ptmax);
4188 Form(
"Number of isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4189 nptbins,ptmin,ptmax);
4195 Form(
"Number of not isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4196 nptbins,ptmin,ptmax);
4202 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4203 nptbins,ptmin,ptmax);
4209 fhTimeENoCut =
new TH2F (
"hTimeE_NoCut",
"time of cluster vs E of clusters, no cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4214 fhTimeESPD =
new TH2F (
"hTimeE_SPD",
"time of cluster vs E of clusters, SPD cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4219 fhTimeESPDMulti =
new TH2F (
"hTimeE_SPDMulti",
"time of cluster vs E of clusters, SPD multi cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4224 fhTimeNPileUpVertSPD =
new TH2F (
"hTime_NPileUpVertSPD",
"time of cluster vs N pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4229 fhTimeNPileUpVertTrack =
new TH2F (
"hTime_NPileUpVertTracks",
"time of cluster vs N pile-up Tracks vertex", ntimebins,timemin,timemax, 50,0,50 );
4234 fhTimeNPileUpVertContributors =
new TH2F (
"hTime_NPileUpVertContributors",
"time of cluster vs N constributors to pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4239 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);
4244 fhTimePileUpMainVertexZDiamond =
new TH2F (
"hTime_PileUpMainVertexZDiamond",
"time of cluster vs distance in Z pile-up SPD vertex - z diamond",ntimebins,timemin,timemax,100,0,50);
4250 return outputContainer ;
4310 AliInfo(
"Open default isolation cuts for multiple Isolation analysis");
4317 AliFatal(
"STOP!: You want to use CTS tracks in analysis but not read!! \n!!Check the configuration file!!");
4382 AliAODPWG4ParticleCorrelation* pLeading = 0;
4388 AliAODPWG4ParticleCorrelation* particle = (AliAODPWG4ParticleCorrelation*) (
GetInputAODBranch()->At(iaod));
4389 particle->SetLeadingParticle(kFALSE);
4395 if(check == 0)
continue;
4396 if(check == -1)
return kFALSE;
4408 if (particle->Pt() > ptTrig)
4410 ptTrig = particle->Pt() ;
4411 phiTrig = particle->Phi();
4413 pLeading = particle ;
4417 if(index < 0)
return kFALSE;
4423 if(phiTrig < 0 ) phiTrig += TMath::TwoPi();
4429 AliVTrack * track = (AliVTrack *) (
GetCTSTracks()->At(ipr)) ;
4437 Bool_t contained = kFALSE;
4439 for(
Int_t i = 0; i < 4; i++)
4441 if( trackID == pLeading->GetTrackLabel(i) ) contained = kTRUE;
4444 if ( contained ) continue ;
4447 fTrackVector.SetXYZ(track->Px(),track->Py(),track->Pz());
4450 if(phi < 0) phi+=TMath::TwoPi();
4454 Float_t deltaPhi = phiTrig-phi;
4460 if(deltaPhi <= -TMath::PiOver2()) deltaPhi+=TMath::TwoPi();
4461 if(deltaPhi > 3*TMath::PiOver2()) deltaPhi-=TMath::TwoPi();
4463 if(pt > ptTrig && deltaPhi < TMath::PiOver2())
return kFALSE;
4473 if (pLeading->GetDetectorTag() ==
kPHOS )
4478 if(!nePl)
return kTRUE;
4480 for(
Int_t ipr = 0;ipr < nePl->GetEntriesFast() ; ipr ++ )
4482 AliVCluster * cluster = (AliVCluster *) (nePl->At(ipr)) ;
4484 if(cluster->GetID() == pLeading->GetCaloLabel(0) || cluster->GetID() == pLeading->GetCaloLabel(1) )
continue ;
4490 if(phi < 0) phi+=TMath::TwoPi();
4497 Float_t deltaPhi = phiTrig-phi;
4498 if(deltaPhi <= -TMath::PiOver2()) deltaPhi+=TMath::TwoPi();
4499 if(deltaPhi > 3*TMath::PiOver2()) deltaPhi-=TMath::TwoPi();
4501 if(pt > ptTrig && deltaPhi < TMath::PiOver2())
return kFALSE ;
4507 pLeading->SetLeadingParticle(kTRUE);
4509 AliDebug(1,Form(
"Particle AOD with index %d is leading with pT %2.2f",idLeading, pLeading->Pt()));
4522 AliFatal(Form(
"No input particles in AOD with name branch < %s >, STOP",
GetInputAODName().
Data()));