| 1c5acb87 |
1 | #ifndef ALIANAPARTICLEJETLEADINGCONECORRELATION_H |
| 2 | #define ALIANAPARTICLEJETLEADINGCONECORRELATION_H |
| 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
| 4 | * See cxx source for full Copyright notice */ |
| 5 | /* $Id: $ */ |
| 6 | |
| 7 | //_________________________________________________________________________ |
| 8 | // Class that contains the algorithm for the reconstruction of jet, cone around leading particle |
| 9 | // The seed is a backward particle (direct photon) |
| 78219bac |
10 | // 1) Take the a trigger particle found stored in AliAODPWG4ParticleCorrelation, |
| 1c5acb87 |
11 | // 2) Search for the highest pt leading particle opposite to the trigger within a phi, pt window |
| 12 | // 3) Take all particles around leading in a cone R with pt larger than threshold and construct the jet |
| 13 | // |
| 14 | // Class created from old AliPHOSGammaJet |
| 15 | // (see AliRoot versions previous Release 4-09) |
| 16 | // |
| 17 | //-- Author: Gustavo Conesa (INFN-LNF) |
| 18 | |
| 19 | // --- ROOT system --- |
| 20 | class TH2F; |
| 21 | |
| 22 | //---- Analysis system ---- |
| 23 | class AliAODTrack; |
| 0ae57829 |
24 | class AliVCluster; |
| 1c5acb87 |
25 | class AliCaloTrackReader; |
| 26 | class AliNeutralMesonSelection; |
| 1c5acb87 |
27 | |
| 28 | #include "AliAnaPartCorrBaseClass.h" |
| 29 | |
| 30 | class AliAnaParticleJetLeadingConeCorrelation : public AliAnaPartCorrBaseClass { |
| 31 | |
| 32 | public: |
| 1c5acb87 |
33 | AliAnaParticleJetLeadingConeCorrelation() ; // default ctor |
| 78219bac |
34 | virtual ~AliAnaParticleJetLeadingConeCorrelation() ; //virtual dtor |
| 78219bac |
35 | |
| 1c5acb87 |
36 | TList * GetCreateOutputObjects(); |
| 37 | |
| 38 | void InitParameters(); |
| 39 | |
| 40 | void Print(const Option_t * opt) const; |
| 41 | |
| 42 | Bool_t AreJetsRecalculated() const {return fReMakeJet ; } |
| 43 | void SwitchOnJetsRecalculation(){fReMakeJet = kTRUE; } |
| 44 | void SwitchOffJetsRecalculation(){fReMakeJet = kFALSE; } |
| 45 | |
| 46 | Bool_t AreJetsOnlyInCTS() const {return fJetsOnlyInCTS ; } |
| 47 | void SwitchOnJetsOnlyInCTS(){fJetsOnlyInCTS = kTRUE; } |
| 48 | void SwitchOffJetsOnlyInCTS(){fJetsOnlyInCTS = kFALSE; } |
| 49 | |
| 50 | Bool_t AreSeveralConeAndPtCuts() const {return fSeveralConeAndPtCuts ; } |
| 51 | void SwitchOnSeveralConeAndPtCuts(){fSeveralConeAndPtCuts = kTRUE ;} |
| 52 | void SwitchOffSeveralConeAndPtCuts(){fSeveralConeAndPtCuts = kFALSE ;} |
| 53 | |
| 54 | Bool_t IsPbPb() const {return fPbPb ; } |
| 55 | void SetppCollisions(){fPbPb = kFALSE; } |
| 56 | void SetPbPbCollisions(){fPbPb = kTRUE; } |
| 57 | |
| 58 | Double_t GetDeltaPhiMaxCut() const {return fDeltaPhiMaxCut ; } |
| 59 | Double_t GetDeltaPhiMinCut() const {return fDeltaPhiMinCut ; } |
| 60 | Double_t GetLeadingRatioMaxCut() const {return fLeadingRatioMaxCut ; } |
| 61 | Double_t GetLeadingRatioMinCut() const {return fLeadingRatioMinCut ; } |
| 62 | |
| 63 | Double_t GetPtTriggerSelectionCut() const {return fPtTriggerSelectionCut ; } |
| 64 | Double_t GetJetRatioMaxCut() const {return fJetRatioMaxCut ; } |
| 65 | Double_t GetJetRatioMinCut() const {return fJetRatioMinCut ; } |
| 66 | |
| 67 | void SetPtTriggerSelectionCut(Double_t cut){fPtTriggerSelectionCut = cut; } |
| 68 | void SetJetSelectionMode(UInt_t select){ fSelect= select ; } |
| 69 | |
| 70 | Int_t GetJetNCones() const {return fJetNCone ; } |
| 71 | Int_t GetJetNPtThres() const {return fJetNPt ; } |
| 72 | Float_t GetJetCone() const {return fJetCone ; } |
| 73 | Float_t GetJetPtThreshold() const {return fJetPtThreshold ; } |
| 74 | Float_t GetJetPtThresPbPb() const {return fJetPtThresPbPb ; } |
| 75 | Float_t GetJetCones(Int_t i) const {return fJetCones[i] ; } |
| 76 | Float_t GetJetPtThreshold(Int_t i) const {return fJetPtThres[i] ; } |
| 77 | TString GetJetConeName(Int_t i) const {return fJetNameCones[i] ; } |
| 78 | TString GetJetPtThresName(Int_t i) const {return fJetNamePtThres[i] ; } |
| 79 | |
| 80 | |
| 81 | void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax) |
| 82 | {fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin;} |
| 83 | void SetLeadingRatioCutRange(Double_t ratiomin, Double_t ratiomax) |
| 84 | {fLeadingRatioMaxCut =ratiomax; fLeadingRatioMinCut = ratiomin ; } |
| 85 | |
| 86 | void SetJetNCones(Int_t n){fJetNCone = n ; } |
| 87 | void SetJetNPtThresholds(Int_t n){fJetNPt = n ; } |
| 88 | void SetJetCones(Int_t i, Float_t cone, TString sc) {fJetCones[i] = cone ; fJetNameCones[i] = sc; }; |
| 89 | void SetCone(Float_t cone) {fJetCone = cone; } |
| 90 | void SetJetPtThreshold(Float_t pt){fJetPtThreshold = pt; }; |
| 91 | void SetJetPtThresPbPb(Float_t pt){fJetPtThresPbPb = pt; }; |
| 92 | void SetJetPtThresholds(Int_t i,Float_t pt, TString spt){fJetPtThres[i] = pt ; fJetNamePtThres[i] = spt; }; |
| 93 | |
| 94 | void SetJetRatioCutRange(Double_t ratiomin, Double_t ratiomax) |
| 95 | {fJetRatioMaxCut =ratiomax; fJetRatioMinCut = ratiomin ; } |
| 96 | void SetJetCTSRatioCutRange(Double_t ratiomin, Double_t ratiomax) |
| 97 | {fJetCTSRatioMaxCut =ratiomax; fJetCTSRatioMinCut = ratiomin ; } |
| 98 | |
| 477d6cee |
99 | Bool_t OnlyIsolated() const {return fSelectIsolated ; } |
| 100 | void SelectIsolated(Bool_t select) {fSelectIsolated = select ; } |
| 101 | |
| 1c5acb87 |
102 | private: |
| 103 | |
| 104 | Double_t CalculateJetRatioLimit(const Double_t ptTrig, const Double_t *param, const Double_t *x) const ; |
| 105 | |
| 106 | void FillJetHistos(AliAODPWG4ParticleCorrelation * particle, const TLorentzVector leading, const TLorentzVector jet, const TString type, const TString lastname); |
| 107 | |
| 108 | TList * GetOutputContainer() const {return fOutCont; } |
| 109 | |
| 110 | Bool_t IsJetSelected(const Double_t ptTrig, const Double_t ptjet) const ; |
| 111 | Bool_t IsParticleInJetCone(const Double_t eta, Double_t phi, const Double_t etal, Double_t phil) const ; |
| 112 | |
| 233e0df8 |
113 | void GetLeadingCharge(AliAODPWG4ParticleCorrelation* const particle, TLorentzVector & pLeading) const ; |
| decca433 |
114 | void GetLeadingPi0 (AliAODPWG4ParticleCorrelation* const particle, TLorentzVector & pLeading) ; |
| 115 | Bool_t GetLeadingParticle(AliAODPWG4ParticleCorrelation *particle, TLorentzVector & pLeading) ; |
| 1c5acb87 |
116 | |
| 117 | void MakeAnalysisFillAOD(); |
| 118 | void MakeAnalysisFillHistograms(); |
| 119 | void MakeAODJet(AliAODPWG4ParticleCorrelation * particle, const TLorentzVector pLeading) const ; |
| 120 | void MakeJetFromAOD(AliAODPWG4ParticleCorrelation * particle, const TLorentzVector pLeading, |
| 121 | TLorentzVector & jet, TLorentzVector & bkg) const ; |
| 122 | |
| 0ae57829 |
123 | Bool_t SelectCluster(AliVCluster * calo, Double_t *vertex, TLorentzVector & mom, Int_t & pdg) ; |
| 1c5acb87 |
124 | |
| 125 | private: |
| 126 | |
| 127 | Bool_t fJetsOnlyInCTS ; // Jets measured only in TPC+ITS. |
| 128 | Bool_t fPbPb; // PbPb event |
| 129 | Bool_t fSeveralConeAndPtCuts; // To play with the jet cone size and pt th. |
| 130 | Bool_t fReMakeJet ; //Re make the jet reconstruction from AODParticleCorrelation input |
| 131 | |
| 132 | //Leading particle selection parameters |
| 133 | Double_t fDeltaPhiMaxCut ; // Minimum Delta Phi Gamma-Leading |
| 134 | Double_t fDeltaPhiMinCut ; // Maximum Delta Phi Gamma-Leading |
| 135 | Double_t fLeadingRatioMaxCut ; // Leading /gamma Ratio cut maximum |
| 136 | Double_t fLeadingRatioMinCut ; // Leading/gamma Ratio cut minimum |
| 137 | |
| 138 | //Jet selection parameters |
| 139 | //Fixed cuts (old) |
| 140 | Double_t fJetCTSRatioMaxCut ; // Jet(CTS) /gamma Ratio cut maximum |
| 141 | Double_t fJetCTSRatioMinCut ; // Jet(CTS) /gamma Ratio cut maximum |
| 142 | Double_t fJetRatioMaxCut ; // Jet(EMCAL+CTS)/gamma Ratio cut maximum |
| 143 | Double_t fJetRatioMinCut ; // Jet(EMCAL+CTS)/gamma Ratio cut minimum |
| 144 | |
| 145 | //Cuts depending on jet pt |
| 146 | Double_t fJetE1[2]; //Rec. jet energy parameters |
| 147 | Double_t fJetE2[2]; //Rec. jet energy parameters |
| 148 | Double_t fJetSigma1[2];//Rec. sigma of jet energy parameters |
| 149 | Double_t fJetSigma2[2];//Rec. sigma of jet energy parameters |
| 150 | Double_t fBkgMean[6]; //Background mean energy |
| 151 | Double_t fBkgRMS[6]; //Background RMS |
| 152 | Double_t fJetXMin1[6]; //X Factor to set jet min limit for pp |
| 153 | Double_t fJetXMin2[6]; //X Factor to set jet min limit for PbPb |
| 154 | Double_t fJetXMax1[6]; //X Factor to set jet max limit for pp |
| 155 | Double_t fJetXMax2[6]; //X Factor to set jet max limit for PbPb |
| 156 | |
| 157 | Int_t fJetNCone ; // Number of jet cones sizes, maximum 5 |
| 158 | Int_t fJetNPt ; // Number of jet particle pT threshold, maximum 5 |
| 159 | Double_t fJetCone ; // Jet cone sizes under study (!fSeveralConeAndPtCuts) |
| 160 | Double_t fJetCones[5]; // Jet cone sizes under study (fSeveralConeAndPtCuts) |
| 161 | TString fJetNameCones[5]; // String name of cone to append to histos |
| 162 | Double_t fJetPtThreshold; // Jet pT threshold under study(!fSeveralConeAndPtCuts) |
| 163 | Double_t fJetPtThresPbPb; // Jet pT threshold under study(!fSeveralConeAndPtCuts) |
| 164 | Double_t fJetPtThres[5]; // Jet pT threshold under study(fSeveralConeAndPtCuts) |
| 165 | TString fJetNamePtThres[5]; // String name of pt th to append to histos |
| 166 | Double_t fPtTriggerSelectionCut; // Jet pt to change to low pt jets analysis |
| 167 | UInt_t fSelect ; //kTRUE: Selects all jets, no limits. |
| 477d6cee |
168 | Bool_t fSelectIsolated ; // Select only trigger particles isolated |
| 169 | |
| 1c5acb87 |
170 | //Histograms |
| 171 | //Leading particle distributions |
| 172 | TList * fOutCont ; //! Container for histograms |
| 173 | |
| 174 | TH2F * fhChargedLeadingPt ; //! Pt(Pt trigger) distribution of charged hadrons |
| 175 | TH2F * fhChargedLeadingPhi ; //! Phi(Pt trigger) distribution of charged hadrons |
| 176 | TH2F * fhChargedLeadingEta ; //! Eta(Pt trigger) distribution of charged hadrons |
| 177 | TH2F * fhChargedLeadingDeltaPt ; //! Difference of charged hadron and trigger pT as function of trigger p |
| 178 | TH2F * fhChargedLeadingDeltaPhi ; //! Difference of charged hadron and trigger phi as function of trigger pT |
| 179 | TH2F * fhChargedLeadingDeltaEta ; //! Difference of charged particle and trigger eta as function of trigger pT |
| 180 | TH2F * fhChargedLeadingRatioPt ; //! Ratio of Pt leading charge and trigger |
| 181 | |
| 182 | TH2F * fhNeutralLeadingPt ; //! Pt(Pt trigger) distribution of neutral hadrons |
| 183 | TH2F * fhNeutralLeadingPhi ; //! Phi(Pt trigger) distribution of neutral hadrons |
| 184 | TH2F * fhNeutralLeadingEta ; //! Eta(Pt trigger) distribution of neutral hadrons |
| 185 | TH2F * fhNeutralLeadingDeltaPt ; //! Difference of neutral hadron and trigger pT as function of trigger pT |
| 186 | TH2F * fhNeutralLeadingDeltaPhi ; //! Difference of neutral hadron and trigger phi as function of trigger pT |
| 187 | TH2F * fhNeutralLeadingDeltaEta ; //! Difference of charged particle and trigger eta as function of trigger pT |
| 188 | TH2F * fhNeutralLeadingRatioPt ; //! Ratio of Pt leading neutral and trigger |
| 9415d854 |
189 | |
| 190 | TH2F * fhChargedLeadingXi ; //! Ln (pt leading charge / pt trigger) |
| 191 | TH2F * fhNeutralLeadingXi ; //! Ln (pt leading neutral / pt trigger) |
| dde5a268 |
192 | |
| 193 | TH2F * fhChargedLeadingDeltaPhiRatioPt30 ; //! Difference of charged hadron and trigger phi as function of pT leading / trigger pT, pT Trigger > 30 GeV |
| 194 | TH2F * fhNeutralLeadingDeltaPhiRatioPt30 ; //! Difference of neutral hadron and trigger phi as function of pT leading / trigger pT, pT Trigger > 30 GeV |
| 195 | TH2F * fhChargedLeadingDeltaPhiRatioPt50 ; //! Difference of charged hadron and trigger phi as function of pT leading / trigger pT, pT Trigger > 50 GeV |
| 196 | TH2F * fhNeutralLeadingDeltaPhiRatioPt50 ; //! Difference of neutral hadron and trigger phi as function of pT leading / trigger pT, pT Trigger > 50 GeV |
| 9415d854 |
197 | |
| 1c5acb87 |
198 | // Jet distributions |
| 199 | // Fixed cone and pt threshold |
| 200 | TH2F * fhJetPt ; //! leading pt jet vs pt trigger |
| 201 | TH2F * fhJetRatioPt ; //! Ratio of pt jet and pt trigger |
| 202 | TH2F * fhJetDeltaPhi ; //! Delta phi jet-trigger |
| 203 | TH2F * fhJetDeltaEta ; //! Delta eta jet-trigger |
| 204 | TH2F * fhJetLeadingRatioPt ; //! Ratio of pt leading and pt jet |
| 205 | TH2F * fhJetLeadingDeltaPhi ; //! Delta phi jet-leading |
| 206 | TH2F * fhJetLeadingDeltaEta ; //! Delta eta jet-leading |
| 207 | TH2F * fhJetFFz; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig |
| 208 | TH2F * fhJetFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet) |
| 209 | TH2F * fhJetFFpt; //! Jet particle pt distribution in cone |
| 210 | TH2F * fhJetNTracksInCone ; //! jet multiplicity in cone |
| 211 | |
| 212 | TH2F * fhBkgPt ; //! leading pt bakground vs pt trigger |
| 213 | TH2F * fhBkgRatioPt ; //! Ratio of pt background and pt trigger |
| 214 | TH2F * fhBkgDeltaPhi ; //! Delta phi background-trigger |
| 215 | TH2F * fhBkgDeltaEta ; //! Delta eta background-trigger |
| 216 | TH2F * fhBkgLeadingRatioPt ; //! Ratio of pt leading and pt background |
| 217 | TH2F * fhBkgLeadingDeltaPhi ; //! Delta phi background-leading |
| 218 | TH2F * fhBkgLeadingDeltaEta ; //! Delta eta background-leading |
| 219 | TH2F * fhBkgFFz; //! Accepted reconstructed background fragmentation function, z=ptjet/pttrig |
| 220 | TH2F * fhBkgFFxi; //! Accepted reconstructed background fragmentation function, xsi = ln(pttrig/ptjet) |
| 221 | TH2F * fhBkgFFpt; //! Background particle pt distribution in cone |
| 222 | TH2F * fhBkgNTracksInCone ; //! Background multiplicity in cone |
| 223 | |
| 224 | // Variable cone and pt threshold |
| 225 | |
| 226 | TH2F * fhJetPts[5][5]; //! leading pt jet vs pt trigger |
| 227 | TH2F * fhJetRatioPts[5][5]; //! Ratio of pt jet and pt trigger |
| 228 | TH2F * fhJetDeltaPhis[5][5]; //! Delta phi jet-trigger |
| 229 | TH2F * fhJetDeltaEtas[5][5]; //! Delta eta jet-trigger |
| 230 | TH2F * fhJetLeadingRatioPts[5][5]; //! Ratio of pt leading and pt jet |
| 231 | TH2F * fhJetLeadingDeltaPhis[5][5]; //! Delta phi jet-leading |
| 232 | TH2F * fhJetLeadingDeltaEtas[5][5]; //! Delta eta jet-leading |
| 233 | TH2F * fhJetFFzs[5][5]; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig |
| 234 | TH2F * fhJetFFxis[5][5]; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet) |
| 235 | TH2F * fhJetFFpts[5][5]; //! Jet particle pt distribution in cone |
| 236 | TH2F * fhJetNTracksInCones[5][5]; //! jet multiplicity in cone |
| 237 | |
| 238 | TH2F * fhBkgPts[5][5]; //! leading pt bakground vs pt trigger |
| 239 | TH2F * fhBkgRatioPts[5][5]; //! Ratio of pt background and pt trigger |
| 240 | TH2F * fhBkgDeltaPhis[5][5]; //! Delta phi background-trigger |
| 241 | TH2F * fhBkgDeltaEtas[5][5]; //! Delta eta background-trigger |
| 242 | TH2F * fhBkgLeadingRatioPts[5][5]; //! Ratio of pt leading and pt background |
| 243 | TH2F * fhBkgLeadingDeltaPhis[5][5]; //! Delta phi background-leading |
| 244 | TH2F * fhBkgLeadingDeltaEtas[5][5]; //! Delta eta background-leading |
| 245 | TH2F * fhBkgFFzs[5][5]; //! Accepted reconstructed background fragmentation function, z=ptjet/pttrig |
| 246 | TH2F * fhBkgFFxis[5][5]; //! Accepted reconstructed background fragmentation function, xsi = ln(pttrig/ptjet) |
| 247 | TH2F * fhBkgFFpts[5][5]; //! Background particle pt distribution in cone |
| 248 | TH2F * fhBkgNTracksInCones[5][5]; //! Background multiplicity in cone |
| 249 | |
| c5693f62 |
250 | AliAnaParticleJetLeadingConeCorrelation(const AliAnaParticleJetLeadingConeCorrelation & g) ; // cpy ctor |
| 251 | AliAnaParticleJetLeadingConeCorrelation & operator = (const AliAnaParticleJetLeadingConeCorrelation & g) ;//cpy assignment |
| 1c5acb87 |
252 | |
| 253 | ClassDef(AliAnaParticleJetLeadingConeCorrelation,1) |
| 254 | } ; |
| 255 | |
| 256 | |
| 257 | #endif //ALIANAPARTICLEJETLEADINGCONECORRELATION_H |
| 258 | |
| 259 | |
| 260 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame