The revisited SPD cluster analysis in the new AOD framework and the corrections for...
[u/mrichter/AliRoot.git] / PWG2 / FORWARD / analysis2 / AliAODCentralMult.cxx
1 //
2 // Class that contains the central multiplicity data per event 
3 //
4 // This class contains a histogram of 
5 // @f[
6 //   \frac{d^2N_{ch}}{d\eta d\phi}\quad,
7 // @f]
8 // as well as a trigger mask for each analysed event.  
9 // 
10 // The eta acceptance of the event is stored in the underflow bins of
11 // the histogram.  So to build the final histogram, one needs to
12 // correct for this acceptance (properly weighted by the events), and
13 // the vertex efficiency.  This simply boils down to defining a 2D
14 // histogram and summing the event histograms in that histogram.  One
15 // should of course also do proper book-keeping of the accepted event.
16 //
17 #include "AliAODCentralMult.h"
18 #include <TBrowser.h>
19 #include <iostream>
20 #include <TMath.h>
21 #include <TObjString.h>
22
23 ClassImp(AliAODCentralMult)
24 #if 0 
25 ; // For Emacs 
26 #endif
27
28 //____________________________________________________________________
29 AliAODCentralMult::AliAODCentralMult()
30   : fIsMC(false),
31     fHist()
32 {
33   // 
34   // Constructor 
35   // 
36 }
37
38 //____________________________________________________________________
39 AliAODCentralMult::AliAODCentralMult(Bool_t isMC) 
40   : fIsMC(isMC),
41     fHist("centralMult", "d^{2}N_{ch}/d#etad#varphi in the central regions", 
42           200, -4, 6, 20, 0, 2*TMath::Pi())
43 {
44   // 
45   // Constructor 
46   // 
47   // Parameters: 
48   //  isMC   If set to true this is for MC data (effects branch name)
49   // 
50   fHist.SetXTitle("#eta");
51   fHist.SetYTitle("#varphi [radians]");
52   fHist.SetZTitle("#frac{d^{2}N_{ch}}{d#etad#varphi}");
53   fHist.SetDirectory(0);
54   fHist.Sumw2();
55 }
56 //____________________________________________________________________
57 void
58 AliAODCentralMult::Clear(Option_t*) {
59   
60   fHist.Reset();
61   
62 }
63 //____________________________________________________________________
64 void
65 AliAODCentralMult::Init(const TAxis& etaAxis)
66 {
67   // Initialize the histogram with an eta axis 
68   // 
69   // Parameters: 
70   //   etaAxis       Eta axis to use 
71   // 
72   fHist.SetBins(etaAxis.GetNbins(), etaAxis.GetXmin(), etaAxis.GetXmax(), 
73                 20, 0, 2*TMath::Pi());
74 }
75
76 //____________________________________________________________________
77 void
78 AliAODCentralMult::Browse(TBrowser* b)
79 {
80   // Browse this object 
81   // 
82   // Parameters: 
83   //   b   Browser to use 
84
85   b->Add(&fHist);
86
87 }
88 //____________________________________________________________________
89 void
90 AliAODCentralMult::Print(Option_t* option) const
91 {
92   // Print this object 
93   // 
94   // Parameters: 
95   //  option   Passed to TH1::Print 
96   fHist.Print(option);
97 }
98
99 //____________________________________________________________________
100 //
101 // EOF
102 //