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

/* 
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. 
 * See cxx source for full Copyright notice 
 * $Id$ 
 */

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

#ifndef ALIFLOWANALYSISWITHCUMULANTS_H
#define ALIFLOWANALYSISWITHCUMULANTS_H

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

class TObjArray;
class TList;
class TFile;

class TH1D;
class TProfile;
class TProfile2D;
class TProfile3D;
class TDirectoryFile;

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

//================================================================================================================

class AliFlowAnalysisWithCumulants{
 public:
  AliFlowAnalysisWithCumulants();
  virtual ~AliFlowAnalysisWithCumulants(); 
  
  virtual void Init();
  virtual void Make(AliFlowEventSimple* anEvent);
  virtual void GetOutputHistograms(TList *outputListHistos); //get pointers to all output histograms (called before Finish()) 
  virtual void Finish();
  virtual void WriteHistograms(TString* outputFileName);
  virtual void WriteHistograms(TString outputFileName);
  virtual void WriteHistograms(TDirectoryFile *outputFileName);
  
//----------------------------------------------------------------------------------------------------------------
//                                            setters and getters                                                 
//----------------------------------------------------------------------------------------------------------------
  TList* GetHistList() const {return this->fHistList;}     
  
  void SetWeightsList(TList* wlist) {this->fWeightsList = wlist;}
  TList* GetWeightsList() const {return this->fWeightsList;}
  
  void SetIntFlowResults(TH1D* const ifr)  {this->fIntFlowResultsGFC = ifr;};
  TH1D* GetIntFlowResults() const    {return this->fIntFlowResultsGFC;};
  
  void SetDiffFlowResults2nd(TH1D* const diff2nd)  {this->fDiffFlowResults2ndOrderGFC = diff2nd;};
  TH1D* GetDiffFlowResults2nd() const        {return this->fDiffFlowResults2ndOrderGFC;};
  
  void SetDiffFlowResults4th(TH1D* const diff4th)  {this->fDiffFlowResults4thOrderGFC = diff4th;};
  TH1D* GetDiffFlowResults4th() const        {return this->fDiffFlowResults4thOrderGFC;};
  
  void SetDiffFlowResults6th(TH1D* const diff6th)  {this->fDiffFlowResults6thOrderGFC = diff6th;};
  TH1D* GetDiffFlowResults6th() const        {return this->fDiffFlowResults6thOrderGFC;};
  
  void SetDiffFlowResults8th(TH1D* const diff8th)  {this->fDiffFlowResults8thOrderGFC = diff8th;};
  TH1D* GetDiffFlowResults8th() const        {return this->fDiffFlowResults8thOrderGFC;};
  
  void SetCommonHistsResults2nd(AliFlowCommonHistResults* const chr2nd)  {this->fCommonHistsResults2nd = chr2nd;};
  AliFlowCommonHistResults* GetCommonHistsResults2nd() const       {return this->fCommonHistsResults2nd;};
  
  void SetCommonHistsResults4th(AliFlowCommonHistResults* const chr4th)  {this->fCommonHistsResults4th = chr4th;};
  AliFlowCommonHistResults* GetCommonHistsResults4th() const       {return this->fCommonHistsResults4th;};
  
  void SetCommonHistsResults6th(AliFlowCommonHistResults* const chr6th)  {this->fCommonHistsResults6th = chr6th;};
  AliFlowCommonHistResults* GetCommonHistsResults6th() const       {return this->fCommonHistsResults6th;};
  
  void SetCommonHistsResults8th(AliFlowCommonHistResults* const chr8th)  {this->fCommonHistsResults8th = chr8th;};
  AliFlowCommonHistResults* GetCommonHistsResults8th() const       {return this->fCommonHistsResults8th;};
  
  void SetIntFlowGenFun(TProfile2D* const ifgf)  {this->fIntFlowGenFun = ifgf;};
  TProfile2D* GetIntFlowGenFun() const     {return this->fIntFlowGenFun;};

  void SetIntFlowGenFun4(TProfile2D* const ifgf4)  {this->fIntFlowGenFun4 = ifgf4;}; //(only for other system of Eq.)
  TProfile2D* GetIntFlowGenFun4() const      {return this->fIntFlowGenFun4;};  //(only for other system of Eq.)
  
  void SetIntFlowGenFun6(TProfile2D* const ifgf6)  {this->fIntFlowGenFun6 = ifgf6;}; //(only for other system of Eq.)
  TProfile2D* GetIntFlowGenFun6() const      {return this->fIntFlowGenFun6;};  //(only for other system of Eq.)
  
  void SetIntFlowGenFun8(TProfile2D* const ifgf8)  {this->fIntFlowGenFun8 = ifgf8;}; //(only for other system of Eq.)
  TProfile2D* GetIntFlowGenFun8() const      {return this->fIntFlowGenFun8;};  //(only for other system of Eq.)
  
  void SetIntFlowGenFun16(TProfile2D* const ifgf16)  {this->fIntFlowGenFun16 = ifgf16;}; //(only for other system of Eq.)
  TProfile2D* GetIntFlowGenFun16() const       {return this->fIntFlowGenFun16;};   //(only for other system of Eq.)
 
  void SetAverageMultiplicity4(TProfile* const am4)     {this->fAvMultIntFlow4GFC = am4;};  //(only for other system of Eq.)
  TProfile* GetAverageMultiplicity4() const       {return this->fAvMultIntFlow4GFC;}; //(only for other system of Eq.)
  
  void SetAverageMultiplicity6(TProfile* const am6)     {this->fAvMultIntFlow6GFC = am6;};  //(only for other system of Eq.)
  TProfile* GetAverageMultiplicity6() const       {return this->fAvMultIntFlow6GFC;}; //(only for other system of Eq.)
  
  void SetAverageMultiplicity8(TProfile* const am8)     {this->fAvMultIntFlow8GFC = am8;};  //(only for other system of Eq.)
  TProfile* GetAverageMultiplicity8() const       {return this->fAvMultIntFlow8GFC;}; //(only for other system of Eq.)
  
  void SetAverageMultiplicity16(TProfile* const am16)    {this->fAvMultIntFlow16GFC = am16;}; //(only for other system of Eq.)
  TProfile* GetAverageMultiplicity16() const       {return this->fAvMultIntFlow16GFC;}; //(only for other system of Eq.)
  
  //----------------------------------------------------------------------------------------------------------
  //RP = Reaction Plane (RP denotes particles used to determine the reaction plane)
  void SetDiffFlowPtRPGenFunRe(TProfile3D* const dfRPPtgfRe)  {this->fDiffFlowPtRPGenFunRe = dfRPPtgfRe;};
  TProfile3D* GetDiffFlowPtRPGenFunRe() const            {return this->fDiffFlowPtRPGenFunRe;};
  
  void SetDiffFlowPtRPGenFunIm(TProfile3D* const dfRPPtgfIm)  {this->fDiffFlowPtRPGenFunIm = dfRPPtgfIm;};
  TProfile3D* GetDiffPtRPFlowGenFunIm() const            {return this->fDiffFlowPtRPGenFunIm;};
  
  void SetNumberOfParticlesPerPtBinRP(TProfile* const nopppbRP) {this->fPtBinRPNoOfParticles = nopppbRP;};
  TProfile* GetNumberOfParticlesPerPtBinRP() const         {return this->fPtBinRPNoOfParticles;};  
  
  void SetDiffFlowEtaRPGenFunRe(TProfile3D* const dfEtaRPgfRe)   {this->fDiffFlowEtaRPGenFunRe = dfEtaRPgfRe;};
  TProfile3D* GetDiffFlowEtaRPGenFunRe() const              {return this->fDiffFlowEtaRPGenFunRe;};
  
  void SetDiffFlowEtaRPGenFunIm(TProfile3D* const dfEtaRPgfIm)   {this->fDiffFlowEtaRPGenFunIm = dfEtaRPgfIm;};
  TProfile3D* GetDiffFlowEtaRPGenFunIm() const              {return this->fDiffFlowEtaRPGenFunIm;};
  
  void SetNumberOfParticlesPerEtaBinRP(TProfile* const noppebRP)    {this->fEtaBinRPNoOfParticles = noppebRP;};
  TProfile* GetNumberOfParticlesPerEtaBinRP() const            {return this->fEtaBinRPNoOfParticles;}; 
  //----------------------------------------------------------------------------------------------------------
   
  //----------------------------------------------------------------------------------------------------------
  //POI = Particles Of Interest (POI denotes particles of interest for the final physical results for int. and diff. flow)
  void SetDiffFlowPtPOIGenFunRe(TProfile3D* const dfPOIPtgfRe)  {this->fDiffFlowPtPOIGenFunRe = dfPOIPtgfRe;};
  TProfile3D* GetDiffFlowPtPOIGenFunRe() const            {return this->fDiffFlowPtPOIGenFunRe;};
  
  void SetDiffFlowPtPOIGenFunIm(TProfile3D* const dfPOIPtgfIm)  {this->fDiffFlowPtPOIGenFunIm = dfPOIPtgfIm;};
  TProfile3D* GetDiffPtPOIFlowGenFunIm() const            {return this->fDiffFlowPtPOIGenFunIm;};
  
  void SetNumberOfParticlesPerPtBinPOI(TProfile* const nopppbPOI) {this->fPtBinPOINoOfParticles = nopppbPOI;};
  TProfile* GetNumberOfParticlesPerPtBinPOI() const         {return this->fPtBinPOINoOfParticles;};  
  
  void SetDiffFlowEtaPOIGenFunRe(TProfile3D* const dfEtaPOIgfRe)   {this->fDiffFlowEtaPOIGenFunRe = dfEtaPOIgfRe;};
  TProfile3D* GetDiffFlowEtaPOIGenFunRe() const              {return this->fDiffFlowEtaPOIGenFunRe;};
  
  void SetDiffFlowEtaPOIGenFunIm(TProfile3D* const dfEtaPOIgfIm)   {this->fDiffFlowEtaPOIGenFunIm = dfEtaPOIgfIm;};
  TProfile3D* GetDiffFlowEtaPOIGenFunIm() const              {return this->fDiffFlowEtaPOIGenFunIm;};
  
  void SetNumberOfParticlesPerEtaBinPOI(TProfile* const noppebPOI)    {this->fEtaBinPOINoOfParticles = noppebPOI;};
  TProfile* GetNumberOfParticlesPerEtaBinPOI() const            {return this->fEtaBinPOINoOfParticles;};  
  //----------------------------------------------------------------------------------------------------------
    
  void SetAverageMultiplicity(TProfile* const am)      {this->fAvMultIntFlowGFC = am;};
  TProfile* GetAverageMultiplicity() const       {return this->fAvMultIntFlowGFC;};
  
