1 #ifndef ALIANAGAMMAHADRON_H
2 #define ALIANAGAMMAHADRON_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 /* History of cvs commits:
10 * Revision 1.2 2007/02/09 18:40:40 schutz
11 * B
\bNew version from Gustavo
13 * Revision 1.1 2007/01/23 17:17:29 schutz
19 //_________________________________________________________________________
20 // Class for the analysis of gamma-jet correlations.
21 // Basically it seaches for a prompt photon in the Calorimeters acceptance,
22 // if so we construct a jet around the highest pt particle in the opposite
23 // side in azimuth. This jet has to fullfill several conditions to be
24 // accepted. Then the fragmentation function of this jet is constructed
25 // Class created from old AliPHOSGammaPion
27 //*-- Author: Gustavo Conesa (INFN-LNF)
29 // --- ROOT system ---
37 #include "AliAnaGammaDirect.h"
41 class AliAnaGammaHadron : public AliAnaGammaDirect {
45 AliAnaGammaHadron(const char *name) ; // default ctor
46 AliAnaGammaHadron(const AliAnaGammaHadron & gj) ; // cpy ctor
47 AliAnaGammaHadron & operator = (const AliAnaGammaHadron & g) ;//cpy assignment
48 virtual ~AliAnaGammaHadron() ; //virtual dtor
50 virtual void Exec(Option_t * opt = "") ;
51 virtual void ConnectInputData(Option_t *);
52 virtual void CreateOutputObjects();
53 virtual void Terminate(Option_t * opt = "");
55 void InitParameters();
56 Double_t GetAngleMaxParam(Int_t i) const {return fAngleMaxParam.At(i) ; }
57 Double_t GetInvMassMaxCut() const {return fInvMassMaxCut ; }
58 Double_t GetInvMassMinCut() const {return fInvMassMinCut ; }
59 Double_t GetPhiMaxCut() const {return fPhiMaxCut ; }
60 Double_t GetPhiMinCut() const {return fPhiMinCut ; }
61 Float_t GetMinPtPion() const {return fMinPtPion ; }
63 void Print(const Option_t * opt)const;
65 void SetAngleMaxParam(Int_t i, Double_t par)
66 {fAngleMaxParam.AddAt(par,i) ; }
67 void SetMinPtPion(Float_t pt){fMinPtPion = pt; };
68 void SetInvMassCutRange(Double_t invmassmin, Double_t invmassmax)
69 {fInvMassMaxCut =invmassmax; fInvMassMinCut =invmassmin;}
70 void SetPhiCutRange(Double_t phimin, Double_t phimax)
71 {fPhiMaxCut =phimax; fPhiMinCut =phimin;}
75 Bool_t IsAngleInWindow(const Float_t angle, const Float_t e);
76 void MakeGammaChargedCorrelation(TClonesArray * pl, TParticle *pGamma) const ;
77 void MakeGammaNeutralCorrelation(TClonesArray * pl, TParticle *pGamma) ;
81 Double_t fPhiMaxCut ; //
82 Double_t fPhiMinCut ; //
83 Double_t fInvMassMaxCut ; // Invariant Mass cut maximum
84 Double_t fInvMassMinCut ; // Invariant Masscut minimun
86 Double_t fMinPtPion; // Minimum pt of pion
87 TObjArray *fOutputContainer ; //! output data container
88 TArrayD fAngleMaxParam ; //Max opening angle selection parameters
96 TH2F * fhDeltaPhiGammaCharged ;
97 TH2F * fhDeltaPhiGammaNeutral ;
98 TH2F * fhDeltaEtaGammaCharged ;
99 TH2F * fhDeltaEtaGammaNeutral ;
101 TH2F * fhCorrelationGammaNeutral ;
102 TH2F * fhCorrelationGammaCharged ;
104 TH2F * fhAnglePairAccepted ;
105 TH2F * fhAnglePairNoCut ;
106 TH2F * fhAnglePairAzimuthCut ;
107 TH2F * fhAnglePairOpeningAngleCut ;
108 TH2F * fhAnglePairAllCut ;
109 TH2F * fhInvMassPairNoCut ;
110 TH2F * fhInvMassPairAzimuthCut ;
111 TH2F * fhInvMassPairOpeningAngleCut ;
112 TH2F * fhInvMassPairAllCut ;
114 ClassDef(AliAnaGammaHadron,0)
118 #endif //ALIANAGAMMAHADRON_H