AliFMDMCTrackDensity.h

Go to the documentation of this file.

#ifndef ALIFMDMCTRACKDENSITY_MC
#define ALIFMDMCTRACKDENSITY_MC
#include "AliForwardUtil.h"
#include "AliBaseMCTrackDensity.h"
class TH1D;
class AliESDFMD;

class AliFMDMCTrackDensity : public AliBaseMCTrackDensity 
{
public:
  AliFMDMCTrackDensity();
  AliFMDMCTrackDensity(const char* name);
  AliFMDMCTrackDensity(const AliFMDMCTrackDensity& o);
  AliFMDMCTrackDensity& operator=(const AliFMDMCTrackDensity& o);
  virtual ~AliFMDMCTrackDensity() {}

  void SetMaxConsequtiveStrips(UShort_t n) { fMaxConsequtiveStrips = n; }
  void SetLowCutvalue(Double_t v) { fLowCutvalue = v; }
  Bool_t Calculate(const AliESDFMD&    esd, 
       const AliMCEvent&   event, 
       const TVector3&     ip, 
       AliESDFMD&          output,
       TH2D*               primary);
  void CreateOutputObjects(TList* list);
  
  void Print(Option_t* option="") const;
protected:
  Int_t GetDetectorId() const;
  AliTrackReference* ProcessRef(AliMCParticle* particle, 
        const AliMCParticle* mother, 
        AliTrackReference* ref);
  void BeginTrackRefs();
  void EndTrackRefs(Int_t nRefs);
  Double_t StoreParticle(AliMCParticle*       particle, 
       const AliMCParticle* mother,
       AliTrackReference*   ref) const;
  mutable struct State 
  {
    Double_t angle;            // Angle 
    UShort_t oldDetector;      // Last detector
    Char_t   oldRing;          // Last ring
    UShort_t oldSector;        // Last sector
    UShort_t oldStrip;         // Last strip 
    UShort_t startStrip;       // First strip 
    UShort_t nRefs;            // Number of references
    UShort_t nStrips;          // Number of strips 
    UShort_t count;            // Count of hit strips 
    AliTrackReference* longest; 

    void Clear(Bool_t alsoCount=false);
    State& operator=(const State& o);
  } fState; 
  
    
  UShort_t   fMaxConsequtiveStrips; // Max 'cluster' size
  Double_t   fLowCutvalue;      // cut on minimal value of the dE to be part of cluster 
  TH1D*      fNr;                   // Number of track-refs per cluster
  TH1D*      fNt;                   // Size of cluster in strips 
  TH1D*      fNc;                   // Number of clusters per track
  TH2D*      fNcr;                  // Number of clusters per track
  AliESDFMD* fOutput;               

  ClassDef(AliFMDMCTrackDensity,6); // Calculate track-ref density
};

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