[u/mrichter/AliRoot.git] / PWG2 / FLOW / AliFlowAnalysisWithCumulants.h

#ifndef AliFlowAnalysisWithCumulants_H
#define AliFlowAnalysisWithCumulants_H

//******************************* 
// flow analysis with cumulants *   
// author: Ante Bilandzic       * 
// email: anteb@nikhef.nl       *
//******************************* 

#include "AliFlowCommonConstants.h"
#include "AliFlowCumuConstants.h"

class TH1;
class TObjArray;
class TList;
class TFile;

class AliFlowEventSimple;
class AliFlowTrackSimple;
class AliFlowCommonHist;
class AliFlowCommonHistResults;
class AliFlowVector;

class AliFlowAnalysisWithCumulants {
 public:
  AliFlowAnalysisWithCumulants();
  virtual ~AliFlowAnalysisWithCumulants();
  
  virtual void CreateOutputObjects();
  virtual void Exec(AliFlowEventSimple* anEvent);
  virtual void Terminate(Int_t nEvents);

 private:
  AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants&);
  AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants&);


  AliFlowTrackSimple* fTrack;//track
  static const Int_t fgkQmax=AliFlowCumuConstants::kQmax;//needed for numerics
  static const Int_t fgkPmax=AliFlowCumuConstants::kPmax;//needed for numerics  
  static const Int_t fgkFlow=AliFlowCumuConstants::kFlow;//integrated flow coefficient to be calculated
  static const Int_t fgkMltpl=AliFlowCumuConstants::kMltpl;//the multiple in p=m*n (diff. flow) 
  static const Int_t fgknBins=100;//number of pt bins
      
  Double_t fAvM;//avarage SELECTED multiplicity

  Double_t fR0;//needed for numerics
  Double_t fPtMax;//maximum pt
  Double_t fPtMin;//minimum pt
  Double_t fBinWidth;//width of pt bin (in GeV)
      
  Double_t fAvQx;//<Q_x>
  Double_t fAvQy;//<Q_y>
  Double_t fAvQ2x;//<(Q_x)^2>
  Double_t fAvQ2y;//<(Q_y)^2>
 
  AliFlowCommonHist* fCommonHists;//control histograms
  AliFlowCommonHistResults *fCommonHistsRes2, *fCommonHistsRes4, *fCommonHistsRes6, *fCommonHistsRes8;//histograms with various order final results 
  
  Double_t fAvG[fgkPmax][fgkQmax];//avarage of the generating function used for integrated flow
  Int_t fBinEventEntries[fgknBins];//counts how many events have at least 1 particle in particular bin
  Int_t fBinNoOfParticles[fgknBins];//number of particles per bin
  Double_t fBinMeanPt[fgknBins];//mean pt per bin
  Double_t fBinEventDRe[fgknBins][fgkPmax][fgkQmax];//real part of the generating function used for differential flow
  Double_t fBinEventDIm[fgknBins][fgkPmax][fgkQmax];//imaginary part of the generating function used for differential flow
  
  ClassDef(AliFlowAnalysisWithCumulants, 0);
};
#endif
Commit	Line	Data
f1d945a1	1	#ifndef AliFlowAnalysisWithCumulants_H
	2	#define AliFlowAnalysisWithCumulants_H
	3
	4	//*******************************
	5	// flow analysis with cumulants *
	6	// author: Ante Bilandzic *
	7	// email: anteb@nikhef.nl *
	8	//*******************************
	9
	10	#include "AliFlowCommonConstants.h"
	11	#include "AliFlowCumuConstants.h"
	12
	13	class TH1;
	14	class TObjArray;
	15	class TList;
	16	class TFile;
	17
	18	class AliFlowEventSimple;
	19	class AliFlowTrackSimple;
	20	class AliFlowCommonHist;
	21	class AliFlowCommonHistResults;
	22	class AliFlowVector;
	23
	24	class AliFlowAnalysisWithCumulants {
	25	public:
	26	AliFlowAnalysisWithCumulants();
	27	virtual ~AliFlowAnalysisWithCumulants();
729ec982	28
f1d945a1	29	virtual void CreateOutputObjects();
fa92c1a0	30	virtual void Exec(AliFlowEventSimple* anEvent);
fa92c1a0	31	virtual void Terminate(Int_t nEvents);
729ec982	32
f1d945a1	33	private:
fa92c1a0	34	AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants&);
	35	AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants&);
	36
	37
f1d945a1	38	AliFlowTrackSimple* fTrack;//track
	39	static const Int_t fgkQmax=AliFlowCumuConstants::kQmax;//needed for numerics
	40	static const Int_t fgkPmax=AliFlowCumuConstants::kPmax;//needed for numerics
	41	static const Int_t fgkFlow=AliFlowCumuConstants::kFlow;//integrated flow coefficient to be calculated
	42	static const Int_t fgkMltpl=AliFlowCumuConstants::kMltpl;//the multiple in p=m*n (diff. flow)
	43	static const Int_t fgknBins=100;//number of pt bins
	44
f1d945a1	45	Double_t fAvM;//avarage SELECTED multiplicity
	46
	47	Double_t fR0;//needed for numerics
	48	Double_t fPtMax;//maximum pt
	49	Double_t fPtMin;//minimum pt
	50	Double_t fBinWidth;//width of pt bin (in GeV)
	51
	52	Double_t fAvQx;//<Q_x>
	53	Double_t fAvQy;//<Q_y>
	54	Double_t fAvQ2x;//<(Q_x)^2>
	55	Double_t fAvQ2y;//<(Q_y)^2>
	56
f1d945a1	57	AliFlowCommonHist* fCommonHists;//control histograms
	58	AliFlowCommonHistResults fCommonHistsRes2, fCommonHistsRes4, fCommonHistsRes6, fCommonHistsRes8;//histograms with various order final results
	59
	60	Double_t fAvG[fgkPmax][fgkQmax];//avarage of the generating function used for integrated flow
	61	Int_t fBinEventEntries[fgknBins];//counts how many events have at least 1 particle in particular bin
	62	Int_t fBinNoOfParticles[fgknBins];//number of particles per bin
	63	Double_t fBinMeanPt[fgknBins];//mean pt per bin
	64	Double_t fBinEventDRe[fgknBins][fgkPmax][fgkQmax];//real part of the generating function used for differential flow
	65	Double_t fBinEventDIm[fgknBins][fgkPmax][fgkQmax];//imaginary part of the generating function used for differential flow
	66
	67	ClassDef(AliFlowAnalysisWithCumulants, 0);
	68	};
	69	#endif
	70
	71
	72