PWG4/PartCorrBase/AliAnaPartCorrBaseClass.h

   1 #ifndef ALIANAPARTCORRBASECLASS_H
   2 #define ALIANAPARTCORRBASECLASS_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 // Base class for analysis algorithms
   9 //-- Author: Gustavo Conesa (INFN-LNF)
  10
  11 #include <cstdlib>
  12
  13 //ROOT
  14 class TClonesArray ;
  15 class TObjArray ;
  16 #include <TList.h>
  17 #include <TObject.h>
  18
  19 //Analysis
  20 class AliAODCaloCluster;
  21 class AliAODCaloCells;
  22 class AliCaloTrackReader ;
  23 class AliCaloPID ;
  24 class AliFidutialCut ;
  25 class AliIsolationCut ;
  26 class AliMCAnalysisUtils ;
  27 class AliNeutralMesonSelection ;
  28 class AliStack ;
  29 class AliHeader ;
  30 class AliGenEventHeader ;
  31 #include "AliAODPWG4ParticleCorrelation.h"
  32
  33 class AliAnaPartCorrBaseClass : public TObject {
  34
  35 public:
  36
  37   AliAnaPartCorrBaseClass() ; // default ctor
  38   AliAnaPartCorrBaseClass(const AliAnaPartCorrBaseClass & g) ; // cpy ctor
  39   AliAnaPartCorrBaseClass & operator = (const AliAnaPartCorrBaseClass & g) ;//cpy assignment
  40   virtual ~AliAnaPartCorrBaseClass() ; //virtual dtor
  41
  42 //      virtual void AddAODCaloCluster(AliAODCaloCluster calo) ;
  43   virtual void AddAODParticle(AliAODPWG4Particle part) ;
  44
  45 //      virtual void ConnectAODCaloClusters();
  46   virtual void ConnectAODPHOSCells();
  47   virtual void ConnectAODEMCALCells();
  48   virtual void ConnectInputOutputAODBranches();
  49
  50   virtual TList * GetCreateOutputObjects() { return (new TList) ;}
  51
  52   virtual void AddToHistogramsName(TString add) { fAddToHistogramsName = add; }
  53   virtual TString GetAddedHistogramsStringToName() {return fAddToHistogramsName ;}
  54
  55   virtual void Init() {;}
  56   virtual void InitParameters() ;
  57
  58   virtual void Print(const Option_t * ) const ;
  59
  60   virtual void MakeAnalysisFillAOD()  {;}
  61
  62   virtual void MakeAnalysisFillHistograms() {;}
  63
  64   virtual Int_t GetDebug() const  { return fDebug ; }
  65   virtual void SetDebug(Int_t d)   { fDebug = d ; }
  66
  67   virtual Int_t GetEventNumber() const ;
  68
  69   virtual Int_t GetLabelShift() const  { return fLabelShift ; }
  70   virtual void SetLabelShift(Int_t shift)   { fLabelShift = shift ; }
  71
  72   virtual AliCaloTrackReader * GetReader() const {return fReader ; }
  73   virtual void SetReader(AliCaloTrackReader * reader) { fReader = reader ; }
  74
  75   virtual void Terminate(TList * /*outputList*/) {;}
  76
  77   //analysis AOD branch
  78   virtual TClonesArray * GetCreateOutputAODBranch() ;
  79   virtual TString GetInputAODName() const {return fInputAODName ; }
  80   virtual void SetInputAODName(TString name)   { fInputAODName = name; }
  81   virtual TString GetOutputAODName()  const {return fOutputAODName ; }
  82   virtual void SetOutputAODName(TString name)   { fNewAOD = kTRUE ; fOutputAODName = name; }
  83   virtual Bool_t NewOutputAOD() const {return fNewAOD;}
  84   virtual TString GetOutputAODClassName() const {return fOutputAODClassName;}
  85   virtual void SetOutputAODClassName(TString name) {fOutputAODClassName = name; }
  86
  87   virtual TString GetAODObjArrayName() const {return fAODObjArrayName;}
  88   virtual void SetAODObjArrayName(TString name) {fAODObjArrayName = name; }
  89
  90   virtual TClonesArray* GetInputAODBranch() const {return fInputAODBranch ;}
  91   virtual TClonesArray* GetOutputAODBranch() const {return fOutputAODBranch ;}
  92
  93 //      virtual TClonesArray* GetAODCaloClusters() const {return fAODCaloClusters ;}
  94   virtual TClonesArray* GetAODCaloClusters() const ;
  95   virtual TClonesArray* GetAODTracks() const ;
  96   virtual AliAODCaloCells* GetAODCaloCells() const {return fAODCaloCells ;}
  97
  98   virtual TObjArray* GetAODCTS() const ;
  99   virtual TObjArray* GetAODEMCAL() const ;
 100   virtual TObjArray* GetAODPHOS() const ;
 101
 102   virtual TString       GetBaseParametersList();
 103
 104   virtual TNamed * GetEMCALCells() const ;
 105   virtual TNamed * GetPHOSCells() const ;
 106
 107   virtual AliStack * GetMCStack() const ;
 108   virtual AliHeader* GetMCHeader() const ;
 109   virtual AliGenEventHeader* GetMCGenEventHeader() const ;
 110
 111   //Analysis helpers classes pointers setters and getters
 112   virtual AliCaloPID * GetCaloPID() const {return  fCaloPID ;}
 113   virtual void SetCaloPID(AliCaloPID * pid) { fCaloPID = pid ;}
 114
 115   virtual AliFidutialCut * GetFidutialCut() const {return  fFidCut ;}
 116   virtual void SetFidutialCut(AliFidutialCut * fc) { fFidCut = fc ;}
 117
 118   virtual AliIsolationCut * GetIsolationCut() const {return  fIC ;}
 119   virtual void SetIsolationCut(AliIsolationCut * fc) { fIC = fc ;}
 120
 121   virtual AliMCAnalysisUtils * GetMCAnalysisUtils() const {return  fMCUtils ;}
 122   virtual void SetMCAnalysisUtils(AliMCAnalysisUtils * mcutils) { fMCUtils = mcutils ;}
 123
 124   virtual AliNeutralMesonSelection * GetNeutralMesonSelection() const {return  fNMS ;}
 125   virtual void SetNeutralMesonSelection(AliNeutralMesonSelection * nms) { fNMS = nms ;}
 126
 127   virtual Bool_t     IsDataMC() const {return fDataMC ; }
 128   virtual void SwitchOnDataMC()    {fDataMC = kTRUE ; }
 129   virtual void SwitchOffDataMC()    {fDataMC = kFALSE ; }
 130
 131   virtual Bool_t IsFidutialCutOn() const {return fCheckFidCut ; }
 132   virtual void SwitchOnFidutialCut() { fCheckFidCut = kTRUE;}
 133   virtual void SwitchOffFidutialCut() { fCheckFidCut = kFALSE;}
 134
 135   virtual Bool_t IsCaloPIDOn() const {return fCheckCaloPID ; }
 136   virtual void SwitchOnCaloPID() { fCheckCaloPID = kTRUE;}
 137   virtual void SwitchOffCaloPID() { fCheckCaloPID = kFALSE;}
 138
 139   virtual Bool_t IsCaloPIDRecalculationOn() const {return fRecalculateCaloPID ; }
 140   virtual void SwitchOnCaloPIDRecalculation() { fRecalculateCaloPID  = kTRUE;}
 141   virtual void SwitchOffCaloPIDRecalculation() { fRecalculateCaloPID  = kFALSE;}
 142
 143   virtual Float_t    GetMaxPt()         const {return fMaxPt ; }
 144   virtual Float_t    GetMinPt()         const {return fMinPt ; }
 145   virtual void SetMaxPt(Float_t pt)              {fMaxPt = pt ; }
 146   virtual void SetMinPt(Float_t pt)              {fMinPt = pt ; }
 147   void SetPtCutRange(Double_t ptmin, Double_t ptmax)
 148   {  fMaxPt=ptmax;   fMinPt=ptmin;}
 149
 150   //Histogrammes setters and getters
 151   virtual void SetHistoPtRangeAndNBins(Float_t min, Float_t max, Int_t n) {
 152     fHistoNPtBins = n ;
 153     fHistoPtMax = max ;
 154     fHistoPtMin = min ;
 155   }
 156
 157   Int_t   GetHistoNPtBins() const { return fHistoNPtBins ; }
 158   Float_t GetHistoPtMin()   const { return fHistoPtMin ; }
 159   Float_t GetHistoPtMax()   const { return fHistoPtMax ; }
 160
 161   virtual void SetHistoPhiRangeAndNBins(Float_t min, Float_t max, Int_t n) {
 162     fHistoNPhiBins = n ;
 163     fHistoPhiMax = max ;
 164     fHistoPhiMin = min ;
 165   }
 166
 167   Int_t   GetHistoNPhiBins() const { return fHistoNPhiBins ; }
 168   Float_t GetHistoPhiMin()   const { return fHistoPhiMin ; }
 169   Float_t GetHistoPhiMax()   const { return fHistoPhiMax ; }
 170
 171   virtual void SetHistoEtaRangeAndNBins(Float_t min, Float_t max, Int_t n) {
 172     fHistoNEtaBins = n ;
 173     fHistoEtaMax = max ;
 174     fHistoEtaMin = min ;
 175   }
 176
 177   Int_t   GetHistoNEtaBins() const { return fHistoNEtaBins ; }
 178   Float_t GetHistoEtaMin()   const { return fHistoEtaMin ; }
 179   Float_t GetHistoEtaMax()   const { return fHistoEtaMax ; }
 180
 181  private:
 182
 183   Bool_t  fDataMC ;             // Flag to access MC data when using ESD or AOD
 184   Int_t   fDebug ;              // Debug level
 185   Bool_t  fCheckFidCut ;        // Do analysis for clusters in defined region
 186   Bool_t  fCheckCaloPID ;       // Do analysis for calorimeters
 187   Bool_t  fRecalculateCaloPID ; // Recalculate PID or use PID weights in calorimeters
 188   Float_t fMinPt ;              // Maximum pt of (trigger) particles in the analysis
 189   Float_t fMaxPt ;              // Minimum pt of (trigger) particles in the analysis
 190   Int_t   fLabelShift;          // In case AOD analysis with 2 input, shift in lable of second input.
 191
 192   AliCaloTrackReader * fReader; // Acces to ESD/AOD/MC data
 193
 194   TClonesArray* fInputAODBranch ;    //! Selected input particles branch
 195   TString       fInputAODName ;      //  Name of input AOD branch;
 196   TClonesArray* fOutputAODBranch ;   //! Selected output particles branch
 197   Bool_t        fNewAOD ;            //  Flag, new aod branch added to the analysis or not.
 198   TString       fOutputAODName ;     //  Name of output AOD branch;
 199   TString       fOutputAODClassName; //  Type of aod objects to be stored in the TClonesArray (AliAODPWG4Particle, AliAODPWG4ParticleCorrelation ...)
 200   TString       fAODObjArrayName ;   // Name of ref array kept in a TList in AliAODParticleCorrelation with clusters or track references.
 201   TString       fAddToHistogramsName;// Add this string to histograms name
 202
 203   //TClonesArray* fAODCaloClusters ;  //! selected PHOS/EMCAL CaloClusters
 204   AliAODCaloCells * fAODCaloCells ; //! selected PHOS/EMCAL CaloCells
 205
 206   //Analysis helper classes access pointers
 207   AliCaloPID               * fCaloPID; // PID calculation
 208   AliFidutialCut           * fFidCut;  // Acceptance cuts
 209   AliIsolationCut          * fIC;      // Isolation cut
 210   AliMCAnalysisUtils       * fMCUtils; // MonteCarlo Analysis utils
 211   AliNeutralMesonSelection * fNMS;     // Neutral Meson Selection
 212
 213   //Histograms binning and range
 214   Int_t   fHistoNPtBins ;  // Number of bins in pt axis
 215   Float_t fHistoPtMax ;    // Maximum value of pt histogram range
 216   Float_t fHistoPtMin ;    // Minimum value of pt histogram range
 217   Int_t   fHistoNPhiBins ; // Number of bins in phi axis
 218   Float_t fHistoPhiMax ;   // Maximum value of phi histogram range
 219   Float_t fHistoPhiMin ;   // Minimum value of phi histogram range
 220   Int_t   fHistoNEtaBins ; // Number of bins in eta axis
 221   Float_t fHistoEtaMax ;   // Maximum value of eta histogram range
 222   Float_t fHistoEtaMin ;   // Minimum value of eta histogram range
 223
 224   ClassDef(AliAnaPartCorrBaseClass,4)
 225     } ;
 226
 227
 228 #endif //ALIANAPARTCORRBASECLASS_H
 229
 230
 231