PWG4/PartCorrDep/AliAnaPi0.h

   1 #ifndef ALIANAPI0_H
   2 #define ALIANAPI0_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice     */
   5 /* $Id: $ */
   6
   7 //_________________________________________________________________________
   8 // Class to fill two-photon invariant mass hisograms
   9 // to be used to extract pi0 raw yield.
  10 //
  11 //-- Author: Dmitri Peressounko (RRC "KI")
  12 //-- Adapted to PartCorr frame by Lamia Benhabib (SUBATECH)
  13 //-- and Gustavo Conesa (INFN-Frascati)
  14
  15 //Root
  16 class TList;
  17 class TH3D ;
  18
  19 //Analysis
  20 class AliAODEvent ;
  21 class AliESDEvent ;
  22 #include "AliAnaPartCorrBaseClass.h"
  23
  24 class AliAnaPi0 : public AliAnaPartCorrBaseClass {
  25
  26   public:
  27
  28   AliAnaPi0() ; // default ctor
  29   AliAnaPi0(const char *name) ; // default ctor
  30   AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
  31   AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
  32   virtual ~AliAnaPi0() ;//virtual dtor
  33
  34   TList * GetCreateOutputObjects();
  35
  36   void Print(const Option_t * opt) const;
  37
  38   //void Init();
  39   void InitParameters();
  40
  41   //void MakeAnalysisFillAOD() {;} //Not needed
  42   void MakeAnalysisFillHistograms();
  43
  44   //    void SetBadRunsList(){;} ;     //Set list of runs which can be used for this analysis
  45   //To be defined in future.
  46
  47   //void SetEtalonHisto(TH3D * h);//Provide etalon of binning for histograms
  48
  49   //Setters for parameters of event buffers
  50   void SetNCentrBin(Int_t n=5){fNCentrBin=n ;} //number of bins in centrality
  51   void SetNZvertBin(Int_t n=5){fNZvertBin=n ;} //number of bins for vertex position
  52   void SetNRPBin(Int_t n=6)   {fNrpBin=n ;}    //number of bins in reaction plain
  53   void SetNMaxEvMix(Int_t n=20){fNmaxMixEv=n ;}//Maximal number of events for mixing
  54
  55   //Setters for event selection
  56   void SetZvertexCut(Float_t zcut=40.){fZvtxCut=zcut ;} //cut on vertex position
  57
  58   TString GetCalorimeter()   const {return fCalorimeter ; }
  59   void SetCalorimeter(TString det)    {fCalorimeter = det ; }
  60
  61   void Terminate(TList* outputList);
  62   void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
  63
  64   Int_t GetModuleNumber(AliAODPWG4Particle * particle);
  65   void SetNumberOfModules(Int_t nmod) {fNModules = nmod;}
  66
  67   Int_t GetNPID()   const {return fNPID ; }
  68   void SetNPID(Int_t n)    {fNPID = n ; }
  69
  70   private:
  71   Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list
  72
  73   private:
  74   Int_t fNCentrBin ;      // Number of bins in event container for centrality
  75   Int_t fNZvertBin ;      // Number of bins in event container for vertex position
  76   Int_t fNrpBin ;                 // Number of bins in event container for reaction plain
  77   Int_t fNPID ;                   // Number of possible PID combinations
  78   Int_t fNmaxMixEv ;      // Maximal number of events stored in buffer for mixing
  79   Float_t fZvtxCut ;      // Cut on vertex position
  80   TString fCalorimeter ;  // Select Calorimeter for IM
  81   Int_t fNModules ;       // Number of EMCAL/PHOS modules, set as many histogras as modules
  82
  83   TList ** fEventsList ;  //! containers for photons in stored events
  84
  85   //Histograms
  86
  87   //TH3D * fhEtalon ; //Etalon histo, all distributions will have same binning as this one
  88
  89   TH3D ** fhReMod ;  //!REAL two-photon invariant mass distribution for different calorimeter modules.
  90
  91   TH3D ** fhRe1 ;  //!REAL two-photon invariant mass distribution for different centralities and PID
  92   TH3D ** fhMi1 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  93   TH3D ** fhRe2 ;  //!REAL two-photon invariant mass distribution for different centralities and PID
  94   TH3D ** fhMi2 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  95   TH3D ** fhRe3 ;  //!REAL two-photon invariant mass distribution for different centralities and PID
  96   TH3D ** fhMi3 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  97   TH3D * fhEvents;  //!Number of events per centrality, RP, zbin
  98
  99   //Acceptance
 100   TH1D * fhPrimPt ;    //! Spectrum of Primary
 101   TH1D * fhPrimAccPt ; //! Spectrum of primary with accepted daughters
 102   TH1D * fhPrimY ;     //! Rapidity distribution of primary particles
 103   TH1D * fhPrimAccY ;  //! Rapidity distribution of primary with accepted daughters
 104   TH1D * fhPrimPhi ;   //! Azimutal distribution of primary particles
 105   TH1D * fhPrimAccPhi; //! Azimutal distribution of primary with accepted daughters
 106
 107   ClassDef(AliAnaPi0,7)
 108 } ;
 109
 110
 111 #endif //ALIANAPI0_H
 112
 113
 114