AliFMDEventPlaneFinder.h

Go to the documentation of this file.

// This task finds the eventplane
// using the FMD
// 
#ifndef ALIFMDEVENTPLANEFINDER_H 
#define ALIFMDEVENTPLANEFINDER_H 

#include <TNamed.h>
#include <TVector2.h>
#include "AliForwardUtil.h"
class AliVEvent;
class TH1D;
class TH2F;
class TH2D;
class TString;
class AliOADBContainer;
class AliAODForwardEP;

class AliFMDEventPlaneFinder : public TNamed
{
public:
  AliFMDEventPlaneFinder();
  AliFMDEventPlaneFinder(const char* name);
  AliFMDEventPlaneFinder(const AliFMDEventPlaneFinder& o);
  virtual ~AliFMDEventPlaneFinder();
  AliFMDEventPlaneFinder& operator=(const AliFMDEventPlaneFinder& o);
  virtual void SetupForData(const TAxis& etaAxis);
  Bool_t FindEventplane(AliVEvent* esd,
                        AliAODForwardEP& aodEp,
                        TH2D* h,
                        AliForwardUtil::Histos* hists);
  virtual void CreateOutputObjects(TList* dir);
  void Print(Option_t* option="") const;
  void SetDebug(Int_t dbg=1) { fDebug = dbg; }
  void CalcQVectors(TH2D* h, TH1D* eHist);
  Double_t CalcEventplane(const TVector2& v) const;
  void SetRunNumber(Int_t run);
  Int_t GetRunNumber() { return fRunNumber; }
  void GetPhiDist();
  void SetUsePhiWeights(Bool_t use = kTRUE) { fUsePhiWeights = use; }
  void FillHists(AliAODForwardEP* fmdEP);
  void SetOADBPath(Char_t* fname) { fOADBFileName = fname; }

protected:
  Double_t GetPhiWeight(Int_t etaBin, Int_t phiBin) const;
  Double_t CalcDifference(Double_t a1, Double_t a2) const;
  TH1D* MakePsiRHist(const char* name, 
         const char* title,
         Int_t color);
  TH1D* MakeDiffHist(const char* name, 
         const char* first,
         const char* second,
         Int_t color);
  TH2F* MakeCorrHist(const char* name, 
         const char* first,
         const char* second);
  TList*            fList;               // List for diag. hists.
  AliVEvent*        fEvent;              // Current event
  TVector2          fQt;                 // Q vector for total ep
  TVector2          fQa;                 // Q vector for sub-ep A
  TVector2          fQc;                 // Q vector for sub-ep C
  TVector2          fQ1;                 // Q vector for sub-ep 1
  TVector2          fQ2;                 // Q vector for sub-ep 2
  TVector2          fQeta;               // Q vector for psi eta-dependence
  TH1D*             fHepFMD;             // Diagnostics histogram
  TH1D*             fHepFMDA;            // Diagnostics histogram
  TH1D*             fHepFMDC;            // Diagnostics histogram
  TH1D*             fHepFMDQC1;          // Diagnostics histogram
  TH1D*             fHepFMDQC2;          // Diagnostics histogram
  TH1D*             fHdiffFMDAC;         // Diagnostics histogram
  TH1D*             fHdiffFMDTPC;        // Diagnostics histogram
  TH1D*             fHdiffFMDVZERO;      // Diagnostics histogram
  TH2F*             fHcorrFMDAC;         // Diagnostics histogram
  TH2F*             fHcorrFMDTPC;        // Diagnostics histogram
  TH2F*             fHcorrFMDVZERO;      // Diagnostics histogram
  TH2D*             fHPhi;               // Diagnostics histogram
  Int_t             fDebug;              // Debug flag
  TString           fOADBFileName;       // Path to OADB container
  AliOADBContainer* fOADBContainer;      // OADBContainer object
  TH2D*             fPhiDist;            // Phi dist. for phi weights
  Int_t             fRunNumber;          // Run number supplied
  Bool_t            fUsePhiWeights;      // Flag for phi weights

  ClassDef(AliFMDEventPlaneFinder,2); //  
};

#endif
// Local Variables:
//   mode: C++
// End: