]>
Commit | Line | Data |
---|---|---|
cf6522d1 | 1 | //_________________________________________________________________________ |
2 | //Created by Christine Nattrass, Rebecca Scott, Irakli Martashvili | |
641e1e0c | 3 | //University of Tennessee at Knoxville |
cf6522d1 | 4 | // |
5 | // This class is designed for the analysis of the hadronic component of | |
6 | // transverse energy. It is used by AliAnalysisTaskHadEt. | |
7 | //_________________________________________________________________________ | |
641e1e0c | 8 | #ifndef ALIANALYSISHADET_H |
9 | #define ALIANALYSISHADET_H | |
10 | ||
11 | #include "TString.h" | |
7d2d1773 | 12 | #include "AliAnalysisEtCommon.h" |
641e1e0c | 13 | |
14 | class TH2F; | |
15 | class TH1F; | |
16 | class AliVEvent; | |
17 | class TList; | |
020abff0 | 18 | class AliESDtrackCuts; |
19 | class Rtypes; | |
20 | class TParticle; | |
21 | class TDatabasePDG; | |
4998becf | 22 | class AliAnalysisEtCuts; |
9a365626 | 23 | class AliPIDResponse; |
c692b181 | 24 | class AliCentrality; |
641e1e0c | 25 | |
7d2d1773 | 26 | class AliAnalysisHadEt : public AliAnalysisEtCommon |
641e1e0c | 27 | { |
28 | public: | |
29 | ||
30 | AliAnalysisHadEt(); | |
31 | virtual ~AliAnalysisHadEt(); | |
32 | ||
33 | /** Analyse the event! */ | |
34 | virtual Int_t AnalyseEvent(AliVEvent *event); | |
35 | ||
641e1e0c | 36 | |
37 | /** Initialise the analysis, must be overloaded. */ | |
38 | virtual void Init(); | |
39 | ||
641e1e0c | 40 | /** Reset event specific values (Et etc.) */ |
41 | virtual void ResetEventValues(); | |
4998becf | 42 | |
4998becf | 43 | |
44 | /** Cuts info */ | |
45 | AliAnalysisEtCuts * GetCuts() const { return fCuts; } | |
46 | virtual void SetCuts(const AliAnalysisEtCuts *cuts) | |
47 | { fCuts = (AliAnalysisEtCuts *) cuts; } | |
641e1e0c | 48 | |
49 | /** Sum of the total Et for all events */ | |
020abff0 | 50 | Double_t GetSumEt() const { return fSumEt; } |
641e1e0c | 51 | |
52 | /** Sum of the total Et within our acceptance for all events */ | |
020abff0 | 53 | Double_t GetSumEtAcc() const { return fSumEtAcc; } |
641e1e0c | 54 | |
55 | /** Total Et in the event (without acceptance cuts) */ | |
020abff0 | 56 | Double_t GetTotEt() const { return fTotEt; } |
641e1e0c | 57 | |
58 | /** Total Et in the event within the acceptance cuts */ | |
020abff0 | 59 | Double_t GetTotEtAcc() const { return fTotEtAcc; } |
641e1e0c | 60 | |
61 | /** Total neutral Et in the event (without acceptance cuts) */ | |
020abff0 | 62 | Double_t GetTotNeutralEt() const { return fTotNeutralEt; } |
641e1e0c | 63 | |
64 | /** Total neutral Et in the event within the acceptance cuts */ | |
020abff0 | 65 | Double_t GetTotNeutralEtAcc() const { return fTotNeutralEtAcc; } |
641e1e0c | 66 | |
67 | /** Total charged Et in the event (without acceptance cuts) */ | |
020abff0 | 68 | Double_t GetTotChargedEt() const { return fTotChargedEt; } |
641e1e0c | 69 | |
70 | /** Total charged Et in the event within the acceptance cuts */ | |
020abff0 | 71 | Double_t GetTotChargedEtAcc() const { return fTotChargedEtAcc; } |
641e1e0c | 72 | |
73 | ||
cf6522d1 | 74 | void SetHistoList(const TList *mylist){fhistoList = (TList *) mylist;} |
641e1e0c | 75 | |
9a365626 | 76 | AliPIDResponse *fPIDResponse; |
c692b181 | 77 | Int_t GetCentralityBin(Int_t numberofbins,AliCentrality *centrality);//number of bins is 11 (10%), 21 (5%), or 41 (2.5%). Number of bins includes a bin for those events with no centrality bin determined |
641e1e0c | 78 | |
cf6522d1 | 79 | protected: |
020abff0 | 80 | |
020abff0 | 81 | Double_t fSumEt;/** Sum of the total Et for all events */ |
020abff0 | 82 | Double_t fSumEtAcc;/** Sum of the total Et within our acceptance for all events */ |
020abff0 | 83 | Double_t fTotEt;/** Total Et in the event (without acceptance cuts) */ |
020abff0 | 84 | Double_t fTotEtAcc;/** Total Et in the event within the acceptance cuts */ |
85 | ||
86 | Double_t fTotNeutralEt;/** Total neutral Et in the event */ | |
020abff0 | 87 | Double_t fTotNeutralEtAcc;/** Total neutral Et in the event within the acceptance cuts */ |
020abff0 | 88 | Double_t fTotChargedEt;/** Total charged Et in the event */ |
020abff0 | 89 | Double_t fTotChargedEtAcc;/** Total charged Et in the event within the acceptance cuts */ |
90 | ||
91 | Int_t fMultiplicity;/** Multiplicity of particles in the event */ | |
020abff0 | 92 | Int_t fChargedMultiplicity;/** Multiplicity of charged particles in the event */ |
020abff0 | 93 | Int_t fNeutralMultiplicity; /** Multiplicity of neutral particles in the event */ |
94 | ||
641e1e0c | 95 | void CreateEtaPtHisto2D(TString name, TString title); |
66e4ff36 | 96 | void CreateResolutionPtHisto2D(TString name, TString title, TString xtitle, TString ytitle); |
97 | void CreatePtHisto1D(TString name, TString title, TString xtitle, TString ytitle); | |
641e1e0c | 98 | void CreateEtaHisto1D(TString name, TString title); |
99 | void CreateHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh,Int_t ybins,Float_t ylow,Float_t yhigh); | |
100 | void CreateHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Float_t xlow,Float_t xhigh); | |
101 | void CreateIntHisto1D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh); | |
102 | void CreateIntHisto2D(TString name, TString title, TString xtitle, TString ytitle,Int_t xbins, Int_t xlow,Int_t xhigh,Int_t ybins,Int_t ylow,Int_t yhigh); | |
103 | void FillHisto1D(TString histname, Float_t x, Float_t weight); | |
104 | void FillHisto2D(TString histname, Float_t x, Float_t y, Float_t weight); | |
66e4ff36 | 105 | Bool_t GoodEvent() const {return fGoodEvent;} |
d263b18f | 106 | Float_t TrueP(float pTrec) const; |
641e1e0c | 107 | |
108 | Float_t Et(TParticle *part, float mass = -1000); | |
464aa50c | 109 | Float_t Et(Float_t p, Float_t theta, Int_t pid, Short_t charge) const; |
020abff0 | 110 | |
4575e9f5 | 111 | TList *fhistoList;//->list of histograms saved out to file |
4b40b2b1 | 112 | //static Float_t fgEtaAxis[47];//bins for eta axis of histograms |
113 | static Float_t fgEtaAxis[17];//bins for eta axis of histograms | |
020abff0 | 114 | static Int_t fgnumOfEtaBins;//number of eta bins |
115 | static Float_t fgPtAxis[117];//bins for pt axis of histograms | |
116 | static Int_t fgNumOfPtBins;//number of pt bins | |
66e4ff36 | 117 | static Float_t fgResAxis[81];//axis for resolution histograms |
118 | static Int_t fgNumOfResBins;//number of bins for resolution axis | |
641e1e0c | 119 | |
3ce6b879 | 120 | |
66e4ff36 | 121 | Bool_t fGoodEvent;//boolean to keep track of whether or not this is a good event. |
3ce6b879 | 122 | |
641e1e0c | 123 | private: |
4998becf | 124 | //Declare it private to avoid compilation warning |
125 | AliAnalysisHadEt & operator = (const AliAnalysisHadEt & g) ;//cpy assignment | |
126 | AliAnalysisHadEt(const AliAnalysisHadEt & g) ; // cpy ctor | |
641e1e0c | 127 | |
e573e46d | 128 | ClassDef(AliAnalysisHadEt, 1); |
641e1e0c | 129 | }; |
130 | ||
131 | #endif // ALIANALYSISHADET_H |