[u/mrichter/AliRoot.git] / PWGPP / VZERO / AliAnaVZEROEPFlatenning.h

#ifndef AliAnaVZEROEPFlatenning_cxx
#define AliAnaVZEROEPFlatenning_cxx

class AliVEvent;

#include "AliAnalysisTaskSE.h"

class AliAnaVZEROEPFlatenning : public AliAnalysisTaskSE {
 public:
  AliAnaVZEROEPFlatenning();
  AliAnaVZEROEPFlatenning(const char *name);
  virtual ~AliAnaVZEROEPFlatenning() {}
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);

  virtual void Init();

  void SetMBTrigName(const char *name = "CPBI") {fMBTrigName = name;}
  void SetUsePhysSel(Bool_t usePhysSel) {fUsePhysSel = usePhysSel;}
  void SetInput(const char *filename);

  Double_t CalculateVZEROEventPlane(const AliVEvent *  event, Int_t ring, Float_t centrality, Double_t &qxTierce, Double_t &qyTierce) const;

 private:
  AliVEvent   *fEvent;    //! ESD ot AOD object
  TList       *fOutputList; //! Output list

  TString fMBTrigName; // MB trigger name (for evt sel)
  Bool_t  fUsePhysSel; // Use or not phys sel

  TProfile *fX2[11]; //! Profile histogram for Q^2_x
  TProfile *fY2[11]; //! Profile histogram for Q^2_y
  TProfile *fX2Y2[11]; //! Profile histogram for Q^2_x*Q^2_y
  TProfile *fCos8Psi[11]; //! Profile histogram for Cos(8*Psi)
  TProfile *fC2[8]; //! Profile histogram for Cos(2*phi)
  TProfile *fS2[8]; //! Profile histogram for Sin(2*phi)
  TProfile *fC4[8]; //! Profile histogram for Cos(4*phi)
  TProfile *fS4[8]; //! Profile histogram for Sin(4*phi)
  TProfile *fX2Corr[11]; //! Profile histogram for Q^2_x
  TProfile *fY2Corr[11]; //! Profile histogram for Q^2_y
  TProfile *fX2Y2Corr[11]; //! Profile histogram for Q^2_x*Q^2_y

  TProfile *fX2In[8]; // Profile histogram for Q^2_x (read from input file)
  TProfile *fY2In[8]; // Profile histogram for Q^2_y (read from input file)
  TProfile *fX2Y2In[8]; // Profile histogram for Q^2_x*Q^2_y (read from input file)
  TProfile *fCos8PsiIn[8]; // Profile histogram for Cos(8*Psi) (read from input file)

  TH2F *fPsiRingRawCentr[8]; //! Raw VZERO event plane vs centrality (ring-by-ring)
  TH2F *fPsiRingFlatCentr[8]; //! Flatenned with corrections on cumulants VZERO event plane vs centrality (ring-by-ring)
  TH2F *fPsiRingFlatFinalCentr[8]; //! Flatenned with corrections on cumulants and fourier VZERO event plane vs centrality (ring-by-ring)
  TH2F *fPsiARawCentr; //! Raw VZEROA event plane vs centrality
  TH2F *fPsiAFlatCentr; //! Flatenned with corrections on cumulants VZEROA event plane vs centrality 
  TH2F *fPsiCRawCentr; //! Raw VZEROC event plane vs centrality
  TH2F *fPsiCFlatCentr; //! Flatenned with corrections on cumulants VZEROC event plane vs centrality 
  TH2F *fPsiACRawCentr; //! Raw VZERO event plane vs centrality
  TH2F *fPsiACFlatCentr; //! Flatenned with corrections on cumulants VZERO event plane vs centrality 

  AliAnaVZEROEPFlatenning(const AliAnaVZEROEPFlatenning&); // not implemented
  AliAnaVZEROEPFlatenning& operator=(const AliAnaVZEROEPFlatenning&); // not implemented
  
  ClassDef(AliAnaVZEROEPFlatenning, 1); // VZERO analysis task for extraction EP flatenning params
};

#endif
Commit	Line	Data
345e2385	1	#ifndef AliAnaVZEROEPFlatenning_cxx
	2	#define AliAnaVZEROEPFlatenning_cxx
	3
