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

#ifndef AliAnaVZEROEPFlatenning_cxx
#define AliAnaVZEROEPFlatenning_cxx

class AliESDEvent;

#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) const;

 private:
  AliESDEvent *fESD;    //! ESD object
  TList       *fOutputList; //! Output list

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

  TProfile *fX2[8]; //! Profile histogram for Q^2_x
  TProfile *fY2[8]; //! Profile histogram for Q^2_y
  TProfile *fX2Y2[8]; //! Profile histogram for Q^2_x*Q^2_y
  TProfile *fCos8Psi[8]; //! Profile histogram for Cos(8*Psi)

  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)

  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
	4	class AliESDEvent;
	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
	24	Double_t CalculateVZEROEventPlane(const AliVEvent * event, Int_t ring, Float_t centrality) const;
	25
	26	private:
	27	AliESDEvent *fESD; //! ESD object
	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
	33	TProfile *fX2[8]; //! Profile histogram for Q^2_x
	34	TProfile *fY2[8]; //! Profile histogram for Q^2_y
	35	TProfile fX2Y2[8]; //! Profile histogram for Q^2_xQ^2_y
	36	TProfile fCos8Psi[8]; //! Profile histogram for Cos(8Psi)
	37
	38	TProfile *fX2In[8]; // Profile histogram for Q^2_x (read from input file)
	39	TProfile *fY2In[8]; // Profile histogram for Q^2_y (read from input file)
	40	TProfile fX2Y2In[8]; // Profile histogram for Q^2_xQ^2_y (read from input file)
	41	TProfile fCos8PsiIn[8]; // Profile histogram for Cos(8Psi) (read from input file)
	42
	43	TH2F *fPsiRingRawCentr[8]; //! Raw VZERO event plane vs centrality (ring-by-ring)
	44	TH2F *fPsiRingFlatCentr[8]; //! Flatenned with corrections on cumulants VZERO event plane vs centrality (ring-by-ring)
	45	TH2F *fPsiRingFlatFinalCentr[8]; //! Flatenned with corrections on cumulants and fourier VZERO event plane vs centrality (ring-by-ring)
	46
	47	AliAnaVZEROEPFlatenning(const AliAnaVZEROEPFlatenning&); // not implemented
	48	AliAnaVZEROEPFlatenning& operator=(const AliAnaVZEROEPFlatenning&); // not implemented
	49
	50	ClassDef(AliAnaVZEROEPFlatenning, 1); // VZERO analysis task for extraction EP flatenning params
	51	};
	52
	53	#endif