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];
1594 Int_t overlab[nlabels];
1624 if ( noverlaps == 1 )
1651 else if (noverlaps > 1 )
1737 if(TMath::Abs(dZ) < 0.05 && TMath::Abs(dR) < 0.05)
1743 Float_t dEdx = track->GetTPCsignal();
1783 const Int_t buffersize = 255;
1784 char onePar[buffersize] ;
1786 snprintf(onePar, buffersize,
"--- AliAnaParticleIsolation ---:") ;
1790 snprintf(onePar, buffersize,
"Isolation Cand Detector: %s;",
fIsoDetectorString.Data()) ;
1792 snprintf(onePar, buffersize,
"fReMakeIC =%d (Flag for reisolation during histogram filling);",
fReMakeIC) ;
1794 snprintf(onePar, buffersize,
"fMakeSeveralIC=%d (Flag for isolation with several cuts at the same time );",
fMakeSeveralIC) ;
1796 snprintf(onePar, buffersize,
"fFillTMHisto=%d (Flag for track matching histograms);",
fFillTMHisto) ;
1798 snprintf(onePar, buffersize,
"fFillSSHisto=%d (Flag for shower shape histograms);",
fFillSSHisto) ;
1803 snprintf(onePar, buffersize,
"fNCones =%d (Number of cone sizes);",
fNCones) ;
1805 snprintf(onePar, buffersize,
"fNPtThresFrac=%d (Flag for isolation with several cuts at the same time);",
fNPtThresFrac) ;
1810 snprintf(onePar, buffersize,
"fConeSizes[%d]=%1.2f (isolation cone size);",icone,
fConeSizes[icone]) ;
1815 snprintf(onePar, buffersize,
"fPtThresholds[%d]=%1.2f (isolation pt threshold);",ipt,
fPtThresholds[ipt]) ;
1820 snprintf(onePar, buffersize,
"fPtFractions[%d]=%1.2f (isolation pt fraction threshold);",ipt,
fPtFractions[ipt]) ;
1825 snprintf(onePar, buffersize,
"fSumPtThresholds[%d]=%1.2f (isolation sum pt threshold);",ipt,
fSumPtThresholds[ipt]) ;
1836 return new TObjString(parList) ;
1846 outputContainer->SetName(
"IsolatedParticleHistos") ;
1895 sThreshold = Form(
", %2.2f < #Sigma #it{p}_{T}^{in cone} < %2.2f GeV/#it{c}",
1898 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone} = %2.2f GeV/#it{c}",
1903 sThreshold = Form(
", %2.2f < #it{p}_{T}^{th} < %2.2f GeV/#it{c}",
1906 sThreshold = Form(
", #it{p}_{T}^{th} = %2.2f GeV/#it{c}",
1910 sThreshold = Form(
", #Sigma #it{p}_{T}^{in cone}/#it{p}_{T}^{trig} = %2.2f" ,
1913 TString sParticle =
", x^{0,#pm}";
1917 TString parTitle = Form(
"#it{R} = %2.2f%s%s",
GetIsolationCut()->GetConeSize(), sThreshold.Data(),sParticle.Data());
1919 TString pileUpName[] = {
"SPD",
"EMCAL",
"SPDOrEMCAL",
"SPDAndEMCAL",
"SPDAndNotEMCAL",
"EMCALAndNotSPD",
"NotSPDAndNotEMCAL"} ;
1922 TString mcPartType[] = {
"#gamma" ,
"#gamma_{prompt}" ,
"#gamma_{fragmentation}",
1923 "#pi^{0} (merged #gamma)",
"#gamma_{#pi decay}" ,
"#gamma_{#pi decay} lost companion",
1924 "#eta (merged #gamma)" ,
"#gamma_{#eta decay}",
"#gamma_{#eta decay} lost companion",
1925 "#gamma_{other decay}" ,
"e^{#pm}" ,
"hadrons?"} ;
1927 TString mcPartName[] = {
"Photon",
"PhotonPrompt",
"PhotonFrag",
1928 "Pi0" ,
"Pi0Decay" ,
"Pi0DecayLostPair",
1929 "Eta" ,
"EtaDecay" ,
"EtaDecayLostPair",
1930 "OtherDecay",
"Electron",
"Hadron"} ;
1933 TString pptype[] = {
"#gamma" ,
"#gamma_{#pi decay}" ,
"#gamma_{#eta decay}",
"#gamma_{other decay}",
1934 "#gamma_{prompt}",
"#gamma_{fragmentation}",
"#gamma_{ISR}" ,
1935 "#pi^{0}" ,
"#eta"} ;
1937 TString ppname[] = {
"Photon" ,
"PhotonPi0Decay",
"PhotonEtaDecay",
"PhotonOtherDecay",
1938 "PhotonPrompt",
"PhotonFrag" ,
"PhotonISR" ,
1944 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
1945 nptbins,ptmin,ptmax);
1946 fhENoIso->SetYTitle(
"#it{counts}");
1947 fhENoIso->SetXTitle(
"E (GeV/#it{c})");
1951 Form(
"Number of not isolated leading particles vs #it{p}_{T}, %s",parTitle.Data()),
1952 nptbins,ptmin,ptmax);
1954 fhPtNoIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
1958 Form(
"Number of not isolated leading particles #eta vs #phi, %s",parTitle.Data()),
1959 netabins,etamin,etamax,nphibins,phimin,phimax);
1970 fhPtNoIsoMC[imc] =
new TH1F(Form(
"hPtNoIsoMC%s",mcPartName[imc].
Data()),
1971 Form(
"#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1972 nptbins,ptmin,ptmax);
1974 fhPtNoIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1977 fhPtIsoMC[imc] =
new TH1F(Form(
"hPtMC%s",mcPartName[imc].
Data()),
1978 Form(
"#it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1979 nptbins,ptmin,ptmax);
1980 fhPtIsoMC[imc]->SetYTitle(
"#it{counts}");
1981 fhPtIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1985 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1986 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
1988 fhPhiIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
1992 Form(
"#phi vs #it{p}_{T} of isolated %s, %s",mcPartType[imc].
Data(),parTitle.Data()),
1993 nptbins,ptmin,ptmax,netabins,etamin,etamax);
1995 fhEtaIsoMC[imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2003 TString isoName [] = {
"NoIso",
"Iso"};
2004 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2008 for(
Int_t iso = 0; iso < 2; iso++)
2012 new TH1F(Form(
"hPtDecay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2013 Form(
"Number of %s leading pi0 decay particles vs #it{p}_{T}, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2014 nptbins,ptmin,ptmax);
2015 fhPtDecay[iso][ibit]->SetYTitle(
"#it{counts}");
2016 fhPtDecay[iso][ibit]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2017 outputContainer->Add(
fhPtDecay[iso][ibit]) ;
2021 Form(
"Number of %s leading Pi0 decay particles #eta vs #phi, bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],parTitle.Data()),
2022 netabins,etamin,etamax,nphibins,phimin,phimax);
2030 (Form(
"hPtLambda0Decay%s_bit%d",isoName[iso].
Data(),
fDecayBits[ibit]),
2031 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, decay bit %d, %s",isoTitle[iso].
Data(),
fDecayBits[ibit], parTitle.Data()),
2032 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
2043 new TH1F(Form(
"hPtDecay%s_bit%d_MC%s",isoName[iso].
Data(),
fDecayBits[ibit],mcPartName[imc].
Data()),
2044 Form(
"#it{p}_{T} of %s, decay bit %d, %s, %s",isoTitle[iso].
Data(),
fDecayBits[ibit],mcPartType[imc].
Data(),parTitle.Data()),
2045 nptbins,ptmin,ptmax);
2046 fhPtDecayMC[iso][ibit][imc]->SetYTitle(
"#it{counts}");
2047 fhPtDecayMC[iso][ibit][imc]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
2048 outputContainer->Add(
fhPtDecayMC[iso][ibit][imc]) ;
2057 TString isoName [] = {
"NoIso",
"Iso"};
2058 TString isoTitle[] = {
"Not isolated" ,
"isolated"};
2061 Form(
"Number of isolated particles vs E, %s",parTitle.Data()),
2062 nptbins,ptmin,ptmax);
2063 fhEIso->SetYTitle(
"d#it{N} / d#it{E}");
2064 fhEIso->SetXTitle(
"#it{E} (GeV/#it{c})");
2065 outputContainer->Add(
fhEIso) ;
2068 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2069 nptbins,ptmin,ptmax);
2070 fhPtIso->SetYTitle(
"d#it{N} / #it{p}_{T}");
2071 fhPtIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2072 outputContainer->Add(
fhPtIso) ;
2075 Form(
"Number of isolated particles vs #phi, %s",parTitle.Data()),
2076 nptbins,ptmin,ptmax,nphibins,phimin,phimax);
2078 fhPhiIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2082 Form(
"Number of isolated particles vs #eta, %s",parTitle.Data()),
2083 nptbins,ptmin,ptmax,netabins,etamin,etamax);
2085 fhEtaIso->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2089 Form(
"Number of isolated particles #eta vs #phi, %s",parTitle.Data()),
2090 netabins,etamin,etamax,nphibins,phimin,phimax);
2098 Form(
"centrality vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2099 nptbins,ptmin,ptmax, 100,0,100);
2105 Form(
"event plane angle vs #it{p}_{T} for isolated particles, %s",parTitle.Data()),
2106 nptbins,ptmin,ptmax, 100,0,TMath::Pi());
2115 Form(
"Number of isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2116 nptbins,ptmin,ptmax,10,0,10);
2121 Form(
"Number of not isolated particles vs #it{p}_{T}, %s",parTitle.Data()),
2122 nptbins,ptmin,ptmax,10,0,10);
2129 Form(
"Track or Cluster leading #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2130 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2131 fhConePtLead->SetYTitle(
"#it{p}_{T, leading} (GeV/#it{c})");
2132 fhConePtLead->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2136 Form(
"Track and Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2137 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2139 fhConeSumPt->SetXTitle(
"#it{p}_{T, trigger} (GeV/#it{c})");
2143 (Form(
"hPtLambda0Eiso"),
2144 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.);
2151 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2152 netabins,etamin,etamax,nphibins,phimin,phimax);
2159 Form(
"#it{p}_{T} of clusters and tracks in isolation cone for #it{R} = %2.2f",r),
2160 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2161 fhPtInCone->SetYTitle(
"#it{p}_{T in cone} (GeV/#it{c})");
2162 fhPtInCone->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
2196 (Form(
"hPtLeadCone_Bin%d",ibin),
2197 Form(
"cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2204 (Form(
"hSumPtCone_Bin%d",ibin),
2205 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2218 (Form(
"hPtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2219 Form(
"Decay bit %d, cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2226 (Form(
"hSumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2227 Form(
"Decay bit %d, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2241 (Form(
"hPtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2242 Form(
"in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2249 (Form(
"hSumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2250 Form(
"in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2261 (Form(
"hPtLeadConeLambda0_Bin%d",ibin),
2262 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, %s",
2269 (Form(
"hSumPtConeLambda0_Bin%d",ibin),
2270 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, %s",
2282 (Form(
"hPtLeadConeLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2283 Form(
"#lambda_{0}, in cone %2.2f<#it{p}_{T}^{leading}<%2.2f GeV/#it{c}, MC %s, %s",
2290 (Form(
"hSumPtConeLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2291 Form(
"#lambda_{0}, in cone %2.2f <#Sigma #it{p}_{T}< %2.2f GeV/#it{c}, MC %s, %s",
2331 (Form(
"hPtTrigBin_PtLeadCone_Bin%d",ibin),
2332 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2339 (Form(
"hPtTrigBin_SumPtCone_Bin%d",ibin),
2340 Form(
"#Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2353 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2354 Form(
"Decay bit %d, #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2361 (Form(
"hPtTrigBin_SumPtCone_Bin%d_DecayBit%d",ibin,
fDecayBits[idecay]),
2362 Form(
"Decay bit %d, #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2376 (Form(
"hPtTrigBin_PtLeadCone_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2377 Form(
"#it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2384 (Form(
"hPtTrigBin_SumPtCone_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2385 Form(
"#Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2396 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d",ibin),
2397 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, %s",
2404 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d",ibin),
2405 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone} %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, %s",
2417 (Form(
"hPtTrigBin_PtLeadConeVSLambda0_Bin%d_MC%s",ibin, mcPartName[imc].
Data()),
2418 Form(
"#lambda_{0} vs #it{p}_{T}^{lead. in cone}, %2.2f<#it{p}_{T}^{cand}<%2.2f GeV/#it{c}, MC %s, %s",
2425 (Form(
"hPtTrigBin_SumPtConeVSLambda0_Bin%d_MC%s",ibin,mcPartName[imc].
Data()),
2426 Form(
"#lambda_{0} vs #Sigma #it{p}_{T}^{in cone}, %2.2f <#it{p}_{T}^{cand}< %2.2f GeV/#it{c}, MC %s, %s",
2440 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2441 100,0,100,nptinconebins,ptinconemin,ptinconemax);
2451 Form(
"Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2452 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2458 Form(
"Cluster leading in isolation cone for #it{R} = %2.2f",r),
2459 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2468 Form(
"Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2469 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2478 "#Sigma cluster #it{p}_{T} in UE Eta Band",
2479 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2485 "#Sigma cluster #it{p}_{T} UE Phi Band",
2486 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2492 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} in UE Eta Band",
2493 netabins,etamin,etamax,nphibins,phimin,phimax);
2500 "Trigger #eta vs #phi, #Sigma cluster #it{p}_{T} UE Phi Band",
2501 netabins,etamin,etamax,nphibins,phimin,phimax);
2510 "#Sigma cell #it{p}_{T} in UE Eta Band",
2511 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2517 "#Sigma cell #it{p}_{T} UE Phi Band",
2518 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2524 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} in UE Eta Band",
2525 netabins,etamin,etamax,nphibins,phimin,phimax);
2532 "Trigger #eta vs #phi, #Sigma cell #it{p}_{T} UE Phi Band",
2533 netabins,etamin,etamax,nphibins,phimin,phimax);
2541 Form(
"#eta vs #phi of clusters in #eta band isolation cone for #it{R} = %2.2f",r),
2542 netabins,-1,1,nphibins,0,TMath::TwoPi());
2548 Form(
"#eta vs #phi of clusters in #phi band isolation cone for #it{R} = %2.2f",r),
2549 netabins,-1,1,nphibins,0,TMath::TwoPi());
2555 Form(
"#eta vs #phi of clusters in cone for #it{R} = %2.2f",r),
2556 netabins,-1,1,nphibins,0,TMath::TwoPi());
2562 Form(
"#eta vs #phi of all clusters"),
2563 netabins,-1,1,nphibins,0,TMath::TwoPi());
2571 Form(
"#it{p}_{T} of clusters in isolation cone for #it{R} = %2.2f",r),
2572 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2580 Form(
"#it{p}_{T} of cells in isolation cone for #it{R} = %2.2f",r),
2581 nptbins,ptmin,ptmax,1000,0,50);
2587 Form(
"#col vs #row of cells in #eta band isolation cone for #it{R} = %2.2f",r),
2594 Form(
"#col vs #row of cells in #phi band isolation cone for #it{R} = %2.2f",r),
2604 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2605 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2611 Form(
"Clusters #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2612 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2618 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),
2619 netabins,etamin,etamax,nphibins,phimin,phimax);
2626 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),
2627 netabins,etamin,etamax,nphibins,phimin,phimax);
2636 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2637 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2643 Form(
"Cells #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2644 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2650 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),
2651 netabins,etamin,etamax,nphibins,phimin,phimax);
2658 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),
2659 netabins,etamin,etamax,nphibins,phimin,phimax);
2667 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance",r),
2668 nptbins,ptmin,ptmax,100,0,1);
2674 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #eta acceptance, in trigger #eta-#phi ",r),
2675 netabins,etamin,etamax,nphibins,phimin,phimax);
2682 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance",r),
2683 nptbins,ptmin,ptmax,100,0,1);
2689 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of clusters #phi acceptance, in trigger #eta-#phi ",r),
2690 netabins,etamin,etamax,nphibins,phimin,phimax);
2697 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),
2698 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2704 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),
2705 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2711 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),
2712 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2718 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),
2719 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2725 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for cluster (before normalization), R=%2.2f",r),
2726 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
2732 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for cluster (before normalization), R=%2.2f",r),
2733 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
2741 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance",r),
2742 nptbins,ptmin,ptmax,100,0,1);
2748 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #eta acceptance, in trigger #eta-#phi ",r),
2749 netabins,etamin,etamax,nphibins,phimin,phimax);
2756 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance",r),
2757 nptbins,ptmin,ptmax,100,0,1);
2763 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of cells #phi acceptance, in trigger #eta-#phi ",r),
2764 netabins,etamin,etamax,nphibins,phimin,phimax);
2772 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),
2773 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2779 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),
2780 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2786 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),
2787 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2793 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),
2794 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2806 Form(
"Track #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2807 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2813 Form(
"Track leading in isolation cone for #it{R} = %2.2f",r),
2814 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
2820 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f",r),
2821 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2830 "#Sigma track #it{p}_{T} in UE Eta Band",
2831 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2837 "#Sigma track #it{p}_{T} in UE Phi Band",
2838 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax*8);
2845 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Eta Band",
2846 netabins,etamin,etamax,nphibins,phimin,phimax);
2853 "Trigger #eta vs #phi, #Sigma track #it{p}_{T} in UE Phi Band",
2854 netabins,etamin,etamax,nphibins,phimin,phimax);
2861 Form(
"#eta vs #phi of tracks in #eta band isolation cone for #it{R} = %2.2f",r),
2862 netabins,-1,1,nphibins,0,TMath::TwoPi());
2868 Form(
"#eta vs #phi of tracks in #phi band isolation cone for #it{R} = %2.2f",r),
2869 netabins,-1,1,nphibins,0,TMath::TwoPi());
2875 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
2876 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2882 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
2883 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
2889 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),
2890 netabins,etamin,etamax,nphibins,phimin,phimax);
2897 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),
2898 netabins,etamin,etamax,nphibins,phimin,phimax);
2905 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance",r),
2906 nptbins,ptmin,ptmax,100,0,1);
2912 Form(
"Fraction of the isolation cone #it{R} = %2.2f, out of tracks #eta acceptance, in trigger #eta-#phi ",r),
2913 netabins,etamin,etamax,nphibins,phimin,phimax);
2920 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),
2921 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2927 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),
2928 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2934 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),
2935 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2941 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),
2942 nptsumbins,ptsummin,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
2950 Form(
"#Sigma #it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
2951 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
2957 Form(
"#it{p}_{T} in isolation cone at #pm 45 degree phi from trigger particle, #it{R} = %2.2f",r),
2958 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
2964 Form(
"#eta vs #phi of all Tracks"),
2965 netabins,-1,1,nphibins,0,TMath::TwoPi());
2971 Form(
"#eta vs #phi of Tracks in cone for #it{R} = %2.2f",r),
2972 netabins,-1,1,nphibins,0,TMath::TwoPi());
2978 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in eta band for tracks (before normalization), R=%2.2f",r),
2979 nptsumbins,ptsummin,ptsummax,2*nptsumbins,ptsummin,2*ptsummax);
2985 Form(
"#Sigma #it{p}_{T} in cone versus #Sigma #it{p}_{T} in phi band for tracks (before normalization), R=%2.2f",r),
2986 nptsumbins,ptsummin,ptsummax,8*nptsumbins,ptsummin,8*ptsummax);
2996 Form(
"Track vs Cluster #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
2997 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3003 Form(
"#Sigma #it{p}_{T}^{cluster}/#Sigma #it{p}_{T}^{track} in isolation cone for #it{R} = %2.2f",r),
3004 nptbins,ptmin,ptmax,200,0,5);
3011 Form(
"Track vs Cluster lead #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3012 nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);
3018 Form(
" #it{p}^{leading cluster}_{T}/#it{p}^{leading track}_{T} in isolation cone for #it{R} = %2.2f",r),
3019 nptbins,ptmin,ptmax,200,0,5);
3028 Form(
"Track vs cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3029 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3035 Form(
"Track and Cell #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3036 nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3042 Form(
"Trigger #eta vs #phi, #Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",r),
3043 netabins,etamin,etamax,nphibins,phimin,phimax);
3053 Form(
"#Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3054 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3060 Form(
"#Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3061 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3067 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),
3068 netabins,etamin,etamax,nphibins,phimin,phimax);
3075 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),
3076 netabins,etamin,etamax,nphibins,phimin,phimax);
3083 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3084 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3090 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3091 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3097 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3098 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3104 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3105 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3111 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3112 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3118 Form(
"Track vs Cluster #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3119 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3125 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3126 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3132 Form(
"Track vs Cluster #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3133 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3142 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub eta band in isolation cone for #it{R} = %2.2f",r),
3143 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3149 Form(
"Track vs Cell #Sigma #it{p}_{T} UE sub phi band in isolation cone for #it{R} = %2.2f",r),
3150 2*nptsumbins,-ptsummax,ptsummax,2*nptsumbins,-ptsummax,ptsummax);
3156 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3157 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3163 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3164 nptsumbins,ptsummin,ptsummax*4,nptsumbins,ptsummin,ptsummax*8);
3170 Form(
"Track vs Cell #Sigma #it{p}_{T} in Eta band in isolation cone for #it{R} = %2.2f",r),
3171 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3177 Form(
"Track vs Cell #Sigma #it{p}_{T} in Phi band in isolation cone for #it{R} = %2.2f",r),
3178 nptsumbins,ptsummin,ptsummax,nptsumbins,ptsummin,ptsummax);
3184 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from eta band in the isolation cone for #it{R} = %2.2f",r),
3185 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3191 Form(
"Tracks #Sigma #it{p}_{T} after bkg subtraction from phi band in the isolation cone for #it{R} = %2.2f",r),
3192 nptbins,ptmin,ptmax,2*nptsumbins,-ptsummax,ptsummax);
3198 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),
3199 netabins,etamin,etamax,nphibins,phimin,phimax);
3206 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),
3207 netabins,etamin,etamax,nphibins,phimin,phimax);
3216 TString region[] = {
"ITS",
"TPC",
"TRD",
"TOF",
"Top EMCal",
"In EMCal"};
3218 for(
Int_t iso = 0; iso < 2; iso++)
3223 (Form(
"hTrackMatchedDEta%s",isoName[iso].
Data()),
3224 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3225 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3230 (Form(
"hTrackMatchedDPhi%s",isoName[iso].
Data()),
3231 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3232 nptbins,ptmin,ptmax,nresphibins,resphimin,resphimax);
3237 (Form(
"hTrackMatchedDEtaDPhi%s",isoName[iso].
Data()),
3238 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3239 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3252 Form(
"%s - d#eta of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3253 nptbins,ptmin,ptmax,nresetabins,resetamin,resetamax);
3258 Form(
"%s - d#phi of cluster-track vs cluster energy, %s",isoTitle[iso].
Data(),parTitle.Data()),
3259 nptbins,ptmin,ptmax,nresetabins,resphimin,resphimax);
3264 (Form(
"hTrackMatchedDEtaDPhi%s_MC%s",isoName[iso].
Data(),mcPartName[imc].
Data()),
3265 Form(
"%s - d#eta vs d#phi of cluster-track, %s",isoTitle[iso].
Data(),parTitle.Data()),
3266 nresetabins,resetamin,resetamax,nresphibins,resphimin,resphimax);
3278 (Form(
"hdEdx%s",isoName[iso].
Data()),
3279 Form(
"%s - Matched track <d#it{E}/d#it{x}> vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3280 nptbins,ptmin,ptmax,ndedxbins, dedxmin, dedxmax);
3281 fhdEdx[iso]->SetXTitle(
"#it{E} (GeV)");
3282 fhdEdx[iso]->SetYTitle(
"<d#it{E}/d#it{x}>");
3283 outputContainer->Add(
fhdEdx[iso]);
3286 (Form(
"hEOverP%s",isoName[iso].
Data()),
3287 Form(
"%s - Matched track #it{E}/#it{p} vs cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3288 nptbins,ptmin,ptmax,nPoverEbins,pOverEmin,pOverEmax);
3289 fhEOverP[iso]->SetXTitle(
"#it{E} (GeV)");
3290 fhEOverP[iso]->SetYTitle(
"#it{E}/#it{p}");
3291 outputContainer->Add(
fhEOverP[iso]);
3296 (Form(
"hTrackMatchedMCParticle%s",isoName[iso].
Data()),
3297 Form(
"%s - Origin of particle vs cluster #it{E}, %s",isoTitle[iso].
Data(),parTitle.Data()),
3298 nptbins,ptmin,ptmax,8,0,8);
3318 (Form(
"hELambda0%s",isoName[iso].
Data()),
3319 Form(
"%s cluster : #it{E} vs #lambda_{0}, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3320 fhELambda0[iso]->SetYTitle(
"#lambda_{0}^{2}");
3332 (Form(
"hPtLambda0%s",isoName[iso].
Data()),
3333 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}, %s",isoTitle[iso].
Data(), parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3335 fhPtLambda0[iso]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3343 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3344 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3346 fhPtLambda0MC[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3350 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3351 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3359 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3360 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3366 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3367 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3373 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3374 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3380 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3381 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3388 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3389 nptbins,ptmin,ptmax,10,0,10);
3391 fhPtNOverlap[imc][iso]->SetXTitle(
"#it{p}_{T}(GeV/#it{c})");
3395 Form(
"%s cluster : #it{p}_{T} vs #lambda_{0}: %s %s, from conversion, 1 overlap",isoTitle[iso].
Data(),mcPartType[imc].
Data(),parTitle.Data()),
3396 nptbins,ptmin,ptmax,10,0,10);
3408 (Form(
"hPtLambda0TRD%s",isoName[iso].
Data()),
3409 Form(
"%s cluster: #it{p}_{T} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3415 (Form(
"hELambda0TRD%s",isoName[iso].
Data()),
3416 Form(
"%s cluster: #it{E} vs #lambda_{0}, SM behind TRD, %s",isoTitle[iso].
Data(),parTitle.Data()),nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3432 (Form(
"hNLocMax%s",isoName[iso].
Data()),
3433 Form(
"%s - Number of local maxima in cluster, %s",isoTitle[iso].
Data(),parTitle.Data()),
3434 nptbins,ptmin,ptmax,10,0,10);
3436 fhNLocMax[iso]->SetXTitle(
"#it{E} (GeV)");
3440 (Form(
"hELambda0LocMax1%s",isoName[iso].
Data()),
3441 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);
3447 (Form(
"hELambda1LocMax1%s",isoName[iso].
Data()),
3448 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);
3454 (Form(
"hELambda0LocMax2%s",isoName[iso].
Data()),
3455 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);
3461 (Form(
"hELambda1LocMax2%s",isoName[iso].
Data()),
3462 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);
3468 ( Form(
"hELambda0LocMaxN%s",isoName[iso].
Data()),
3469 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);
3475 (Form(
"hELambda1LocMaxN%s",isoName[iso].
Data()),
3476 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);
3486 for(
Int_t ieta = 0; ieta < 4; ieta++)
3488 for(
Int_t iphi = 0; iphi < 3; iphi++)
3514 new TH2F(Form(
"hLam0_%s_eta%d_phi%d_sm%d",isoName[iso].
Data(),ieta,iphi,ism),
3515 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, sm %d, region eta %d, phi %d",
3516 isoTitle[iso].
Data(),ism,ieta,iphi),
3517 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3525 Float_t ptLimit[] = {2,3,4,5,6,8,10,12};
3526 for(
Int_t ipt = 0; ipt < 7; ipt++)
3529 (Form(
"hEtaPhiLam0BinPtBin%d%s",ipt,isoName[iso].
Data()),
3530 Form(
"%s, #eta vs #phi in #it{p}_{T}=[%2.1f,%2.1f] GeV/#it{c} and #lambda^{2}_{0}=[0.3,0.4]",
3531 isoTitle[iso].
Data(),ptLimit[ipt],ptLimit[ipt+1]),
3532 netabins,etamin,etamax,nphibins,phimin,phimax);
3542 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s",
3543 isoTitle[iso].
Data(),parTitle.Data()),
3544 nptbins,ptmin,ptmax,500,0,500);
3552 Form(
"%s, cluster from converted photon, #it{p}_{T} vs vertex distance, %s, SM covered by TRD",
3553 isoTitle[iso].
Data(),parTitle.Data()),
3554 nptbins,ptmin,ptmax,500,0,500);
3562 for(
Int_t iR = 0; iR < 6; iR++)
3565 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s",
3566 isoTitle[iso].
Data(),region[iR].
Data()),
3567 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3575 Form(
"%s, cluster from converted photon, #it{p}_{T} vs #lambda_{0}^{2}, conversion in %s, SM covered by TRD",
3576 isoTitle[iso].
Data(),region[iR].
Data()),
3577 nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
3627 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0",r),
3628 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3634 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC!=0, pile-up from SPD",r),
3635 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3641 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3642 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3648 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0",r),
3649 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3656 Form(
"#it{p}_{T} of tracks in isolation cone for #it{R} = %2.2f, TOF from BC==0, pile-up from SPD",r),
3657 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3663 for (
Int_t i = 0; i < 7 ; i++)
3666 Form(
"#it{p}_{T} in isolation cone for #it{R} = %2.2f, from pile-up (%s)",r,pileUpName[i].
Data()),
3667 nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3680 fhEPrimMC[i] =
new TH1F(Form(
"hEPrim_MC%s",ppname[i].
Data()),
3681 Form(
"primary photon %s : #it{E}, %s",pptype[i].
Data(),parTitle.Data()),
3682 nptbins,ptmin,ptmax);
3683 fhEPrimMC[i]->SetXTitle(
"#it{E} (GeV)");
3687 Form(
"primary photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3688 nptbins,ptmin,ptmax);
3689 fhPtPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3693 Form(
"primary isolated photon %s : #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3694 nptbins,ptmin,ptmax);
3699 Form(
"primary photon %s : #eta vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3700 nptbins,ptmin,ptmax,200,-2,2);
3702 fhEtaPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3706 Form(
"primary photon %s : #phi vs #it{p}_{T}, %s",pptype[i].
Data(),parTitle.Data()),
3707 nptbins,ptmin,ptmax,200,0.,TMath::TwoPi());
3709 fhPhiPrimMC[i]->SetXTitle(
"#it{p}_{T} (GeV/#it{c})");
3716 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3717 nptbins,ptmin,ptmax);
3722 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
3724 nptbins,ptmin,ptmax);
3729 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
3731 nptbins,ptmin,ptmax);
3736 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
3738 nptbins,ptmin,ptmax);
3743 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
3745 nptbins,ptmin,ptmax);
3750 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
3752 nptbins,ptmin,ptmax);
3758 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
3760 nptbins,ptmin,ptmax);
3765 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
3767 nptbins,ptmin,ptmax);
3772 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3773 nptbins,ptmin,ptmax);
3778 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
3780 nptbins,ptmin,ptmax);
3785 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3786 nptbins,ptmin,ptmax);
3791 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimPi0Decay].
Data(),parTitle.Data()),
3792 nptbins,ptmin,ptmax);
3797 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
3799 nptbins,ptmin,ptmax);
3804 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
3806 nptbins,ptmin,ptmax);
3811 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3813 nptbins,ptmin,ptmax);
3818 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3820 nptbins,ptmin,ptmax);
3832 Form(
"primary photon %s : #it{p}_{T}, pair in cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3833 nptbins,ptmin,ptmax);
3838 Form(
"isolated primary photon %s, pair in cone : #it{p}_{T}, %s",
3840 nptbins,ptmin,ptmax);
3845 Form(
"primary photon %s, no overlap, pair in cone : #it{p}_{T}, %s",
3847 nptbins,ptmin,ptmax);
3852 Form(
"isolated primary photon %s, pair in cone,no overlap : #it{p}_{T}, %s",
3854 nptbins,ptmin,ptmax);
3859 Form(
"primary photon %s, no overlap, pair in cone, E > calo min: #it{p}_{T}, %s",
3861 nptbins,ptmin,ptmax);
3866 Form(
"isolated primary photon %s, pair in cone,no overlap, E > calo min: #it{p}_{T}, %s",
3868 nptbins,ptmin,ptmax);
3874 Form(
"primary photon %s, no overlap: #it{p}_{T}, %s",
3876 nptbins,ptmin,ptmax);
3881 Form(
"isolated primary photon %s, no overlap: #it{p}_{T}, %s",
3883 nptbins,ptmin,ptmax);
3888 Form(
"primary photon %s : #it{p}_{T}, pair out of cone, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3889 nptbins,ptmin,ptmax);
3894 Form(
"isolated primary photon %s, pair out of cone : #it{p}_{T}, %s",
3896 nptbins,ptmin,ptmax);
3901 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3902 nptbins,ptmin,ptmax);
3907 Form(
"primary photon %s : #it{p}_{T}, pair out of acceptance, no overlap, %s",pptype[
kmcPrimEtaDecay].
Data(),parTitle.Data()),
3908 nptbins,ptmin,ptmax);
3913 Form(
"isolated primary photon %s, pair out of acceptance : #it{p}_{T}, %s",
3915 nptbins,ptmin,ptmax);
3920 Form(
"isolated primary photon %s, pair out of acceptance, no overlap : #it{p}_{T}, %s",
3922 nptbins,ptmin,ptmax);
3927 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3929 nptbins,ptmin,ptmax);
3934 Form(
"primary %s, overlap: #it{p}_{T}, %s",
3936 nptbins,ptmin,ptmax);
3948 const Int_t buffersize = 255;
3949 char name[buffersize];
3950 char title[buffersize];
3954 snprintf(name, buffersize,
"hSumPtLeadingPt_Cone_%d",icone);
3955 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3956 fhSumPtLeadingPt[icone] =
new TH2F(name, title,nptbins,ptmin,ptmax,nptsumbins,ptsummin,ptsummax);
3962 snprintf(name, buffersize,
"hPtLeadingPt_Cone_%d",icone);
3963 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3964 fhPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3965 fhPtLeadingPt[icone] ->SetYTitle(
"#it{p}_{T}^{cone} (GeV/#it{c})");
3966 fhPtLeadingPt[icone] ->SetXTitle(
"#it{p}_{T}^{leading} (GeV/#it{c})");
3970 snprintf(name, buffersize,
"hPerpSumPtLeadingPt_Cone_%d",icone);
3971 snprintf(title, buffersize,
"#Sigma #it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3978 snprintf(name, buffersize,
"hPerpPtLeadingPt_Cone_%d",icone);
3979 snprintf(title, buffersize,
"#it{p}_{T} in isolation cone for #it{R} = %2.2f",
fConeSizes[icone]);
3980 fhPerpPtLeadingPt[icone] =
new TH2F(name, title, nptbins,ptmin,ptmax,nptinconebins,ptinconemin,ptinconemax);
3989 snprintf(name , buffersize,
"hSumPtLeadingPt_MC%s_Cone_%d",mcPartName[imc].
Data(),icone);
3990 snprintf(title, buffersize,
"Candidate %s #it{p}_{T} vs cone #Sigma #it{p}_{T} for #it{R}=%2.2f",mcPartType[imc].
Data(),
fConeSizes[icone]);
4000 snprintf(name, buffersize,
"hPtThres_Cone_%d_Pt%d",icone,ipt);
4001 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]);
4006 snprintf(name, buffersize,
"hPtFrac_Cone_%d_Pt%d",icone,ipt);
4007 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]);
4012 snprintf(name, buffersize,
"hSumPt_Cone_%d_Pt%d",icone,ipt);
4013 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]);
4014 fhSumPtIsolated[icone][ipt] =
new TH1F(name, title,nptbins,ptmin,ptmax);
4019 snprintf(name, buffersize,
"hPtSumDensity_Cone_%d_Pt%d",icone,ipt);
4020 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]);
4026 snprintf(name, buffersize,
"hPtFracPtSum_Cone_%d_Pt%d",icone,ipt);
4027 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]);
4034 snprintf(name, buffersize,
"hEtaPhiPtThres_Cone_%d_Pt%d",icone,ipt);
4035 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]);
4041 snprintf(name, buffersize,
"hEtaPhiPtFrac_Cone_%d_Pt%d",icone,ipt);
4042 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]);
4043 fhEtaPhiPtFracIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4048 snprintf(name, buffersize,
"hEtaPhiPtSum_Cone_%d_Pt%d",icone,ipt);
4049 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]);
4050 fhEtaPhiPtSumIso[icone][ipt] =
new TH2F(name, title,netabins,etamin,etamax,nphibins,phimin,phimax);
4055 snprintf(name, buffersize,
"hEtaPhiSumDensity_Cone_%d_Pt%d",icone,ipt);
4056 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]);
4062 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Cone_%d_Pt%d",icone,ipt);
4063 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]);
4072 snprintf(name, buffersize,
"hPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4073 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]);
4078 snprintf(name, buffersize,
"hPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4079 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]);
4084 snprintf(name, buffersize,
"hPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4085 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]);
4091 snprintf(name, buffersize,
"hPtSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4092 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]);
4098 snprintf(name, buffersize,
"hPtFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4099 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]);
4106 snprintf(name, buffersize,
"hEtaPhiPtThres_Decay_Cone_%d_Pt%d",icone,ipt);
4107 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]);
4113 snprintf(name, buffersize,
"hEtaPhiPtFrac_Decay_Cone_%d_Pt%d",icone,ipt);
4114 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]);
4121 snprintf(name, buffersize,
"hEtaPhiPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4122 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]);
4128 snprintf(name, buffersize,
"hEtaPhiSumDensity_Decay_Cone_%d_Pt%d",icone,ipt);
4129 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]);
4135 snprintf(name, buffersize,
"hEtaPhiFracPtSum_Decay_Cone_%d_Pt%d",icone,ipt);
4136 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]);
4148 snprintf(name , buffersize,
"hPtThreshMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4149 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #it{p}_{T}^{th}=%2.2f",
4157 snprintf(name , buffersize,
"hPtFracMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4158 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",
4165 snprintf(name , buffersize,
"hSumPtMC%s_Cone_%d_Pt%d",mcPartName[imc].
Data(),icone,ipt);
4166 snprintf(title, buffersize,
"Isolated %s #it{p}_{T} for #it{R}=%2.2f and #Sigma #it{p}_{T}^{in cone}=%2.2f",
4180 for (
Int_t i = 0; i < 7 ; i++)
4183 Form(
"Number of isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4184 nptbins,ptmin,ptmax);
4190 Form(
"Number of isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4191 nptbins,ptmin,ptmax);
4197 Form(
"Number of not isolated particles vs E, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4198 nptbins,ptmin,ptmax);
4204 Form(
"Number of not isolated particles vs #it{p}_{T}, %s, pile-up event by %s",parTitle.Data(),pileUpName[i].Data()),
4205 nptbins,ptmin,ptmax);
4211 fhTimeENoCut =
new TH2F (
"hTimeE_NoCut",
"time of cluster vs E of clusters, no cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4216 fhTimeESPD =
new TH2F (
"hTimeE_SPD",
"time of cluster vs E of clusters, SPD cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4221 fhTimeESPDMulti =
new TH2F (
"hTimeE_SPDMulti",
"time of cluster vs E of clusters, SPD multi cut", nptbins,ptmin,ptmax, ntimebins,timemin,timemax);
4226 fhTimeNPileUpVertSPD =
new TH2F (
"hTime_NPileUpVertSPD",
"time of cluster vs N pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4231 fhTimeNPileUpVertTrack =
new TH2F (
"hTime_NPileUpVertTracks",
"time of cluster vs N pile-up Tracks vertex", ntimebins,timemin,timemax, 50,0,50 );
4236 fhTimeNPileUpVertContributors =
new TH2F (
"hTime_NPileUpVertContributors",
"time of cluster vs N constributors to pile-up SPD vertex", ntimebins,timemin,timemax,50,0,50);
4241 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);
4246 fhTimePileUpMainVertexZDiamond =
new TH2F (
"hTime_PileUpMainVertexZDiamond",
"time of cluster vs distance in Z pile-up SPD vertex - z diamond",ntimebins,timemin,timemax,100,0,50);
4252 return outputContainer ;
4312 AliInfo(
"Open default isolation cuts for multiple Isolation analysis");
4319 AliFatal(
"STOP!: You want to use CTS tracks in analysis but not read!! \n!!Check the configuration file!!");
4384 AliAODPWG4ParticleCorrelation* pLeading = 0;
4390 AliAODPWG4ParticleCorrelation* particle = (AliAODPWG4ParticleCorrelation*) (
GetInputAODBranch()->At(iaod));
4391 particle->SetLeadingParticle(kFALSE);
4397 if(check == 0)
continue;
4398 if(check == -1)
return kFALSE;
4410 if (particle->Pt() > ptTrig)
4412 ptTrig = particle->Pt() ;
4413 phiTrig = particle->Phi();
4415 pLeading = particle ;
4419 if(index < 0)
return kFALSE;
4425 if(phiTrig < 0 ) phiTrig += TMath::TwoPi();
4431 AliVTrack * track = (AliVTrack *) (
GetCTSTracks()->At(ipr)) ;
4439 Bool_t contained = kFALSE;
4441 for(
Int_t i = 0; i < 4; i++)
4443 if( trackID == pLeading->GetTrackLabel(i) ) contained = kTRUE;
4446 if ( contained ) continue ;
4449 fTrackVector.SetXYZ(track->Px(),track->Py(),track->Pz());
4452 if(phi < 0) phi+=TMath::TwoPi();
4456 Float_t deltaPhi = phiTrig-phi;
4462 if(deltaPhi <= -TMath::PiOver2()) deltaPhi+=TMath::TwoPi();
4463 if(deltaPhi > 3*TMath::PiOver2()) deltaPhi-=TMath::TwoPi();
4465 if(pt > ptTrig && deltaPhi < TMath::PiOver2())
return kFALSE;
4475 if (pLeading->GetDetectorTag() ==
kPHOS )
4480 if(!nePl)
return kTRUE;
4482 for(
Int_t ipr = 0;ipr < nePl->GetEntriesFast() ; ipr ++ )
4484 AliVCluster * cluster = (AliVCluster *) (nePl->At(ipr)) ;
4486 if(cluster->GetID() == pLeading->GetCaloLabel(0) || cluster->GetID() == pLeading->GetCaloLabel(1) )
continue ;
4492 if(phi < 0) phi+=TMath::TwoPi();
4499 Float_t deltaPhi = phiTrig-phi;
4500 if(deltaPhi <= -TMath::PiOver2()) deltaPhi+=TMath::TwoPi();
4501 if(deltaPhi > 3*TMath::PiOver2()) deltaPhi-=TMath::TwoPi();
4503 if(pt > ptTrig && deltaPhi < TMath::PiOver2())
return kFALSE ;
4509 pLeading->SetLeadingParticle(kTRUE);
4511 AliDebug(1,Form(
"Particle AOD with index %d is leading with pT %2.2f",idLeading, pLeading->Pt()));
4524 AliFatal(Form(
"No input particles