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

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

#ifndef AliFlowAnalysisWithLYZEventPlane_H
#define AliFlowAnalysisWithLYZEventPlane_H

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

class TString;
class TFile;
class TProfile;
class TH1F;
class TH1D;
class TList;
class Riostream;

// AliFlowAnalysisWithLYZEventPlane:
// Class to do flow analysis with the event plane from the LYZ method
// author: N. van der Kolk (kolk@nikhef.nl)


class AliFlowAnalysisWithLYZEventPlane {

 public:

  AliFlowAnalysisWithLYZEventPlane();                 //default constructor
  virtual ~AliFlowAnalysisWithLYZEventPlane();        //destructor
  
  virtual void   Init();
  virtual void   Make(AliFlowEventSimple* fEvent, AliFlowLYZEventPlane* fLYZEP);
  virtual void   Finish();
  void           WriteHistograms(TString* outputFileName);
  void           WriteHistograms(TString outputFileName);

  void      SetEventNumber(Int_t n)      { this->fEventNumber = n; }
  Int_t     GetEventNumber() const       { return this->fEventNumber; }
  void      SetQ2sum(Double_t d)         { this->fQ2sum = d; }
  Double_t  GetQ2sum()                   { return this->fQ2sum; }

  //output
  TList*             GetHistList() const     {return this->fHistList; }
  AliFlowCommonHist* GetCommonHists() const  { return this->fCommonHists; }
  void               SetCommonHists(AliFlowCommonHist* aCommonHist)  
     { this->fCommonHists = aCommonHist; }
  AliFlowCommonHistResults* GetCommonHistsRes() const  
     { return this->fCommonHistsRes; }
  void               SetCommonHistsRes(AliFlowCommonHistResults* aCommonHistResult) 
     { this->fCommonHistsRes = aCommonHistResult; }

  // !!!!! make getters and setters for all histograms
  TProfile*  GetSecondReDtheta() {return this->fSecondReDtheta; } 
  void       SetSecondReDtheta(TProfile* aSecondReDtheta) 
    {this->fSecondReDtheta = aSecondReDtheta; }
  TProfile*  GetSecondImDtheta() {return this->fSecondImDtheta; }
  void       SetSecondImDtheta(TProfile* aSecondImDtheta)
    {this->fSecondImDtheta = aSecondImDtheta; }
  TProfile*  GetFirstr0theta()   {return this->fFirstr0theta; }
  void       SetFirstr0theta(TProfile* aFirstr0theta)
    {this->fFirstr0theta = aFirstr0theta; }
  
  TProfile*  GetHistProVetaRP()    {return this->fHistProVetaRP;}
  void       SetHistProVetaRP(TProfile* aHistProVetaRP)
    {this->fHistProVetaRP =aHistProVetaRP; }
  TProfile*  GetHistProVetaPOI()   {return this->fHistProVetaPOI;} 
  void       SetHistProVetaPOI(TProfile* aHistProVetaPOI)
    {this->fHistProVetaPOI = aHistProVetaPOI; } 
  TProfile*  GetHistProVPtRP()    {return this->fHistProVPtRP;}
  void       SetHistProVPtRP(TProfile* aHistProVPtRP)
    {this->fHistProVPtRP =aHistProVPtRP; } 
  TProfile*  GetHistProVPtPOI()   {return this->fHistProVPtPOI;} 
  void       SetHistProVPtPOI(TProfile* aHistProVPtPOI)
    {this->fHistProVPtPOI = aHistProVPtPOI; }      
  TProfile*  GetHistProWr()      {return this->fHistProWr; }
  void       SetHistProWr(TProfile* aHistProWr)
    {this->fHistProWr = aHistProWr; }
  TProfile*  GetHistProWrCorr()  {return this->fHistProWrCorr; }
  void       SetHistProWrCorr(TProfile* aHistProWrCorr)
    {this->fHistProWrCorr = aHistProWrCorr; }
  TH1F*      GetHistQsumforChi() {return this->fHistQsumforChi; }
  void       SetHistQsumforChi(TH1F* aHistQsumforChi) 
    {this->fHistQsumforChi = aHistQsumforChi; }
  TH1F*      GetHistDeltaPhi()   {return this->fHistDeltaPhi; }  
  void       SetHistDeltaPhi(TH1F* aHistDeltaPhi)
    {this->fHistDeltaPhi = aHistDeltaPhi; }
  TH1F*      GetHistDeltaPhi2()  {return this->fHistDeltaPhi2; } 
  void       SetHistDeltaPhi2(TH1F* aHistDeltaPhi2)
    {this->fHistDeltaPhi2 = aHistDeltaPhi2; }
  TH1F*      GetHistDeltaPhihere() {return this->fHistDeltaPhihere; }
  void       SetHistDeltaPhihere(TH1F* aHistDeltaPhihere)
    {this->fHistDeltaPhihere = aHistDeltaPhihere; }
  TH1F*      GetHistPhiEP()      {return this->fHistPhiEP; }   
  void       SetHistPhiEP(TH1F* aHistPhiEP)
    {this->fHistPhiEP = aHistPhiEP; }  
  TH1F*      GetHistPhiEPhere()  {return this->fHistPhiEPhere; }
  void       SetHistPhiEPhere(TH1F* aHistPhiEPhere)
    {this->fHistPhiEPhere = aHistPhiEPhere; }       
  TH1F*      GetHistPhiLYZ()     {return this->fHistPhiLYZ; }  
  void       SetHistPhiLYZ(TH1F* aHistPhiLYZ)
    {this->fHistPhiLYZ = aHistPhiLYZ; }       
  TH1F*      GetHistPhiLYZ2()    {return this->fHistPhiLYZ2;}               
  void       SetHistPhiLYZ2(TH1F* aHistPhiLYZ2)
    {this->fHistPhiLYZ2 = aHistPhiLYZ2; }

