]>
Commit | Line | Data |
---|---|---|
2817d3e2 | 1 | // $Id$ |
2 | // Category: run | |
3 | // | |
e5967ab3 | 4 | // Author: I. Hrivnacova |
5 | // | |
6 | // Class TGeant4 | |
7 | // ------------- | |
8 | // Geant4 implementation of the MonteCarlo interface. | |
2817d3e2 | 9 | |
10 | #ifndef TGEANT4_H | |
11 | #define TGEANT4_H | |
12 | ||
13 | #include "AliMC.h" | |
9ee4729b | 14 | #include "AliMCProcess.h" |
2817d3e2 | 15 | |
277e1919 | 16 | #include <TArrayI.h> |
17 | ||
2817d3e2 | 18 | class TG4VRunConfiguration; |
19 | class TG4GeometryManager; | |
2da53674 | 20 | class TG4SDManager; |
2817d3e2 | 21 | class TG4PhysicsManager; |
22 | class TG4StepManager; | |
23 | class TG4VisManager; | |
24 | class TG4RunManager; | |
25 | class TG4Messenger; | |
26 | ||
89381a6c | 27 | class AliDecayer; |
28 | ||
2817d3e2 | 29 | class TGeant4: public AliMC |
30 | { | |
31 | public: | |
32 | TGeant4(const char* name, const char* title, | |
33 | TG4VRunConfiguration* configuration, int argc, char** argv); | |
34 | TGeant4(const char* name, const char* title, | |
35 | TG4VRunConfiguration* configuration); | |
36 | // --> protected | |
37 | // TGeant4(); | |
38 | // TGeant4(const TGeant4& right); | |
39 | virtual ~TGeant4(); | |
40 | ||
41 | // | |
42 | // methods for building/management of geometry | |
43 | // ------------------------------------------------ | |
44 | // | |
45 | ||
46 | // functions from GBASE | |
47 | virtual void FinishGeometry(); | |
48 | //Ggclos(); | |
49 | ||
50 | // functions from GCONS | |
51 | virtual void Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z, | |
52 | Float_t &dens, Float_t &radl, Float_t &absl, | |
53 | Float_t* ubuf, Int_t& nbuf); | |
54 | ||
55 | // detector composition | |
56 | virtual void Material(Int_t& kmat, const char* name, Float_t a, | |
57 | Float_t z, Float_t dens, Float_t radl, Float_t absl, | |
58 | Float_t* buf, Int_t nwbuf); | |
59 | virtual void Mixture(Int_t& kmat, const char *name, Float_t *a, | |
60 | Float_t *z, Float_t dens, Int_t nlmat, Float_t *wmat); | |
61 | virtual void Medium(Int_t& kmed, const char *name, Int_t nmat, | |
62 | Int_t isvol, Int_t ifield, Float_t fieldm, Float_t tmaxfd, | |
63 | Float_t stemax, Float_t deemax, Float_t epsil, | |
64 | Float_t stmin, Float_t* ubuf, Int_t nbuf); | |
3859747a | 65 | virtual void Matrix(Int_t& krot, Double_t thetaX, Double_t phiX, |
66 | Double_t thetaY, Double_t phiY, Double_t thetaZ, | |
67 | Double_t phiZ); | |
2817d3e2 | 68 | virtual void Matrix(Int_t& krot, Float_t thetaX, Float_t phiX, |
69 | Float_t thetaY, Float_t phiY, Float_t thetaZ, | |
70 | Float_t phiZ); | |
75578018 | 71 | virtual void SetCerenkov(Int_t itmed, Int_t npckov, Float_t *ppckov, |
72 | Float_t *absco, Float_t *effic, Float_t *rindex); | |
2817d3e2 | 73 | |
74 | // functions from GGEOM | |
3859747a | 75 | virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed, |
76 | Double_t *upar, Int_t np); | |
2817d3e2 | 77 | virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed, |
78 | Float_t *upar, Int_t np); | |
79 | virtual void Gsdvn(const char *name, const char *mother, Int_t ndiv, | |
80 | Int_t iaxis); | |
3859747a | 81 | virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv, |
82 | Int_t iaxis, Double_t c0i, Int_t numed); | |
2817d3e2 | 83 | virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv, |
84 | Int_t iaxis, Float_t c0i, Int_t numed); | |
3859747a | 85 | virtual void Gsdvt(const char *name, const char *mother, Double_t step, |
86 | Int_t iaxis, Int_t numed, Int_t ndvmx); | |
2817d3e2 | 87 | virtual void Gsdvt(const char *name, const char *mother, Float_t step, |
88 | Int_t iaxis, Int_t numed, Int_t ndvmx); | |
3859747a | 89 | virtual void Gsdvt2(const char *name, const char *mother, Double_t step, |
90 | Int_t iaxis, Double_t c0, Int_t numed, Int_t ndvmx); | |
2817d3e2 | 91 | virtual void Gsdvt2(const char *name, const char *mother, Float_t step, |
92 | Int_t iaxis, Float_t c0, Int_t numed, Int_t ndvmx); | |
93 | virtual void Gsord(const char *name, Int_t iax); | |
3859747a | 94 | virtual void Gspos(const char *name, Int_t nr, const char *mother, |
95 | Double_t x, Double_t y, Double_t z, Int_t irot, | |
96 | const char *konly="ONLY"); | |
2817d3e2 | 97 | virtual void Gspos(const char *name, Int_t nr, const char *mother, |
98 | Float_t x, Float_t y, Float_t z, Int_t irot, | |
99 | const char *konly); | |
3859747a | 100 | virtual void Gsposp(const char *name, Int_t nr, const char *mother, |
101 | Double_t x, Double_t y, Double_t z, Int_t irot, | |
102 | const char *konly, Double_t *upar, Int_t np); | |
2817d3e2 | 103 | virtual void Gsposp(const char *name, Int_t nr, const char *mother, |
104 | Float_t x, Float_t y, Float_t z, Int_t irot, | |
105 | const char *konly, Float_t *upar, Int_t np); | |
106 | ||
107 | // Euclid | |
108 | virtual void WriteEuclid(const char* fileName, const char* topVol, | |
109 | Int_t number, Int_t nlevel); | |
110 | ||
111 | // get methods | |
112 | virtual Int_t VolId(const Text_t* volName) const; | |
113 | virtual const char* VolName(Int_t id) const; | |
114 | virtual Int_t NofVolumes() const; | |
277e1919 | 115 | virtual Int_t VolId2Mate(Int_t id) const; |
2817d3e2 | 116 | |
117 | // | |
118 | // methods for physics management | |
119 | // ------------------------------------------------ | |
120 | // | |
121 | ||
122 | virtual void BuildPhysics(); | |
2da53674 | 123 | virtual void Gstpar(Int_t itmed, const char *param, Float_t parval); |
2817d3e2 | 124 | |
125 | // set methods | |
126 | virtual void SetCut(const char* cutName, Float_t cutValue); | |
127 | virtual void SetProcess(const char* flagName, Int_t flagValue); | |
128 | virtual Float_t Xsec(char* reac, Float_t energy, Int_t part, Int_t mate); | |
89381a6c | 129 | virtual void SetExternalDecayer(AliDecayer* decayer); //NEW |
130 | ||
131 | // get methods | |
132 | virtual AliDecayer* Decayer() const; //NEW | |
2817d3e2 | 133 | |
134 | // particle table usage | |
135 | virtual Int_t IdFromPDG(Int_t pdgID) const; | |
136 | virtual Int_t PDGFromId(Int_t mcID) const; | |
137 | virtual void DefineParticles(); | |
138 | ||
139 | // | |
140 | // methods for step management | |
141 | // ------------------------------------------------ | |
142 | // | |
143 | ||
144 | // action methods | |
145 | virtual void StopTrack(); | |
146 | virtual void StopEvent(); | |
147 | ||
148 | // set methods | |
149 | virtual void SetMaxStep(Float_t); | |
150 | virtual void SetMaxNStep(Int_t); | |
151 | virtual void SetUserDecay(Int_t); //NEW | |
152 | ||
153 | // get methods | |
154 | // tracking volume(s) | |
155 | virtual Int_t CurrentVolID(Int_t& copyNo) const; | |
156 | virtual Int_t CurrentVolOffID(Int_t off, Int_t& copyNo) const; | |
157 | virtual const char* CurrentVolName() const; | |
158 | virtual const char* CurrentVolOffName(Int_t off) const; | |
159 | virtual Int_t CurrentMaterial(Float_t &a, Float_t &z, | |
160 | Float_t &dens, Float_t &radl, Float_t &absl) const; | |
161 | virtual void Gmtod(Float_t* xm, Float_t* xd, Int_t iflag); | |
162 | virtual void Gdtom(Float_t* xd, Float_t* xm, Int_t iflag); | |
163 | virtual Float_t MaxStep() const; | |
164 | virtual Int_t GetMaxNStep() const; | |
165 | virtual Int_t GetMedium() const; | |
166 | ||
167 | // tracking particle | |
168 | // dynamic properties | |
169 | virtual void TrackPosition(TLorentzVector& position) const; | |
170 | virtual void TrackMomentum(TLorentzVector& momentum) const; | |
171 | virtual void TrackVertexPosition(TLorentzVector& position) const; | |
172 | virtual void TrackVertexMomentum(TLorentzVector& momentum) const; | |
173 | virtual Float_t TrackStep() const; | |
174 | virtual Float_t TrackLength() const; | |
175 | virtual Float_t TrackTime() const; | |
176 | virtual Float_t Edep() const; | |
177 | // static properties | |
178 | virtual Int_t TrackPid() const; | |
179 | virtual Float_t TrackCharge() const; | |
180 | virtual Float_t TrackMass() const; | |
181 | virtual Float_t Etot() const; | |
182 | ||
183 | // track status | |
184 | virtual Bool_t IsTrackInside() const; | |
185 | virtual Bool_t IsTrackEntering() const; | |
186 | virtual Bool_t IsTrackExiting() const; | |
187 | virtual Bool_t IsTrackOut() const; | |
188 | virtual Bool_t IsTrackDisappeared() const; | |
189 | virtual Bool_t IsTrackStop() const; | |
190 | virtual Bool_t IsTrackAlive() const; | |
191 | virtual Bool_t IsNewTrack() const; | |
192 | ||
193 | // secondaries | |
194 | virtual Int_t NSecondaries() const; | |
195 | virtual void GetSecondary(Int_t isec, Int_t& particleId, | |
196 | TLorentzVector& position, TLorentzVector& momentum); | |
277e1919 | 197 | virtual AliMCProcess ProdProcess(Int_t isec) const; |
198 | virtual Int_t StepProcesses(TArrayI &proc) const; | |
2817d3e2 | 199 | |
200 | // | |
201 | // methods for visualization | |
202 | // ------------------------------------------------ | |
203 | // | |
204 | // functions for drawing | |
205 | virtual void DrawOneSpec(const char* name); | |
206 | virtual void Gsatt(const char* name, const char* att, Int_t val); | |
207 | virtual void Gdraw(const char* name, Float_t theta, Float_t phi, | |
208 | Float_t psi, Float_t u0, Float_t v0, | |
209 | Float_t ul, Float_t vl); | |
210 | ||
211 | // | |
212 | // NEW | |
213 | // Geant3 specific methods | |
214 | // !!! need to be transformed to common interface | |
215 | // | |
216 | virtual void Gdopt(const char* name , const char* value); | |
217 | virtual void SetClipBox(const char *name, Float_t xmin, Float_t xmax, | |
218 | Float_t ymin, Float_t ymax, Float_t zmin, Float_t zmax); | |
219 | virtual void DefaultRange(); | |
220 | virtual void Gdhead(Int_t isel, const char* name, Float_t chrsiz); | |
221 | virtual void Gdman(Float_t u, Float_t v, const char* type); | |
222 | virtual void SetColors(); | |
223 | virtual void Gtreve(); | |
aee8290b | 224 | virtual void GtreveRoot(); |
2817d3e2 | 225 | virtual void Gckmat(Int_t itmed, char* natmed); |
226 | virtual void InitLego(); | |
227 | virtual void Gfpart(Int_t ipart, char *name, Int_t& itrtyp, | |
228 | Float_t& amass, Float_t& charge, Float_t& tlife); | |
229 | virtual void Gspart(Int_t ipart, const char *name, Int_t itrtyp, | |
230 | Float_t amass, Float_t charge, Float_t tlife); | |
231 | ||
232 | // | |
233 | // methods for run control | |
234 | // ------------------------------------------------ | |
235 | // | |
236 | ||
237 | virtual void Init(); | |
238 | virtual void ProcessEvent(); | |
239 | virtual void ProcessRun(Int_t nofEvents); | |
240 | ||
241 | // UI control methods | |
242 | void StartGeantUI(); | |
243 | void StartRootUI(); | |
244 | void ProcessGeantMacro(const char* macroName); | |
245 | void ProcessGeantCommand(const char* commandPath); | |
246 | ||
247 | // get methods | |
248 | virtual Int_t CurrentEvent() const; | |
249 | ||
250 | protected: | |
251 | TGeant4(); | |
252 | TGeant4(const TGeant4& right); | |
253 | ||
254 | // operators | |
255 | TGeant4& operator=(const TGeant4& right); | |
256 | ||
257 | private: | |
258 | // data members | |
259 | TG4GeometryManager* fGeometryManager; //geometry manager | |
2da53674 | 260 | TG4SDManager* fSDManager; //sensitive detectors manager |
2817d3e2 | 261 | TG4PhysicsManager* fPhysicsManager; //physics manager |
262 | TG4StepManager* fStepManager; //step manager | |
263 | TG4VisManager* fVisManager; //visualization manager | |
264 | TG4RunManager* fRunManager; //run manager | |
265 | TG4Messenger* fMessenger; //messenger | |
266 | }; | |
267 | ||
268 | #ifndef __CINT__ | |
269 | ||
270 | // inline methods | |
271 | #include "TGeant4.icc" | |
272 | ||
273 | #endif | |
274 | #endif // TGEANT4_H | |
275 |