1 #ifndef AliFlowAnalysisWithCumulants_H
2 #define AliFlowAnalysisWithCumulants_H
4 //*******************************
5 // flow analysis with cumulants *
6 // author: Ante Bilandzic *
7 // email: anteb@nikhef.nl *
8 //*******************************
10 #include "AliFlowCommonConstants.h"
11 #include "AliFlowCumuConstants.h"
27 class AliFlowEventSimple;
28 class AliFlowTrackSimple;
29 class AliFlowCommonHist;
30 class AliFlowCommonHistResults;
33 class AliFlowAnalysisWithCumulants {
35 AliFlowAnalysisWithCumulants();
36 virtual ~AliFlowAnalysisWithCumulants();
38 virtual void CreateOutputObjects();
39 virtual void Make(AliFlowEventSimple* anEvent);
40 virtual void Finish();
45 TList* GetHistList() const { return this->fHistList ; } // Gets output histogram list //NEW
48 AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants& aAnalysis);
49 AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants& aAnalysis);
51 AliFlowTrackSimple* fTrack;//track
52 static const Int_t fgkQmax=AliFlowCumuConstants::kQmax;//needed for numerics
53 static const Int_t fgkPmax=AliFlowCumuConstants::kPmax;//needed for numerics
54 static const Int_t fgkFlow=AliFlowCumuConstants::kFlow;//integrated flow coefficient to be calculated
55 static const Int_t fgkMltpl=AliFlowCumuConstants::kMltpl;//the multiple in p=m*n (diff. flow)
56 static const Int_t fgknBins=100;//number of pt bins
58 //TString fHistFileName; //! The output file name // NEW: In the old version this was not commented out
59 //TFile* fHistFile; // histogram file for Cumulants // NEW: In the old version this was not commented out
61 TList* fHistList; //list to hold all output histograms //NEW
63 Double_t fAvM;//avarage SELECTED multiplicity
65 Double_t fR0;//needed for numerics
66 Double_t fPtMax;//maximum pt
67 Double_t fPtMin;//minimum pt
68 Double_t fBinWidth;//width of pt bin (in GeV)
70 Double_t fAvQx;//<Q_x>
71 Double_t fAvQy;//<Q_y>
72 Double_t fAvQ2x;//<(Q_x)^2>
73 Double_t fAvQ2y;//<(Q_y)^2>
75 Int_t fNumberOfEvents;//number of events
78 TProfile* fHistProAvM; //avarage multiplicity
80 TH1D* fIntFlowResults; //integrated flow final results
82 TH1D* fDiffFlowResults2; //differential flow final results //to be improved
83 TH1D* fDiffFlowResults4; //differential flow final results //to be improved
84 TH1D* fDiffFlowResults6; //differential flow final results //to be improved
85 TH1D* fDiffFlowResults8; //differential flow final results //to be improved
87 TProfile2D* fIntFlowGenFun; // avarage of the generating function for integrated flow
88 TProfile3D* fDiffFlowGenFunRe; // avarage of the generating function for differential flow (real part)
89 TProfile3D* fDiffFlowGenFunIm; // avarage of the generating function for differential flow (imaginary part)
92 AliFlowCommonHist* fCommonHists;//control histograms
93 //AliFlowCommonHistResults *fCommonHistsRes2, *fCommonHistsRes4, *fCommonHistsRes6, *fCommonHistsRes8;//final results//to be improved
95 Double_t fAvG[fgkPmax][fgkQmax];//avarage of the generating function used for integrated flow
96 Int_t fBinEventEntries[fgknBins];//counts how many events have at least 1 particle in particular bin
97 Int_t fBinNoOfParticles[fgknBins];//number of particles per bin
98 Double_t fBinMeanPt[fgknBins];//mean pt per bin
99 Double_t fBinEventDRe[fgknBins][fgkPmax][fgkQmax];//real part of the generating function used for differential flow
100 Double_t fBinEventDIm[fgknBins][fgkPmax][fgkQmax];//imaginary part of the generating function used for differential flow
102 ClassDef(AliFlowAnalysisWithCumulants, 0);