  //input
  void       SetSecondRunList(TList* list) { this->fSecondRunList = list; }
  TList*     GetSecondRunList()            { return this->fSecondRunList; }

 private:

  AliFlowAnalysisWithLYZEventPlane(const AliFlowAnalysisWithLYZEventPlane& aAnalysis);             // copy constructor
  AliFlowAnalysisWithLYZEventPlane& operator=(const AliFlowAnalysisWithLYZEventPlane& aAnalysis);  // assignment operator

  //histograms
  TList*     fHistList;                       //list ro hold all histograms
  TList*     fSecondRunList;                  //list from Second LYZ run output
  
  //input
  TProfile*  fSecondReDtheta;                 // input profile
  TProfile*  fSecondImDtheta;                 // input profile
  TProfile*  fFirstr0theta;                   // input profile
  
  //output
  TProfile*  fHistProVetaRP;                  //
  TProfile*  fHistProVetaPOI;                 //
  TProfile*  fHistProVPtRP;                   //
  TProfile*  fHistProVPtPOI;                  //
  TProfile*  fHistProWr;                      //
  TProfile*  fHistProWrCorr;                  //
  TH1F*      fHistQsumforChi;                 //
  TH1F*      fHistDeltaPhi;                   //
  TH1F*      fHistDeltaPhi2;                  //
  TH1F*      fHistDeltaPhihere;               //
  TH1F*      fHistPhiEP;                      //
  TH1F*      fHistPhiEPhere;                  //
  TH1F*      fHistPhiLYZ;                     //
  TH1F*      fHistPhiLYZ2;                    //
  
  AliFlowCommonHist* fCommonHists;            //
  AliFlowCommonHistResults* fCommonHistsRes;  //

  Int_t     fEventNumber;                     // event counter

  TVector2  *fQsum;                           // flow vector sum
  Double_t  fQ2sum;                           // flow vector sum squared
     

  ClassDef(AliFlowAnalysisWithLYZEventPlane, 1);          // lyz analysis 
};

 #endif
Commit	Line	Data
f1d945a1	1	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	2	* See cxx source for full Copyright notice */
	3	/* $Id$ */
	4
	5	#ifndef AliFlowAnalysisWithLYZEventPlane_H
	6	#define AliFlowAnalysisWithLYZEventPlane_H
	7
	8	class AliFlowVector;
	9	class AliFlowTrackSimple;
	10	class AliFlowEventSimple;
	11	class AliFlowCommonHist;
	12	class AliFlowCommonHistResults;
	13	class AliFlowLYZEventPlane;
	14
	15	class TString;
	16	class TFile;
	17	class TProfile;
	18	class TH1F;
	19	class TH1D;
