PWG0/AliPWG0Helper.h

   1 /* $Id$ */
   2
   3 #ifndef ALIPWG0HELPER_H
   4 #define ALIPWG0HELPER_H
   5
   6 #include <TObject.h>
   7 #include <AliTriggerAnalysis.h>
   8
   9 // static helper functions
  10
  11 class AliESDEvent;
  12 class AliESDVertex;
  13 class TParticle;
  14 class TH1;
  15 class TH2;
  16 class TH3;
  17 class AliHeader;
  18 class AliGenEventHeader;
  19 class AliStack;
  20 class TTree;
  21 class AliOfflineTrigger;
  22 class AliMultiplicity;
  23
  24 class AliPWG0Helper : public TObject
  25 {
  26   public:
  27     enum AnalysisMode { kInvalid = -1, kSPD = 0x1, kTPC = 0x2, kTPCITS = 0x4, kFieldOn = 0x8, kSPDOnlyL0 = 0x10 };
  28     // in case we want to use bitmaps...
  29     enum MCProcessType { kInvalidProcess = -1, kND = 0x1, kDD = 0x2, kSD = 0x4, kOnePart = 0x8 };
  30     enum DiffTreatment { kMCFlags = 0, kUA5Cuts = 1, kE710Cuts, kALICEHadronLevel };
  31
  32     static const AliESDVertex* GetVertex(const AliESDEvent* aEsd, AnalysisMode analysisMethod, Bool_t debug = kFALSE);
  33     static Bool_t TestVertex(const AliESDVertex* vertex, AnalysisMode analysisMode, Bool_t debug = kFALSE);
  34
  35     static Bool_t IsPrimaryCharged(TParticle* aParticle, Int_t aTotalPrimaries, Bool_t adebug = kFALSE);
  36
  37     static MCProcessType GetEventProcessType(AliESDEvent* esd, AliHeader* header, AliStack* stack, DiffTreatment diffTreatment);
  38     static MCProcessType GetEventProcessType(AliHeader* aHeader, Bool_t adebug = kFALSE);
  39     static MCProcessType GetPythiaEventProcessType(AliGenEventHeader* aHeader, Bool_t adebug = kFALSE);
  40     static MCProcessType GetDPMjetEventProcessType(AliGenEventHeader* aHeader, Bool_t adebug = kFALSE);
  41     static Int_t GetLastProcessType() { return fgLastProcessType; }
  42
  43     static TParticle* FindPrimaryMother(AliStack* stack, Int_t label);
  44     static Int_t FindPrimaryMotherLabel(AliStack* stack, Int_t label);
  45
  46     static void CreateProjections(TH3* hist, Bool_t save = kFALSE);
  47     static void CreateDividedProjections(TH3* hist, TH3* hist2, const char* axis = 0, Bool_t putErrors = kFALSE, Bool_t save = kFALSE);
  48     static const char* GetAxisTitle(TH3* hist, const char axis);
  49
  50     static void NormalizeToBinWidth(TH1* hist);
  51     static void NormalizeToBinWidth(TH2* hist);
  52
  53     static void PrintConf(AnalysisMode analysisMode, AliTriggerAnalysis::Trigger trigger, DiffTreatment diffTreatment);
  54
  55     static Double_t Rapidity(Double_t pt, Double_t pz, Double_t m);
  56     static Bool_t IsHadronLevelSingleDiffractive(AliStack* stack, Float_t cms, Float_t xiMin, Float_t xiMax);
  57
  58   protected:
  59     static Int_t fgLastProcessType;    // stores the raw value of the last process type extracted
  60
  61     ClassDef(AliPWG0Helper, 0)
  62
  63   private:
  64     AliPWG0Helper(const AliPWG0Helper&);
  65     AliPWG0Helper& operator=(const AliPWG0Helper&);
  66 };
  67
  68 #endif
  69