AliFMDMCDensityCalculator.h

Go to the documentation of this file.

// 
// This class calculates the inclusive charged particle density
// in each for the 5 FMD rings based on the MC truth.
//
#ifndef ALIFMDMCDENSITYCALCULATOR_H
#define ALIFMDMCDENSITYCALCULATOR_H

#include "AliFMDDensityCalculator.h"
#include <TList.h>
#include "AliForwardUtil.h"
class AliMCEvent;
class TH2;
class TH2D;
class TProfile2D;

class AliFMDMCDensityCalculator : public AliFMDDensityCalculator
{
public:
  AliFMDMCDensityCalculator() 
    : AliFMDDensityCalculator(),
      fFMD1i(0), 
      fFMD2i(0),
      fFMD2o(0),
      fFMD3i(0),
      fFMD3o(0),
      fFMD1iC(0), 
      fFMD2iC(0),
      fFMD2oC(0),
      fFMD3iC(0),
      fFMD3oC(0),
      fFMD1iD(0), 
      fFMD2iD(0),
      fFMD2oD(0),
      fFMD3iD(0),
      fFMD3oD(0),
      fComps(0)
  {}
  AliFMDMCDensityCalculator(const char* name) 
   : AliFMDDensityCalculator(name),
      fFMD1i(0), 
      fFMD2i(0),
      fFMD2o(0),
      fFMD3i(0),
      fFMD3o(0),
      fFMD1iC(0), 
      fFMD2iC(0),
      fFMD2oC(0),
      fFMD3iC(0),
      fFMD3oC(0),
      fFMD1iD(0), 
      fFMD2iD(0),
      fFMD2oD(0),
      fFMD3iD(0),
      fFMD3oD(0),
      fComps(0)
  {}
  AliFMDMCDensityCalculator(const AliFMDMCDensityCalculator& o)
   : AliFMDDensityCalculator(o) ,
      fFMD1i(o.fFMD1i), 
      fFMD2i(o.fFMD2i),
      fFMD2o(o.fFMD2o),
      fFMD3i(o.fFMD3i),
      fFMD3o(o.fFMD3o),
      fFMD1iC(o.fFMD1iC), 
      fFMD2iC(o.fFMD2iC),
      fFMD2oC(o.fFMD2oC),
      fFMD3iC(o.fFMD3iC),
      fFMD3oC(o.fFMD3oC),
      fFMD1iD(o.fFMD1iD), 
      fFMD2iD(o.fFMD2iD),
      fFMD2oD(o.fFMD2oD),
      fFMD3iD(o.fFMD3iD),
      fFMD3oD(o.fFMD3oD),
      fComps(0)
  {}
  virtual ~AliFMDMCDensityCalculator();
  AliFMDMCDensityCalculator& operator=(const AliFMDMCDensityCalculator& o);
  void SetupForData(const TAxis& etaAxis);
  virtual Bool_t CalculateMC(const AliESDFMD&        fmd,
           AliForwardUtil::Histos& hists);
           
  virtual Bool_t CompareResults(AliForwardUtil::Histos& esd, 
        AliForwardUtil::Histos& mc);
  void CreateOutputObjects(TList* dir);
protected:
  TProfile2D* Make(UShort_t d, Char_t r, const TAxis& axis) const;
  TH2D* Make(UShort_t d, Char_t r) const;
  TH1D* Make(UShort_t d, Char_t r, Int_t max) const;
  void Fill(UShort_t d, Char_t r, TH2* esd, TH2* mc);

  TProfile2D* fFMD1i; // Comparison
  TProfile2D* fFMD2i; // Comparison
  TProfile2D* fFMD2o; // Comparison
  TProfile2D* fFMD3i; // Comparison
  TProfile2D* fFMD3o; // Comparison
  TH2D*       fFMD1iC; // Correlation in FMD1i
  TH2D*       fFMD2iC; // Correlation in FMD2i
  TH2D*       fFMD2oC; // Correlation in FMD2o
  TH2D*       fFMD3iC; // Correlation in FMD3i
  TH2D*       fFMD3oC; // Correlation in FMD3o
  TH1D*       fFMD1iD; // Correlation in FMD1i
  TH1D*       fFMD2iD; // Correlation in FMD2i
  TH1D*       fFMD2oD; // Correlation in FMD2o
  TH1D*       fFMD3iD; // Correlation in FMD3i
  TH1D*       fFMD3oD; // Correlation in FMD3o
  TList*      fComps; // List of comparisons 

  ClassDef(AliFMDMCDensityCalculator,2); // Calculate Nch density 
};

#endif
// Local Variables:
//   mode: C++
// End: