PWGLF/FORWARD/analysis2/AliAODCentralMult.h

   1 //
   2 // See implementation or Doxygen comments for more information
   3 //
   4 #ifndef ALIAODCENTRALMULT_H
   5 #define ALIAODCENTRALMULT_H
   6 /**
   7  * @file   AliAODCentralMult.h
   8  * @author Hans Hjersing Dalsgaard
   9  * @date   Wed Mar 23 13:57:24 2011
  10  *
  11  * @brief
  12  *
  13  * @ingroup pwg2_forward_aod
  14  *
  15  */
  16 #include <TObject.h>
  17 #include <TH2D.h>
  18 class TBrowser;
  19 /**
  20  * Class that contains the central multiplicity data per event
  21  *
  22  * This class contains a histogram of
  23  * @f[
  24  *   \frac{d^2N_{ch}}{d\eta d\phi}\quad,
  25  * @f]
  26  * as well as a trigger mask for each analysed event.
  27  *
  28  * The eta acceptance of the event is stored in the underflow bins of
  29  * the histogram.  So to build the final histogram, one needs to
  30  * correct for this acceptance (properly weighted by the events), and
  31  * the vertex efficiency.  This simply boils down to defining a 2D
  32  * histogram and summing the event histograms in that histogram.  One
  33  * should of course also do proper book-keeping of the accepted event.
  34  *
  35  * @ingroup pwg2_forward
  36  * @ingroup pwg2_forward_aod
  37  */
  38 class AliAODCentralMult : public TObject
  39 {
  40 public:
  41
  42   /**
  43    * Default constructor
  44    *
  45    * Used by ROOT I/O sub-system - do not use
  46    */
  47   AliAODCentralMult();
  48   /**
  49    * Constructor
  50    *
  51    * @param isMC Whether this was from MC or not
  52    */
  53   AliAODCentralMult(Bool_t isMC);
  54   /**
  55    * Destructor
  56    */
  57   virtual ~AliAODCentralMult() {} // Destructor
  58   /**
  59    * Initialize
  60    *
  61    * @param etaAxis  Pseudo-rapidity axis
  62    */
  63   void Init(const TAxis& etaAxis);
  64   /**
  65    * Get the @f$ d^2N_{ch}/d\eta d\phi@f$ histogram,
  66    *
  67    * @return @f$ d^2N_{ch}/d\eta d\phi@f$ histogram,
  68    */
  69   const TH2D& GetHistogram() const { return fHist; } // Get histogram
  70   /**
  71    * Get the @f$ d^2N_{ch}/d\eta d\phi@f$ histogram,
  72    *
  73    * @return @f$ d^2N_{ch}/d\eta d\phi@f$ histogram,
  74    */
  75   TH2D& GetHistogram() { return fHist; } // Get histogram
  76   /**
  77    * Clear Object between events
  78    *
  79    * @param opt option (not used)
  80    *
  81    */
  82   void Clear(Option_t* opt);
  83   /**
  84    * browse this object
  85    *
  86    * @param b Browser
  87    */
  88   void Browse(TBrowser* b);
  89   /**
  90    * This is a folder
  91    *
  92    * @return Always true
  93    */
  94   Bool_t IsFolder() const { return kTRUE; } // Always true
  95   /**
  96    * Print content
  97    *
  98    * @param option Passed verbatim to TH2::Print
  99    */
 100   void Print(Option_t* option="") const;
 101
 102   /**
 103    * Get the name of the object
 104    *
 105    * @return Name of object
 106    */
 107   const Char_t* GetName() const { return (fIsMC ? "CentralClustersMC" : "CentralClusters"); }
 108
 109 protected:
 110   Bool_t  fIsMC;     // Whether this is from MC
 111   TH2D    fHist;     // Histogram of d^2N_{ch}/(deta dphi) for this event
 112
 113   ClassDef(AliAODCentralMult,1); // AOD forward multiplicity
 114 };
 115
 116 #endif
 117 // Local Variables:
 118 //  mode: C++
 119 // End:
 120