]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TGeant4/TGeant4.cxx
updated to renaming in global category
[u/mrichter/AliRoot.git] / TGeant4 / TGeant4.cxx
CommitLineData
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 17TGeant4::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
53TGeant4::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
89TGeant4::TGeant4() {
90//
91}
92
93TGeant4::TGeant4(const TGeant4& right) {
94//
95 TG4Globals::Exception("TGeant4 is protected from copying.");
96}
97
98TGeant4::~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
110TGeant4& 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
124void TGeant4::FinishGeometry() {
125//
126 fGeometryManager->Ggclos();
127}
128
129void 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
137void 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
145void 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
152void 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
162void 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
170void TGeant4::Gstpar(Int_t itmed, const char *param, Float_t parval) {
171//
172 fGeometryManager->Gstpar(itmed, param, parval);
173}
174
2817d3e2 175Int_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
181void 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
187void 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
193void 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
199void 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
205void TGeant4::Gsord(const char *name, Int_t iax) {
206//
207 fGeometryManager->Gsord(name, iax);
208}
209
210void 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
217void 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 224void 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 230void 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
236Int_t TGeant4::VolId(const Text_t* volName) const {
237//
238 return fGeometryManager->VolId(volName);
239}
240
241const char* TGeant4::VolName(Int_t id) const {
242//
243 return fGeometryManager->VolName(id);
244}
245
246Int_t TGeant4::NofVolumes() const {
247//
248 return fGeometryManager->NofVolumes();
249}
250
fdd271a5 251Int_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
261void TGeant4::BuildPhysics() {
262//
263 fPhysicsManager->BuildPhysics();
264}
265
266void TGeant4::SetCut(const char* cutName, Float_t cutValue) {
267//
268 fPhysicsManager->SetCut(cutName, cutValue);
269}
270
271void TGeant4::SetProcess(const char* flagName, Int_t flagValue) {
272//
273 fPhysicsManager->SetProcess(flagName, flagValue);
274}
275
276Float_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 281void TGeant4::SetExternalDecayer(AliDecayer* decayer) {
282//
283 TG4Globals:: Warning("TGeant4::SetExternalDecayer() is not implemented.");
284}
285
286AliDecayer* TGeant4::Decayer() const {
287//
288 TG4Globals:: Warning("TGeant4::Decayer() is not implemented.");
289 return 0;
290}
291
292
2817d3e2 293Int_t TGeant4::IdFromPDG(Int_t pdgID) const {
294//
295 return fPhysicsManager->IdFromPDG(pdgID);
296}
297
298Int_t TGeant4::PDGFromId(Int_t mcID) const {
299//
300 return fPhysicsManager->PDGFromId(mcID);
301}
302
303void 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
316void TGeant4::DrawOneSpec(const char* name) {
317//
318 fVisManager->DrawOneSpec(name);
319}
320
321void TGeant4::Gsatt(const char* name, const char* att, Int_t val) {
322//
323 fVisManager->Gsatt(name, att, val);
324}
325
326void 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
334void TGeant4::DrawOneSpec(const char* name) {
335//
336 TG4Globals:: Warning("TGeant4::DrawOneSpec(): no visualization available.");
337}
338
339void TGeant4::Gsatt(const char* name, const char* att, Int_t val) {
340//
341 TG4Globals:: Warning("TGeant4::Gsatt(): no visualization available.");
342}
343
344void 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
356void TGeant4::Init() {
357//
358 fRunManager->Initialize();
359}
360
361void TGeant4::ProcessEvent() {
362//
363 fRunManager->ProcessEvent();
364}
365
366void TGeant4::ProcessRun(Int_t nofEvents) {
367//
368 fRunManager->ProcessRun(nofEvents);
369}
370
371void TGeant4::StartGeantUI() {
372//
373 fRunManager->StartGeantUI();
374}
375
376void TGeant4::StartRootUI() {
377//
562476b1 378 fRunManager->StartRootUI();
2817d3e2 379}
380
381void TGeant4::ProcessGeantMacro(const char* macroName) {
382//
383 fRunManager->ProcessGeantMacro(macroName);
384}
385
386void TGeant4::ProcessGeantCommand(const char* command) {
387//
388 fRunManager->ProcessGeantCommand(command);
389}
390
391Int_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
400void TGeant4::Gdopt(const char* name, const char* value) {
401//
402 TG4Globals:: Warning("TGeant4::Gdopt(..) is not implemented.");
403}
404
405void 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
411void TGeant4::DefaultRange() {
412//
413 TG4Globals:: Warning("TGeant4::DefaultRange() is not implemented.");
414}
415
416void TGeant4::Gdhead(Int_t isel, const char* name, Float_t chrsiz) {
417//
418 TG4Globals:: Warning("TGeant4::Gdhead(..) is not implemented.");
419}
420
421void TGeant4::Gdman(Float_t u, Float_t v, const char* type) {
422//
423 TG4Globals:: Warning("TGeant4::Gdman(..) is not implemented.");
424}
425
426void TGeant4::SetColors() {
427//
428 TG4Globals:: Warning("TGeant4::SetColours() is not implemented.");
429}
430
431void TGeant4::Gtreve() {
432//
433 TG4Globals:: Warning("TGeant4::Gtreve() is not implemented.");
434}
435
aee8290b 436void TGeant4::GtreveRoot() {
2817d3e2 437//
aee8290b 438 TG4Globals:: Warning("TGeant4::GtreveRoot() is not implemented.");
2817d3e2 439}
440
441void TGeant4::Gckmat(Int_t itmed, char* natmed) {
442//
443 TG4Globals:: Warning("TGeant4::Gckmat(..) is not implemented.");
444}
445
446void TGeant4::InitLego() {
447//
448 TG4Globals:: Warning("TGeant4::InitLego() is not implemented.");
449}
450
451void 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
457void 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}