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