TGeant4/TGeant4.cxx

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