PWG1/AliComparisonRes.h

   1 #ifndef ALICOMPARISONRES_H
   2 #define ALICOMPARISONRES_H
   3
   4 //------------------------------------------------------------------------------
   5 // Class to keep information from comparison of
   6 // reconstructed and MC particle tracks (TPC resolution).
   7 //
   8 // Author: J.Otwinowski 04/02/2008
   9 //------------------------------------------------------------------------------
  10
  11 class TString;
  12 class TNamed;
  13 class TCanvas;
  14 class TH1F;
  15 class TH2F;
  16
  17 class AliESDVertex;
  18 class AliESDtrack;
  19 class AliMCInfo;
  20 class AliESDRecInfo;
  21 class AliESDEvent;
  22 class AliMCInfoCuts;
  23 class AliRecInfoCuts;
  24
  25 #include "THnSparse.h"
  26 #include "AliComparisonObject.h"
  27
  28 class AliComparisonRes : public AliComparisonObject {
  29 public :
  30   AliComparisonRes();
  31   AliComparisonRes(Char_t* name, Char_t* title, Int_t analysisMode, Bool_t hptGenerator);
  32   virtual ~AliComparisonRes();
  33
  34   // Init data members
  35   virtual void  Init();
  36
  37   // Execute analysis
  38   virtual void  Exec(AliMCInfo* const infoMC, AliESDRecInfo *const infoRC);
  39
  40   // Merge output objects (needed by PROOF)
  41   virtual Long64_t Merge(TCollection* const list);
  42
  43   // Analyse output histograms
  44   virtual void Analyse();
  45
  46   // Get analysis folder
  47   virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
  48
  49   // Process events
  50   void ProcessConstrained(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
  51   void ProcessTPC(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
  52   void ProcessTPCITS(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
  53
  54   // Create folder for analysed histograms
  55   TFolder *CreateFolder(TString folder = "folderRes",TString title = "Analysed Resolution histograms");
  56
  57   // Export objects to folder
  58   TFolder *ExportToFolder(TObjArray * array=0);
  59
  60   // Selection cuts
  61   void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
  62   void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
  63
  64   AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}
  65   AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}
  66
  67   static TH1F*     MakeResol(TH2F * his, Int_t integ, Bool_t type);
  68
  69   // getters
  70   //
  71   THnSparse *GetResolHisto() const  { return fResolHisto; }
  72   THnSparse *GetPullHisto()  const  { return fPullHisto; }
  73
  74 private:
  75   //
  76   // Control histograms
  77   // 5 track parameters (details in STEER/AliExternalTrackParam.h)
  78   //
  79
  80   // resolution histogram
  81   THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_1pt:y:z:eta:phi:pt
  82
  83   // pull histogram
  84   THnSparseF *fPullHisto;  //-> pull_y:pull_z:pull_phi:pull_lambda:pull_1pt:y:z:eta:phi:pt
  85
  86   // Global cuts objects
  87   AliRecInfoCuts*  fCutsRC;      // selection cuts for reconstructed tracks
  88   AliMCInfoCuts*  fCutsMC;       // selection cuts for MC tracks
  89
  90   // analysis folder
  91   TFolder *fAnalysisFolder; // folder for analysed histograms
  92
  93   AliComparisonRes(const AliComparisonRes&); // not implemented
  94   AliComparisonRes& operator=(const AliComparisonRes&); // not implemented
  95
  96   ClassDef(AliComparisonRes,1);
  97 };
  98
  99 #endif