CRT/AliGenCRT.h

   1 #ifndef ALIGENCRT_H
   2 #define ALIGENCRT_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 #include "AliGenerator.h"
   9 #include "AliCRTConstants.h"
  10
  11 class TClonesArray;
  12 class TF1;
  13
  14 class AliGenCRT : public AliGenerator {
  15  public:
  16   AliGenCRT();
  17   AliGenCRT(Int_t npart);
  18   AliGenCRT(const AliGenCRT& gen);
  19   virtual ~AliGenCRT();
  20
  21   AliGenCRT& operator= (const AliGenCRT& gen);
  22
  23   virtual void Init();
  24   virtual void Generate();
  25   virtual void SetPart(Int_t part) {fIpart = part;}
  26
  27   void SetMode(ECRMode mode) {fCRMode = mode;}
  28   const TString*  GetMode() const {return fCRModeName;}
  29
  30   void SetZenithalAngleRange(Float_t min,Float_t max=0) {fZenithMin=min;fZenithMax=max;}
  31   void SetAzimuthalAngleRange(Float_t min, Float_t max=0) {fAzimuthMin=min;fAzimuthMax=max;}
  32
  33   void SetGridRange(Int_t nx,Float_t xwidth, Int_t nz, Float_t zwidth);
  34   const Float_t GetMomentumResolution() const {return fPResolution;}
  35
  36   void SetMomentumDistrubutionFunction(TF1 *func) {fMomentumDist=func;}
  37   void SetZenithalDistributionFunction(TF1 *func) {fZenithDist = func;}
  38   void SetMomentumResolution(Float_t res=1.) {fPResolution=res;}
  39
  40   const Float_t GetMomentum() const;
  41   const Float_t GetZenithAngle(Float_t mom) const;
  42
  43   // The following methods are for testing pourpuses
  44   TF1* GetMomentumDistibution() const {return fMomentumDist;}
  45   TF1* GetUnfoldedDistribution() const {return fUnfoldedMomentumDist;}
  46
  47  protected:
  48   void InitApWeightFactors();
  49   void InitMomentumGeneration();
  50   void InitZenithalAngleGeneration();
  51   void GenerateOneMuonBundle();
  52   void GenerateOneSingleMuon(Bool_t withFlatMomentum=kFALSE);
  53
  54  private:
  55   Int_t    fIpart;              //! Particle type.
  56   ECRMode  fCRMode;             //! Cosmic muons generation method flag
  57   TString* fCRModeName;         //! Cosmic muons generation mode name
  58
  59   Float_t  fXwidth;             //! X width of the grid
  60   Int_t    fNx;                 //! Number of divisions in x
  61   Float_t  fZwidth;             //! Z widht of the  grid
  62   Int_t    fNz;                 //! Number of divisions in z
  63   Bool_t   fMuonGrid;           //! Flag for method (Muon-bundles) checkout
  64
  65   Float_t  fZenithMin;          //! Minimum zenithal angle.
  66   Float_t  fZenithMax;          //! Maximum zenithal angle.
  67
  68   Float_t  fAzimuthMin;         //! Minimum azimuthal angle.
  69   Float_t  fAzimuthMax;         //! Maximum azimuthal angle.
  70
  71   Float_t  fPRange;             //! Cosmic muon momentum range width in GeVs.
  72   Float_t  fPResolution;        //! Momentum resolution in GeVs.
  73
  74   TArrayF* fAp;                 //! a(p) correction factors for the ang. dist.
  75
  76   TF1*     fMomentumDist;       //! Function to generate the momentum dist.
  77   TF1*     fUnfoldedMomentumDist; //!
  78   TF1*     fZenithDist;         //! Function to generate the zenith angle dist.
  79
  80   TClonesArray* fPDist;         //! Array of fZenithDist, to be used by a(p).
  81
  82   ClassDef(AliGenCRT, 1) // Generator for AliCRT class
  83 };
  84 #endif // ALIGENCRT_H