28ca24ad	20	class TList;
28ca24ad	21	class Riostream;
f1d945a1	22
	23	// AliFlowAnalysisWithLYZEventPlane:
	24	// Class to do flow analysis with the event plane from the LYZ method
	25	// author: N. van der Kolk (kolk@nikhef.nl)
	26
	27
	28	class AliFlowAnalysisWithLYZEventPlane {
28ca24ad	29
f1d945a1	30	public:
28ca24ad	31
	32	AliFlowAnalysisWithLYZEventPlane(); //default constructor
	33	virtual ~AliFlowAnalysisWithLYZEventPlane(); //destructor
f1d945a1	34
	35	virtual void Init();
	36	virtual void Make(AliFlowEventSimple* fEvent, AliFlowLYZEventPlane* fLYZEP);
	37	virtual void Finish();
1dfa3c16	38	void WriteHistograms(TString* outputFileName);
b0fda271	39	void WriteHistograms(TString outputFileName);
f1d945a1	40
9d062fe3	41	void SetEventNumber(Int_t n) { this->fEventNumber = n; }
	42	Int_t GetEventNumber() const { return this->fEventNumber; }
	43	void SetQ2sum(Double_t d) { this->fQ2sum = d; }
	44	Double_t GetQ2sum() { return this->fQ2sum; }
28ca24ad	45
28ca24ad	46	//output
9d062fe3	47	TList* GetHistList() const {return this->fHistList; }
	48	AliFlowCommonHist* GetCommonHists() const { return this->fCommonHists; }
	49	void SetCommonHists(AliFlowCommonHist* aCommonHist)
	50	{ this->fCommonHists = aCommonHist; }
	51	AliFlowCommonHistResults* GetCommonHistsRes() const
	52	{ return this->fCommonHistsRes; }
	53	void SetCommonHistsRes(AliFlowCommonHistResults* aCommonHistResult)
	54	{ this->fCommonHistsRes = aCommonHistResult; }
	55
	56	// !!!!! make getters and setters for all histograms
	57	TProfile* GetSecondReDtheta() {return this->fSecondReDtheta; }
	58	void SetSecondReDtheta(TProfile* aSecondReDtheta)
	59	{this->fSecondReDtheta = aSecondReDtheta; }
	60	TProfile* GetSecondImDtheta() {return this->fSecondImDtheta; }
	61	void SetSecondImDtheta(TProfile* aSecondImDtheta)
	62	{this->fSecondImDtheta = aSecondImDtheta; }
	63	TProfile* GetFirstr0theta() {return this->fFirstr0theta; }
	64	void SetFirstr0theta(TProfile* aFirstr0theta)
	65	{this->fFirstr0theta = aFirstr0theta; }
80e93fef	66
	67	TProfile* GetHistProVetaRP() {return this->fHistProVetaRP;}
	68	void SetHistProVetaRP(TProfile* aHistProVetaRP)
	69	{this->fHistProVetaRP =aHistProVetaRP; }
	70	TProfile* GetHistProVetaPOI() {return this->fHistProVetaPOI;}
	71	void SetHistProVetaPOI(TProfile* aHistProVetaPOI)
	72	{this->fHistProVetaPOI = aHistProVetaPOI; }
	73	TProfile* GetHistProVPtRP() {return this->fHistProVPtRP;}
	74	void SetHistProVPtRP(TProfile* aHistProVPtRP)
	75	{this->fHistProVPtRP =aHistProVPtRP; }
	76	TProfile* GetHistProVPtPOI() {return this->fHistProVPtPOI;}
	77	void SetHistProVPtPOI(TProfile* aHistProVPtPOI)
	78	{this->fHistProVPtPOI = aHistProVPtPOI; }
9d062fe3	79	TProfile* GetHistProWr() {return this->fHistProWr; }
	80	void SetHistProWr(TProfile* aHistProWr)
	81	{this->fHistProWr = aHistProWr; }
	82	TProfile* GetHistProWrCorr() {return this->fHistProWrCorr; }
	83	void SetHistProWrCorr(TProfile* aHistProWrCorr)
	84	{this->fHistProWrCorr = aHistProWrCorr; }
	85	TH1F* GetHistQsumforChi() {return this->fHistQsumforChi; }
	86	void SetHistQsumforChi(TH1F* aHistQsumforChi)
	87	{this->fHistQsumforChi = aHistQsumforChi; }
	88	TH1F* GetHistDeltaPhi() {return this->fHistDeltaPhi; }
	89	void SetHistDeltaPhi(TH1F* aHistDeltaPhi)
	90	{this->fHistDeltaPhi = aHistDeltaPhi; }
	91	TH1F* GetHistDeltaPhi2() {return this->fHistDeltaPhi2; }
	92	void SetHistDeltaPhi2(TH1F* aHistDeltaPhi2)
	93	{this->fHistDeltaPhi2 = aHistDeltaPhi2; }
	94	TH1F* GetHistDeltaPhihere() {return this->fHistDeltaPhihere; }
	95	void SetHistDeltaPhihere(TH1F* aHistDeltaPhihere)
	96	{this->fHistDeltaPhihere = aHistDeltaPhihere; }
	97	TH1F* GetHistPhiEP() {return this->fHistPhiEP; }
	98	void SetHistPhiEP(TH1F* aHistPhiEP)
	99	{this->fHistPhiEP = aHistPhiEP; }
	100	TH1F* GetHistPhiEPhere() {return this->fHistPhiEPhere; }
	101	void SetHistPhiEPhere(TH1F* aHistPhiEPhere)
	102	{this->fHistPhiEPhere = aHistPhiEPhere; }
	103	TH1F* GetHistPhiLYZ() {return this->fHistPhiLYZ; }
	104	void SetHistPhiLYZ(TH1F* aHistPhiLYZ)
	105	{this->fHistPhiLYZ = aHistPhiLYZ; }
	106	TH1F* GetHistPhiLYZ2() {return this->fHistPhiLYZ2;}
	107	void SetHistPhiLYZ2(TH1F* aHistPhiLYZ2)
	108	{this->fHistPhiLYZ2 = aHistPhiLYZ2; }
f1d945a1	109
9d062fe3	110	//input
	111	void SetSecondRunList(TList* list) { this->fSecondRunList = list; }
	112	TList* GetSecondRunList() { return this->fSecondRunList; }
f1d945a1	113
	114	private:
	115
882ffd6a	116	AliFlowAnalysisWithLYZEventPlane(const AliFlowAnalysisWithLYZEventPlane& aAnalysis); // copy constructor
882ffd6a	117	AliFlowAnalysisWithLYZEventPlane& operator=(const AliFlowAnalysisWithLYZEventPlane& aAnalysis); // assignment operator
af795c87	118
f1d945a1	119	//histograms
28ca24ad	120	TList* fHistList; //list ro hold all histograms
9d062fe3	121	TList* fSecondRunList; //list from Second LYZ run output
80e93fef	122
f1d945a1	123	//input
882ffd6a	124	TProfile* fSecondReDtheta; // input profile
	125	TProfile* fSecondImDtheta; // input profile
	126	TProfile* fFirstr0theta; // input profile
80e93fef	127
f1d945a1	128	//output
80e93fef	129	TProfile* fHistProVetaRP; //
	130	TProfile* fHistProVetaPOI; //
	131	TProfile* fHistProVPtRP; //
	132	TProfile* fHistProVPtPOI; //
882ffd6a	133	TProfile* fHistProWr; //
882ffd6a	134	TProfile* fHistProWrCorr; //
9d062fe3	135	TH1F* fHistQsumforChi; //
882ffd6a	136	TH1F* fHistDeltaPhi; //
	137	TH1F* fHistDeltaPhi2; //
	138	TH1F* fHistDeltaPhihere; //
	139	TH1F* fHistPhiEP; //
	140	TH1F* fHistPhiEPhere; //
	141	TH1F* fHistPhiLYZ; //
	142	TH1F* fHistPhiLYZ2; //
f1d945a1	143
882ffd6a	144	AliFlowCommonHist* fCommonHists; //
882ffd6a	145	AliFlowCommonHistResults* fCommonHistsRes; //
f1d945a1	146
882ffd6a	147	Int_t fEventNumber; // event counter
af795c87	148
28ca24ad	149	TVector2 *fQsum; // flow vector sum
28ca24ad	150	Double_t fQ2sum; // flow vector sum squared
882ffd6a	151
f1d945a1	152
9d062fe3	153	ClassDef(AliFlowAnalysisWithLYZEventPlane, 1); // lyz analysis
f1d945a1	154	};
	155
	156	#endif
882ffd6a	157