fixing compilation problem

[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaNeutralMeson.h

#ifndef ALIANANEUTRALMESON_H\r
#define ALIANANEUTRALMESON_H\r
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
 * See cxx source for full Copyright notice     */\r
/* $Id: $ */\r
\r
//_________________________________________________________________________\r
//a general class to fill two-photon and pi0+gamma invariant mass hisograms \r
// to be used to extract pi0, eta and omega raw yield.\r
// also for PHOS acceptance \r
//--\r
//-- by Renzhuo Wan (Iopp-wuhan) May 13,2009\r
//_________________________________________________________________________\r
\r
//Root\r
class TList;\r
class TH3F ;\r
class TH2F ;\r
class TLorentzVector;\r
//Analysis\r
class AliAODEvent ;\r
class AliESDEvent ;\r
#include "AliAnaPartCorrBaseClass.h"\r
class TParticle;\r
\r
#ifdef __PHOSUTIL__\r
  class AliPHOSGeoUtils;\r
#endif\r

#ifdef __EMCALUTIL__\r
  class AliEMCALGeoUtils;\r
#endif\r

\r
class AliAnaNeutralMeson : public AliAnaPartCorrBaseClass {\r
  \r
  public: \r
  \r
  AliAnaNeutralMeson() ; // default ctor\r
  AliAnaNeutralMeson(const char *name) ; // default ctor\r
  AliAnaNeutralMeson(const AliAnaNeutralMeson & g) ; // cpy ctor\r
  AliAnaNeutralMeson & operator = (const AliAnaNeutralMeson & api0) ;//cpy assignment\r
  virtual ~AliAnaNeutralMeson() ;//virtual dtor\r
  \r
  //......\r
  TList * GetCreateOutputObjects(); \r
  void Print(const Option_t * opt) const;\r
  \r
  void Init();\r
  void InitParameters();\r
  void MakeAnalysisFillHistograms();\r
  //void Terminate(TList * outList);\r
  //......\r
  void Get2GammaAsyDistPtM(AliAODPWG4Particle *p1, AliAODPWG4Particle *p2, Double_t &asy, Double_t &dist, \r
                         Double_t &pt, Double_t &mass); //get dist,asy, pair pt, pairmass btween two ptc\r
  void GetAsyDistPtM(TLorentzVector photon1, TLorentzVector photon2, Double_t &asy, Double_t &dist, \r
                   Double_t &pt, Double_t &mass); //get dist,asy, pair pt, pairmass btween two ptc\r
  Bool_t IsPhotonSelected(AliAODPWG4Particle *p1, Int_t ipid, Int_t nmod); //select the photon based on pid and nmod to be analyzed\r
\r
  void RealPi0Eta(TClonesArray * array) ; //perform the 2gamma IVM analysis for one event to extract pi0 and eta\r
  void MixPi0Eta(TClonesArray * array1, TClonesArray *array2) ; //perform the 2gamma IVM analysis for mixed eventsto reconstruct the background\r
  void RealOmega(TClonesArray *array) ; //perform the pi0+gamma->3gamma IVM analysis for one event to extract omega\r
  void MixOmegaAB(TClonesArray * array1, TClonesArray * array2); ////perform the 2gamma IVM analysis for mixed events to reconstruct the bacground (r1_event1+r2_event1)+r3_event2 and (r1_event1+r2_event2)+r3_event2\r
  void MixOmegaC(TClonesArray * array1, TClonesArray * array2, TClonesArray * array3); //perform the 2gamma IVM analysis for mixed events to reconstruct the bacground (r1_event1+r2_event2)+r3_event3\r
      \r
  void PhotonAcceptance(AliStack *stack); //to get MC photon pt and photon can be hit on dectector\r
  void Pi0EtaAcceptance(AliStack *stack); //to get MC pi0 and eta pt and photon can be hit on dectector\r
  void OmegaAcceptance(AliStack * stack); //to get MC omega pt and photon can be hit on dectector\r
      \r
  void GetMCDistAsy(TParticle * p1, TParticle * p2, Double_t & dist, Double_t &asy); //to get MC information between two particle\r
  \r
  void SetNAsyBinsMinMax(Int_t bins, Double_t min, Double_t max) {fNbinsAsy=bins; fMinAsy=min; fMaxAsy=max; }   //set Asy bins, min and max\r
  void SetNPtBinsMinMax(Int_t bins, Double_t min, Double_t max) {fNbinsPt=bins; fMinPt=min; fMaxPt=max; } //set pt bins, min and max\r
  void SetNMassBinsMinMas(Int_t bins, Double_t min, Double_t max) {fNbinsM=bins; fMinM=min; fMaxM=max; } //set mass pt bins, min and max\r
  void SetNEventsMixed(Int_t nevents) { fNmaxMixEv=nevents;} //events to be mixed \r
  \r
  TString GetCalorimeter()   const { return fCalorimeter ; } //PHOS or EMCAL\r
  void SetCalorimeter(TString det);// { fCalorimeter = det ; } \r
 \r
  void SetNPID(Int_t pid) {fNPID=pid;}  //number of PID cut\r
  void SetInvMassCut(Int_t mass) {fInvMassCut = mass;}  // \r
  void SetPi0MassPeakWidthCut(Double_t cut) {fPi0MassPeakWidthCut = cut;}\r
 \r
  Bool_t SetAnaPi0Eta(Bool_t ana)  {fAnaPi0Eta = ana; return fAnaPi0Eta;} //analysis pi0 and eta ????\r
  Bool_t SetAnaOmega(Bool_t ana)   {fAnaOmega = ana; return fAnaOmega ;}  //analysis omega ????\r
 \r
  Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list\r
  \r
  private:\r
\r
  Bool_t fAnaPi0Eta; //whether ana the 2gamma IVM to extract pi0 and eta spectrum\r
  Bool_t fAnaOmega;  //whether ana the 3gamma IVM to extract omega spectrum\r
\r
  Int_t fNPID;       //PID bin\r
  Int_t fNmaxMixEv ; //number events to be mixed\r
  Int_t fNAsy;       //number of Asy cut\r
  Int_t fNMod;       //number of PHOS or EMCAL modules to be analyzed\r
  Int_t fNHistos;    //number of histograms\r
  Double_t fPerPtBin;//plot the IVM distribution per pt bin\r
 \r
  Double_t *fAsyCut; // asy cut(< 0.7, 0.8, 0.9, 1)\r
  Int_t    *fModCut; // mod cut(1,3,5) for PHOS, and 0 (default) for EMCAL\r
\r
  Double_t fDist; //distance between two particle in (eta, phi) plane\r
  Double_t fAsy;  //asysmetry |(e1-e2)|/(e1+e2)\r
  Double_t fPt;   //pair pt\r
  Double_t fMass; //pair mass\r
  Double_t fInvMassCut; // cut of invmass\r
\r
  Int_t fNbinsAsy; //Asy bin number, min and max\r
  Double_t fMinAsy;\r
  Double_t fMaxAsy;\r
  Int_t fNbinsPt;  //Pt bin number, min and max\r
  Double_t fMinPt;\r
  Double_t fMaxPt;\r
  Int_t fNbinsM;  //mass bin number, min and max\r
  Double_t fMinM;\r
  Double_t fMaxM;\r
\r
  //! containers for photons in stored events\r
  TList * fEventsListPi0Eta ; //! events with cluster larger than 2 for pi0 and eta analysis\r
  TList * fEventsListOmega ; //! events with cluster larger than 3 for omega analysis\r
\r
  Double_t fPi0Mass;    //nominal pi0 mass\r
  Double_t fPi0MassPeakWidthCut; //pi0 candidate\r
  \r
  TH1I *  fhNClusters;  //cluster multiplicity distribution per event\r
  TH1D ** fhRecPhoton; //reconstructed photon pt distribution\r
  TH2F ** fhRecPhotonEtaPhi; //photon (eta, phi) distruction\r
\r
  TH1F ** fReal2Gamma; //real 2gamma IVM in per pt bin (ipid, imod, iasy)\r
  TH1F ** fMix2Gamma;  //mix 2gamma  IVM in per pt bin (ipid, imod, iasy, pt)\r
  TH1F ** fRealOmega;  //real omega  IVM in per pt bin (ipid, imod, iasy, pt)\r
  TH1F ** fMixOmegaA;  //mixA omega  IVM in per pt bin (ipid, imod, iasy, pt)  (r1_event1+r2_event1)+r3_event2\r
  TH1F ** fMixOmegaB;  //mixB omega  IVM in per pt bin (ipid, imod, iasy, pt)  (r1_event1+r2_event2)+r3_event2\r
  TH1F ** fMixOmegaC;  //mixC omega  IVM in per pt bin (ipid, imod, iasy, pt)  (r1_event1+r2_event2)+r3_event3 \r
\r
  TH3F ** fRealTwoGammaAsyPtM; //real 2gamma IVM(asy, pt, m) without PID, mod cut\r
  TH3F ** fMixTwoGammaAsyPtM;  //mix  2gamma IVM(asy, pt, m) without PID, mod cut\r
  TH3F ** fRealPi0GammaAsyPtM; //real omega  IVM(asy, pt, m) without PID, mod cut\r
  TH3F ** fMixAPi0GammaAsyPtM; //mixA omega  IVM(asy, pt, m) without PID, mod cut\r
  TH3F ** fMixBPi0GammaAsyPtM; //mixB omega  IVM(asy, pt, m) without PID, mod cut\r
  TH3F ** fMixCPi0GammaAsyPtM; //mixC omega  IVM(asy, pt, m) without PID, mod cut\r
\r
  //Histograms filled only if MC data is requested       \r
  //pt distribution\r
  //currently only PHOS included \r
  TH1D * fhPrimPhotonPt ;    //! Spectrum of Primary \r
  TH1D * fhPrimPhotonAccPt ; //! Spectrum of primary with accepted daughters \r
  TH1D * fhPrimPhotonY ;     //! Rapidity distribution of primary particles\r
  TH1D * fhPrimPhotonAccY ;  //! Rapidity distribution of primary with accepted daughters\r
  TH1D * fhPrimPhotonPhi ;   //! Azimutal distribution of primary particles\r
  TH1D * fhPrimPhotonAccPhi; //! Azimutal distribution of primary with accepted daughters \r
 \r
  TH1D * fhPrimPi0Pt ;    //! Spectrum of Primary \r
  TH1D * fhPrimPi0AccPt ; //! Spectrum of primary with accepted daughters \r
  TH1D * fhPrimPi0Y ;     //! Rapidity distribution of primary particles\r
  TH1D * fhPrimPi0AccY ;  //! Rapidity distribution of primary with accepted daughters\r
  TH1D * fhPrimPi0Phi ;   //! Azimutal distribution of primary particles\r
  TH1D * fhPrimPi0AccPhi; //! Azimutal distribution of primary with accepted daughters  \r
 \r
  TH1D * fhPrimEtaPt ;    //! Spectrum of Primary \r
  TH1D * fhPrimEtaAccPt ; //! Spectrum of primary with accepted daughters \r
  TH1D * fhPrimEtaY ;     //! Rapidity distribution of primary particles\r
  TH1D * fhPrimEtaAccY ;  //! Rapidity distribution of primary with accepted daughters\r
  TH1D * fhPrimEtaPhi ;   //! Azimutal distribution of primary particles\r
  TH1D * fhPrimEtaAccPhi; //! Azimutal distribution of primary with accepted daughters\r
 \r
  TH1D * fhPrimOmegaPt ;    //! Spectrum of Primary \r
  TH1D * fhPrimOmegaAccPt ; //! Spectrum of primary with accepted daughters \r
  TH1D * fhPrimOmegaY ;     //! Rapidity distribution of primary particles\r
  TH1D * fhPrimOmegaAccY ;  //! Rapidity distribution of primary with accepted daughters\r
  TH1D * fhPrimOmegaPhi ;   //! Azimutal distribution of primary particles\r
  TH1D * fhPrimOmegaAccPhi; //! Azimutal distribution of primary with accepted daughters\r
\r
  TString fCalorimeter ;     //Select Calorimeter       \r
    \r
#ifdef __PHOSUTIL__\r
  AliPHOSGeoUtils * fPHOSGeo ; //! PHOS geometry pointer\r
#endif	\r
#ifdef __EMCALUTIL__\r
  AliEMCALGeoUtils * fEMCALGeo ; //! EMCAL geometry pointer\r
#endif	\r
  \r
  ClassDef(AliAnaNeutralMeson,1)\r
} ;\r
\r
\r
#endif //ALIANANEUTRALMESON\r
\r
\r
\r