AliPhysics  f05a842 (f05a842)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnaParticleJetFinderCorrelation.h
Go to the documentation of this file.
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 */
5 
6 //_________________________________________________________________________
21 
22 // --- ROOT system ---
23 class TH2F;
24 class TTree;
25 class TRandom2;
26 
27 //---- Analysis system ----
29 
31 
32  public:
33 
35 
37 
38  // General methods
39 
40  void InitParameters();
41 
43 
44  void MakeAnalysisFillAOD() ;
45 
47 
48  // To access non standard branch
49  Int_t SelectJet(AliAODPWG4Particle * particle, TClonesArray * aodRecJets) ;
50 
51  void Print(const Option_t * opt) const;
52 
53  // Settings
54 
55  Bool_t OnlyIsolated() const { return fSelectIsolated ; }
56  void SelectIsolated(Bool_t select) { fSelectIsolated = select ; }
57 
58  Float_t GetConeSize() const { return fConeSize ; }
62  Double_t GetRatioMaxCut() const { return fRatioMaxCut ; }
63  Double_t GetRatioMinCut() const { return fRatioMinCut ; }
66  Double_t GetJetConeSize() const { return fJetConeSize ; }
67  Double_t GetJetMinPt() const { return fJetMinPt ; }
71 
72  void SetConeSize(Float_t cone) { fConeSize = cone ; }
74 
76  { fDeltaPhiMaxCut = phimax ; fDeltaPhiMinCut = phimin ; }
77 
78  void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
79  { fRatioMaxCut = ratiomax; fRatioMinCut = ratiomin ; }
80 
81  void UseJetRefTracks(Bool_t use) { fUseJetRefTracks = use ; }
83  void SetJetConeSize(Double_t cone) { fJetConeSize = cone ; }
84  void SetJetMinPt(Double_t minpt) { fJetMinPt = minpt ; }
85  void SetJetMinPtBkgSub(Double_t minpt) { fJetMinPtBkgSub = minpt ; }
86  void SetJetAreaFraction(Double_t areafr) { fJetAreaFraction = areafr ; }
87  void SetGammaConeSize(Float_t cone) { fGammaConeSize = cone ; }
88 
89  // Settings for non standard jet branch
91  void SetJetBranchName(const char *name) { fJetBranchName = name ; }
92 //void SwitchOnNonStandardJetFromReader() { fNonStandardJetFromReader = kTRUE ; }
93 //void SwitchOffNonStandardJetFromReader() { fNonStandardJetFromReader = kFALSE ; }
94 //Bool_t IsNonStandardJetFromReader() { return fNonStandardJetFromReader ; }
95 
97  void SetBkgJetBranchName(const char *name) { fBkgJetBranchName = name ; }
101 
102  //switches for photons
106 
107  void CalculateBkg(TVector3 gamma, TVector3 jet,Double_t *vector,Int_t type);
108 
109  void FindMCgenInfo();//gives information on generated level
110 
111  void SwitchOnSaveGJTree() { fSaveGJTree = kTRUE ; }
112  void SwitchOffSaveGJTree() { fSaveGJTree = kFALSE ; }
113  Bool_t IsSaveGJTree() const { return fSaveGJTree ; }
114 
115  void SwitchOnMostEnergetic() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
116  void SwitchOffMostEnergetic() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
117  void SwitchOffMostOpposite() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
118  void SwitchOnMostOpposite() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
119  Bool_t IsMostEnergetic() const { return fMostEnergetic ; }
120  Bool_t IsMostOpposite() const { return fMostOpposite ; }
121 
122  //switches for histograms
126 
130 
134 
135  void SwitchOnMCStudies() { fMCStudies = kTRUE ; }
136  void SwitchOffMCStudies() { fMCStudies = kFALSE ; }
137  Bool_t IsMCStudies() const { return fMCStudies ; }
138 
139 private:
140 
141  // selection parameters
146 
152 
157 //Bool_t fNonStandardJetFromReader; ///< use non standard jet from reader //new
161 
167 
171 
173 
174  TRandom2 * fGenerator;
175 
176  TLorentzVector fMomentum;
177 
178  // Histograms
180 //TH2F * fhDeltaPhi; //!<! Difference of jet phi and trigger particle phi as function of trigger particle pT
183 
186  TH2F * fhPt;
187 
188  TH2F * fhFFz ;
192 
198 
199  TH1F * fhGamPtPerTrig ;
201 
202  // background from RC
203  TH2F * fhBkgFFz[5] ;
209 
210  // temporary histograms
212  TH1F * fhCuts;
213 
220 
221  // temporary jet histograms
222  TH1F * fhJetPtBefore;
224  TH1F * fhJetPt;
225  TH1F * fhJetPtMostEne;
226  TH1F * fhJetPhi;
227  TH1F * fhJetEta;
236  TH1F * fhJetNjetOverPtCut[10];
243 
250 
251  // temporary background jet histograms
252  TH1F * fhBkgJetBackground[4];
253  TH1F * fhBkgJetSigma[4];
254  TH1F * fhBkgJetArea[4];
255 
256  // temporary photon histograms
269  TH1F * fhPhotonPt;
272  TH1F * fhPhotonPtDiff;
277 
281 
282 
283  // temporary jet histograms after selection
290 
291  TH1F * fhCuts2;
292 
293  // temporary photon histogram after selection
297 
298  // MC generated histograms
299  TH1F * fhMCPhotonCuts;
300  TH1F * fhMCPhotonPt;
302  TH1F * fhMCJetOrigin;
315 
316  // tree with data gamma and jet
345 
366 
371 
374 
377 
381 
382  } ;
383 
384 #endif //ALIANAPARTICLEJETFINDERCORRELATION_H
385 
386 
387 
Int_t fMCJetNPart
MC gen number of full jet particles.
TH2F * fhJetNtrackRatioMostEne[5]
! the same for most energetic jet
Int_t fGamNtracks
number of tracks in iso cone
Bool_t fSaveGJTree
flag to save gamma-jet tree
TH2F * fhNjetsNgammas
! Number of jets vs number of photons in the event
Bool_t fMakeCorrelationInHistoMaker
Make particle-jet correlation in histogram maker.
void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
double Double_t
Definition: External.C:58
Bool_t fSelectIsolated
Select only trigger particles isolated.
TH2F * fhMCJetChNPart150VsPt
! generated N parts (pt>150 MeV/c) vs pt charged jet
Double_t fMCJet150Eta
MC gen full jet eta (pt>150MeV/c)
Int_t fMCJet150NPart
MC gen number of full jet particles (pt>150MeV/c)
Definition: External.C:236
TH2F * fhPhotonPtDiffVsNtracks
! correction vs Ntracks
Bool_t fMostEnergetic
flag to choose gamma-jet pairs most energetic
Double_t fJetConeSize
Reconstructed jet cone size.
Double_t fGamRho
background energy for photons per cell in EMCal
void InitParameters()
Initialize the parameters of the analysis.
TH2F * fhMCJetEtaPhi
! generated jet eta vs phi for full jet
TH2F * fhDeltaPhiBefore
! Difference of jet phi and trigger particle phi as function of trigger particle pT ...
TH1F * fhJetPtMostEne
! Pt of the most energetic jet
TH1F * fhBkgJetSigma[4]
! sigma of jet in backgroud branch
TH2F * fhSelectedPhotonNLMVsPt
! nlm vs pt for selected photons
Double_t fMCJetCh150Eta
MC gen charged jet eta (pt>150MeV/c)
TH2F * fhRandomPhiEta[5]
! eta and phi from random generator
TH2F * fhSelectedTrackPhiVsEta
! Phi vs eta of all chosen tracks in selected events
TH1F * fhPhotonAverageEnergyMinus1
! average energy of photon w/o most ene photon
TH2F * fhBkgFFz[5]
! Background fragmentation function, z=ptjet/pttrig
TH2F * fhJetChAreaVsPt
! area of each charged jet vs jet pt
Int_t fJetNtracks2
number of jet tracks with pt>2 GeV/c
TH2F * fhPtRatio
! Ratio of jet pT and trigger particle pT as function of trigger particle pT
void FindMCgenInfo()
Find information about photon and (quark or gluon) on generated level.
Double_t fMCJet150Phi
MC gen full jet phi (pt>150MeV/c)
TH2F * fhSelectedPhotonLambda0VsPt
! lambda0 vs pt for selected photons
TH2F * fhDeltaPhi0PiCorrect
! Difference of jet phi and trigger particle phi as function of trigger particle pT ...
TH1F * fhJetPt
! Pt of all jets after bkg correction
Int_t fMCJetCh150NPart
MC gen number of charged jet particles (pt>150MeV/c)
TH1F * fhPhotonSumPtInCone
! sum pt in cone before correction
TH2F * fhFFz
! Accepted reconstructed jet fragmentation function, z=pt^particle,jet/pttrig
TH2F * fhJetChBkgEnergyVsPt
! background energy of each charged jet vs jet pt
TH1F * fhPhotonSumPtChargedInCone
! sum pt of charged tracks in the cone before correction
Int_t fMCJetCh150ConeNPart
MC gen number of charged jet particles (pt>150MeV/c),R=0.4.
Double_t fGamSumPtNeu
energy in isolation cone neutral
TH2F * fhMCJetChEtaPhi
! generated jet eta vs phi for charged jet
TH1F * fhPhotonPt
! pt of gamma before bkg correction
TH1F * fhJetNjetOverPtCut[10]
! number of reconstructed jets in event over pT threshold
Int_t fMCPartonType
MC gen parton type origin of jet.
TH2F * fhBkgSumPtInCone[5]
! Background sum pt in cone
TH2F * fhPtRatioBefore
! Ratio of jet pT and trigger particle pT as function of trigger particle pT
TH2F * fhMCJetNPartVsPt
! generated N parts vs pt full jet
TH2F * fhMCJetNPart150VsPt
! generated N parts (pt>150 MeV/c) vs pt full jet
TH2F * fhJetNtrackRatioJet5GeV[5]
! the same for pt jet above 5 GeV
TH2F * fhJetRatioNTrkAboveToNTrk[5]
! ratio tracks in jet with pt above 1,2,3,4,5GeV to ntracks
TH2F * fhJetNtrackRatioLead5GeV[5]
! the same for jet with leading particle pt>5GeV
TH1F * fhJetPtBeforeCut
! Pt of all jets after bkg correction, raw jet pt>fJetMinPt
TString fJetBranchName
name of jet branch not set in reader part //new
TH2F * fhDeltaPt
! Difference of jet pT and trigger particle pT as function of trigger particle pT ...
TRandom2 * fGenerator
! pointer to random generator object
Base class for CaloTrackCorr analysis algorithms.
Double_t fMCJet150Pt
MC gen full jet (pt^particles>150MeV/c) pt.
AliAnaParticleJetFinderCorrelation()
Default constructor. Initialize parameters.
TH2F * fhDeltaPhi0PiCorrectBefore
! Difference of jet phi and trigger particle phi (0,pi) as function of trigger particle pT ...
TH1F * fhMCJetRatioChFull
! generated ratio pt charged/full jet
int Int_t
Definition: External.C:63
Bool_t fUseHistogramJetTracks
flag to save jet tracks features
void MakeAnalysisFillHistograms()
Particle-Jet Correlation Analysis, fill histograms.
TH2F * fhPhotonPtDiffVsCentrality
! correction vs centrality
TH2F * fhJetDeltaEtaDeltaPhi
! delta eta vs delta phi for (jet-track) <-0.8,0.8>
float Float_t
Definition: External.C:68
TH1F * fhGamPtPerTrig
! per trigger normalisation
TH2F * fhJetEtaVsNpartInJetBkg
! Eta vs number of particles in jet for background subtracted jets
void CalculateBkg(TVector3 gamma, TVector3 jet, Double_t *vector, Int_t type)
TH2F * fhBkgSumPtnTracksInCone[5]
! Background sum pt over ntracks in cone
TH2F * fhSelectedJetPhiVsEta
! phi vs eta of selected jet
TH2F * fhDeltaEta
! Difference of jet eta and trigger particle eta as function of trigger particle pT ...
TH1F * fhMCJetRatioCh150Ch
! generated ratio pt charged(pt>150MeV/c)/charged jet
Double_t fMCJetCh150Pt
MC gen charged jet (pt^particles>150MeV/c) pt.
Double_t fRatioMaxCut
Jet/particle Ratio cut maximum.
TH1F * fhCuts
! Number of events after cuts
Bool_t fUseHistogramJetBkg
flag to save bkg jet histograms
TH1F * fhPhotonNgammaMoreAverageMinus1ToNgamma
! number of gammas with ene. more than average ene (w/o most ene gamma) divided by no...
TH1F * fhPhotonPtCorrectedZoom
! pt of gamma after background correction in +-5 GeV/c
TH2F * fhJetFFz
! Accepted reconstructed jet fragmentation function, z=pt^particle,jet/ptjet
TH2F * fhFFpt
! Jet particle pt distribution in cone
TH2F * fhJetChBkgEnergyVsArea
! area of each charged jet vs jet background
TH1F * fhSelectedJetNjet
! number of jets in selected event
TH1F * fhBkgJetArea[4]
! area of jet in bkg branch
TH2F * fhMCJet150EtaPhi
! generated jet eta vs phi full jet (pt>150 MeV/c)
Double_t fJetMinPt
Minumum jet pt, default 5GeV/c.
Bool_t fMostOpposite
flag to choose gamma-jet pairs most opposite
Bool_t fBackgroundJetFromReader
use background jet from reader //new
TH1F * fhSelectedNtracks
! number of tracks in selected event
TH2F * fhDeltaPtBefore
! Difference of jet pT and trigger particle pT as function of trigger particle pT ...
Double_t fGamAvEne
average energy of photons (without most ene)
TH2F * fhPhotonBkgRhoVsNclusters
! average energy in one cell vs n clusters
Double_t fMCJetCh150ConePhi
MC gen charged jet phi (pt>150MeV/c),R=0.4.
TH2F * fhMCJetCh150EtaPhi
! generated jet eta vs phi charged jet (pt>150 MeV/c)
Double_t fConeSize
Jet cone size to calculate fragmentation function.
TH2F * fhFFxi
! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/pt^particle,jet) ...
TH2F * fhJetNtracksInJetAboveThr[6]
! number of tracks in jet with pt above 0,1,2,3,4,5GeV
TH2F * fhSelectedJetChBkgEnergyVsPtJet
! background energy of selected charged jet vs jet pt
TH1F * fhTrackAveTrackPt
! average track pt in event
Double_t fMCJetCh150Phi
MC gen charged jet phi (pt>150MeV/c)
Int_t fGamNclusters
number of clusters in iso cone
TH2F * fhJetFFzCor
! Accepted reconstructed jet fragmentation function, z=pt^particle,jet*-cos(jet,trig)/ptjet ...
TH1F * fhPhotonPtDiff
! bkg correction = n_cells * median_rho
void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
TH2F * fhJetDeltaEtaDeltaPhiAllTracks
! delta eta vs delta phi for (jet-track) <-pi,pi>
TH1F * fhBkgJetBackground[4]
! background from jet bkg branch
TH2F * fhNTracksInCone
! jet multiplicity in cone
TH2F * fhPt
! jet pT vs trigger particle pT
TH2F * fhPhotonBkgRhoVsNtracks
! average energy in one cell vs n tracks
TH2F * fhPhotonPtDiffVsNclusters
! correction vs Nclustres
TH2F * fhJetFFxi
! Accepted reconstructed jet fragmentation function, xsi = ln(ptjet/pt^particle,jet) ...
TH2F * fhMCPhotonEtaPhi
! generated direct photon eta vs phi
TH1F * fhMCPhotonPt
! generated direct photon pt
TH2F * fhSelectedJetChAreaVsPtJet
! area of selected charged jet vs jet pt
TH2F * fhPtGamPtJet
! gamma jet correlation filling
Bool_t fUseJetRefTracks
Use track references from JETAN not the AOD tracks to calculate fragmentation function.
Int_t fJetNtracks1
number of jet tracks with pt>1 GeV/c
Bool_t fUseHistogramTracks
flag to save CTS tracks features
Double_t fDeltaPhiMaxCut
Minimum Delta Phi Gamma-Leading.
TH2F * fhJetFFxiCor
! Accepted reconstructed jet fragmentation function, xsi = ln(ptjet/pt^particle*-cos(jet,trig),jet)
TH1F * fhPhotonAverageEnergy
! average energy of photon
TH2F * fhPtBefore
! jet pT vs trigger particle pT
TH2F * fhDeltaPhiCorrect
! Difference of jet phi and trigger particle phi as function of trigger particle pT ...
TH1F * fhJetAveTrackPt
! average track from jets pt in event
TH1F * fhJetNparticlesInJet
! number of particles in jets
TH2F * fhBkgFFxi[5]
! Background fragmentation function, xsi = ln(pttrig/ptjet)
TString fBkgJetBranchName
name of background jet branch not set in reader part //new
TH2F * fhBkgFFpt[5]
! Background particle pt distribution in cone
Double_t fJetMinPtBkgSub
Minumum jet pt after bkg subtraction, default -100 GeV/c.
TH2F * fhBkgNTracksInCone[5]
! Background multiplicity in cone
const char Option_t
Definition: External.C:48
Double_t fMCJetCh150ConeEta
MC gen charged jet eta (pt>150MeV/c),R=0.4.
Double_t fGamSumPtCh
energy in isolation cone charged
TH1F * fhPhotonSumPtCorrectInCone
! sum pt in cone afrer correction
Double_t fRatioMinCut
Jet/particle Ratio cut minimum.
AliAnaParticleJetFinderCorrelation & operator=(const AliAnaParticleJetFinderCorrelation &g)
Assignment operator not implemented.
Bool_t fUseBackgroundSubtractionGamma
flag to use backgrouind subtraction for photons or not
Double_t fPtThresholdInCone
Jet pT threshold in jet cone.
TH1F * fhPhotonNgammaOverPtCut[10]
! number of photons in event over pT threshold
bool Bool_t
Definition: External.C:53
TH2F * fhJetRhoVsCentrality
! jet energy density vs centrality
TH2F * fhPhotonBkgRhoVsCentrality
! average energy in one cell vs centrality
Int_t SelectJet(AliAODPWG4Particle *particle, TClonesArray *aodRecJets)
TH1F * fhPhotonPtCorrected
! pt of gamma after background correction
TH1F * fhPhotonRatioAveEneMinus1ToMostEne
! ratio average energy of photon w/o most ene photon to most energetic photon
TH2F * fhJetFFpt
! Jet particle pt distribution in jet cone
Double_t fDeltaPhiMinCut
Maximum Delta Phi Gamma-Leading.
Int_t fMCJetChNPart
MC gen number of charged jet particles.
TH2F * fhTrackPhiVsEta
! Phi vs eta of all chosen tracks in all events
TH2F * fhJetRhoVsPt
! jet energy density vs jet pt
TH2F * fhDeltaEtaBefore
! Difference of jet eta and trigger particle eta as function of trigger particle pT ...
TH1F * fhPhotonNgammaMoreAverageToNgamma
! number of gammas with ene. more than average ene divided by no. of gammas
TH2F * fhMCJetChNPartVsPt
! generated N parts vs pt charged jet
TH1F * fhPhotonRatioAveEneToMostEne
! ratio average energy to most energetic photon
TH2F * fhJetEtaVsNpartInJet
! Eta vs number of particles in jet for all jets
Double_t fJetAreaFraction
Jet area fraction X in X*pi*R^2, default 0.6.
TH2F * fhPhotonBkgRhoVsNcells
! average energy in one cell vs n cells
Double_t fMCJetCh150ConePt
MC gen charged jet (pt^particles>150MeV/c),R=0.4 pt.
TH2F * fhMCJetCh150ConeEtaPhi
! generated jet eta vs phi charged jet (pt>150 MeV/c) R=0.4
void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
const Double_t phimin
void MakeAnalysisFillAOD()
Particle-Jet Correlation Analysis, fill AODs.
TH2F * fhMCJetChNPart150ConeVsPt
! generated N parts (pt>150 MeV/c) vs pt charged jet R=0.4