CbmGeant4Settings.cxx

Go to the documentation of this file.

#include "CbmGeant4Settings.h"
 
#include <FairLogger.h>
#include <FairRunSim.h>
 
#include <TGeant4.h>
#include <TPythia6Decayer.h>
#include <TVirtualMC.h>
 
#include <G4NuclearLevelData.hh>
 
void CbmGeant4Settings::Init(TVirtualMC* vmc) {
 
  CbmVMCSettings::Init(vmc);
 
  TGeant4* vmcg4 = dynamic_cast<TGeant4*>(vmc);
 
  assert(vmcg4);
 
  // TODO: These settings were taken over from g4Config.C. To be documented.
 
  // --- Set external decayer (Pythia) if required
  if (FairRunSim::Instance()->IsExtDecayer()) {
    TVirtualMCDecayer* decayer = TPythia6Decayer::Instance();
    vmcg4->SetExternalDecayer(decayer);
    LOG(info) << ": Using Phythia6 decayer";
  }
 
  // --- Random seed and maximum number of steps
  Text_t buffer[50];
  sprintf(
    buffer, "/random/setSeeds %i  %i ", gRandom->GetSeed(), gRandom->GetSeed());
  vmcg4->ProcessGeantCommand(buffer);
 
  vmcg4->SetMaxNStep(fMaxNumSteps);
 
  for (std::string command : fG4Commands) {
    LOG(info) << "Passing the following command to Geant4: " << command;
    vmcg4->ProcessGeantCommand(command.c_str());
  }
 
  LOG(warning) << "Fix for Geant4 bug applied (see "
                  "https://bugzilla-geant4.kek.jp/show_bug.cgi?id=2226)";
  G4NuclearLevelData::GetInstance()->GetParameters()->SetMaxLifeTime(
    1 * CLHEP::microsecond);
}
 
ClassImp(CbmGeant4Settings);