AliPhysics  3337bb0 (3337bb0)
 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 //_________________________________________________________________________
22 
23 // --- ROOT system ---
24 class TH2F;
25 class TTree;
26 class TRandom2;
27 
28 //---- Analysis system ----
30 
32 
33  public:
34 
36 
38 
39  // General methods
40 
41  void InitParameters();
42 
44 
45  void MakeAnalysisFillAOD() ;
46 
48 
49  // To access non standard branch
50  Int_t SelectJet(AliCaloTrackParticle * particle, TClonesArray * aodRecJets) ;
51 
52  void Print(const Option_t * opt) const;
53 
54  // Settings
55 
56  Bool_t OnlyIsolated() const { return fSelectIsolated ; }
57  void SelectIsolated(Bool_t select) { fSelectIsolated = select ; }
58 
59  Float_t GetConeSize() const { return fConeSize ; }
63  Double_t GetRatioMaxCut() const { return fRatioMaxCut ; }
64  Double_t GetRatioMinCut() const { return fRatioMinCut ; }
67  Double_t GetJetConeSize() const { return fJetConeSize ; }
68  Double_t GetJetMinPt() const { return fJetMinPt ; }
72 
73  void SetConeSize(Float_t cone) { fConeSize = cone ; }
75 
77  { fDeltaPhiMaxCut = phimax ; fDeltaPhiMinCut = phimin ; }
78 
79  void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
80  { fRatioMaxCut = ratiomax; fRatioMinCut = ratiomin ; }
81 
82  void UseJetRefTracks(Bool_t use) { fUseJetRefTracks = use ; }
84  void SetJetConeSize(Double_t cone) { fJetConeSize = cone ; }
85  void SetJetMinPt(Double_t minpt) { fJetMinPt = minpt ; }
86  void SetJetMinPtBkgSub(Double_t minpt) { fJetMinPtBkgSub = minpt ; }
87  void SetJetAreaFraction(Double_t areafr) { fJetAreaFraction = areafr ; }
88  void SetGammaConeSize(Float_t cone) { fGammaConeSize = cone ; }
89 
90  // Settings for non standard jet branch
92  void SetJetBranchName(const char *name) { fJetBranchName = name ; }
93 //void SwitchOnNonStandardJetFromReader() { fNonStandardJetFromReader = kTRUE ; }
94 //void SwitchOffNonStandardJetFromReader() { fNonStandardJetFromReader = kFALSE ; }
95 //Bool_t IsNonStandardJetFromReader() { return fNonStandardJetFromReader ; }
96 
98  void SetBkgJetBranchName(const char *name) { fBkgJetBranchName = name ; }
102 
103  //switches for photons
107 
108  void CalculateBkg(TVector3 gamma, TVector3 jet,Double_t *vector,Int_t type);
109 
110  void FindMCgenInfo();//gives information on generated level
111 
112  void SwitchOnSaveGJTree() { fSaveGJTree = kTRUE ; }
113  void SwitchOffSaveGJTree() { fSaveGJTree = kFALSE ; }
114  Bool_t IsSaveGJTree() const { return fSaveGJTree ; }
115 
116  void SwitchOnMostEnergetic() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
117  void SwitchOffMostEnergetic() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
118  void SwitchOffMostOpposite() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
119  void SwitchOnMostOpposite() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
120  Bool_t IsMostEnergetic() const { return fMostEnergetic ; }
121  Bool_t IsMostOpposite() const { return fMostOpposite ; }
122 
123  //switches for histograms
127 
131 
135 
136  void SwitchOnMCStudies() { fMCStudies = kTRUE ; }
137  void SwitchOffMCStudies() { fMCStudies = kFALSE ; }
138  Bool_t IsMCStudies() const { return fMCStudies ; }
139 
140 private:
141 
142  // selection parameters
147 
153 
158 //Bool_t fNonStandardJetFromReader; ///< use non standard jet from reader //new
162 
168 
172 
174 
175  TRandom2 * fGenerator;
176 
177  TLorentzVector fMomentum;
178 
179  // Histograms
181 //TH2F * fhDeltaPhi; //!<! Difference of jet phi and trigger particle phi as function of trigger particle pT
184 
187  TH2F * fhPt;
188 
189  TH2F * fhFFz ;
193 
199 
200  TH1F * fhGamPtPerTrig ;
202 
203  // background from RC
204  TH2F * fhBkgFFz[5] ;
210 
211  // temporary histograms
213  TH1F * fhCuts;
214 
221 
222  // temporary jet histograms
223  TH1F * fhJetPtBefore;
225  TH1F * fhJetPt;
226  TH1F * fhJetPtMostEne;
227  TH1F * fhJetPhi;
228  TH1F * fhJetEta;
237  TH1F * fhJetNjetOverPtCut[10];
244 
251 
252  // temporary background jet histograms
253  TH1F * fhBkgJetBackground[4];
254  TH1F * fhBkgJetSigma[4];
255  TH1F * fhBkgJetArea[4];
256 
257  // temporary photon histograms
270  TH1F * fhPhotonPt;
273  TH1F * fhPhotonPtDiff;
278 
282 
283 
284  // temporary jet histograms after selection
291 
292  TH1F * fhCuts2;
293 
294  // temporary photon histogram after selection
298 
299  // MC generated histograms
300  TH1F * fhMCPhotonCuts;
301  TH1F * fhMCPhotonPt;
303  TH1F * fhMCJetOrigin;
316 
317  // tree with data gamma and jet
346 
367 
372 
375 
378 
382 
383  } ;
384 
385 #endif //ALIANAPARTICLEJETFINDERCORRELATION_H
386 
387 
388 
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.
Container for input particle information on CaloTrackCorr package.
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
Int_t SelectJet(AliCaloTrackParticle *particle, TClonesArray *aodRecJets)
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
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