  void SetQVectorComponents(TProfile* const sqvc)    {this->fQVectorComponentsGFC = sqvc;};
  TProfile* GetQVectorComponents() const       {return this->fQVectorComponentsGFC;};
  
  void SetCommonHists(AliFlowCommonHist* const ch)  {this->fCommonHists = ch;};
  AliFlowCommonHist* GetCommonHists() const   {return this->fCommonHists;};
  
  void SetUsePhiWeights(Bool_t const uPhiW) {this->fUsePhiWeights = uPhiW;};
  Bool_t GetUsePhiWeights() const {return this->fUsePhiWeights;};
  
  void SetUsePtWeights(Bool_t const uPtW) {this->fUsePtWeights = uPtW;};
  Bool_t GetUsePtWeights() const {return this->fUsePtWeights;};
  
  void SetUseEtaWeights(Bool_t const uEtaW) {this->fUseEtaWeights = uEtaW;};
  Bool_t GetUseEtaWeights() const {return this->fUseEtaWeights;};
  
  void SetAverageOfSquaredWeight(TProfile* const aosw) {this->fAverageOfSquaredWeight = aosw;};
  TProfile* GetSumOfSquaredWeight() const {return this->fAverageOfSquaredWeight;}; 
//----------------------------------------------------------------------------------------------------------------
 
 private:
  AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants& afawc);
  AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants& afawc);
  AliFlowTrackSimple* fTrack;                                   //track
  static const Int_t fgkQmax = AliFlowCumuConstants::fgkQmax;     //needed for numerics
  static const Int_t fgkPmax = AliFlowCumuConstants::fgkPmax;     //needed for numerics  
  static const Int_t fgkQmax4 = AliFlowCumuConstants::fgkQmax4;   //needed for numerics (only for different system of Eq.)
  static const Int_t fgkPmax4 = AliFlowCumuConstants::fgkPmax4;   //needed for numerics (only for different system of Eq.) 
  static const Int_t fgkQmax6 = AliFlowCumuConstants::fgkQmax6;   //needed for numerics (only for different system of Eq.)
  static const Int_t fgkPmax6 = AliFlowCumuConstants::fgkPmax6;   //needed for numerics (only for different system of Eq.) 
  static const Int_t fgkQmax8 = AliFlowCumuConstants::fgkQmax8;   //needed for numerics (only for different system of Eq.)
  static const Int_t fgkPmax8 = AliFlowCumuConstants::fgkPmax8;   //needed for numerics (only for different system of Eq.)  
  static const Int_t fgkQmax16 = AliFlowCumuConstants::fgkQmax16; //needed for numerics (only for different system of Eq.)
  static const Int_t fgkPmax16 = AliFlowCumuConstants::fgkPmax16; //needed for numerics (only for different system of Eq.)
  static const Int_t fgkFlow = AliFlowCumuConstants::fgkFlow;     //integrated flow coefficient to be calculated
  static const Int_t fgkMltpl = AliFlowCumuConstants::fgkMltpl;   //the multiple in p=m*n (diff. flow) 
  TList* fHistList;                                             //list to hold all output histograms
  TList* fWeightsList;                                          //list to hold all histograms with weights

  Double_t fR0;         //needed for numerics
  Double_t fPtMax;      //maximum pt
  Double_t fPtMin;      //minimum pt
  Double_t fBinWidthPt; //width of pt bin (in GeV)
  Int_t    fgknBinsPt;  //number of pt bins
  
  Double_t fEtaMax;      //maximum pt
  Double_t fEtaMin;      //minimum pt
  Double_t fBinWidthEta; //width of pt bin (in GeV)
  Int_t    fgknBinsEta;  //number of pt bins
      
  Double_t fAvQx;  //<Q_x>
  Double_t fAvQy;  //<Q_y>
  Double_t fAvQ2x; //<(Q_x)^2>
  Double_t fAvQ2y; //<(Q_y)^2>

  TProfile*          fAvMultIntFlowGFC; //average selected multiplicity
 
  TProfile*          fQVectorComponentsGFC; //averages of Q-vector components (1st bin: <Q_x>, 2nd bin: <Q_y>, 3rd bin: <(Q_x)^2>, 4th bin: <(Q_y)^2>)
    
  TH1D*              fIntFlowResultsGFC; //integrated flow final results
  
  TH1D*              fDiffFlowResults2ndOrderGFC; //differential flow final results (2nd order estimate) 
  TH1D*              fDiffFlowResults4thOrderGFC; //differential flow final results (4th order estimate)
  TH1D*              fDiffFlowResults6thOrderGFC; //differential flow final results (6th order estimate)
  TH1D*              fDiffFlowResults8thOrderGFC; //differential flow final results (8th order estimate)
  
  AliFlowCommonHistResults*  fCommonHistsResults2nd; //final results for 2nd order int. and diff. flow stored in the common histograms 
  AliFlowCommonHistResults*  fCommonHistsResults4th; //final results for 4th order int. and diff. flow stored in the common histograms 
  AliFlowCommonHistResults*  fCommonHistsResults6th; //final results for 6th order int. and diff. flow stored in the common histograms
  AliFlowCommonHistResults*  fCommonHistsResults8th; //final results for 8th order int. and diff. flow stored in the common histograms
  
  TProfile2D*        fIntFlowGenFun;   //avarage of the generating function for integrated flow 
  
  TProfile2D*        fIntFlowGenFun4;  //avarage of the generating function for integrated flow (only for different system of Eq.)
  TProfile2D*        fIntFlowGenFun6;  //avarage of the generating function for integrated flow (only for different system of Eq.)
  TProfile2D*        fIntFlowGenFun8;  //avarage of the generating function for integrated flow (only for different system of Eq.)
  TProfile2D*        fIntFlowGenFun16; //avarage of the generating function for integrated flow (only for different system of Eq.)

  TProfile*          fAvMultIntFlow4GFC;  //average selected multiplicity (only for different system of Eq.)
  TProfile*          fAvMultIntFlow6GFC;  //average selected multiplicity (only for different system of Eq.) 
  TProfile*          fAvMultIntFlow8GFC;  //average selected multiplicity (only for different system of Eq.)
  TProfile*          fAvMultIntFlow16GFC; //average selected multiplicity (only for different system of Eq.)

  //RP = Reaction Plane (RP denotes particles used to determine the reaction plane)                                                       
  TProfile3D*        fDiffFlowPtRPGenFunRe; //avarage of the generating function for differential flow in Pt (real part)
  TProfile3D*        fDiffFlowPtRPGenFunIm; //avarage of the generating function for differential flow in Pt (imaginary part)
  TProfile*          fPtBinRPNoOfParticles; //number of particles per pt bin
  TProfile3D*        fDiffFlowEtaRPGenFunRe; //avarage of the generating function for differential flow in Eta (real part)
  TProfile3D*        fDiffFlowEtaRPGenFunIm; //avarage of the generating function for differential flow in Eta (imaginary part)
  TProfile*          fEtaBinRPNoOfParticles; //number of particles per eta bin
                                    
  //POI = Particles Of Interest (POI denotes particles of interest for the final physical results for int. and diff. flow)                                                        
  TProfile3D*        fDiffFlowPtPOIGenFunRe; //avarage of the generating function for differential flow in Pt (real part)
  TProfile3D*        fDiffFlowPtPOIGenFunIm; //avarage of the generating function for differential flow in Pt (imaginary part)
  TProfile*          fPtBinPOINoOfParticles; //number of particles per pt bin
  TProfile3D*        fDiffFlowEtaPOIGenFunRe; //avarage of the generating function for differential flow in Eta (real part)
  TProfile3D*        fDiffFlowEtaPOIGenFunIm; //avarage of the generating function for differential flow in Eta (imaginary part)
  TProfile*          fEtaBinPOINoOfParticles; //number of particles per eta bin
  
  /*
  TProfile2D *fDiffFlowGenFunRe0,*fDiffFlowGenFunRe1,*fDiffFlowGenFunRe2,*fDiffFlowGenFunRe3;//differential flow 
  TProfile2D *fDiffFlowGenFunRe4,*fDiffFlowGenFunRe5,*fDiffFlowGenFunRe6,*fDiffFlowGenFunRe7;//differential flow
  TProfile2D *fDiffFlowGenFunIm0,*fDiffFlowGenFunIm1,*fDiffFlowGenFunIm2,*fDiffFlowGenFunIm3;//differential flow
  TProfile2D *fDiffFlowGenFunIm4,*fDiffFlowGenFunIm5,*fDiffFlowGenFunIm6,*fDiffFlowGenFunIm7;//differential flow
  */
 
  AliFlowCommonHist* fCommonHists;      //common control histograms
  
  Bool_t fOtherEquations; //numerical equations for cumulants solved up to different highest order 
  
  Bool_t fUsePhiWeights;  //phi weights
  Bool_t fUsePtWeights;   //v_2(pt) weights
  Bool_t fUseEtaWeights;  //v_2(eta) weights
  
  TProfile* fAverageOfSquaredWeight; //<w^2>      
  
  ClassDef(AliFlowAnalysisWithCumulants, 0);
};