6b4456fd	4	class AliVEvent;
345e2385	5
	6	#include "AliAnalysisTaskSE.h"
	7
	8	class AliAnaVZEROEPFlatenning : public AliAnalysisTaskSE {
	9	public:
	10	AliAnaVZEROEPFlatenning();
	11	AliAnaVZEROEPFlatenning(const char *name);
	12	virtual ~AliAnaVZEROEPFlatenning() {}
	13
	14	virtual void UserCreateOutputObjects();
	15	virtual void UserExec(Option_t *option);
	16	virtual void Terminate(Option_t *);
	17
	18	virtual void Init();
	19
	20	void SetMBTrigName(const char *name = "CPBI") {fMBTrigName = name;}
	21	void SetUsePhysSel(Bool_t usePhysSel) {fUsePhysSel = usePhysSel;}
	22	void SetInput(const char *filename);
	23
fe509094	24	Double_t CalculateVZEROEventPlane(const AliVEvent * event, Int_t ring, Float_t centrality, Double_t &qxTierce, Double_t &qyTierce) const;
345e2385	25
345e2385	26	private:
6b4456fd	27	AliVEvent *fEvent; //! ESD ot AOD object
345e2385	28	TList *fOutputList; //! Output list
	29
	30	TString fMBTrigName; // MB trigger name (for evt sel)
	31	Bool_t fUsePhysSel; // Use or not phys sel
	32
ed806549	33	TProfile *fX2[11]; //! Profile histogram for Q^2_x
	34	TProfile *fY2[11]; //! Profile histogram for Q^2_y
	35	TProfile fX2Y2[11]; //! Profile histogram for Q^2_xQ^2_y
	36	TProfile fCos8Psi[11]; //! Profile histogram for Cos(8Psi)
fe509094	37	TProfile fC2[8]; //! Profile histogram for Cos(2phi)
	38	TProfile fS2[8]; //! Profile histogram for Sin(2phi)
	39	TProfile fC4[8]; //! Profile histogram for Cos(4phi)
	40	TProfile fS4[8]; //! Profile histogram for Sin(4phi)
ed806549	41	TProfile *fX2Corr[11]; //! Profile histogram for Q^2_x
	42	TProfile *fY2Corr[11]; //! Profile histogram for Q^2_y
	43	TProfile fX2Y2Corr[11]; //! Profile histogram for Q^2_xQ^2_y
345e2385	44
	45	TProfile *fX2In[8]; // Profile histogram for Q^2_x (read from input file)
	46	TProfile *fY2In[8]; // Profile histogram for Q^2_y (read from input file)
	47	TProfile fX2Y2In[8]; // Profile histogram for Q^2_xQ^2_y (read from input file)
	48	TProfile fCos8PsiIn[8]; // Profile histogram for Cos(8Psi) (read from input file)
	49
	50	TH2F *fPsiRingRawCentr[8]; //! Raw VZERO event plane vs centrality (ring-by-ring)
	51	TH2F *fPsiRingFlatCentr[8]; //! Flatenned with corrections on cumulants VZERO event plane vs centrality (ring-by-ring)
	52	TH2F *fPsiRingFlatFinalCentr[8]; //! Flatenned with corrections on cumulants and fourier VZERO event plane vs centrality (ring-by-ring)
ed806549	53	TH2F *fPsiARawCentr; //! Raw VZEROA event plane vs centrality
	54	TH2F *fPsiAFlatCentr; //! Flatenned with corrections on cumulants VZEROA event plane vs centrality
	55	TH2F *fPsiCRawCentr; //! Raw VZEROC event plane vs centrality
	56	TH2F *fPsiCFlatCentr; //! Flatenned with corrections on cumulants VZEROC event plane vs centrality
	57	TH2F *fPsiACRawCentr; //! Raw VZERO event plane vs centrality
	58	TH2F *fPsiACFlatCentr; //! Flatenned with corrections on cumulants VZERO event plane vs centrality
fe509094	59
345e2385	60	AliAnaVZEROEPFlatenning(const AliAnaVZEROEPFlatenning&); // not implemented
	61	AliAnaVZEROEPFlatenning& operator=(const AliAnaVZEROEPFlatenning&); // not implemented
	62
	63	ClassDef(AliAnaVZEROEPFlatenning, 1); // VZERO analysis task for extraction EP flatenning params
	64	};
	65
	66	#endif