]>
Commit | Line | Data |
---|---|---|
1c5acb87 | 1 | #ifndef ALIANAPARTICLEJETFINDERCORRELATION_H |
2 | #define ALIANAPARTICLEJETFINDERCORRELATION_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
1c5acb87 | 5 | |
6 | //_________________________________________________________________________ | |
7 | // Class that contains the algorithm for the analysis of particle (direct gamma) - jet | |
8 | // (standard jet found with JETAN) correlation | |
9 | // Particle and jet for correlation found by independent algorithms. | |
10 | // For Example direct isolated photon found in AliAnaGammaDirect and the jet with JETAN | |
11 | // | |
12 | //-- Author: Gustavo Conesa (INFN-LNF) | |
290c1e1b | 13 | //-- Modified by Adam Matyja (INP PAN, Krakow) |
1c5acb87 | 14 | |
15 | // --- ROOT system --- | |
16 | class TH2F; | |
290c1e1b | 17 | class TTree; |
18 | class TRandom2; | |
1c5acb87 | 19 | //---- Analysis system ---- |
745913ae | 20 | #include "AliAnaCaloTrackCorrBaseClass.h" |
1c5acb87 | 21 | |
745913ae | 22 | class AliAnaParticleJetFinderCorrelation : public AliAnaCaloTrackCorrBaseClass { |
1c5acb87 | 23 | |
78219bac | 24 | public: |
745913ae | 25 | AliAnaParticleJetFinderCorrelation() ; // default ctor |
290c1e1b | 26 | virtual ~AliAnaParticleJetFinderCorrelation() ; // virtual dtor |
78219bac | 27 | |
c5693f62 | 28 | // General methods |
477d6cee | 29 | |
c5693f62 | 30 | void InitParameters(); |
477d6cee | 31 | |
c5693f62 | 32 | TList * GetCreateOutputObjects(); |
33 | ||
34 | void MakeAnalysisFillAOD() ; | |
477d6cee | 35 | |
c5693f62 | 36 | void MakeAnalysisFillHistograms() ; |
477d6cee | 37 | |
290c1e1b | 38 | Int_t SelectJet(AliAODPWG4Particle * particle, TClonesArray * aodRecJets) ;//to access non standard branch |
c5693f62 | 39 | |
40 | void Print(const Option_t * opt) const; | |
477d6cee | 41 | |
c5693f62 | 42 | // Settings |
477d6cee | 43 | |
c5693f62 | 44 | Bool_t OnlyIsolated() const { return fSelectIsolated ; } |
45 | void SelectIsolated(Bool_t select) { fSelectIsolated = select ; } | |
477d6cee | 46 | |
c5693f62 | 47 | Float_t GetConeSize() const { return fConeSize ; } |
48 | Float_t GetPtThresholdInCone() const { return fPtThresholdInCone ; } | |
49 | Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; } | |
50 | Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; } | |
51 | Double_t GetRatioMaxCut() const { return fRatioMaxCut ; } | |
52 | Double_t GetRatioMinCut() const { return fRatioMinCut ; } | |
53 | Bool_t AreJetRefTracks() const { return fUseJetRefTracks ; } | |
54 | Bool_t IsCorrelationMadeInHistoMaker() const { return fMakeCorrelationInHistoMaker ; } | |
290c1e1b | 55 | Double_t GetJetConeSize() const { return fJetConeSize ; } |
56 | Double_t GetJetMinPt() const { return fJetMinPt ; } | |
57 | Double_t GetJetAreaFraction() const { return fJetAreaFraction ; } | |
58 | Double_t GetGammaConeSize() const { return fGammaConeSize ; } | |
59 | ||
c5693f62 | 60 | void SetConeSize(Float_t cone) { fConeSize = cone ; } |
61 | void SetPtThresholdInCone(Float_t pt) { fPtThresholdInCone = pt ; } | |
62 | void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax) | |
63 | { fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin ; } | |
64 | void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax) | |
65 | { fRatioMaxCut =ratiomax; fRatioMinCut = ratiomin ; } | |
66 | void UseJetRefTracks(Bool_t use) { fUseJetRefTracks = use ; } | |
290c1e1b | 67 | void SetMakeCorrelationInHistoMaker(Bool_t make) { fMakeCorrelationInHistoMaker = make ; } |
68 | void SetJetConeSize(Double_t cone) { fJetConeSize = cone ; } | |
69 | void SetJetMinPt(Double_t minpt) { fJetMinPt = minpt ; } | |
70 | void SetJetAreaFraction(Double_t areafr) { fJetAreaFraction = areafr ; } | |
71 | void SetGammaConeSize(Float_t cone) { fGammaConeSize = cone ; } | |
72 | ||
73 | // Settings for non standard jet branch | |
74 | TString GetJetBranchName() const { return fJetBranchName ; } | |
75 | void SetJetBranchName(const char *name) { fJetBranchName = name ; } | |
9717eff6 | 76 | //void SwitchOnNonStandardJetFromReader() { fNonStandardJetFromReader = kTRUE ; } |
77 | //void SwitchOffNonStandardJetFromReader() { fNonStandardJetFromReader = kFALSE; } | |
78 | //Bool_t IsNonStandardJetFromReader() { return fNonStandardJetFromReader ; } | |
290c1e1b | 79 | |
80 | TString GetBkgJetBranchName() const { return fBkgJetBranchName ; } | |
81 | void SetBkgJetBranchName(const char *name) { fBkgJetBranchName = name ; } | |
82 | void SwitchOnBackgroundJetFromReader() { fBackgroundJetFromReader = kTRUE ; } | |
83 | void SwitchOffBackgroundJetFromReader() { fBackgroundJetFromReader = kFALSE ; } | |
84 | Bool_t IsBackgroundJetFromReader() { return fBackgroundJetFromReader ; } | |
85 | ||
86 | //switches for photons | |
87 | void SwitchOnBackgroundSubtractionGamma() { fUseBackgroundSubtractionGamma = kTRUE ; } | |
88 | void SwitchOffBackgroundSubtractionGamma() { fUseBackgroundSubtractionGamma = kFALSE; } | |
89 | Bool_t IsBackgroundSubtractionGamma() { return fUseBackgroundSubtractionGamma ; } | |
90 | ||
91 | void CalculateBkg(TVector3 gamma, TVector3 jet,Double_t *vector,Int_t type); | |
92 | ||
2ac4a57a | 93 | void FindMCgenInfo();//gives information on generated level |
94 | ||
290c1e1b | 95 | void SwitchOnSaveGJTree() { fSaveGJTree = kTRUE ; } |
96 | void SwitchOffSaveGJTree() { fSaveGJTree = kFALSE; } | |
97 | Bool_t IsSaveGJTree() { return fSaveGJTree ; } | |
98 | ||
99 | void SwitchOnMostEnergetic() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE; } | |
100 | void SwitchOffMostEnergetic() { fMostEnergetic = kFALSE; fMostOpposite = kTRUE ; } | |
101 | void SwitchOffMostOpposite() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE; } | |
102 | void SwitchOnMostOpposite() { fMostEnergetic = kFALSE; fMostOpposite = kTRUE ; } | |
103 | Bool_t IsMostEnergetic() { return fMostEnergetic ; } | |
104 | Bool_t IsMostOpposite() { return fMostOpposite; } | |
105 | ||
106 | //switches for histograms | |
107 | void SwitchOnHistogramJetBkg() { fUseHistogramJetBkg = kTRUE ; } | |
108 | void SwitchOffHistogramJetBkg() { fUseHistogramJetBkg = kFALSE; } | |
109 | Bool_t IsHistogramJetBkg() { return fUseHistogramJetBkg ; } | |
110 | ||
111 | void SwitchOnHistogramTracks() { fUseHistogramTracks = kTRUE ; } | |
112 | void SwitchOffHistogramTracks() { fUseHistogramTracks = kFALSE; } | |
113 | Bool_t IsHistogramTracks() { return fUseHistogramTracks ; } | |
114 | ||
115 | void SwitchOnHistogramJetTracks() { fUseHistogramJetTracks = kTRUE ; } | |
116 | void SwitchOffHistogramJetTracks() { fUseHistogramJetTracks = kFALSE; } | |
117 | Bool_t IsHistogramJetTracks() { return fUseHistogramJetTracks ; } | |
118 | ||
2ac4a57a | 119 | void SwitchOnMCStudies() { fMCStudies = kTRUE ; } |
120 | void SwitchOffMCStudies() { fMCStudies = kFALSE; } | |
121 | Bool_t IsMCStudies() { return fMCStudies ; } | |
122 | ||
c5693f62 | 123 | private: |
124 | ||
477d6cee | 125 | //selection parameters |
c5693f62 | 126 | Double_t fDeltaPhiMaxCut ; //! Minimum Delta Phi Gamma-Leading |
290c1e1b | 127 | Double_t fDeltaPhiMinCut ; //! Maximum Delta Phi Gamma-Leading |
128 | Double_t fRatioMaxCut ; //! Jet/particle Ratio cut maximum | |
c5693f62 | 129 | Double_t fRatioMinCut ; //! Jet/particle Ratio cut minimum |
477d6cee | 130 | |
290c1e1b | 131 | Double_t fConeSize ; //! Jet cone size to calculate fragmentation function |
477d6cee | 132 | Double_t fPtThresholdInCone ; //! Jet pT threshold in jet cone |
290c1e1b | 133 | Bool_t fUseJetRefTracks ; //! Use track references from JETAN not the AOD tracks to calculate fragmentation function |
134 | Bool_t fMakeCorrelationInHistoMaker ; //!Make particle-jet correlation in histogram maker | |
135 | Bool_t fSelectIsolated ; //! Select only trigger particles isolated | |
c5693f62 | 136 | |
290c1e1b | 137 | Double_t fJetConeSize ; //! Reconstructed jet cone size |
138 | Double_t fJetMinPt ; //! Minumum jet pt, default 5GeV/c | |
139 | Double_t fJetAreaFraction ; //! Jet area fraction X in X*pi*R^2, default 0.6 | |
9717eff6 | 140 | //Bool_t fNonStandardJetFromReader; //! use non standard jet from reader //new |
290c1e1b | 141 | TString fJetBranchName ;//! name of jet branch not set in reader part //new |
142 | Bool_t fBackgroundJetFromReader; //! use background jet from reader //new | |
143 | TString fBkgJetBranchName ;//! name of background jet branch not set in reader part //new | |
144 | ||
145 | Double_t fGammaConeSize ; //! Isolation cone radius | |
146 | Bool_t fUseBackgroundSubtractionGamma;//! flag to use backgrouind subtraction for photons or not | |
147 | Bool_t fSaveGJTree;//! flag to save gamma-jet tree | |
148 | Bool_t fMostEnergetic;//! flag to choose gamma-jet pairs most energetic | |
149 | Bool_t fMostOpposite;//! flag to choose gamma-jet pairs most opposite | |
150 | ||
151 | Bool_t fUseHistogramJetBkg;//! flag to save bkg jet histograms | |
152 | Bool_t fUseHistogramTracks;//! flag to save CTS tracks features | |
153 | Bool_t fUseHistogramJetTracks;//! flag to save jet tracks features | |
154 | ||
2ac4a57a | 155 | Bool_t fMCStudies; //! flag to use MC methods |
156 | ||
290c1e1b | 157 | TRandom2 * fGenerator;//! pointer to random generator object |
158 | ||
1a8c88c1 | 159 | TLorentzVector fMomentum; //! momentum |
160 | ||
c5693f62 | 161 | // Histograms |
290c1e1b | 162 | TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT |
163 | //TH2F * fhDeltaPhi; //! Difference of jet phi and trigger particle phi as function of trigger particle pT | |
164 | TH2F * fhDeltaPhiCorrect; //! Difference of jet phi and trigger particle phi as function of trigger particle pT | |
165 | TH2F * fhDeltaPhi0PiCorrect; //! Difference of jet phi and trigger particle phi as function of trigger particle pT | |
166 | ||
c5693f62 | 167 | TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT |
168 | TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT | |
169 | TH2F * fhPt; //! jet pT vs trigger particle pT | |
477d6cee | 170 | |
290c1e1b | 171 | TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=pt^particle,jet/pttrig |
172 | TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/pt^particle,jet) | |
c5693f62 | 173 | TH2F * fhFFpt; //! Jet particle pt distribution in cone |
174 | TH2F * fhNTracksInCone; //! jet multiplicity in cone | |
477d6cee | 175 | |
290c1e1b | 176 | TH2F * fhJetFFz ; //! Accepted reconstructed jet fragmentation function, z=pt^particle,jet/ptjet |
177 | TH2F * fhJetFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(ptjet/pt^particle,jet) | |
178 | TH2F * fhJetFFpt; //! Jet particle pt distribution in jet cone | |
179 | TH2F * fhJetFFzCor ; //! Accepted reconstructed jet fragmentation function, z=pt^particle,jet*-cos(jet,trig)/ptjet | |
180 | TH2F * fhJetFFxiCor; //! Accepted reconstructed jet fragmentation function, xsi = ln(ptjet/pt^particle*-cos(jet,trig),jet) | |
181 | ||
182 | //background from RC | |
183 | TH2F * fhBkgFFz[5] ; //! Background fragmentation function, z=ptjet/pttrig | |
184 | TH2F * fhBkgFFxi[5]; //! Background fragmentation function, xsi = ln(pttrig/ptjet) | |
185 | TH2F * fhBkgFFpt[5]; //! Background particle pt distribution in cone | |
186 | TH2F * fhBkgNTracksInCone[5]; //! Background multiplicity in cone | |
187 | TH2F * fhBkgSumPtInCone[5]; //! Background sum pt in cone | |
188 | TH2F * fhBkgSumPtnTracksInCone[5];//! Background sum pt over ntracks in cone | |
189 | ||
190 | //temporary histograms | |
191 | TH2F * fhNjetsNgammas; //! Number of jets vs number of photons in the event | |
192 | TH1F * fhCuts; //! Number of events after cuts | |
193 | ||
194 | TH2F * fhDeltaEtaBefore; //! Difference of jet eta and trigger particle eta as function of trigger particle pT | |
195 | TH2F * fhDeltaPhiBefore; //! Difference of jet phi and trigger particle phi as function of trigger particle pT | |
196 | TH2F * fhDeltaPtBefore; //! Difference of jet pT and trigger particle pT as function of trigger particle pT | |
197 | TH2F * fhPtRatioBefore; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT | |
198 | TH2F * fhPtBefore; //! jet pT vs trigger particle pT | |
199 | TH2F * fhDeltaPhi0PiCorrectBefore; //! Difference of jet phi and trigger particle phi (0,pi) as function of trigger particle pT | |
200 | ||
201 | //temporary jet histograms | |
202 | TH1F * fhJetPtBefore; //! Pt of all jets | |
203 | TH1F * fhJetPt; //! Pt of all jets after bkg correction | |
204 | TH1F * fhJetPtMostEne; //! Pt of the most energetic jet | |
205 | TH1F * fhJetPhi; //! Phi of all jets | |
206 | TH1F * fhJetEta; //! Eta of all jets | |
207 | TH2F * fhJetEtaVsPt; //! Eta of all jets vs pt | |
208 | TH2F * fhJetPhiVsEta; //! Phi vs eta of all jets | |
209 | TH2F * fhJetEtaVsNpartInJet; //! Eta vs number of particles in jet for all jets | |
210 | TH2F * fhJetEtaVsNpartInJetBkg; //! Eta vs number of particles in jet for background subtracted jets | |
211 | TH2F * fhJetChBkgEnergyVsPt; //! background energy of each charged jet vs jet pt | |
212 | TH2F * fhJetChAreaVsPt; //! area of each charged jet vs jet pt | |
213 | TH2F * fhTrackPhiVsEta; //! Phi vs eta of all chosen tracks in all events | |
214 | TH1F * fhTrackAveTrackPt; //! average track pt in event | |
215 | TH1F * fhJetNjetOverPtCut[10]; //! number of reconstructed jets in event over pT threshold | |
216 | TH2F * fhJetChBkgEnergyVsArea; //! area of each charged jet vs jet background | |
217 | TH2F * fhJetRhoVsPt; //! jet energy density vs jet pt | |
218 | TH2F * fhJetRhoVsCentrality; //! jet energy density vs centrality | |
219 | TH1F * fhJetNparticlesInJet; //! number of particles in jets | |
220 | TH2F * fhJetDeltaEtaDeltaPhi; //! delta eta vs delta phi for (jet-track) <-0.8,0.8> | |
221 | TH2F * fhJetDeltaEtaDeltaPhiAllTracks;//! delta eta vs delta phi for (jet-track) <-pi,pi> | |
222 | ||
223 | ||
224 | ||
225 | TH1F * fhJetAveTrackPt; //! average track from jets pt in event | |
226 | TH2F * fhJetNtracksInJetAboveThr[6]; //! number of tracks in jet with pt above 0,1,2,3,4,5GeV | |
227 | TH2F * fhJetRatioNTrkAboveToNTrk[5]; //! ratio tracks in jet with pt above 1,2,3,4,5GeV to ntracks | |
228 | TH2F * fhJetNtrackRatioMostEne[5]; //! the same for most energetic jet | |
229 | TH2F * fhJetNtrackRatioJet5GeV[5]; //! the same for pt jet above 5 GeV | |
230 | TH2F * fhJetNtrackRatioLead5GeV[5]; //! the same for jet with leading particle pt>5GeV | |
231 | ||
232 | //temporary background jet histograms | |
233 | TH1F * fhBkgJetBackground[4]; //! background from jet bkg branch | |
234 | TH1F * fhBkgJetSigma[4]; //! sigma of jet in backgroud branch | |
235 | TH1F * fhBkgJetArea[4]; //! area of jet in bkg branch | |
236 | ||
237 | //temporary photon histograms | |
238 | TH1F * fhPhotonPtMostEne; //! most pt photon | |
239 | TH1F * fhPhotonAverageEnergy; //! average energy of photon | |
240 | TH1F * fhPhotonRatioAveEneToMostEne; //! ratio average energy to most energetic photon | |
241 | TH1F * fhPhotonAverageEnergyMinus1; //! average energy of photon w/o most ene photon | |
242 | TH1F * fhPhotonRatioAveEneMinus1ToMostEne; //! ratio average energy of photon w/o most ene photon to most energetic photon | |
243 | TH1F * fhPhotonNgammaMoreAverageToNgamma; //! number of gammas with ene. more than average ene divided by no. of gammas | |
244 | TH1F * fhPhotonNgammaMoreAverageMinus1ToNgamma; //! number of gammas with ene. more than average ene (w/o most ene gamma) divided by no. of gammas | |
245 | TH1F * fhPhotonNgammaOverPtCut[10]; //! number of photons in event over pT threshold | |
246 | TH2F * fhPhotonBkgRhoVsNtracks; //! average energy in one cell vs n tracks | |
247 | TH2F * fhPhotonBkgRhoVsNclusters; //! average energy in one cell vs n clusters | |
248 | TH2F * fhPhotonBkgRhoVsCentrality; //! average energy in one cell vs centrality | |
249 | TH2F * fhPhotonBkgRhoVsNcells; //! average energy in one cell vs n cells | |
250 | TH1F * fhPhotonPt; //! pt of gamma before bkg correction | |
251 | TH1F * fhPhotonPtCorrected; //! pt of gamma after background correction | |
252 | TH1F * fhPhotonPtCorrectedZoom; //! pt of gamma after background correction in +-5 GeV/c | |
253 | TH1F * fhPhotonPtDiff; //! bkg correction = n_cells * median_rho | |
254 | TH2F * fhPhotonPtDiffVsCentrality; //! correction vs centrality | |
255 | TH2F * fhPhotonPtDiffVsNcells; //! correction vs Ncells | |
256 | TH2F * fhPhotonPtDiffVsNtracks; //! correction vs Ntracks | |
257 | TH2F * fhPhotonPtDiffVsNclusters; //! correction vs Nclustres | |
258 | ||
259 | TH1F * fhPhotonSumPtInCone; //! sum pt in cone before correction | |
260 | TH1F * fhPhotonSumPtCorrectInCone; //! sum pt in cone afrer correction | |
261 | TH1F * fhPhotonSumPtChargedInCone; //! sum pt of charged tracks in the cone before correction | |
262 | ||
263 | ||
264 | //temporary jet histograms after selection | |
265 | TH2F * fhSelectedJetPhiVsEta; //! phi vs eta of selected jet | |
266 | TH2F * fhSelectedJetChBkgEnergyVsPtJet; //! background energy of selected charged jet vs jet pt | |
267 | TH2F * fhSelectedJetChAreaVsPtJet; //! area of selected charged jet vs jet pt | |
268 | TH1F * fhSelectedJetNjet; //! number of jets in selected event | |
269 | TH1F * fhSelectedNtracks; //! number of tracks in selected event | |
270 | TH2F * fhSelectedTrackPhiVsEta; //! Phi vs eta of all chosen tracks in selected events | |
271 | ||
272 | TH1F * fhCuts2; //! efficienct cuts | |
273 | ||
274 | //temporary photon histogram after selection | |
275 | TH2F * fhSelectedPhotonNLMVsPt; //! nlm vs pt for selected photons | |
276 | TH2F * fhSelectedPhotonLambda0VsPt; //! lambda0 vs pt for selected photons | |
277 | TH2F * fhRandomPhiEta[5]; //! eta and phi from random generator | |
278 | ||
2ac4a57a | 279 | //MC generated histograms |
280 | TH1F * fhMCPhotonCuts; //! generated photon cuts | |
281 | TH1F * fhMCPhotonPt; //! generated direct photon pt | |
282 | TH2F * fhMCPhotonEtaPhi;//! generated direct photon eta vs phi | |
283 | TH1F * fhMCJetOrigin; //! generated origin of jet | |
284 | TH2F * fhMCJetNPartVsPt; //! generated N parts vs pt full jet | |
285 | TH2F * fhMCJetChNPartVsPt; //! generated N parts vs pt charged jet | |
286 | TH2F * fhMCJetNPart150VsPt; //! generated N parts (pt>150 MeV/c) vs pt full jet | |
287 | TH2F * fhMCJetChNPart150VsPt;//! generated N parts (pt>150 MeV/c) vs pt charged jet | |
288 | TH2F * fhMCJetChNPart150ConeVsPt;//! generated N parts (pt>150 MeV/c) vs pt charged jet R=0.4 | |
289 | TH1F * fhMCJetRatioChFull; //! generated ratio pt charged/full jet | |
290 | TH1F * fhMCJetRatioCh150Ch; //! generated ratio pt charged(pt>150MeV/c)/charged jet | |
291 | TH2F * fhMCJetEtaPhi; //! generated jet eta vs phi for full jet | |
292 | TH2F * fhMCJetChEtaPhi; //! generated jet eta vs phi for charged jet | |
293 | TH2F * fhMCJet150EtaPhi; //! generated jet eta vs phi full jet (pt>150 MeV/c) | |
294 | TH2F * fhMCJetCh150EtaPhi;//! generated jet eta vs phi charged jet (pt>150 MeV/c) | |
295 | TH2F * fhMCJetCh150ConeEtaPhi;//! generated jet eta vs phi charged jet (pt>150 MeV/c) R=0.4 | |
296 | ||
290c1e1b | 297 | //tree with data gamma and jet |
298 | TTree * fTreeGJ; //! gamma-jet tree | |
299 | Double_t fGamPt; //! pt | |
300 | Double_t fGamLambda0; //! lambda 0 | |
301 | Int_t fGamNLM; //! NLM | |
302 | Double_t fGamSumPtCh; //! energy in isolation cone charged | |
303 | Double_t fGamTime; //! time | |
304 | Int_t fGamNcells; //! ncells | |
305 | Double_t fGamEta; //! eta photon | |
306 | Double_t fGamPhi; //! phi photon | |
307 | Double_t fGamSumPtNeu; //! energy in isolation cone neutral | |
308 | Int_t fGamNtracks; //! number of tracks in iso cone | |
309 | Int_t fGamNclusters; //! number of clusters in iso cone | |
310 | Double_t fGamAvEne; //! average energy of photons (without most ene) | |
311 | Double_t fJetPhi; //! jet phi | |
312 | Double_t fJetEta; //! eta phi | |
313 | Double_t fJetPt; //! jet pt | |
314 | Double_t fJetBkgChEne; //! bkg energy of jet | |
315 | Double_t fJetArea; //! jet area | |
316 | Int_t fJetNtracks; //! number of jet tracks | |
317 | Int_t fJetNtracks1; //! number of jet tracks with pt>1 GeV/c | |
318 | Int_t fJetNtracks2; //! number of jet tracks with pt>2 GeV/c | |
319 | Double_t fJetRho; //! jet rho in event | |
320 | Int_t fEventNumber; //! event number | |
321 | Int_t fNtracks; //! n tracks in event | |
322 | Double_t fZvertex; //! z vertex | |
323 | Double_t fCentrality; //! centrality | |
324 | Bool_t fIso; //! flag isolated or not | |
325 | Double_t fGamRho; //! background energy for photons per cell in EMCal | |
326 | ||
2ac4a57a | 327 | Double_t fMCGamPt; //! MC gen pt photon |
328 | Double_t fMCGamEta; //! MC gen eta photon | |
329 | Double_t fMCGamPhi; //! MC gen phi photon | |
330 | Int_t fMCPartonType; //! MC gen parton type origin of jet | |
331 | Double_t fMCJetPt; //! MC gen full jet pt | |
332 | Double_t fMCJetChPt; //! MC gen charged jet pt | |
333 | Double_t fMCJet150Pt; //! MC gen full jet (pt^particles>150MeV/c) pt | |
334 | Double_t fMCJetCh150Pt; //! MC gen charged jet (pt^particles>150MeV/c) pt | |
335 | Int_t fMCJetNPart; //! MC gen number of full jet particles | |
336 | Int_t fMCJetChNPart; //! MC gen number of charged jet particles | |
337 | Int_t fMCJet150NPart; //! MC gen number of full jet particles (pt>150MeV/c) | |
338 | Int_t fMCJetCh150NPart;//! MC gen number of charged jet particles (pt>150MeV/c) | |
339 | Double_t fMCJetEta; //! MC gen full jet eta | |
340 | Double_t fMCJetPhi; //! MC gen full jet phi | |
341 | Double_t fMCJetChEta; //! MC gen charged jet eta | |
342 | Double_t fMCJetChPhi; //! MC gen charged jet phi | |
343 | Double_t fMCJet150Eta; //! MC gen full jet eta (pt>150MeV/c) | |
344 | Double_t fMCJet150Phi; //! MC gen full jet phi (pt>150MeV/c) | |
345 | Double_t fMCJetCh150Eta; //! MC gen charged jet eta (pt>150MeV/c) | |
346 | Double_t fMCJetCh150Phi; //! MC gen charged jet phi (pt>150MeV/c) | |
347 | ||
348 | Double_t fMCJetCh150ConePt;//! MC gen charged jet (pt^particles>150MeV/c),R=0.4 pt | |
349 | Int_t fMCJetCh150ConeNPart;//! MC gen number of charged jet particles (pt>150MeV/c),R=0.4 | |
350 | Double_t fMCJetCh150ConeEta; //! MC gen charged jet eta (pt>150MeV/c),R=0.4 | |
351 | Double_t fMCJetCh150ConePhi; //! MC gen charged jet phi (pt>150MeV/c),R=0.4 | |
352 | ||
745913ae | 353 | AliAnaParticleJetFinderCorrelation( const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor |
c5693f62 | 354 | AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ; // cpy assignment |
477d6cee | 355 | |
290c1e1b | 356 | ClassDef(AliAnaParticleJetFinderCorrelation,3) |
c5693f62 | 357 | |
1c5acb87 | 358 | } ; |
359 | ||
1c5acb87 | 360 | #endif //ALIANAPARTICLEJETFINDERCORRELATION_H |
361 | ||
362 | ||
363 |