AliPhysics  master (3d17d9d)
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 ; }
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 SelectIsolated(Bool_t select) { fSelectIsolated = select ; }
74  void SetIsolationMatters(Bool_t select) { fIsolationMatters = select ; }
75 
76  void SetConeSize(Float_t cone) { fConeSize = cone ; }
78 
80  { fDeltaPhiMaxCut = phimax ; fDeltaPhiMinCut = phimin ; }
81 
82  void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
83  { fRatioMaxCut = ratiomax; fRatioMinCut = ratiomin ; }
84 
85  void UseJetRefTracks(Bool_t use) { fUseJetRefTracks = use ; }
87  void SetJetConeSize(Double_t cone) { fJetConeSize = cone ; }
88  void SetJetMinPt(Double_t minpt) { fJetMinPt = minpt ; }
89  void SetJetMinPtBkgSub(Double_t minpt) { fJetMinPtBkgSub = minpt ; }
90  void SetJetAreaFraction(Double_t areafr) { fJetAreaFraction = areafr ; }
91  void SetGammaConeSize(Float_t cone) { fGammaConeSize = cone ; }
92 
93  // Settings for non standard jet branch
94 //TString GetJetBranchName() const { return fJetBranchName ; }
95 //void SetJetBranchName(const char *name) { fJetBranchName = name ; }
96 //void SwitchOnNonStandardJetFromReader() { fNonStandardJetFromReader = kTRUE ; }
97 //void SwitchOffNonStandardJetFromReader() { fNonStandardJetFromReader = kFALSE ; }
98 //Bool_t IsNonStandardJetFromReader() { return fNonStandardJetFromReader ; }
99 
100 // TString GetBkgJetBranchName() const { return fBkgJetBranchName ; }
101 // void SetBkgJetBranchName(const char *name) { fBkgJetBranchName = name ; }
105 
106  //switches for photons
110 
111  void CalculateBkg(TVector3 gamma, TVector3 jet,Double_t *vector,Int_t type);
112 
113  void FindMCgenInfo();//gives information on generated level
114 
115  void SwitchOnSaveGJTree() { fSaveGJTree = kTRUE ; }
116  void SwitchOffSaveGJTree() { fSaveGJTree = kFALSE ; }
117  Bool_t IsSaveGJTree() const { return fSaveGJTree ; }
118 
119  void SwitchOnMostEnergetic() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
120  void SwitchOffMostEnergetic() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
121  void SwitchOffMostOpposite() { fMostEnergetic = kTRUE ; fMostOpposite = kFALSE ; }
122  void SwitchOnMostOpposite() { fMostEnergetic = kFALSE ; fMostOpposite = kTRUE ; }
123  Bool_t IsMostEnergetic() const { return fMostEnergetic ; }
124  Bool_t IsMostOpposite() const { return fMostOpposite ; }
125 
126  //switches for histograms
130 
134 
138 
139  void SwitchOnMCStudies() { fMCStudies = kTRUE ; }
140  void SwitchOffMCStudies() { fMCStudies = kFALSE ; }
141  Bool_t IsMCStudies() const { return fMCStudies ; }
142 
143 private:
144 
145  // selection parameters
150 
157 
162 //Bool_t fNonStandardJetFromReader; ///< use non standard jet from reader //new
163 // TString fJetBranchName ; ///< name of jet branch not set in reader part //new
165  // TString fBkgJetBranchName ; ///< name of background jet branch not set in reader part //new
166 
172 
176 
178 
179  TRandom2 * fGenerator;
180 
181  TLorentzVector fMomentum;
182 
183  // Histograms
185 //TH2F * fhDeltaPhi; //!<! Difference of jet phi and trigger particle phi as function of trigger particle pT
188 
191  TH2F * fhPt;
192 
193  TH2F * fhFFz ;
197 
203 
204  TH1F * fhGamPtPerTrig ;
206 
207  // background from RC
208  TH2F * fhBkgFFz[5] ;
214 
215  // temporary histograms
217  TH1F * fhCuts;
218 
225 
226  // jet histograms
227  TH1F * fhJetPtBefore;
229  TH1F * fhJetPt;
230  TH1F * fhJetPtMostEne;
231  TH1F * fhJetPhi;
232  TH1F * fhJetEta;
241  TH1F * fhNoCTSTracks;
243  TH1F * fhJetNjetOverPtCut[10];
250 
257 
258  // background jet histograms
259  TH1F * fhBkgJetBackground[3];
260  TH1F * fhBkgJetSigma[3];
261  TH1F * fhBkgJetArea[3];
262 
263  // temporary photon histograms
276  TH1F * fhPhotonPt;
279  TH1F * fhPhotonPtDiff;
284 
288 
294 
295 
296  // temporary jet histograms after selection
303 
304  TH1F * fhCuts2;
305 
306  // temporary photon histogram after selection
310 
311  // MC generated histograms
312  TH1F * fhMCPhotonCuts;
313  TH1F * fhMCPhotonPt;
315  TH1F * fhMCJetOrigin;
328 
329  // tree with data gamma and jet
358 
379 
384 
387 
390 
394 
395  } ;
396 
397 #endif //ALIANAPARTICLEJETFINDERCORRELATION_H
398 
399 
400 
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.
TH1F * fhPhotonRatioSumPtChargedInCone
! ratio (sum pt/gamma) of charged tracks in the cone before correction
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
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 neutral 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
TH1F * fhNoCTSTracksCut
! number of CTS tracks with pT>cut
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
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
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>
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
TH1F * fhPhotonAvePtChargedInCone
! average pt of charged tracks in the cone before correction
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
TH1F * fhBkgJetArea[3]
! area of jet in bkg branch
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)
TH1F * fhPhotonNChargedInCone
! number of charged tracks in the cone before correction
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 * fhPhotonRatioMostPtChargedInCone
! ratio (most pt/gamma) of charged tracks in the cone before correction
Bool_t fIsolationMatters
flag if isolation matters
TH1F * fhPhotonSumPtCorrectInCone
! sum pt neutral 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
TH1F * fhPhotonRatioAvePtChargedInCone
! ratio (average pt/gamma) of charged tracks in the cone before correction
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
TH1F * fhBkgJetSigma[3]
! sigma of jet in backgroud branch
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
TH1F * fhBkgJetBackground[3]
! background from jet bkg branch
void MakeAnalysisFillAOD()
Particle-Jet Correlation Analysis, fill AODs.
TH2F * fhMCJetChNPart150ConeVsPt
! generated N parts (pt>150 MeV/c) vs pt charged jet R=0.4