AliFMDMCTrackInspector.h

Go to the documentation of this file.

#ifndef ALIFMDMCTRACKINSPECTOR_H
#define ALIFMDMCTRACKINSPECTOR_H
#include "AliFMDEnergyFitter.h"
#include "AliFMDMCTrackELoss.h"
#include <TArrayF.h>
class AliMCEvent;
class AliESDEvent;

class AliFMDMCTrackInspector : public AliFMDEnergyFitter
{
public:
  AliFMDMCTrackInspector();
  AliFMDMCTrackInspector(const char* title);
  virtual ~AliFMDMCTrackInspector();

  virtual void CreateOutputObjects(TList* dir);
  virtual Bool_t PreEvent(const AliMCEvent& mcInput);
  virtual Bool_t Event(const AliESDEvent&  esdInput,
           const AliMCEvent&   mcInput,
           Double_t            cent=-1);
  virtual Bool_t PostEvent();

  AliFMDMCTrackELoss& GetTracker() { return fTracker; }
  const AliFMDMCTrackELoss& GetTracker() const { return fTracker; }
protected:
  AliFMDMCTrackInspector(const AliFMDMCTrackInspector&);
  AliFMDMCTrackInspector& operator=(const AliFMDMCTrackInspector&);

public:
  struct RingHistos : public AliFMDEnergyFitter::RingHistos
  {
    RingHistos();
    RingHistos(UShort_t d, Char_t r);
    ~RingHistos() {}
    RingHistos(const RingHistos& o){;}
    RingHistos& operator=(const RingHistos& o){return *this;}
    TArrayD MakeIncreasingAxis(Int_t nBins, 
             Double_t low,
             Double_t high) const;
    virtual void SetupForData(const TAxis& eAxis, 
            const TAxis& cAxis,
            Double_t     maxDE=10, 
            Int_t        nDEbins=300, 
            Bool_t       useIncrBin=true);
    virtual void FillMC(UShort_t flag, Double_t eta, Double_t mult);
    virtual void  Scale(TH1* dist) const;
    TObjArray* Fit(TList*           dir, 
       Double_t         lowCut, 
       UShort_t         nParticles,
       UShort_t         minEntries,
       UShort_t         minusBins, 
       Double_t         relErrorCut, 
       Double_t         chi2nuCut,
       Double_t         minWeight,
       Double_t         regCut,
       EResidualMethod  residuals) const;
    
    TH2* fPrimary;   // @f$\Delta@f$ vs @f$\eta@f$ for primaries
    TH2* fSecondary; // @f$\Delta@f$ vs @f$\eta@f$ for second.
  public:
    TH2* fBetaGammadEdx;
    TH2* fBetaGammaEta;
    TH2* fDedxEta;
    ClassDef(RingHistos,1); // Cache of histograms per ring 
  };
protected:
  AliFMDEnergyFitter::RingHistos* CreateRingHistos(UShort_t d, Char_t r) const;
  AliFMDMCTrackELoss fTracker;
  TArrayF          fIp;
  Int_t            fNTrack;
  Int_t            fNPrimary;

  ClassDef(AliFMDMCTrackInspector,1);
};

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