TGeant4/TGeant4.cxx

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