]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDsimpleMC.h
Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / TRD / AliTRDsimpleMC.h
CommitLineData
16bf9884 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
0a29d0f1 8///////////////////////////////////////////////////////////////////////////////
9// //
10// Simple TRD Monte Carlo class //
11// //
12///////////////////////////////////////////////////////////////////////////////
13
16bf9884 14#include "AliMC.h"
15#include "AliMCProcess.h"
16
17class AliTRDv1;
5443e65e 18class AliTRDparameter;
16bf9884 19
20class AliTRDsimpleMC : public AliMC {
21
22 public:
23
24 enum {
25 kPdgElectron = 11
26 , kPdgPion = 211
27 };
28
29 AliTRDsimpleMC();
30 AliTRDsimpleMC(const char *name, const char *title);
31 AliTRDsimpleMC(const AliTRDsimpleMC &m);
32
33 virtual ~AliTRDsimpleMC();
34 AliTRDsimpleMC &operator=(const AliTRDsimpleMC &m);
35
36 virtual void Copy(TObject &m);
37
38 //
39 // Methods for building / management of geometry
40 //
41
42 // Functions from GCONS
43 virtual void Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z,
44 Float_t &dens, Float_t &radl, Float_t &absl,
45 Float_t* ubuf, Int_t& nbuf) { };
46
47 // Detector composition
48 virtual void Material(Int_t& kmat, const char* name, Float_t a,
49 Float_t z, Float_t dens, Float_t radl, Float_t absl,
50 Float_t* buf, Int_t nwbuf) { };
51 virtual void Mixture(Int_t& kmat, const char *name, Float_t *a,
52 Float_t *z, Float_t dens, Int_t nlmat, Float_t *wmat) { };
53 virtual void Medium(Int_t& kmed, const char *name, Int_t nmat,
54 Int_t isvol, Int_t ifield, Float_t fieldm, Float_t tmaxfd,
55 Float_t stemax, Float_t deemax, Float_t epsil,
56 Float_t stmin, Float_t* ubuf, Int_t nbuf) { };
57 virtual void Matrix(Int_t& krot, Float_t thetaX, Float_t phiX,
58 Float_t thetaY, Float_t phiY, Float_t thetaZ,
59 Float_t phiZ) { };
60 virtual void Gstpar(Int_t itmed, const char *param, Float_t parval) { };
61
62 // Functions from GGEOM
63 virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed,
64 Float_t *upar, Int_t np) { return 0; };
65 virtual void Gsdvn(const char *name, const char *mother, Int_t ndiv,
66 Int_t iaxis) { };
67 virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv,
68 Int_t iaxis, Float_t c0i, Int_t numed) { };
69 virtual void Gsdvt(const char *name, const char *mother, Float_t step,
70 Int_t iaxis, Int_t numed, Int_t ndvmx) { };
71 virtual void Gsdvt2(const char *name, const char *mother, Float_t step,
72 Int_t iaxis, Float_t c0, Int_t numed, Int_t ndvmx) { };
73 virtual void Gsord(const char *name, Int_t iax) { };
74 virtual void Gspos(const char *name, Int_t nr, const char *mother,
75 Float_t x, Float_t y, Float_t z, Int_t irot,
76 const char *konly="ONLY") { };
77 virtual void Gsposp(const char *name, Int_t nr, const char *mother,
78 Float_t x, Float_t y, Float_t z, Int_t irot,
79 const char *konly, Float_t *upar, Int_t np) { };
b39d9131 80 virtual void Gsbool(const char* onlyVolName, const char* manyVolName) {}
16bf9884 81
82 virtual void SetCerenkov(Int_t itmed, Int_t npckov, Float_t *ppckov,
83 Float_t *absco, Float_t *effic, Float_t *rindex) { };
84
85 // Functions for drawing
86 virtual void DrawOneSpec(const char* name) { };
87 virtual void Gsatt(const char* name, const char* att, Int_t val) { };
88 virtual void Gdraw(const char*,Float_t theta = 30, Float_t phi = 30,
89 Float_t psi = 0, Float_t u0 = 10, Float_t v0 = 10,
90 Float_t ul = 0.01, Float_t vl = 0.01) { };
91
92 // Euclid
93 virtual void WriteEuclid(const char *a, const char *b, Int_t c, Int_t d) { };
94
95 // Get methods
96 virtual Int_t VolId(const Text_t* volName) const;
97 virtual const char *VolName(Int_t id) const { return " "; };
98 virtual Int_t NofVolumes() const { return 0; };
99 virtual Int_t VolId2Mate(Int_t id) const { return 0; };
100
101 //
102 // Methods for physics management
103 //
104
105 // Set methods
106 virtual void SetCut(const char* cutName, Float_t cutValue) { };
107 virtual void SetProcess(const char* flagName, Int_t flagValue) { };
108 virtual Float_t Xsec(char *a, Float_t b, Int_t c, Int_t d) { return 0.0; };
109
110 // Particle table usage
111 virtual Int_t IdFromPDG(Int_t id) const { return 0; };
112 virtual Int_t PDGFromId(Int_t pdg) const { return 0; };
113 virtual void DefineParticles() { };
114
115 //
116 // Methods for step management
117 //
118
119 // Action methods
120 virtual void StopTrack() { };
121 virtual void StopEvent() { };
122
123 // Set methods
124 virtual void SetMaxStep(Float_t step) { fMaxStep = step; };
125 virtual void SetMaxNStep(Int_t n) { };
126 virtual void SetUserDecay(Int_t d) { };
127
128 virtual void NewTrack(Int_t iTrack, Int_t pdg, Double_t px, Double_t py, Double_t pz);
129
130 // Tracking volume(s)
131 virtual Int_t CurrentVolID(Int_t& copyNo) const;
132 virtual Int_t CurrentVolOffID(Int_t off, Int_t& copyNo) const;
133 virtual const char* CurrentVolName() const { return ""; };
134 virtual const char* CurrentVolOffName(Int_t off) const { return ""; };
135 virtual Int_t CurrentMaterial(Float_t &a, Float_t &z,
136 Float_t &dens, Float_t &radl,
137 Float_t &absl) const { return 0; };
138 virtual Int_t CurrentEvent() const { return 0; };
139 virtual void Gmtod(Float_t* xm, Float_t* xd, Int_t iflag) { };
140 virtual void Gdtom(Float_t* xd, Float_t* xm, Int_t iflag) { } ;
141 virtual Float_t MaxStep() const { return fMaxStep; };
142 virtual Int_t GetNStep() const { return fNStep; };
143 virtual Int_t GetMaxNStep() const { return 0; };
144 virtual Int_t GetMedium() const { return 0; };
145
146 // Dynamic properties
147 virtual void TrackPosition(TLorentzVector& position) const;
148 virtual void TrackMomentum(TLorentzVector& momentum) const;
149 virtual Float_t TrackStep() const { return fTrackStep; };
150 virtual Float_t TrackLength() const { return 0.0; };
151 virtual Float_t TrackTime() const { return 0.0; };
152 virtual Float_t Edep() const { return 0.0; };
153
154 // Static properties
155 virtual Int_t TrackPid() const { return fTrackPid; };
156 virtual Float_t TrackCharge() const { return fTrackCharge; };
157 virtual Float_t TrackMass() const { return fTrackMass; };
158 virtual Float_t Etot() const { return fTrackEtot; };
159
160 // Track status
161 virtual Bool_t IsNewTrack() const { return kFALSE; };
162 virtual Bool_t IsTrackInside() const { return kFALSE; };
163 virtual Bool_t IsTrackEntering() const { return fTrackEntering; };
164 virtual Bool_t IsTrackExiting() const { return kFALSE; };
165 virtual Bool_t IsTrackOut() const { return kFALSE; };
166 virtual Bool_t IsTrackDisappeared() const { return kFALSE; };
167 virtual Bool_t IsTrackStop() const { return kFALSE; };
168 virtual Bool_t IsTrackAlive() const { return kFALSE; };
169
170 // Secondaries
171 virtual Int_t NSecondaries() const { return 0; };
172 virtual void GetSecondary(Int_t isec, Int_t& particleId,
173 TLorentzVector& position,
174 TLorentzVector& momentum) { };
175 virtual AliMCProcess ProdProcess(Int_t isec) const { return kPNoProcess; };
176 virtual Int_t StepProcesses(TArrayI &proc) const { return 0; };
177
178 //
179 // Other (then geometry/step/run management) methods
180 //
181
182 // Geant3 specific methods
183 virtual void Gdopt(const char *c1,const char*c2) { };
184 virtual void SetClipBox(const char* cc,Float_t a=-9999,Float_t b=0,
185 Float_t c=-9999,Float_t d=0,
186 Float_t e=-9999,Float_t f=0) { };
187 virtual void DefaultRange() { };
188 virtual void Gdhead(Int_t, const char *c, Float_t a=0) { };
189 virtual void Gdman(Float_t a, Float_t b, const char *c) { };
190 virtual void SetColors() { };
191 virtual void Gtreve() { };
192 virtual void GtreveRoot() { };
193 virtual void Gckmat(Int_t, char*) { };
194 virtual void InitLego() { };
195 virtual void Gfpart(Int_t a, char *b, Int_t &c, Float_t &d,
196 Float_t &e, Float_t &g) { };
197 virtual void Gspart(Int_t a, const char *b, Int_t c, Float_t d,
198 Float_t e, Float_t g) { };
199
200 // Control Methods
201 virtual void Init() { };
202 virtual void FinishGeometry() { };
203 virtual void BuildPhysics() { };
204 virtual void ProcessEvent();
205 virtual void ProcessRun(Int_t nevent) { };
b39d9131 206 virtual AliMCGeomType GetMCGeomType() const { return kGeant3; }
16bf9884 207
208 // External Decayer
209 virtual void SetExternalDecayer(AliDecayer* decayer) { };
210 virtual AliDecayer *Decayer() const { return 0; };
211
212 protected:
213
214 enum {
215 kVolDrRg
216 , kVolAmRg
217 , kVolDrCh
218 };
219
5443e65e 220 Float_t fMaxStep; // Maximum step size
221 Int_t fNStep; // Number of steps
222 Int_t fTrack; // Track number
223 Double_t fTrackPx; // Track px
224 Double_t fTrackPy; // Track py
225 Double_t fTrackPz; // Track pz
226 Double_t fTrackPtot; // Track total momentum
227 Double_t fTrackEtot; // Track total energy
228 Double_t fTrackX; // Track x position
229 Double_t fTrackY; // Track y position
230 Double_t fTrackZ; // Track z position
231 Double_t fX0; // X position of the beginning of the chamber
232 Double_t fTrackStep; // Track step size
233 Int_t fTrackPid; // Track PID
234 Float_t fTrackCharge; // Track charge
235 Float_t fTrackMass; // Track particle mass
236 Bool_t fTrackEntering; // Track entering chamber
237
238 AliTRDv1 *fTRD; //! TRD detector object
239 AliTRDparameter *fPar; //! TRD parameter object
240
241 ClassDef(AliTRDsimpleMC,2) // Simple TRD Monte Carlo class
16bf9884 242
243};
244#endif