Code used for pi0 Raa and v2 analysis in LHC10h

[u/mrichter/AliRoot.git] / PWGGA / PHOSTasks / PHOS_PbPb / AliAnalysisTaskPi0Flow.h

#ifndef AliAnalysisTaskPi0Flow_cxx\r
#define AliAnalysisTaskPi0Flow_cxx\r
\r
// example of an analysis task creating a p_t spectrum\r
// Authors: Panos Cristakoglou, Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing\r
\r
class TObjArray;\r
class TH1F;\r
class TH2I;\r
class TH2F;\r
class TH3F;\r
class TF1 ;\r
class AliStack ;\r
class AliESDtrackCuts;\r
class AliPHOSGeometry;\r
class AliESDEvent ;\r
class AliPHOSCalibData;\r
class AliESDtrack ;\r
class AliESDCaloCluster ;\r
\r
#include "AliAnalysisTaskSE.h"\r
\r
class AliAnalysisTaskPi0Flow : public AliAnalysisTaskSE {\r
public:\r
  AliAnalysisTaskPi0Flow(const char *name = "AliAnalysisTaskPi0Flow");\r
  virtual ~AliAnalysisTaskPi0Flow() {}\r
  \r
  virtual void   UserCreateOutputObjects();\r
  virtual void   UserExec(Option_t *option);\r
  virtual void   Terminate(Option_t *);\r
  void SetPHOSBadMap(Int_t mod,TH2I * h)\r
  {\r
    if(fPHOSBadMap[mod]) delete fPHOSBadMap[mod] ;\r
    fPHOSBadMap[mod]=new TH2I(*h) ;\r
    printf("Set %s \n",fPHOSBadMap[mod]->GetName());\r
  }\r
  \r
private:\r
  AliAnalysisTaskPi0Flow(const AliAnalysisTaskPi0Flow&); // not implemented\r
  AliAnalysisTaskPi0Flow& operator=(const AliAnalysisTaskPi0Flow&); // not implemented\r
  Bool_t IsGoodChannel(const char * det, Int_t mod,Int_t ix, Int_t iz); //Use addisional bad map for PHOS\r
\r
  void FillHistogram(const char * key,Double_t x) const ; //Fill 1D histogram witn name key\r
  void FillHistogram(const char * key,Double_t x, Double_t y) const ; //Fill 2D histogram witn name key\r
  void FillHistogram(const char * key,Double_t x, Double_t y, Double_t z) const ; //Fill 3D histogram witn name key\r
\r
  Bool_t TestLambda(Double_t pt,Double_t l1,Double_t l2) ;  //Evaluate Dispersion cuts for photons\r
  Bool_t TestLambda2(Double_t pt,Double_t l1,Double_t l2) ;  //Evaluate Dispersion cuts for photons\r
  Double_t TestCPV(Double_t dx, Double_t dz, Double_t pt, Int_t charge);\r
  Int_t ConvertRunNumber(Int_t run) ; \r
\r
  void OpenInfoCalbration(Int_t run); //V0 calibration\r
  void EvalV0ReactionPlane(AliESDEvent * event) ;\r
  Double_t ApplyFlattening(Double_t phi, Double_t c) ; //Apply centrality-dependent flattening\r
  Double_t ApplyFlatteningV0A(Double_t phi, Double_t c) ; //Apply centrality-dependent flattening\r
  Double_t ApplyFlatteningV0C(Double_t phi, Double_t c) ; //Apply centrality-dependent flattening\r
\r
  Double_t CoreEnergy(AliESDCaloCluster * clu); \r
  void    Reclusterize(AliESDCaloCluster * clu) ;\r
  Bool_t  AreNeibors(Int_t id1,Int_t id2) ;\r
private:\r
  AliESDtrackCuts *fESDtrackCuts; // Track cut\r
  AliStack * fStack ;\r
  TList * fOutputContainer;        //final histogram container\r
  TList * fPHOSEvents[1][10][11] ; //Containers for events with PHOS photons\r
  TClonesArray * fPHOSEvent ;      //PHOS photons in current event\r
  AliPHOSCalibData *fPHOSCalibData; // PHOS calibration object\r
  TF1 *fNonLinCorr;          // Non-linearity correction\r
 \r
  //Reaction plain for v2\r
  Float_t fRP ; //!Reaction plane calculated with full TPC \r
  Float_t fRPV0A ;    //!Reaction plain calculated with A-side TPC: eta>0.15 \r
  Float_t fRPV0C ;    //!Reaction plain calculated with C-side TPC: eta<-0.15\r
  Bool_t fHaveTPCRP ; //! Is TPC RP defined?\r
  \r
  //V0 calibration\r
\r
  static const Int_t nCentrBinV0 = 9; // # cenrality bins\r
\r
  TProfile *fMultV0;                  // object containing VZERO calibration information\r
\r
  Float_t fV0Cpol,fV0Apol;            // loaded by OADB\r
\r
  Float_t fMeanQ[nCentrBinV0][2][2];    // and recentering\r
\r
  Float_t fWidthQ[nCentrBinV0][2][2];   // ...\r
\r
\r
//  TF1 * fRecent[5][12] ;//Recentering corrections\r
\r
  Int_t fRunNumber ;    //Current run number\r
  Float_t fCentrality ; //!Centrality of the currecnt event\r
\r
  Int_t fCenBin ;       //! Current centrality bin\r
\r
  TH2I *fPHOSBadMap[6] ;    //Container for PHOS bad channels map\r
\r
  AliPHOSGeometry  *fPHOSGeo;  //! PHOS geometry\r
  Int_t fEventCounter;         // number of analyzed events\r
\r
  ClassDef(AliAnalysisTaskPi0Flow, 1); // PHOS analysis task\r
};\r
\r
#endif\r