PWGLF/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 //