TRD/AliTRDsimpleMC.h

   1 #ifndef ALITRDSIMPLEMC_H
   2 #define ALITRDSIMPLEMC_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 ///////////////////////////////////////////////////////////////////////////////
   9 //                                                                           //
  10 //  Simple TRD Monte Carlo class                                             //
  11 //                                                                           //
  12 ///////////////////////////////////////////////////////////////////////////////
  13
  14 #include <TMCProcess.h>
  15 #include <TVirtualMC.h>
  16 #include "AliDecayer.h"
  17
  18 class AliTRDv1;
  19 class AliTRDparameter;
  20
  21 class AliTRDsimpleMC : public TVirtualMC {
  22
  23  public:
  24
  25   enum {
  26       kPdgElectron = 11
  27     , kPdgPion     = 211
  28   };
  29
  30   AliTRDsimpleMC();
  31   AliTRDsimpleMC(const char *name, const char *title);
  32   AliTRDsimpleMC(const AliTRDsimpleMC &m);
  33
  34   virtual ~AliTRDsimpleMC();
  35   AliTRDsimpleMC &operator=(const AliTRDsimpleMC &m);
  36
  37   virtual void          Copy(TObject &m);
  38
  39   //
  40   // Methods for building / management of geometry
  41   //
  42
  43   // Functions from GCONS
  44     virtual void  Gfmate(Int_t , char* , Float_t& , Float_t& ,
  45                          Float_t& , Float_t& , Float_t& ,
  46                          Float_t* , Int_t& )  {}
  47     virtual void  Gfmate(Int_t , char* , Double_t& , Double_t& ,
  48                          Double_t& , Double_t& , Double_t& ,
  49                          Double_t* , Int_t& ) {}
  50
  51   // Detector composition
  52     virtual void  Material(Int_t& , const char* , Double_t ,
  53                      Double_t , Double_t , Double_t , Double_t ,
  54                      Float_t* , Int_t ) {}
  55     virtual void  Material(Int_t& , const char* , Double_t ,
  56                      Double_t , Double_t , Double_t , Double_t ,
  57                      Double_t* , Int_t ) {}
  58     virtual void  Mixture(Int_t& , const char* , Float_t* ,
  59                      Float_t *, Double_t , Int_t , Float_t* ) {}
  60     virtual void  Mixture(Int_t& , const char* , Double_t *,
  61                      Double_t *, Double_t , Int_t , Double_t* ) {}
  62     virtual void  Medium(Int_t& , const char* , Int_t ,
  63                      Int_t , Int_t , Double_t , Double_t ,
  64                      Double_t , Double_t , Double_t ,
  65                      Double_t , Float_t* , Int_t ) {}
  66     virtual void  Medium(Int_t& , const char* , Int_t ,
  67                      Int_t , Int_t , Double_t , Double_t ,
  68                      Double_t , Double_t , Double_t ,
  69                      Double_t , Double_t* , Int_t ) {}
  70     virtual void  Matrix(Int_t& , Double_t , Double_t ,
  71                      Double_t , Double_t , Double_t ,
  72                      Double_t ) {}
  73     virtual void  Gstpar(Int_t , const char* , Double_t ) {}
  74
  75   // Functions from GGEOM
  76     virtual Int_t  Gsvolu(const char* , const char* , Int_t ,
  77                           Float_t* , Int_t )  { return 0; }
  78     virtual Int_t  Gsvolu(const char* , const char* , Int_t ,
  79                           Double_t* , Int_t) { return 0; }
  80     virtual void  Gsdvn(const char* , const char* , Int_t ,
  81                          Int_t ) {}
  82     virtual void  Gsdvn2(const char* , const char* , Int_t ,
  83                          Int_t , Double_t , Int_t ) {}
  84     virtual void  Gsdvt(const char* , const char* , Double_t ,
  85                          Int_t , Int_t , Int_t ) {}
  86     virtual void  Gsdvt2(const char* , const char* , Double_t ,
  87                          Int_t , Double_t , Int_t , Int_t ) {}
  88     virtual void  Gsord(const char* , Int_t ) {}
  89     virtual void  Gspos(const char* , Int_t , const char* ,
  90                          Double_t , Double_t , Double_t , Int_t ,
  91                          const char* ) {}
  92     virtual void  Gsposp(const char* , Int_t , const char* ,
  93                          Double_t, Double_t, Double_t, Int_t ,
  94                          const char* , Float_t* , Int_t ) {}
  95     virtual void  Gsposp(const char* , Int_t , const char* ,
  96                          Double_t , Double_t , Double_t , Int_t ,
  97                          const char* , Double_t* , Int_t ) {}
  98     virtual void  Gsbool(const char* , const char* ) {}
  99
 100     virtual void  SetCerenkov(Int_t , Int_t , Float_t* ,
 101                                Float_t* , Float_t* , Float_t* ) {}
 102     virtual void  SetCerenkov(Int_t , Int_t , Double_t* ,
 103                                Double_t* , Double_t* , Double_t* ) {}
 104
 105   // Functions for drawing
 106     virtual void  DrawOneSpec(const char* ) {}
 107     virtual void  Gsatt(const char* , const char* , Int_t ) {}
 108     virtual void  Gdraw(const char* , Double_t , Double_t ,
 109                         Double_t , Double_t , Double_t ,
 110                         Double_t , Double_t ) {}
 111
 112   // Euclid
 113   virtual void          WriteEuclid(const char* , const char* , Int_t , Int_t ) {}
 114
 115   // Get methods
 116     virtual Int_t VolId(const Text_t* volName) const;
 117     virtual const char* VolName(Int_t ) const { return ""; }
 118     virtual Int_t NofVolumes() const { return 0; }
 119     virtual Int_t VolId2Mate(Int_t ) const { return 0; }
 120
 121   //
 122   // Methods for physics management
 123   //
 124
 125   // Set methods
 126     virtual void     SetCut(const char* , Double_t ) {}
 127     virtual void     SetProcess(const char* , Int_t ) {}
 128     virtual Double_t Xsec(char*, Double_t, Int_t, Int_t) { return 0.; }
 129
 130   // Particle table usage
 131     virtual Int_t   IdFromPDG(Int_t ) const { return 0; }
 132     virtual Int_t   PDGFromId(Int_t ) const { return 0; }
 133     virtual void    DefineParticles() {}
 134
 135   //
 136   // Methods for step management
 137   //
 138
 139   // Action methods
 140   virtual void          StopTrack() { };
 141   virtual void          StopEvent() { };
 142
 143   // Set methods
 144   virtual void          SetMaxStep(Double_t step)                                         { fMaxStep = step; };
 145   virtual void          SetMaxNStep(Int_t )                                              { };
 146   virtual void          SetUserDecay(Int_t )                                             { };
 147
 148   virtual void          NewTrack(Int_t iTrack, Int_t pdg, Double_t px, Double_t py, Double_t pz);
 149
 150   // Tracking volume(s)
 151   virtual Int_t         CurrentVolID(Int_t& copyNo) const;
 152   virtual Int_t         CurrentVolOffID(Int_t off, Int_t& copyNo) const;
 153   virtual const char*   CurrentVolName() const;
 154   virtual const char*   CurrentVolOffName(Int_t ) const                                { return ""; };
 155   virtual Int_t         CurrentMaterial(Float_t& , Float_t& ,
 156                                         Float_t& , Float_t& ,
 157                                         Float_t& ) const                               { return 0;  };
 158   virtual Int_t         CurrentEvent() const                                              { return 0;  };
 159   virtual void          Gmtod(Float_t* , Float_t* , Int_t )   {}
 160   virtual void          Gmtod(Double_t* , Double_t* , Int_t ) {}
 161   virtual void          Gdtom(Float_t* , Float_t* , Int_t )   {}
 162   virtual void          Gdtom(Double_t* , Double_t* , Int_t ) {}
 163   virtual Double_t      MaxStep() const                                                   { return fMaxStep; };
 164   virtual Int_t         GetNStep() const                                                  { return fNStep;   };
 165   virtual Int_t         GetMaxNStep() const                                               { return 0;  };
 166   virtual Int_t         GetMedium() const                                                 { return 0;  };
 167
 168   // Dynamic properties
 169   virtual void          TrackPosition(TLorentzVector& position) const;
 170   virtual void          TrackPosition(Double_t &x, Double_t &y, Double_t &z) const;
 171   virtual void          TrackMomentum(TLorentzVector& momentum) const;
 172   virtual void          TrackMomentum(Double_t &px, Double_t &py, Double_t &pz, Double_t &etot) const;
 173   virtual Double_t      TrackStep() const                                                 { return fTrackStep; };
 174   virtual Double_t      TrackLength() const                                               { return 0.0; };
 175   virtual Double_t      TrackTime() const                                                 { return 0.0; };
 176   virtual Double_t      Edep() const                                                      { return 0.0; };
 177
 178   // Static properties
 179   virtual Int_t         TrackPid() const                                                  { return fTrackPid;    };
 180   virtual Double_t      TrackCharge() const                                               { return fTrackCharge; };
 181   virtual Double_t      TrackMass() const                                                 { return fTrackMass;   };
 182   virtual Double_t      Etot() const                                                      { return fTrackEtot;   };
 183
 184   // Track status
 185   virtual Bool_t        IsNewTrack() const                                                { return kFALSE; };
 186   virtual Bool_t        IsTrackInside() const                                             { return kFALSE; };
 187   virtual Bool_t        IsTrackEntering() const                                           { return fTrackEntering; };
 188   virtual Bool_t        IsTrackExiting() const                                            { return kFALSE; };
 189   virtual Bool_t        IsTrackOut() const                                                { return kFALSE; };
 190   virtual Bool_t        IsTrackDisappeared() const                                        { return kFALSE; };
 191   virtual Bool_t        IsTrackStop() const                                               { return kFALSE; };
 192   virtual Bool_t        IsTrackAlive() const                                              { return kFALSE; };
 193
 194   // Secondaries
 195   virtual Int_t         NSecondaries() const                                              { return 0; };
 196   virtual void          GetSecondary(Int_t , Int_t& ,
 197                                      TLorentzVector& ,
 198                                      TLorentzVector& )                                    { };
 199   virtual TMCProcess  ProdProcess(Int_t ) const                                           { return kPNoProcess; };
 200   virtual Int_t         StepProcesses(TArrayI& ) const                                    { return 0; };
 201
 202   //
 203   // Other (then geometry/step/run management) methods
 204   //
 205
 206   // Geant3 specific methods
 207     virtual void Gdopt(const char*,const char*) {}
 208     virtual void SetClipBox(const char*,Double_t=-9999,Double_t=0, Double_t=-9999,
 209                              Double_t=0,Double_t=-9999,Double_t=0) {}
 210     virtual void DefaultRange() {}
 211     virtual void Gdhead(Int_t, const char*, Double_t=0) {}
 212     virtual void Gdman(Double_t, Double_t, const char*) {}
 213     virtual void SetColors() {}
 214     virtual void Gtreve() {}
 215     virtual void GtreveRoot() {}
 216     virtual void Gckmat(Int_t, char*) {}
 217     virtual void InitLego() {}
 218     virtual void Gfpart(Int_t, char*, Int_t&, Float_t&, Float_t&, Float_t&) {}
 219     virtual void Gspart(Int_t, const char*, Int_t, Double_t, Double_t, Double_t) {}
 220     virtual void DefineParticle(int, const char*, TMCParticleType, double, double, double);
 221     virtual void DefineIon(const char*, int, int, int, double, double);
 222     virtual TString  ParticleName(int) const {return "";}
 223     virtual Double_t ParticleMass(int) const {return 0.;}
 224     virtual Double_t ParticleCharge(int) const {return 0.;}
 225     virtual Double_t ParticleLifeTime(int) const {return 0.;}
 226     virtual TMCParticleType ParticleMCType(int) const {return (TMCParticleType) 0;}
 227
 228   // Control Methods
 229   virtual void          Init()                                                            { };
 230   virtual void          FinishGeometry()                                                  { };
 231   virtual void          BuildPhysics()                                                    { };
 232   virtual void          ProcessEvent();
 233   virtual void          ProcessRun(Int_t )                                                { };
 234   //virtual TMCGeomType   GetMCGeomType() const                                             { return kGeant3; }
 235
 236   // External Decayer
 237   virtual void          SetExternalDecayer(AliDecayer* )                                  { };
 238   virtual AliDecayer   *Decayer() const                                                   { return 0; };
 239
 240  protected:
 241
 242   enum {
 243       kVolDrRg
 244     , kVolAmRg
 245     , kVolDrCh
 246   };
 247
 248   Float_t          fMaxStep;            //  Maximum step size
 249   Int_t            fNStep;              //  Number of steps
 250   Int_t            fTrack;              //  Track number
 251   Double_t         fTrackPx;            //  Track px
 252   Double_t         fTrackPy;            //  Track py
 253   Double_t         fTrackPz;            //  Track pz
 254   Double_t         fTrackPtot;          //  Track total momentum
 255   Double_t         fTrackEtot;          //  Track total energy
 256   Double_t         fTrackX;             //  Track x position
 257   Double_t         fTrackY;             //  Track y position
 258   Double_t         fTrackZ;             //  Track z position
 259   Double_t         fX0;                 //  X position of the beginning of the chamber
 260   Double_t         fTrackStep;          //  Track step size
 261   Int_t            fTrackPid;           //  Track PID
 262   Float_t          fTrackCharge;        //  Track charge
 263   Float_t          fTrackMass;          //  Track particle mass
 264   Bool_t           fTrackEntering;      //  Track entering chamber
 265
 266   AliTRDv1        *fTRD;                //! TRD detector object
 267   AliTRDparameter *fPar;                //! TRD parameter object
 268
 269   ClassDef(AliTRDsimpleMC,2)            //  Simple TRD Monte Carlo class
 270
 271 };
 272 #endif
 273