// Configuration macro for Geant4 VirtualMC
void Config()
{
/// Create the run configuration
/// In constructor user has to specify the geometry input
/// and select geometry navigation via the following options:
/// - geomVMCtoGeant4 - geometry defined via VMC, G4 native navigation
/// - geomVMCtoRoot - geometry defined via VMC, Root navigation
/// - geomRoot - geometry defined via Root, Root navigation
/// - geomRootToGeant4 - geometry defined via Root, G4 native navigation
/// - geomGeant4 - geometry defined via Geant4, G4 native navigation
///
/// The second argument in the constructor selects physics list:
/// Available options:
/// EMonly, EMonly+Extra, Hadron_EM, Hadron_EM+Extra
/// where EMonly = emStandard
/// Hadron = FTFP_BERT FTFP_BERT_TRV FTFP_BERT_HP FTFP_INCLXX FTFP_INCLXX_HP FTF_BIC LBE QBBC QGSP_BERT QGSP_BERT_HP QGSP_BIC QGSP_BIC_HP QGSP_FTFP_BERT QGSP_INCLXX QGSP_INCLXX_HP QGS_BIC Shielding ShieldingLEND
/// EM = _EMV _EMX _EMY _EMZ _LIV _PEN
/// Extra = extra optical radDecay
/// The Extra selections are cumulative, while Hadron selections are exlusive.
/// The third argument activates the special processes in the TG4SpecialPhysicsList,
/// which implement VMC features:
/// - stepLimiter - step limiter (default)
/// - specialCuts - VMC cuts
/// - specialControls - VMC controls for activation/inactivation selected processes
/// - stackPopper - stackPopper process
/// When more than one options are selected, they should be separated with '+'
/// character: eg. stepLimit+specialCuts.
TG4RunConfiguration* runConfiguration
= new TG4RunConfiguration("geomRoot", "QGSP_FTFP_BERT", "stepLimiter+specialCuts+specialControls+stackPopper");
/// Create the G4 VMC
TGeant4* geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
cout << "Geant4 has been created." << endl;
/// create the Specific stack
FccStack *stack = new FccStack(1000);
stack->StoreSecondaries(kTRUE);
stack->SetMinPoints(0);
geant4->SetStack(stack);
if(FairRunSim::Instance()->IsExtDecayer()){
TVirtualMCDecayer* decayer = TPythia6Decayer::Instance();
geant4->SetExternalDecayer(decayer);
}
/// Customise Geant4 setting
/// (verbose level, global range cut, ..)
TString configm(gSystem->Getenv("VMCWORKDIR"));
configm1 = configm + "/gconfig/g4config.in";
cout << " -I g4Config() using g4conf macro: " << configm1 << endl;
//set geant4 specific stuff
geant4->SetMaxNStep(10000); // default is 30000
geant4->ProcessGeantMacro(configm1.Data());
}
// Configuration macro for Geant4 VirtualMC
void Config()
{
/// Create the run configuration
/// In constructor user has to specify the geometry input
/// and select geometry navigation via the following options:
/// - geomVMCtoGeant4 - geometry defined via VMC, G4 native navigation
/// - geomVMCtoRoot - geometry defined via VMC, Root navigation
/// - geomRoot - geometry defined via Root, Root navigation
/// - geomRootToGeant4 - geometry defined via Root, G4 native navigation
/// - geomGeant4 - geometry defined via Geant4, G4 native navigation
///
/// The second argument in the constructor selects physics list:
/// - emStandard - standard em physics (default)
/// - emStandard+optical - standard em physics + optical physics
/// - XYZ - selected hadron physics list ( XYZ = LHEP, QGSP, ...)
/// - XYZ+optical - selected hadron physics list + optical physics
///
/// The third argument activates the special processes in the TG4SpecialPhysicsList,
/// which implement VMC features:
/// - stepLimiter - step limiter (default)
/// - specialCuts - VMC cuts
/// - specialControls - VMC controls for activation/inactivation selected processes
/// - stackPopper - stackPopper process
/// When more than one options are selected, they should be separated with '+'
/// character: eg. stepLimit+specialCuts.
TG4RunConfiguration* runConfiguration
= new TG4RunConfiguration("geomRoot", "QGSP_BERT_EMV", "stepLimiter+specialCuts+specialControls");
/// Create the G4 VMC
TGeant4* geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
cout << "Geant4 has been created." << endl;
/// create the Specific stack
FairStack *stack = new FairStack(1000);
stack->StoreSecondaries(kTRUE);
stack->SetMinPoints(0);
geant4->SetStack(stack);
if(FairRunSim::Instance()->IsExtDecayer()){
TVirtualMCDecayer* decayer = TVirtualMC::GetMC()->GetDecayer();
// TVirtualMCDecayer* decayer = TPythia6Decayer::Instance();
geant4->SetExternalDecayer(decayer);
}
/// Customise Geant4 setting
/// (verbose level, global range cut, ..)
TString configm(gSystem->Getenv("VMCWORKDIR"));
configm1 = configm + "/common/gconfig/g4config.in";
cout << " -I g4Config() using g4conf macro: " << configm1 << endl;
//set geant4 specific stuff
geant4->SetMaxNStep(10000); // default is 30000
geant4->ProcessGeantMacro(configm1.Data());
}
void Config()
{
new TFluka("C++ Interface to Fluka", 1/*verbositylevel*/);
cout << "GConfig: Fluka has been created." << endl;
FairStack *st = new FairStack();
st->SetMinPoints(0);
TVirtual::GetMC()->SetStack( st ) ;
TVirtual::GetMC()->SetProcess("CKOV",1);
// set the common cuts
TString configm(gSystem->Getenv("VMCWORKDIR"));
TString cuts = configm + "/gconfig/SetCuts.C";
cout << "Physics cuts with script \n "<< cuts.Data() << endl;
Int_t cut=gROOT->LoadMacro(cuts.Data());
if(cut==0)gInterpreter->ProcessLine("SetCuts()");
}
