PWG2/FLOW/AliFlowAnalysisWithCumulants.h

   1 /*
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
   3  * See cxx source for full Copyright notice
   4  * $Id$
   5  */
   6
   7 /********************************
   8  * flow analysis with cumulants *
   9  *                              *
  10  * author: Ante Bilandzic       *
  11  *          (anteb@nikhef.nl)   *
  12  *******************************/
  13
  14 #ifndef AliFlowAnalysisWithCumulants_H
  15 #define AliFlowAnalysisWithCumulants_H
  16
  17 #include "AliFlowCommonConstants.h"
  18 #include "AliFlowCumuConstants.h"
  19
  20 class TObjArray;
  21 class TList;
  22 class TFile;
  23
  24 class TH1;
  25 class TProfile;
  26 class TProfile2D;
  27 class TProfile3D;
  28
  29 class AliFlowEventSimple;
  30 class AliFlowTrackSimple;
  31 class AliFlowCommonHist;
  32 class AliFlowCommonHistResults;
  33 class AliFlowVector;
  34
  35 //================================================================================================================
  36
  37 class AliFlowAnalysisWithCumulants{
  38  public:
  39   AliFlowAnalysisWithCumulants();
  40   virtual ~AliFlowAnalysisWithCumulants();
  41
  42   virtual void CreateOutputObjects();
  43   virtual void Make(AliFlowEventSimple* anEvent);
  44   virtual void Finish();
  45
  46   TList* GetHistList() const {return this->fHistList;}      //output histogram list
  47
  48  private:
  49   AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants& afawc);
  50   AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants& afawc);
  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
  57   TList* fHistList;                                         //list to hold all output histograms
  58
  59   Double_t fR0;       //needed for numerics
  60   Double_t fPtMax;    //maximum pt
  61   Double_t fPtMin;    //minimum pt
  62   Double_t fBinWidth; //width of pt bin (in GeV)
  63
  64   Double_t fAvQx;  //<Q_x>
  65   Double_t fAvQy;  //<Q_y>
  66   Double_t fAvQ2x; //<(Q_x)^2>
  67   Double_t fAvQ2y; //<(Q_y)^2>
  68
  69   TProfile*          fAvMultIntFlow;     //avarage selected multiplicity
  70
  71   TProfile*          fQVectorComponents; //avarages of Q-vector components (1st bin: <Q_x>, 2nd bin: <Q_y>, 3rd bin: <(Q_x)^2>, 4th bin: <(Q_y)^2>)
  72
  73   TH1D*              fIntFlowResults;   //integrated flow final results
  74
  75   TH1D*              fDiffFlowResults2; //differential flow final results (2nd order estimate)
  76   TH1D*              fDiffFlowResults4; //differential flow final results (4th order estimate)
  77   TH1D*              fDiffFlowResults6; //differential flow final results (6th order estimate)
  78   TH1D*              fDiffFlowResults8; //differential flow final results (8th order estimate)
  79
  80   TProfile2D*        fIntFlowGenFun;    //avarage of the generating function for integrated flow
  81   TProfile3D*        fDiffFlowGenFunRe; //avarage of the generating function for differential flow (real part)
  82   TProfile3D*        fDiffFlowGenFunIm; //avarage of the generating function for differential flow (imaginary part)
  83
  84   TProfile2D *fDiffFlowGenFunRe0,*fDiffFlowGenFunRe1,*fDiffFlowGenFunRe2,*fDiffFlowGenFunRe3;//differential flow
  85   TProfile2D *fDiffFlowGenFunRe4,*fDiffFlowGenFunRe5,*fDiffFlowGenFunRe6,*fDiffFlowGenFunRe7;//differential flow
  86   TProfile2D *fDiffFlowGenFunIm0,*fDiffFlowGenFunIm1,*fDiffFlowGenFunIm2,*fDiffFlowGenFunIm3;//differential flow
  87   TProfile2D *fDiffFlowGenFunIm4,*fDiffFlowGenFunIm5,*fDiffFlowGenFunIm6,*fDiffFlowGenFunIm7;//differential flow
  88
  89   AliFlowCommonHist* fCommonHists;      //common control histograms
  90
  91   TH1D*              fQDist;            //q-distribution
  92
  93   ClassDef(AliFlowAnalysisWithCumulants, 0);
  94 };
  95
  96 //================================================================================================================
  97
  98 #endif
  99
 100
 101
 102
 103