PWGPP/TPC/AliPerformanceMC.h

   1 #ifndef ALIPERFORMANCEMC_H
   2 #define ALIPERFORMANCEMC_H
   3
   4 //------------------------------------------------------------------------------
   5 // Class to keep information for MC  tracks
   6 // to check propagation algorithm, B-field and material budget.
   7 //
   8 // Author: J.Otwinowski 09/06/2009
   9 //------------------------------------------------------------------------------
  10
  11 class TString;
  12 class TNamed;
  13 class TCanvas;
  14 class TH1F;
  15 class TH2F;
  16 class TParticle;
  17
  18 class AliESDVertex;
  19 class AliESDtrack;
  20 class AliMCEvent;
  21 class AliStack;
  22 class AliTrackReference;
  23 class AliESDEvent;
  24 class AliESDfriend;
  25 class AliESDfriendTrack;
  26 class AliMCEvent;
  27 class AliMCParticle;
  28 class AliMCInfoCuts;
  29 class AliRecInfoCuts;
  30 class AliExternalTrackParam;
  31
  32 #include "THnSparse.h"
  33 #include "AliPerformanceObject.h"
  34
  35 class AliPerformanceMC : public AliPerformanceObject {
  36 public :
  37   AliPerformanceMC(const Char_t* name="AliPerformanceMC", const Char_t* title="AliPerformanceMC",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE);
  38   virtual ~AliPerformanceMC();
  39
  40   // Init data members
  41   virtual void  Init();
  42
  43   // Execute analysis
  44   virtual void  Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
  45
  46   // Merge output objects (needed by PROOF)
  47   virtual Long64_t Merge(TCollection* const list);
  48
  49   // Analyse output histograms
  50   virtual void Analyse();
  51
  52   // Get analysis folder
  53   virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
  54
  55   // Process events
  56   void ProcessTPC(AliTrackReference* const refIn, TParticle* const particle);
  57   void ProcessInnerTPC(AliTrackReference* const refIn,  AliTrackReference* const refOut, TParticle* const particle);
  58   void ProcessOuterTPCExt(TParticle *const  part, TClonesArray * const trefs);
  59   AliExternalTrackParam * MakeTrack(const AliTrackReference* ref, TParticle *const part);
  60
  61   // Create folder for analysed histograms
  62   TFolder *CreateFolder(TString folder = "folderMC", TString title = "Analysed MC histograms");
  63
  64   // Export objects to folder
  65   TFolder *ExportToFolder(TObjArray * array=0);
  66
  67   // Selection cuts
  68   void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
  69   void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
  70
  71   AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}
  72   AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}
  73
  74   TH1F*  MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0);
  75
  76   // getters
  77   //
  78   THnSparse *GetResolHisto() const  { return fResolHisto; }
  79   THnSparse *GetPullHisto()  const  { return fPullHisto; }
  80
  81 private:
  82   //
  83   // Control histograms
  84   // 5 track parameters (details in STEER/AliExternalTrackParam.h)
  85   //
  86
  87   // resolution histogram
  88   THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt
  89
  90   // pull histogram
  91   THnSparseF *fPullHisto;  //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt
  92
  93   // Global cuts objects
  94   AliRecInfoCuts*  fCutsRC;      // selection cuts for reconstructed tracks
  95   AliMCInfoCuts*  fCutsMC;       // selection cuts for MC tracks
  96
  97   // analysis folder
  98   TFolder *fAnalysisFolder; // folder for analysed histograms
  99
 100   AliPerformanceMC(const AliPerformanceMC&); // not implemented
 101   AliPerformanceMC& operator=(const AliPerformanceMC&); // not implemented
 102
 103   ClassDef(AliPerformanceMC,1);
 104 };
 105
 106 #endif