777ca9f6e3d80ab9c449c5986fdb9e18d2cbf191
[u/mrichter/AliRoot.git] / 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);
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     virtual void     SetCut(const char* , Double_t ) {}
128     virtual void     SetProcess(const char* , Int_t ) {}
129     virtual Double_t Xsec(char*, Double_t, Int_t, Int_t) { return 0.; } 
130  
131   // Particle table usage         
132     virtual Int_t   IdFromPDG(Int_t ) const { return 0; }  
133     virtual Int_t   PDGFromId(Int_t ) const { return 0; }  
134     virtual void    DefineParticles() {}      
135   
136   //
137   // Methods for step management
138   //
139
140   // Action methods
141   virtual void          StopTrack() { };
142   virtual void          StopEvent() { };   
143 //change is in 4.0.5 (which is not yet tagged) (rdm)
144 #if ROOT_VERSION_CODE >= 262148
145   virtual void          StopRun()   { }
146 #endif
147
148   // Set methods
149   virtual void          SetMaxStep(Double_t step)                                         { fMaxStep = step; };
150   virtual void          SetMaxNStep(Int_t )                                              { };
151   virtual void          SetUserDecay(Int_t )                                             { };  
152
153   virtual void          NewTrack(Int_t iTrack, Int_t pdg, Double_t px, Double_t py, Double_t pz);
154
155   // Tracking volume(s) 
156   virtual Int_t         CurrentVolID(Int_t& copyNo) const;
157   virtual Int_t         CurrentVolOffID(Int_t off, Int_t& copyNo) const;
158   virtual const char*   CurrentVolName() const;
159   virtual const char*   CurrentVolOffName(Int_t ) const                                { return ""; };
160   virtual Int_t         CurrentMaterial(Float_t& , Float_t& , 
161                                         Float_t& , Float_t& , 
162                                         Float_t& ) const                               { return 0;  };  
163   virtual Int_t         CurrentEvent() const                                              { return 0;  }; 
164   virtual void          Gmtod(Float_t* , Float_t* , Int_t )   {}
165   virtual void          Gmtod(Double_t* , Double_t* , Int_t ) {}
166   virtual void          Gdtom(Float_t* , Float_t* , Int_t )   {}
167   virtual void          Gdtom(Double_t* , Double_t* , Int_t ) {}
168   virtual Double_t      MaxStep() const                                                   { return fMaxStep; };
169   virtual Int_t         GetNStep() const                                                  { return fNStep;   };
170   virtual Int_t         GetMaxNStep() const                                               { return 0;  };
171   virtual Int_t         GetMedium() const                                                 { return 0;  };
172
173   // Dynamic properties
174   virtual void          TrackPosition(TLorentzVector& position) const;
175   virtual void          TrackPosition(Double_t &x, Double_t &y, Double_t &z) const;
176   virtual void          TrackMomentum(TLorentzVector& momentum) const;
177   virtual void          TrackMomentum(Double_t &px, Double_t &py, Double_t &pz, Double_t &etot) const;
178   virtual Double_t      TrackStep() const                                                 { return fTrackStep; };
179   virtual Double_t      TrackLength() const                                               { return 0.0; };
180   virtual Double_t      TrackTime() const                                                 { return 0.0; };
181   virtual Double_t      Edep() const                                                      { return 0.0; };
182   
183   // Static properties
184   virtual Int_t         TrackPid() const                                                  { return fTrackPid;    };
185   virtual Double_t      TrackCharge() const                                               { return fTrackCharge; };
186   virtual Double_t      TrackMass() const                                                 { return fTrackMass;   };
187   virtual Double_t      Etot() const                                                      { return fTrackEtot;   };
188
189   // Track status
190   virtual Bool_t        IsNewTrack() const                                                { return kFALSE; };
191   virtual Bool_t        IsTrackInside() const                                             { return kFALSE; };
192   virtual Bool_t        IsTrackEntering() const                                           { return fTrackEntering; };
193   virtual Bool_t        IsTrackExiting() const                                            { return kFALSE; };
194   virtual Bool_t        IsTrackOut() const                                                { return kFALSE; };
195   virtual Bool_t        IsTrackDisappeared() const                                        { return kFALSE; };
196   virtual Bool_t        IsTrackStop() const                                               { return kFALSE; };
197   virtual Bool_t        IsTrackAlive() const                                              { return kFALSE; };
198
199   // Secondaries
200   virtual Int_t         NSecondaries() const                                              { return 0; };
201   virtual void          GetSecondary(Int_t , Int_t& , 
202                                      TLorentzVector& , 
203                                      TLorentzVector& )                                    { };
204   virtual TMCProcess  ProdProcess(Int_t ) const                                           { return kPNoProcess; }; 
205   virtual Int_t         StepProcesses(TArrayI& ) const                                    { return 0; };
206
207   //
208   // Other (then geometry/step/run management) methods
209   //
210     
211   // Geant3 specific methods
212     virtual void Gdopt(const char*,const char*) {}
213     virtual void SetClipBox(const char*,Double_t=-9999,Double_t=0, Double_t=-9999,
214                              Double_t=0,Double_t=-9999,Double_t=0) {}
215     virtual void DefaultRange() {}
216     virtual void Gdhead(Int_t, const char*, Double_t=0) {}   
217     virtual void Gdman(Double_t, Double_t, const char*) {}
218     virtual void SetColors() {}
219     virtual void Gtreve() {}
220     virtual void GtreveRoot() {}
221     virtual void Gckmat(Int_t, char*) {}
222     virtual void InitLego() {}
223     virtual void Gfpart(Int_t, char*, Int_t&, Float_t&, Float_t&, Float_t&) {} 
224     virtual void Gspart(Int_t, const char*, Int_t, Double_t, Double_t, Double_t) {} 
225     // Dummy methods 
226 #if ROOT_VERSION_CODE > 197895
227     virtual void DefineParticle(int, const char*, TMCParticleType, double, double, double){;}
228     virtual void DefineIon(const char*, int, int, int, double, double){;}
229     virtual TString  ParticleName(int) const {return "";}
230     virtual Double_t ParticleMass(int) const {return 0.;}
231     virtual Double_t ParticleCharge(int) const {return 0.;}
232     virtual Double_t ParticleLifeTime(int) const {return 0.;}
233     virtual TMCParticleType ParticleMCType(int) const {return (TMCParticleType) 0;}
234 #endif
235     //
236   // Control Methods
237   virtual void          Init()                                                            { };
238   virtual void          FinishGeometry()                                                  { };
239   virtual void          BuildPhysics()                                                    { };
240   virtual void          ProcessEvent();
241 //change is in 4.0.5 (which is not yet tagged) (rdm)
242 #if ROOT_VERSION_CODE >= 262148
243   virtual Bool_t        ProcessRun(Int_t )                                                { return kTRUE; }
244 #else
245   virtual void          ProcessRun(Int_t )                                                { };
246 #endif
247   //virtual TMCGeomType   GetMCGeomType() const                                             { return kGeant3; }
248
249   // External Decayer
250   virtual void          SetExternalDecayer(AliDecayer* )                                  { };
251   virtual AliDecayer   *Decayer() const                                                   { return 0; };
252
253  protected:
254
255   enum {
256       kVolDrRg
257     , kVolAmRg
258     , kVolDrCh
259   };
260
261   Float_t          fMaxStep;            //  Maximum step size
262   Int_t            fNStep;              //  Number of steps
263   Int_t            fTrack;              //  Track number
264   Double_t         fTrackPx;            //  Track px
265   Double_t         fTrackPy;            //  Track py
266   Double_t         fTrackPz;            //  Track pz
267   Double_t         fTrackPtot;          //  Track total momentum
268   Double_t         fTrackEtot;          //  Track total energy
269   Double_t         fTrackX;             //  Track x position
270   Double_t         fTrackY;             //  Track y position
271   Double_t         fTrackZ;             //  Track z position
272   Double_t         fX0;                 //  X position of the beginning of the chamber
273   Double_t         fTrackStep;          //  Track step size
274   Int_t            fTrackPid;           //  Track PID
275   Float_t          fTrackCharge;        //  Track charge
276   Float_t          fTrackMass;          //  Track particle mass
277   Bool_t           fTrackEntering;      //  Track entering chamber
278
279   AliTRDv1        *fTRD;                //! TRD detector object
280   AliTRDparameter *fPar;                //! TRD parameter object
281
282   ClassDef(AliTRDsimpleMC,2)            //  Simple TRD Monte Carlo class
283  
284 };
285 #endif                                                                          
286