AliBaseMultTask.h

Go to the documentation of this file.

#ifndef ALIBASEMULTTASK_H
#define ALIBASEMULTTASK_H
#include "AliBaseAODTask.h"

class AliBaseMultTask : public AliBaseAODTask
{
public:

  struct Bin : public TNamed
  {
    Bin();
    Bin(Double_t etaLow, Double_t etaHigh);
    Bin(const Bin&){;}
    Bin& operator=(const Bin&) { return *this; }
    virtual ~Bin(){}
    virtual void CreateOutputObjects(TList* cont,  Int_t max);
    virtual void Process(TH1D*              dndetaForward,
       TH1D*              dndetaCentral,
       TH1D*              normForward,
       TH1D*              normCentral,
       TH1D*              mc,
       Double_t           ipZ,
       Bool_t             pileup,
       Bool_t             selectedTrigger,
       Bool_t             isMCClass,
       Bool_t             isESDClass,
       const AliAODEvent& aodevent,
       Double_t           minIPz,
       Double_t           maxIPz);
    virtual Double_t CalcMult(TH1D*              dndetaForward, 
            TH1D*              dndetaCentral,
            TH1D*              normForward,   
            TH1D*              normCentral,
            TH1D*              mc, 
            Double_t           ipZ,
            Double_t&          statErr,
            Double_t&          sysErr,
            Double_t&          mcMult, 
            Double_t&          mcErr);
    static const Char_t* FormBinName(Double_t etaLow, Double_t etaHigh);
    
    Double_t fEtaLow;                  // low eta limit 
    Double_t fEtaHigh;                 // high eta limit
    TH1D*    fHist;                    
    TH1D*    fHistMC;                  
    TH2D*    fAcceptance;              
    TH2D*    fVtxZvsNdataBins;         
    ClassDef(Bin,1); // Base class for eta bins 
  };

  //__________________________________________________________________
  AliBaseMultTask();
  virtual ~AliBaseMultTask() {}
  //__________________________________________________________________
  void SetMaxMult(Int_t max) { fMaxMult = max; }
  void SetMCIsNSD(Bool_t assume) { fMCIsNSD = assume; }
  /* @} */
  //__________________________________________________________________
  void AddBin(Double_t etaLow, Double_t etaHigh)
  {
    fBins.Add(MakeBin(etaLow, etaHigh));
  }
  virtual Bin* MakeBin(Double_t minEta,Double_t maxEta) = 0;
  virtual void DefaultBins();
  /* @} */
  //__________________________________________________________________
  virtual Bool_t Book();
  virtual Bool_t PreEvent() { fIsSelected = false; return true; }
  virtual Bool_t Event(AliAODEvent& aod);
  virtual Bool_t Finalize() {return true;}
protected:
  AliBaseMultTask(const char* name);
  AliBaseMultTask(const AliBaseMultTask&)
    : fBins(),
      fIsSelected(false),
      fMaxMult(500),
      fMCIsNSD(false)
  {}
  AliBaseMultTask& operator=(const AliBaseMultTask&) { return *this; }
  virtual Bool_t CheckEvent(const AliAODForwardMult& fwd);
  virtual Bool_t IsMCClass(AliAODForwardMult* m) const;
  virtual Bool_t IsESDClass(AliAODForwardMult* m) const;
  /* @} */
  virtual void Process(TH1D*              dndetaForward,
           TH1D*              dndetaCentral,
           TH1D*              normForward,
           TH1D*              normCentral,
           TH1D*              mc,
           Double_t           ipZ,
           Bool_t             pileup,
           Bool_t             selectedTrigger,
           Bool_t             isMCClass,
           Bool_t             isESDClass,
           const AliAODEvent& aodevent);
  TList  fBins;        // List of the bins to use 
  Bool_t fIsSelected;  // Event was selected 
  Int_t  fMaxMult;     // Maximum of (and number of bins in) histograms
  Bool_t fMCIsNSD;     // Assume MC event is always NSD 
  ClassDef(AliBaseMultTask,2); // Base class for multiplicity dist tasks,
};
#endif
// Local Variables:
//  mode: C++
// End: