added correction for events with vertex but 0 tracks
[u/mrichter/AliRoot.git] / PWG0 / dNdEta / AlidNdEtaCorrection.h
1 /* $Id$ */
2
3 #ifndef ALIDNDETACORRECTION_H
4 #define ALIDNDETACORRECTION_H
5
6
7 // ------------------------------------------------------
8 //
9 // Class to handle corrections for dN/dEta measurements
10 //
11 // ------------------------------------------------------
12 //
13 // TODO:
14 // - make the ntrack to npart correction 3D
15 // - add documentation
16 // - add status: generate or use maps
17 // - add functionality to set the bin sizes
18 // - update MERge function
19 //
20
21 #include <TCollection.h>
22 #include <TNamed.h>
23 #include "AliPWG0Helper.h"
24
25 class AliCorrection;
26 class TH1;
27
28 class AlidNdEtaCorrection : public TNamed
29 {
30 public:
31   enum CorrectionType {
32     kNone = 0,
33     kTrack2Particle,  // measured events
34     kVertexReco,      // MB sample
35     kINEL,
36     kNSD,
37     kND
38   };
39
40   AlidNdEtaCorrection();
41   AlidNdEtaCorrection(const Char_t* name, const Char_t* title, AliPWG0Helper::AnalysisMode analysis = AliPWG0Helper::kTPC);
42
43   virtual Long64_t Merge(TCollection* list);
44
45   ~AlidNdEtaCorrection();
46
47   void FillMCParticle(Float_t vtx, Float_t eta, Float_t pt, Bool_t trigger, Bool_t vertex, Int_t processType);
48   void FillTrackedParticle(Float_t vtx, Float_t eta, Float_t pt);
49   void FillEvent(Float_t vtx, Float_t n, Bool_t trigger, Bool_t vertex, Int_t processType);
50
51   void Finish();
52
53   AliCorrection* GetTrack2ParticleCorrection()  {return fTrack2ParticleCorrection;}
54   AliCorrection* GetVertexRecoCorrection()      {return fVertexRecoCorrection;}
55   AliCorrection* GetTriggerBiasCorrectionINEL() {return fTriggerBiasCorrectionMBToINEL;}
56   AliCorrection* GetTriggerBiasCorrectionNSD()  {return fTriggerBiasCorrectionMBToNSD;}
57   AliCorrection* GetTriggerBiasCorrectionND()   {return fTriggerBiasCorrectionMBToND;}
58   AliCorrection* GetCorrection(CorrectionType correctionType);
59
60   void    Reset(void);
61   void    Add(AlidNdEtaCorrection* aCorrectionsToAdd, Float_t c=1);
62
63   void    SaveHistograms();
64   Bool_t  LoadHistograms(const Char_t* dir = 0);
65   void    DrawHistograms();
66   void    DrawOverview(const char* canvasName = 0);
67
68   Float_t GetMeasuredFraction(CorrectionType correctionType, Float_t ptCutOff, Float_t eta = -100, Bool_t debug = kFALSE);
69   TH1*    GetMeasuredEventFraction(CorrectionType correctionType, Int_t multCut);
70
71   void ReduceInformation();
72
73 protected:
74   AliCorrection* fTrack2ParticleCorrection;       //-> handles the track-to-particle correction (only track level (vtx_z, eta, pt))
75   AliCorrection* fVertexRecoCorrection;           //-> handles the vertex reconstruction efficiency, (n, vtx_z)
76   AliCorrection* fTriggerBiasCorrectionMBToINEL;  //-> handles the trigger bias MB->INEL, function of n and vtx_z
77   AliCorrection* fTriggerBiasCorrectionMBToNSD;   //-> handles the trigger bias MB->NSD,  function of n and vtx_z
78   AliCorrection* fTriggerBiasCorrectionMBToND;    //-> handles the trigger bias MB->ND,   function of n and vtx_z
79
80 private:
81   AlidNdEtaCorrection(const AlidNdEtaCorrection&);
82   AlidNdEtaCorrection& operator=(const AlidNdEtaCorrection&);
83
84   ClassDef(AlidNdEtaCorrection, 1)
85 };
86
87 #endif