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 /* $Id: AliAnaParticleJetFinderCorrelation.h 21839 2007-10-29 13:49:42Z gustavo $ */
8 //_________________________________________________________________________
9 // Class that contains the algorithm for the analysis of particle (direct gamma) - jet
10 // (standard jet found with JETAN) correlation
11 // Particle and jet for correlation found by independent algorithms.
12 // For Example direct isolated photon found in AliAnaGammaDirect and the jet with JETAN
14 //-- Author: Gustavo Conesa (INFN-LNF)
16 // --- ROOT system ---
19 //---- Analysis system ----
20 #include "AliAnaPartCorrBaseClass.h"
22 class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
25 AliAnaParticleJetFinderCorrelation() ; // default ctor
26 virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
28 AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
29 AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment
33 TList * GetCreateOutputObjects();
35 Bool_t OnlyIsolated() const {return fSelectIsolated ; }
36 void SelectIsolated(Bool_t select) {fSelectIsolated = select ; }
38 void InitParameters();
40 void Print(const Option_t * opt) const;
42 Float_t GetConeSize() const { return fConeSize ; }
43 Float_t GetPtThresholdInCone() const { return fPtThresholdInCone ; }
44 Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; }
45 Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; }
46 Double_t GetRatioMaxCut() const { return fRatioMaxCut ; }
47 Double_t GetRatioMinCut() const { return fRatioMinCut ; }
48 Bool_t AreJetRefTracks() const { return fUseJetRefTracks ; }
49 Bool_t IsCorrelationMadeInHistoMaker() const {return fMakeCorrelationInHistoMaker ;}
51 void SetConeSize(Float_t cone) {fConeSize = cone; }
52 void SetPtThresholdInCone(Float_t pt){fPtThresholdInCone = pt; };
53 void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
54 {fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin;}
55 void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
56 {fRatioMaxCut =ratiomax; fRatioMinCut = ratiomin ; }
57 void UseJetRefTracks(Bool_t use) {fUseJetRefTracks = use ; }
58 void SetMakeCorrelationInHistoMaker(Bool_t make) {fMakeCorrelationInHistoMaker = make ;}
62 Int_t SelectJet(AliAODPWG4Particle * particle, AliAODEvent * event) const ;
64 void MakeAnalysisFillAOD() ;
65 void MakeAnalysisFillHistograms() ;
69 //selection parameters
70 Double_t fDeltaPhiMaxCut ; //! Minimum Delta Phi Gamma-Leading
71 Double_t fDeltaPhiMinCut ; //! Maximum Delta Phi Gamma-Leading
72 Double_t fRatioMaxCut ; //! Jet/ particle Ratio cut maximum
73 Double_t fRatioMinCut ; //! Jet/particle Ratio cut minimum
75 Double_t fConeSize ; //! Jet cone size
76 Double_t fPtThresholdInCone ; //! Jet pT threshold in jet cone
77 Bool_t fUseJetRefTracks ; //! Use track references from JETAN not the AOD tracks
78 Bool_t fMakeCorrelationInHistoMaker ; //!Make particle-jet correlation in histogram maker
79 Bool_t fSelectIsolated ; // Select only trigger particles isolated
81 TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT
82 TH2F * fhDeltaPhi; //! Difference of jet phi and trigger particle phi as function of trigger particle pT
83 TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT
84 TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT
85 TH2F * fhPt; //! jet pT vs trigger particle pT
87 TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
88 TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
89 TH2F * fhFFpt; //! Jet particle pt distribution in cone
90 TH2F * fhNTracksInCone ; //! jet multiplicity in cone
92 ClassDef(AliAnaParticleJetFinderCorrelation,2)
96 #endif //ALIANAPARTICLEJETFINDERCORRELATION_H