g4Config.C

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file g4Config.C

// Configuration macro for Geant4 VirtualMC
void Config()
{

   TG4RunConfiguration* runConfiguration 
     //= new TG4RunConfiguration("geomRoot", "QGSP_BERT_EMV", "stepLimiter+specialCuts+specialControls");
     //= new TG4RunConfiguration("geomRootToGeant4", "emStandard", "stepLimiter+specialCuts+specialControls");
     //= new TG4RunConfiguration("geomRoot", "emStandard", "stepLimiter");
     = new TG4RunConfiguration("geomRoot", "FTFP_BERT", "stepLimiter");
     //= new TG4RunConfiguration("geomRoot", "emStandard", "stepLimiter+specialControls");
     // In case want to activate SetCuts.C ...;
     //= new TG4RunConfiguration("geomRoot", "emStandard", "stepLimiter+specialCuts");

   TGeant4* geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
   cout << "Geant4 has been created." << endl;

   PndStack *stack = new PndStack(1000); 
   stack->StoreSecondaries(kTRUE);
   stack->SetMinPoints(1);
   geant4->SetStack(stack);

#if _RETURN_BACK_
   if(FairRunSim::Instance()->IsExtDecayer()){
      TVirtualMCDecayer* decayer = TPythia6Decayer::Instance();
      geant4->SetExternalDecayer(decayer);
   }
#endif
  

   TString configm(gSystem->Getenv("VMCWORKDIR"));
   TString configm1 = configm + "/gconfig/g4config.in";
   cout << " -I g4Config() using g4conf  macro: " << configm1 << endl;

     //set geant4 specific stuff
   geant4->SetMaxNStep(10000);  // default is 30000

  Text_t buffer[50];
  sprintf(buffer,"/random/setSeeds %i  %i ",gRandom->GetSeed(), gRandom->GetSeed());
  geant4->ProcessGeantCommand(buffer);
 
  geant4->ProcessGeantMacro(configm1.Data());


}