]>
Commit | Line | Data |
---|---|---|
2817d3e2 | 1 | // $Id$ |
2 | // Category: run | |
3 | // | |
4 | // See the class description in the header file. | |
5 | ||
6 | #include "TGeant4.h" | |
7 | #include "TG4Messenger.h" | |
8 | #include "TG4GeometryManager.h" | |
9 | #include "TG4PhysicsManager.h" | |
10 | #include "TG4StepManager.h" | |
11 | #include "TG4VisManager.h" | |
12 | #include "TG4RunManager.h" | |
13 | #include "TG4Globals.h" | |
14 | ||
e59f13eb | 15 | #include "AliDecayer.h" |
16 | ||
2817d3e2 | 17 | TGeant4::TGeant4(const char* name, const char* title, |
18 | TG4VRunConfiguration* configuration, int argc, char** argv) | |
19 | : AliMC(name, title), | |
20 | fVisManager(0) | |
21 | { | |
22 | // create run manager | |
23 | fRunManager = new TG4RunManager(configuration, argc, argv); | |
24 | // add verbose level | |
25 | //G4cout << "TG4RunManager has been created." << endl; | |
26 | ||
27 | // create geometry manager | |
28 | fGeometryManager = new TG4GeometryManager(); | |
29 | // add verbose level | |
30 | //G4cout << "TG4GeometryManager has been created." << endl; | |
31 | ||
32 | // create physics manager | |
33 | fPhysicsManager = new TG4PhysicsManager(); | |
34 | // add verbose level | |
35 | //G4cout << "TG4GeometryManager has been created." << endl; | |
36 | ||
37 | // create step manager | |
38 | fStepManager = new TG4StepManager(); | |
39 | // add verbose level | |
40 | //G4cout << "TG4StepManager has been created." << endl; | |
41 | ||
42 | #ifdef G4VIS_USE | |
43 | // create visualization manager | |
44 | fVisManager = new TG4VisManager(); | |
45 | fVisManager->Initialize(); | |
46 | #endif | |
47 | ||
48 | // create messenger | |
49 | fMessenger = | |
50 | new TG4Messenger(fGeometryManager, fPhysicsManager, fStepManager); | |
51 | } | |
52 | ||
53 | TGeant4::TGeant4(const char* name, const char* title, | |
54 | TG4VRunConfiguration* configuration) | |
55 | : AliMC(name, title), | |
56 | fVisManager(0) | |
57 | { | |
58 | // create run manager | |
59 | fRunManager = new TG4RunManager(configuration); | |
60 | // add verbose level | |
61 | //G4cout << "TG4RunManager has been created." << endl; | |
62 | ||
63 | // create geometry manager | |
64 | fGeometryManager = new TG4GeometryManager(); | |
65 | // add verbose level | |
66 | //G4cout << "TG4GeometryManager has been created." << endl; | |
67 | ||
68 | // create physics manager | |
69 | fPhysicsManager = new TG4PhysicsManager(); | |
70 | // add verbose level | |
71 | //G4cout << "TG4GeometryManager has been created." << endl; | |
72 | ||
73 | // create step manager | |
74 | fStepManager = new TG4StepManager(); | |
75 | // add verbose level | |
76 | //G4cout << "TG4StepManager has been created." << endl; | |
77 | ||
78 | #ifdef G4VIS_USE | |
79 | // create visualization manager | |
80 | fVisManager = new TG4VisManager(); | |
81 | fVisManager->Initialize(); | |
82 | #endif | |
83 | ||
84 | // create messenger | |
85 | fMessenger = | |
86 | new TG4Messenger(fGeometryManager, fPhysicsManager, fStepManager); | |
87 | } | |
88 | ||
89 | TGeant4::TGeant4() { | |
90 | // | |
91 | } | |
92 | ||
93 | TGeant4::TGeant4(const TGeant4& right) { | |
94 | // | |
95 | TG4Globals::Exception("TGeant4 is protected from copying."); | |
96 | } | |
97 | ||
98 | TGeant4::~TGeant4() { | |
99 | // | |
100 | delete fRunManager; | |
101 | delete fGeometryManager; | |
102 | delete fPhysicsManager; | |
103 | delete fStepManager; | |
104 | // fVisManager is deleted with G4RunManager destructor | |
105 | delete fMessenger; | |
106 | } | |
107 | ||
108 | // operators | |
109 | ||
110 | TGeant4& TGeant4::operator=(const TGeant4& right) | |
111 | { | |
112 | // check assignement to self | |
113 | if (this == &right) return *this; | |
114 | ||
115 | TG4Globals::Exception("TGeant4 is protected from assigning."); | |
116 | ||
117 | return *this; | |
118 | } | |
119 | ||
120 | ||
121 | // methods for building/management of geometry | |
122 | // ------------------------------------------------ | |
123 | ||
124 | void TGeant4::FinishGeometry() { | |
125 | // | |
126 | fGeometryManager->Ggclos(); | |
127 | } | |
128 | ||
129 | void TGeant4::Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z, | |
130 | Float_t &dens, Float_t &radl, Float_t &absl, | |
131 | Float_t* ubuf, Int_t& nbuf) { | |
132 | // | |
133 | fGeometryManager | |
134 | ->Gfmate(imat, name, a, z, dens, radl, absl, ubuf, nbuf); | |
135 | } | |
136 | ||
137 | void TGeant4::Material(Int_t& kmat, const char* name, Float_t a, | |
138 | Float_t z, Float_t dens, Float_t radl, Float_t absl, | |
139 | Float_t* buf, Int_t nwbuf) { | |
140 | // | |
141 | fGeometryManager | |
142 | ->Material(kmat, name, a, z, dens, radl, absl, buf, nwbuf); | |
143 | } | |
144 | ||
145 | void TGeant4::Mixture(Int_t& kmat, const char *name, Float_t *a, | |
146 | Float_t *z, Float_t dens, Int_t nlmat, Float_t *wmat) { | |
147 | // | |
148 | fGeometryManager | |
149 | ->Mixture(kmat, name, a, z, dens, nlmat, wmat); | |
150 | } | |
151 | ||
152 | void TGeant4::Medium(Int_t& kmed, const char *name, Int_t nmat, | |
153 | Int_t isvol, Int_t ifield, Float_t fieldm, Float_t tmaxfd, | |
154 | Float_t stemax, Float_t deemax, Float_t epsil, | |
155 | Float_t stmin, Float_t* ubuf, Int_t nbuf) { | |
156 | // | |
157 | fGeometryManager | |
158 | ->Medium(kmed, name, nmat, isvol, ifield, fieldm, tmaxfd, stemax, deemax, | |
159 | epsil, stmin, ubuf, nbuf); | |
160 | } | |
161 | ||
162 | void TGeant4::Matrix(Int_t& krot, Float_t thetaX, Float_t phiX, | |
163 | Float_t thetaY, Float_t phiY, Float_t thetaZ, | |
164 | Float_t phiZ) { | |
165 | // | |
166 | fGeometryManager | |
167 | ->Matrix(krot, thetaX, phiX, thetaY, phiY, thetaZ, phiZ); | |
168 | } | |
169 | ||
170 | void TGeant4::Gstpar(Int_t itmed, const char *param, Float_t parval) { | |
171 | // | |
172 | fGeometryManager->Gstpar(itmed, param, parval); | |
173 | } | |
174 | ||
2817d3e2 | 175 | Int_t TGeant4::Gsvolu(const char *name, const char *shape, Int_t nmed, |
176 | Float_t *upar, Int_t np) { | |
177 | // | |
178 | return fGeometryManager->Gsvolu(name, shape, nmed, upar, np); | |
179 | } | |
180 | ||
181 | void TGeant4::Gsdvn(const char *name, const char *mother, Int_t ndiv, | |
182 | Int_t iaxis) { | |
183 | // | |
184 | fGeometryManager->Gsdvn(name, mother, ndiv, iaxis); | |
185 | } | |
186 | ||
187 | void TGeant4::Gsdvn2(const char *name, const char *mother, Int_t ndiv, | |
188 | Int_t iaxis, Float_t c0i, Int_t numed) { | |
189 | // | |
190 | fGeometryManager->Gsdvn2(name, mother, ndiv, iaxis, c0i, numed); | |
191 | } | |
192 | ||
193 | void TGeant4::Gsdvt(const char *name, const char *mother, Float_t step, | |
194 | Int_t iaxis, Int_t numed, Int_t ndvmx) { | |
195 | // | |
196 | fGeometryManager->Gsdvt(name, mother, step, iaxis, numed, ndvmx); | |
197 | } | |
198 | ||
199 | void TGeant4::Gsdvt2(const char *name, const char *mother, Float_t step, | |
200 | Int_t iaxis, Float_t c0, Int_t numed, Int_t ndvmx) { | |
201 | // | |
202 | fGeometryManager->Gsdvt2(name, mother, step, iaxis, c0, numed, ndvmx); | |
203 | } | |
204 | ||
205 | void TGeant4::Gsord(const char *name, Int_t iax) { | |
206 | // | |
207 | fGeometryManager->Gsord(name, iax); | |
208 | } | |
209 | ||
210 | void TGeant4::Gspos(const char *name, Int_t nr, const char *mother, | |
211 | Float_t x, Float_t y, Float_t z, Int_t irot, | |
212 | const char *konly) { | |
213 | // | |
214 | fGeometryManager->Gspos(name, nr, mother, x, y, z, irot, konly); | |
215 | } | |
216 | ||
217 | void TGeant4::Gsposp(const char *name, Int_t nr, const char *mother, | |
218 | Float_t x, Float_t y, Float_t z, Int_t irot, | |
219 | const char *konly, Float_t *upar, Int_t np) { | |
220 | // | |
221 | fGeometryManager->Gsposp(name, nr, mother, x, y, z, irot, konly, upar, np); | |
222 | } | |
223 | ||
75578018 | 224 | void TGeant4::SetCerenkov(Int_t itmed, Int_t npckov, Float_t *ppckov, |
225 | Float_t *absco, Float_t *effic, Float_t *rindex) { | |
226 | // | |
227 | fGeometryManager->SetCerenkov(itmed, npckov, ppckov, absco, effic, rindex); | |
228 | } | |
229 | ||
2817d3e2 | 230 | void TGeant4::WriteEuclid(const char* fileName, const char* topVol, |
231 | Int_t number, Int_t nlevel) { | |
232 | // | |
233 | fGeometryManager->WriteEuclid(fileName, topVol, number, nlevel); | |
234 | } | |
235 | ||
236 | Int_t TGeant4::VolId(const Text_t* volName) const { | |
237 | // | |
238 | return fGeometryManager->VolId(volName); | |
239 | } | |
240 | ||
241 | const char* TGeant4::VolName(Int_t id) const { | |
242 | // | |
243 | return fGeometryManager->VolName(id); | |
244 | } | |
245 | ||
246 | Int_t TGeant4::NofVolumes() const { | |
247 | // | |
248 | return fGeometryManager->NofVolumes(); | |
249 | } | |
250 | ||
fdd271a5 | 251 | Int_t TGeant4::VolId2Mate(Int_t id) const { |
252 | // | |
253 | return fGeometryManager->VolId2Mate(id); | |
254 | } | |
255 | ||
256 | ||
257 | ||
2817d3e2 | 258 | // methods for physics management |
259 | // ------------------------------------------------ | |
260 | ||
261 | void TGeant4::BuildPhysics() { | |
262 | // | |
263 | fPhysicsManager->BuildPhysics(); | |
264 | } | |
265 | ||
266 | void TGeant4::SetCut(const char* cutName, Float_t cutValue) { | |
267 | // | |
268 | fPhysicsManager->SetCut(cutName, cutValue); | |
269 | } | |
270 | ||
271 | void TGeant4::SetProcess(const char* flagName, Int_t flagValue) { | |
272 | // | |
273 | fPhysicsManager->SetProcess(flagName, flagValue); | |
274 | } | |
275 | ||
276 | Float_t TGeant4::Xsec(char* reac, Float_t energy, Int_t part, Int_t mate) { | |
277 | // | |
278 | return fPhysicsManager->Xsec(reac, energy, part, mate); | |
279 | } | |
280 | ||
e59f13eb | 281 | void TGeant4::SetExternalDecayer(AliDecayer* decayer) { |
282 | // | |
283 | TG4Globals:: Warning("TGeant4::SetExternalDecayer() is not implemented."); | |
284 | } | |
285 | ||
286 | AliDecayer* TGeant4::Decayer() const { | |
287 | // | |
288 | TG4Globals:: Warning("TGeant4::Decayer() is not implemented."); | |
289 | return 0; | |
290 | } | |
291 | ||
292 | ||
2817d3e2 | 293 | Int_t TGeant4::IdFromPDG(Int_t pdgID) const { |
294 | // | |
295 | return fPhysicsManager->IdFromPDG(pdgID); | |
296 | } | |
297 | ||
298 | Int_t TGeant4::PDGFromId(Int_t mcID) const { | |
299 | // | |
300 | return fPhysicsManager->PDGFromId(mcID); | |
301 | } | |
302 | ||
303 | void TGeant4::DefineParticles() { | |
304 | // | |
305 | fPhysicsManager->DefineParticles(); | |
306 | } | |
307 | ||
308 | // methods for step management | |
309 | // ------------------------------------------------ | |
310 | // inlined (in TGeant4.icc) | |
311 | ||
2817d3e2 | 312 | // methods for visualization |
313 | // ------------------------------------------------ | |
314 | ||
315 | #ifdef G4VIS_USE | |
316 | void TGeant4::DrawOneSpec(const char* name) { | |
317 | // | |
318 | fVisManager->DrawOneSpec(name); | |
319 | } | |
320 | ||
321 | void TGeant4::Gsatt(const char* name, const char* att, Int_t val) { | |
322 | // | |
323 | fVisManager->Gsatt(name, att, val); | |
324 | } | |
325 | ||
326 | void TGeant4::Gdraw(const char* name, Float_t theta, Float_t phi, | |
327 | Float_t psi, Float_t u0, Float_t v0, | |
328 | Float_t ul, Float_t vl) { | |
329 | // | |
330 | fVisManager->Gdraw(name, theta, phi, psi, u0, v0, ul, vl); | |
331 | } | |
332 | ||
333 | #else | |
334 | void TGeant4::DrawOneSpec(const char* name) { | |
335 | // | |
336 | TG4Globals:: Warning("TGeant4::DrawOneSpec(): no visualization available."); | |
337 | } | |
338 | ||
339 | void TGeant4::Gsatt(const char* name, const char* att, Int_t val) { | |
340 | // | |
341 | TG4Globals:: Warning("TGeant4::Gsatt(): no visualization available."); | |
342 | } | |
343 | ||
344 | void TGeant4::Gdraw(const char* p1, Float_t theta, Float_t phi, | |
345 | Float_t psi, Float_t u0, Float_t v0, | |
346 | Float_t ul, Float_t vl) { | |
347 | // | |
348 | TG4Globals:: Warning("TGeant4::Gdraw(): no visualization available."); | |
349 | } | |
350 | ||
351 | #endif //G4VIS_USE | |
352 | ||
353 | // methods for run control | |
354 | // ------------------------------------------------ | |
355 | ||
356 | void TGeant4::Init() { | |
357 | // | |
358 | fRunManager->Initialize(); | |
359 | } | |
360 | ||
361 | void TGeant4::ProcessEvent() { | |
362 | // | |
363 | fRunManager->ProcessEvent(); | |
364 | } | |
365 | ||
366 | void TGeant4::ProcessRun(Int_t nofEvents) { | |
367 | // | |
368 | fRunManager->ProcessRun(nofEvents); | |
369 | } | |
370 | ||
371 | void TGeant4::StartGeantUI() { | |
372 | // | |
373 | fRunManager->StartGeantUI(); | |
374 | } | |
375 | ||
376 | void TGeant4::StartRootUI() { | |
377 | // | |
562476b1 | 378 | fRunManager->StartRootUI(); |
2817d3e2 | 379 | } |
380 | ||
381 | void TGeant4::ProcessGeantMacro(const char* macroName) { | |
382 | // | |
383 | fRunManager->ProcessGeantMacro(macroName); | |
384 | } | |
385 | ||
386 | void TGeant4::ProcessGeantCommand(const char* command) { | |
387 | // | |
388 | fRunManager->ProcessGeantCommand(command); | |
389 | } | |
390 | ||
391 | Int_t TGeant4::CurrentEvent() const { | |
392 | // | |
393 | return fRunManager->CurrentEvent(); | |
394 | } | |
395 | ||
396 | // Geant3 specific methods | |
397 | // !!! need to be transformed to common interface | |
398 | // ------------------------------------------------ | |
399 | ||
400 | void TGeant4::Gdopt(const char* name, const char* value) { | |
401 | // | |
402 | TG4Globals:: Warning("TGeant4::Gdopt(..) is not implemented."); | |
403 | } | |
404 | ||
405 | void TGeant4::SetClipBox(const char *name, Float_t xmin, Float_t xmax, | |
406 | Float_t ymin, Float_t ymax, Float_t zmin, Float_t zmax) { | |
407 | // | |
408 | TG4Globals:: Warning("TGeant4::SetClipBox(..) is not implemented."); | |
409 | } | |
410 | ||
411 | void TGeant4::DefaultRange() { | |
412 | // | |
413 | TG4Globals:: Warning("TGeant4::DefaultRange() is not implemented."); | |
414 | } | |
415 | ||
416 | void TGeant4::Gdhead(Int_t isel, const char* name, Float_t chrsiz) { | |
417 | // | |
418 | TG4Globals:: Warning("TGeant4::Gdhead(..) is not implemented."); | |
419 | } | |
420 | ||
421 | void TGeant4::Gdman(Float_t u, Float_t v, const char* type) { | |
422 | // | |
423 | TG4Globals:: Warning("TGeant4::Gdman(..) is not implemented."); | |
424 | } | |
425 | ||
426 | void TGeant4::SetColors() { | |
427 | // | |
428 | TG4Globals:: Warning("TGeant4::SetColours() is not implemented."); | |
429 | } | |
430 | ||
431 | void TGeant4::Gtreve() { | |
432 | // | |
433 | TG4Globals:: Warning("TGeant4::Gtreve() is not implemented."); | |
434 | } | |
435 | ||
aee8290b | 436 | void TGeant4::GtreveRoot() { |
2817d3e2 | 437 | // |
aee8290b | 438 | TG4Globals:: Warning("TGeant4::GtreveRoot() is not implemented."); |
2817d3e2 | 439 | } |
440 | ||
441 | void TGeant4::Gckmat(Int_t itmed, char* natmed) { | |
442 | // | |
443 | TG4Globals:: Warning("TGeant4::Gckmat(..) is not implemented."); | |
444 | } | |
445 | ||
446 | void TGeant4::InitLego() { | |
447 | // | |
448 | TG4Globals:: Warning("TGeant4::InitLego() is not implemented."); | |
449 | } | |
450 | ||
451 | void TGeant4::Gfpart(Int_t ipart, char *name, Int_t& itrtyp, | |
452 | Float_t& amass, Float_t& charge, Float_t& tlife) { | |
453 | // | |
454 | TG4Globals:: Warning("TGeant4::Gfpart(..) is not implemented."); | |
455 | } | |
456 | ||
457 | void TGeant4::Gspart(Int_t ipart, const char *name, Int_t itrtyp, | |
458 | Float_t amass, Float_t charge, Float_t tlife) { | |
459 | // | |
460 | TG4Globals:: Warning("TGeant4::Gspart(..) is not implemented."); | |
461 | } |