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