//================================================================================================================

#endif
Commit	Line	Data
2188af53	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
b6cd16a9	14	#ifndef ALIFLOWANALYSISWITHCUMULANTS_H
b6cd16a9	15	#define ALIFLOWANALYSISWITHCUMULANTS_H
f1d945a1	16
f1d945a1	17	#include "AliFlowCommonConstants.h"
	18	#include "AliFlowCumuConstants.h"
	19
f1d945a1	20	class TObjArray;
2188af53	21	class TList;
2188af53	22	class TFile;
aaebd73d	23
929098e4	24	class TH1D;
aaebd73d	25	class TProfile;
	26	class TProfile2D;
	27	class TProfile3D;
929098e4	28	class TDirectoryFile;
aaebd73d	29
f1d945a1	30	class AliFlowEventSimple;
	31	class AliFlowTrackSimple;
	32	class AliFlowCommonHist;
	33	class AliFlowCommonHistResults;
	34	class AliFlowVector;
	35
2188af53	36	//================================================================================================================
	37
	38	class AliFlowAnalysisWithCumulants{
f1d945a1	39	public:
f1d945a1	40	AliFlowAnalysisWithCumulants();
2188af53	41	virtual ~AliFlowAnalysisWithCumulants();
729ec982	42
e5e75b58	43	virtual void Init();
dce74562	44	virtual void Make(AliFlowEventSimple* anEvent);
fd46c3dd	45	virtual void GetOutputHistograms(TList *outputListHistos); //get pointers to all output histograms (called before Finish())
dce74562	46	virtual void Finish();
1315fe58	47	virtual void WriteHistograms(TString* outputFileName);
b0fda271	48	virtual void WriteHistograms(TString outputFileName);
ad87ae62	49	virtual void WriteHistograms(TDirectoryFile *outputFileName);
aaebd73d	50
52021ae2	51	//----------------------------------------------------------------------------------------------------------------
e5e75b58	52	// setters and getters
52021ae2	53	//----------------------------------------------------------------------------------------------------------------
e5e75b58	54	TList* GetHistList() const {return this->fHistList;}
	55
	56	void SetWeightsList(TList* wlist) {this->fWeightsList = wlist;}
	57	TList* GetWeightsList() const {return this->fWeightsList;}
52021ae2	58
5270f185	59	void SetIntFlowResults(TH1D* const ifr) {this->fIntFlowResultsGFC = ifr;};
52021ae2	60	TH1D* GetIntFlowResults() const {return this->fIntFlowResultsGFC;};
52021ae2	61
5270f185	62	void SetDiffFlowResults2nd(TH1D* const diff2nd) {this->fDiffFlowResults2ndOrderGFC = diff2nd;};
52021ae2	63	TH1D* GetDiffFlowResults2nd() const {return this->fDiffFlowResults2ndOrderGFC;};
52021ae2	64
5270f185	65	void SetDiffFlowResults4th(TH1D* const diff4th) {this->fDiffFlowResults4thOrderGFC = diff4th;};
52021ae2	66	TH1D* GetDiffFlowResults4th() const {return this->fDiffFlowResults4thOrderGFC;};
52021ae2	67
5270f185	68	void SetDiffFlowResults6th(TH1D* const diff6th) {this->fDiffFlowResults6thOrderGFC = diff6th;};
52021ae2	69	TH1D* GetDiffFlowResults6th() const {return this->fDiffFlowResults6thOrderGFC;};
52021ae2	70
5270f185	71	void SetDiffFlowResults8th(TH1D* const diff8th) {this->fDiffFlowResults8thOrderGFC = diff8th;};
52021ae2	72	TH1D* GetDiffFlowResults8th() const {return this->fDiffFlowResults8thOrderGFC;};
52021ae2	73
5270f185	74	void SetCommonHistsResults2nd(AliFlowCommonHistResults* const chr2nd) {this->fCommonHistsResults2nd = chr2nd;};
52021ae2	75	AliFlowCommonHistResults* GetCommonHistsResults2nd() const {return this->fCommonHistsResults2nd;};
52021ae2	76
5270f185	77	void SetCommonHistsResults4th(AliFlowCommonHistResults* const chr4th) {this->fCommonHistsResults4th = chr4th;};
52021ae2	78	AliFlowCommonHistResults* GetCommonHistsResults4th() const {return this->fCommonHistsResults4th;};
52021ae2	79
5270f185	80	void SetCommonHistsResults6th(AliFlowCommonHistResults* const chr6th) {this->fCommonHistsResults6th = chr6th;};
52021ae2	81	AliFlowCommonHistResults* GetCommonHistsResults6th() const {return this->fCommonHistsResults6th;};
52021ae2	82
5270f185	83	void SetCommonHistsResults8th(AliFlowCommonHistResults* const chr8th) {this->fCommonHistsResults8th = chr8th;};
52021ae2	84	AliFlowCommonHistResults* GetCommonHistsResults8th() const {return this->fCommonHistsResults8th;};
52021ae2	85
5270f185	86	void SetIntFlowGenFun(TProfile2D* const ifgf) {this->fIntFlowGenFun = ifgf;};
b6cd16a9	87	TProfile2D* GetIntFlowGenFun() const {return this->fIntFlowGenFun;};
b6cd16a9	88
5270f185	89	void SetIntFlowGenFun4(TProfile2D* const ifgf4) {this->fIntFlowGenFun4 = ifgf4;}; //(only for other system of Eq.)
cb308e83	90	TProfile2D* GetIntFlowGenFun4() const {return this->fIntFlowGenFun4;}; //(only for other system of Eq.)
b6cd16a9	91
5270f185	92	void SetIntFlowGenFun6(TProfile2D* const ifgf6) {this->fIntFlowGenFun6 = ifgf6;}; //(only for other system of Eq.)
cb308e83	93	TProfile2D* GetIntFlowGenFun6() const {return this->fIntFlowGenFun6;}; //(only for other system of Eq.)
b6cd16a9	94
5270f185	95	void SetIntFlowGenFun8(TProfile2D* const ifgf8) {this->fIntFlowGenFun8 = ifgf8;}; //(only for other system of Eq.)
cb308e83	96	TProfile2D* GetIntFlowGenFun8() const {return this->fIntFlowGenFun8;}; //(only for other system of Eq.)
52021ae2	97
5270f185	98	void SetIntFlowGenFun16(TProfile2D* const ifgf16) {this->fIntFlowGenFun16 = ifgf16;}; //(only for other system of Eq.)
cb308e83	99	TProfile2D* GetIntFlowGenFun16() const {return this->fIntFlowGenFun16;}; //(only for other system of Eq.)
b6cd16a9	100
5270f185	101	void SetAverageMultiplicity4(TProfile* const am4) {this->fAvMultIntFlow4GFC = am4;}; //(only for other system of Eq.)
cb308e83	102	TProfile* GetAverageMultiplicity4() const {return this->fAvMultIntFlow4GFC;}; //(only for other system of Eq.)
cb308e83	103
5270f185	104	void SetAverageMultiplicity6(TProfile* const am6) {this->fAvMultIntFlow6GFC = am6;}; //(only for other system of Eq.)
cb308e83	105	TProfile* GetAverageMultiplicity6() const {return this->fAvMultIntFlow6GFC;}; //(only for other system of Eq.)
cb308e83	106
5270f185	107	void SetAverageMultiplicity8(TProfile* const am8) {this->fAvMultIntFlow8GFC = am8;}; //(only for other system of Eq.)
cb308e83	108	TProfile* GetAverageMultiplicity8() const {return this->fAvMultIntFlow8GFC;}; //(only for other system of Eq.)
cb308e83	109
5270f185	110	void SetAverageMultiplicity16(TProfile* const am16) {this->fAvMultIntFlow16GFC = am16;}; //(only for other system of Eq.)
cb308e83	111	TProfile* GetAverageMultiplicity16() const {return this->fAvMultIntFlow16GFC;}; //(only for other system of Eq.)
52021ae2	112
813a4157	113	//----------------------------------------------------------------------------------------------------------
813a4157	114	//RP = Reaction Plane (RP denotes particles used to determine the reaction plane)
5270f185	115	void SetDiffFlowPtRPGenFunRe(TProfile3D* const dfRPPtgfRe) {this->fDiffFlowPtRPGenFunRe = dfRPPtgfRe;};
813a4157	116	TProfile3D* GetDiffFlowPtRPGenFunRe() const {return this->fDiffFlowPtRPGenFunRe;};
813a4157	117
5270f185	118	void SetDiffFlowPtRPGenFunIm(TProfile3D* const dfRPPtgfIm) {this->fDiffFlowPtRPGenFunIm = dfRPPtgfIm;};
813a4157	119	TProfile3D* GetDiffPtRPFlowGenFunIm() const {return this->fDiffFlowPtRPGenFunIm;};
813a4157	120
5270f185	121	void SetNumberOfParticlesPerPtBinRP(TProfile* const nopppbRP) {this->fPtBinRPNoOfParticles = nopppbRP;};
813a4157	122	TProfile* GetNumberOfParticlesPerPtBinRP() const {return this->fPtBinRPNoOfParticles;};
813a4157	123
5270f185	124	void SetDiffFlowEtaRPGenFunRe(TProfile3D* const dfEtaRPgfRe) {this->fDiffFlowEtaRPGenFunRe = dfEtaRPgfRe;};
813a4157	125	TProfile3D* GetDiffFlowEtaRPGenFunRe() const {return this->fDiffFlowEtaRPGenFunRe;};
813a4157	126
5270f185	127	void SetDiffFlowEtaRPGenFunIm(TProfile3D* const dfEtaRPgfIm) {this->fDiffFlowEtaRPGenFunIm = dfEtaRPgfIm;};
813a4157	128	TProfile3D* GetDiffFlowEtaRPGenFunIm() const {return this->fDiffFlowEtaRPGenFunIm;};
813a4157	129
5270f185	130	void SetNumberOfParticlesPerEtaBinRP(TProfile* const noppebRP) {this->fEtaBinRPNoOfParticles = noppebRP;};
813a4157	131	TProfile* GetNumberOfParticlesPerEtaBinRP() const {return this->fEtaBinRPNoOfParticles;};
	132	//----------------------------------------------------------------------------------------------------------
	133
	134	//----------------------------------------------------------------------------------------------------------
	135	//POI = Particles Of Interest (POI denotes particles of interest for the final physical results for int. and diff. flow)
5270f185	136	void SetDiffFlowPtPOIGenFunRe(TProfile3D* const dfPOIPtgfRe) {this->fDiffFlowPtPOIGenFunRe = dfPOIPtgfRe;};
813a4157	137	TProfile3D* GetDiffFlowPtPOIGenFunRe() const {return this->fDiffFlowPtPOIGenFunRe;};
813a4157	138
5270f185	139	void SetDiffFlowPtPOIGenFunIm(TProfile3D* const dfPOIPtgfIm) {this->fDiffFlowPtPOIGenFunIm = dfPOIPtgfIm;};
813a4157	140	TProfile3D* GetDiffPtPOIFlowGenFunIm() const {return this->fDiffFlowPtPOIGenFunIm;};
52021ae2	141
5270f185	142	void SetNumberOfParticlesPerPtBinPOI(TProfile* const nopppbPOI) {this->fPtBinPOINoOfParticles = nopppbPOI;};
813a4157	143	TProfile* GetNumberOfParticlesPerPtBinPOI() const {return this->fPtBinPOINoOfParticles;};
52021ae2	144
5270f185	145	void SetDiffFlowEtaPOIGenFunRe(TProfile3D* const dfEtaPOIgfRe) {this->fDiffFlowEtaPOIGenFunRe = dfEtaPOIgfRe;};
813a4157	146	TProfile3D* GetDiffFlowEtaPOIGenFunRe() const {return this->fDiffFlowEtaPOIGenFunRe;};
813a4157	147
5270f185	148	void SetDiffFlowEtaPOIGenFunIm(TProfile3D* const dfEtaPOIgfIm) {this->fDiffFlowEtaPOIGenFunIm = dfEtaPOIgfIm;};
813a4157	149	TProfile3D* GetDiffFlowEtaPOIGenFunIm() const {return this->fDiffFlowEtaPOIGenFunIm;};
813a4157	150
5270f185	151	void SetNumberOfParticlesPerEtaBinPOI(TProfile* const noppebPOI) {this->fEtaBinPOINoOfParticles = noppebPOI;};
813a4157	152	TProfile* GetNumberOfParticlesPerEtaBinPOI() const {return this->fEtaBinPOINoOfParticles;};
	153	//----------------------------------------------------------------------------------------------------------
	154
5270f185	155	void SetAverageMultiplicity(TProfile* const am) {this->fAvMultIntFlowGFC = am;};
52021ae2	156	TProfile* GetAverageMultiplicity() const {return this->fAvMultIntFlowGFC;};
52021ae2	157
5270f185	158	void SetQVectorComponents(TProfile* const sqvc) {this->fQVectorComponentsGFC = sqvc;};
52021ae2	159	TProfile* GetQVectorComponents() const {return this->fQVectorComponentsGFC;};
813a4157	160
5270f185	161	void SetCommonHists(AliFlowCommonHist* const ch) {this->fCommonHists = ch;};
813a4157	162	AliFlowCommonHist* GetCommonHists() const {return this->fCommonHists;};
e5e75b58	163
	164	void SetUsePhiWeights(Bool_t const uPhiW) {this->fUsePhiWeights = uPhiW;};
	165	Bool_t GetUsePhiWeights() const {return this->fUsePhiWeights;};
	166
	167	void SetUsePtWeights(Bool_t const uPtW) {this->fUsePtWeights = uPtW;};
	168	Bool_t GetUsePtWeights() const {return this->fUsePtWeights;};
	169
	170	void SetUseEtaWeights(Bool_t const uEtaW) {this->fUseEtaWeights = uEtaW;};
	171	Bool_t GetUseEtaWeights() const {return this->fUseEtaWeights;};
9dd53ff2	172
	173	void SetAverageOfSquaredWeight(TProfile* const aosw) {this->fAverageOfSquaredWeight = aosw;};
	174	TProfile* GetSumOfSquaredWeight() const {return this->fAverageOfSquaredWeight;};
52021ae2	175	//----------------------------------------------------------------------------------------------------------------
924fafb7	176
2188af53	177	private:
	178	AliFlowAnalysisWithCumulants(const AliFlowAnalysisWithCumulants& afawc);
	179	AliFlowAnalysisWithCumulants& operator=(const AliFlowAnalysisWithCumulants& afawc);
b6cd16a9	180	AliFlowTrackSimple* fTrack; //track
bee2efdc	181	static const Int_t fgkQmax = AliFlowCumuConstants::fgkQmax; //needed for numerics
	182	static const Int_t fgkPmax = AliFlowCumuConstants::fgkPmax; //needed for numerics
	183	static const Int_t fgkQmax4 = AliFlowCumuConstants::fgkQmax4; //needed for numerics (only for different system of Eq.)
	184	static const Int_t fgkPmax4 = AliFlowCumuConstants::fgkPmax4; //needed for numerics (only for different system of Eq.)
	185	static const Int_t fgkQmax6 = AliFlowCumuConstants::fgkQmax6; //needed for numerics (only for different system of Eq.)
	186	static const Int_t fgkPmax6 = AliFlowCumuConstants::fgkPmax6; //needed for numerics (only for different system of Eq.)
	187	static const Int_t fgkQmax8 = AliFlowCumuConstants::fgkQmax8; //needed for numerics (only for different system of Eq.)
	188	static const Int_t fgkPmax8 = AliFlowCumuConstants::fgkPmax8; //needed for numerics (only for different system of Eq.)
	189	static const Int_t fgkQmax16 = AliFlowCumuConstants::fgkQmax16; //needed for numerics (only for different system of Eq.)
	190	static const Int_t fgkPmax16 = AliFlowCumuConstants::fgkPmax16; //needed for numerics (only for different system of Eq.)
	191	static const Int_t fgkFlow = AliFlowCumuConstants::fgkFlow; //integrated flow coefficient to be calculated
	192	static const Int_t fgkMltpl = AliFlowCumuConstants::fgkMltpl; //the multiple in p=m*n (diff. flow)
b6cd16a9	193	TList* fHistList; //list to hold all output histograms
e5e75b58	194	TList* fWeightsList; //list to hold all histograms with weights
2188af53	195
813a4157	196	Double_t fR0; //needed for numerics
	197	Double_t fPtMax; //maximum pt
	198	Double_t fPtMin; //minimum pt
	199	Double_t fBinWidthPt; //width of pt bin (in GeV)
	200	Int_t fgknBinsPt; //number of pt bins
	201
	202	Double_t fEtaMax; //maximum pt
	203	Double_t fEtaMin; //minimum pt
	204	Double_t fBinWidthEta; //width of pt bin (in GeV)
	205	Int_t fgknBinsEta; //number of pt bins
f1d945a1	206
2188af53	207	Double_t fAvQx; //<Q_x>
	208	Double_t fAvQy; //<Q_y>
	209	Double_t fAvQ2x; //<(Q_x)^2>
	210	Double_t fAvQ2y; //<(Q_y)^2>
	211
b6cd16a9	212	TProfile* fAvMultIntFlowGFC; //average selected multiplicity
2188af53	213
52021ae2	214	TProfile* fQVectorComponentsGFC; //averages of Q-vector components (1st bin: <Q_x>, 2nd bin: <Q_y>, 3rd bin: <(Q_x)^2>, 4th bin: <(Q_y)^2>)
2188af53	215
b6cd16a9	216	TH1D* fIntFlowResultsGFC; //integrated flow final results
aaebd73d	217
52021ae2	218	TH1D* fDiffFlowResults2ndOrderGFC; //differential flow final results (2nd order estimate)
	219	TH1D* fDiffFlowResults4thOrderGFC; //differential flow final results (4th order estimate)
	220	TH1D* fDiffFlowResults6thOrderGFC; //differential flow final results (6th order estimate)
	221	TH1D* fDiffFlowResults8thOrderGFC; //differential flow final results (8th order estimate)
	222
b6cd16a9	223	AliFlowCommonHistResults* fCommonHistsResults2nd; //final results for 2nd order int. and diff. flow stored in the common histograms
	224	AliFlowCommonHistResults* fCommonHistsResults4th; //final results for 4th order int. and diff. flow stored in the common histograms
	225	AliFlowCommonHistResults* fCommonHistsResults6th; //final results for 6th order int. and diff. flow stored in the common histograms
	226	AliFlowCommonHistResults* fCommonHistsResults8th; //final results for 8th order int. and diff. flow stored in the common histograms
	227
	228	TProfile2D* fIntFlowGenFun; //avarage of the generating function for integrated flow
cb308e83	229
	230	TProfile2D* fIntFlowGenFun4; //avarage of the generating function for integrated flow (only for different system of Eq.)
	231	TProfile2D* fIntFlowGenFun6; //avarage of the generating function for integrated flow (only for different system of Eq.)
	232	TProfile2D* fIntFlowGenFun8; //avarage of the generating function for integrated flow (only for different system of Eq.)
	233	TProfile2D* fIntFlowGenFun16; //avarage of the generating function for integrated flow (only for different system of Eq.)
	234
	235	TProfile* fAvMultIntFlow4GFC; //average selected multiplicity (only for different system of Eq.)
	236	TProfile* fAvMultIntFlow6GFC; //average selected multiplicity (only for different system of Eq.)
	237	TProfile* fAvMultIntFlow8GFC; //average selected multiplicity (only for different system of Eq.)
	238	TProfile* fAvMultIntFlow16GFC; //average selected multiplicity (only for different system of Eq.)
813a4157	239
	240	//RP = Reaction Plane (RP denotes particles used to determine the reaction plane)
	241	TProfile3D* fDiffFlowPtRPGenFunRe; //avarage of the generating function for differential flow in Pt (real part)
	242	TProfile3D* fDiffFlowPtRPGenFunIm; //avarage of the generating function for differential flow in Pt (imaginary part)
	243	TProfile* fPtBinRPNoOfParticles; //number of particles per pt bin
	244	TProfile3D* fDiffFlowEtaRPGenFunRe; //avarage of the generating function for differential flow in Eta (real part)
	245	TProfile3D* fDiffFlowEtaRPGenFunIm; //avarage of the generating function for differential flow in Eta (imaginary part)
	246	TProfile* fEtaBinRPNoOfParticles; //number of particles per eta bin
	247
	248	//POI = Particles Of Interest (POI denotes particles of interest for the final physical results for int. and diff. flow)
	249	TProfile3D* fDiffFlowPtPOIGenFunRe; //avarage of the generating function for differential flow in Pt (real part)
	250	TProfile3D* fDiffFlowPtPOIGenFunIm; //avarage of the generating function for differential flow in Pt (imaginary part)
	251	TProfile* fPtBinPOINoOfParticles; //number of particles per pt bin
	252	TProfile3D* fDiffFlowEtaPOIGenFunRe; //avarage of the generating function for differential flow in Eta (real part)
	253	TProfile3D* fDiffFlowEtaPOIGenFunIm; //avarage of the generating function for differential flow in Eta (imaginary part)
	254	TProfile* fEtaBinPOINoOfParticles; //number of particles per eta bin
aaebd73d	255
52021ae2	256	/*
2188af53	257	TProfile2D fDiffFlowGenFunRe0,fDiffFlowGenFunRe1,fDiffFlowGenFunRe2,fDiffFlowGenFunRe3;//differential flow
	258	TProfile2D fDiffFlowGenFunRe4,fDiffFlowGenFunRe5,fDiffFlowGenFunRe6,fDiffFlowGenFunRe7;//differential flow
	259	TProfile2D fDiffFlowGenFunIm0,fDiffFlowGenFunIm1,fDiffFlowGenFunIm2,fDiffFlowGenFunIm3;//differential flow
	260	TProfile2D fDiffFlowGenFunIm4,fDiffFlowGenFunIm5,fDiffFlowGenFunIm6,fDiffFlowGenFunIm7;//differential flow
52021ae2	261	*/
2188af53	262
2188af53	263	AliFlowCommonHist* fCommonHists; //common control histograms
b6cd16a9	264
b6cd16a9	265	Bool_t fOtherEquations; //numerical equations for cumulants solved up to different highest order
e5e75b58	266
	267	Bool_t fUsePhiWeights; //phi weights
	268	Bool_t fUsePtWeights; //v_2(pt) weights
9dd53ff2	269	Bool_t fUseEtaWeights; //v_2(eta) weights
	270
	271	TProfile* fAverageOfSquaredWeight; //<w^2>
	272
f1d945a1	273	ClassDef(AliFlowAnalysisWithCumulants, 0);
f1d945a1	274	};
2188af53	275
	276	//================================================================================================================
	277
f1d945a1	278	#endif
	279
	280
aaebd73d	281
	282
	283