Remove the depracated warning about using gconfig.C

karabowi · karabowi · commit c2f7a813fe3d · 2022-06-23T09:37:22.000+02:00
Use the FairVMCConfig in the simulation macros. Addresses issue #1010.
diff --git a/examples/MQ/pixelDetector/macros/run_sim.C b/examples/MQ/pixelDetector/macros/run_sim.C
@@ -61,6 +61,7 @@ void run_sim(Int_t nEvents = 10, TString mcEngine = "TGeant3", Int_t fileId = 0,
     // -----   Create simulation run   ----------------------------------------
     auto run = std::make_unique<FairRunSim>();
     run->SetName(mcEngine);   // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);       // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/MQ/pixelDetector/src/CMakeLists.txt b/examples/MQ/pixelDetector/src/CMakeLists.txt
@@ -57,6 +57,7 @@ list(APPEND headers
 add_library(${target} SHARED ${sources} ${no_dict_sources} ${headers} ${no_dict_headers})
 add_library(FairRoot::${target} ALIAS ${target})
 set_target_properties(${target} PROPERTIES ${PROJECT_LIBRARY_PROPERTIES})
+set_property(TARGET ${target} PROPERTY POSITION_INDEPENDENT_CODE ON)
 
 target_include_directories(${target} PUBLIC
   $<BUILD_INTERFACE:${CMAKE_CURRENT_SOURCE_DIR}>
@@ -71,6 +72,7 @@ target_link_libraries(${target} PUBLIC
   FairRoot::Gen # FairBoxGenerator
   FairRoot::GeoBase
   FairRoot::ParBase # FairParGenericSet
+  FairRoot::ExGConfig
   FairRoot::ExMCStack
   FairRoot::ExPassive # FairCave
   FairRoot::FairMQ
@@ -84,6 +86,7 @@ target_link_libraries(${target} PUBLIC
   ROOT::Physics
   ROOT::Geom
   ROOT::Hist
+
   ${VMCLIB}
 )
 
diff --git a/examples/MQ/pixelDetector/src/devices/FairMQSimDevice.cxx b/examples/MQ/pixelDetector/src/devices/FairMQSimDevice.cxx
@@ -19,6 +19,7 @@
 #include "FairPrimaryGenerator.h"
 #include "FairRunSim.h"
 #include "FairRuntimeDb.h"
+#include "FairVMCConfig.h"
 
 #include <TCollection.h>
 #include <TList.h>
@@ -60,7 +61,7 @@ void FairMQSimDevice::InitTask()
     }
 
     fRunSim->SetName(fTransportName.data());
-    //  fRunSim->SetSimulationConfig(new FairVMCConfig());
+    fRunSim->SetSimulationConfig(new FairVMCConfig());
     fRunSim->SetIsMT(kFALSE);
 
     if (fUserConfig.Length() > 0)
diff --git a/examples/advanced/Tutorial3/macro/run_sim.C b/examples/advanced/Tutorial3/macro/run_sim.C
@@ -29,6 +29,7 @@ void run_sim(Int_t nEvents = 100, TString mcEngine = "TGeant3")
     // ------------------------
 
     fRun->SetName(mcEngine);
+    fRun->SetSimulationConfig(new FairVMCConfig());
 
     TString outFile = "data/testrun_";
     outFile = outFile + mcEngine + ".root";
diff --git a/examples/advanced/Tutorial3/macro/run_sim_sep.C b/examples/advanced/Tutorial3/macro/run_sim_sep.C
@@ -23,6 +23,7 @@ void run_sim_sep(Int_t fileId, Int_t nEvents = 1000, TString mcEngine = "TGeant3
     // ------------------------
 
     fRun->SetName(mcEngine);
+    fRun->SetSimulationConfig(new FairVMCConfig());
 
     TString outfile = Form("data/testrun_%s_f%d.root", mcEngine.Data(), fileId);
     TString outparam = Form("data/testpar_%s_f%d.root", mcEngine.Data(), fileId);
diff --git a/examples/advanced/propagator/macros/runMC.C b/examples/advanced/propagator/macros/runMC.C
@@ -50,7 +50,7 @@ int runMC(Int_t nEvents = 1000, TString mcEngine = "TGeant4", Bool_t isMT = fals
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);   // Transport engine
-    //  run->SetSimulationConfig(new FairVMCConfig());
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/advanced/propagator/macros/runMM.C b/examples/advanced/propagator/macros/runMM.C
@@ -43,7 +43,7 @@ int runMM(Int_t nEvents = 1000, TString mcEngine = "TGeant4", Bool_t isMT = fals
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);   // Transport engine
-    //  run->SetSimulationConfig(new FairVMCConfig());
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1.C b/examples/simulation/Tutorial1/macros/run_tutorial1.C
@@ -65,7 +65,7 @@ void run_tutorial1(Int_t nEvents = 10,
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);   // Transport engine
-    FairGenericVMCConfig* config = new FairGenericVMCConfig();
+    FairVMCConfig* config = new FairVMCConfig();
     if (loadPostInitConfig)
         config->UsePostInitConfig();
     run->SetSimulationConfig(config);
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_binary.C b/examples/simulation/Tutorial1/macros/run_tutorial1_binary.C
@@ -17,6 +17,7 @@
 #include "FairSimConfig.h"
 #include "FairSystemInfo.h"
 #include "FairTutorialDet1.h"
+#include "FairVMCConfig.h"
 
 #include <TRandom3.h>
 #include <TStopwatch.h>
@@ -78,6 +79,7 @@ void run_tutorial1_main(const FairSimConfig& config)
     // -----   Create simulation run   ----------------------------------------
     FairRunSim run;
     run.SetName(config.GetEngine());              // Transport engine
+    run.SetSimulationConfig(new FairVMCConfig());
     run.SetIsMT(config.IsMultiThreaded());        // Multi-threading mode (Geant4 only)
     run.SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run.GetRuntimeDb();
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_fastsim.C b/examples/simulation/Tutorial1/macros/run_tutorial1_fastsim.C
@@ -62,6 +62,7 @@ void run_tutorial1_fastsim(Int_t nEvents = 10, TString mcEngine = "TGeant3", Boo
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_mesh.C b/examples/simulation/Tutorial1/macros/run_tutorial1_mesh.C
@@ -60,6 +60,7 @@ void run_tutorial1_mesh(Int_t nEvents = 10, TString mcEngine = "TGeant3")
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
     // ------------------------------------------------------------------------
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_pythia6.C b/examples/simulation/Tutorial1/macros/run_tutorial1_pythia6.C
@@ -64,6 +64,7 @@ void run_tutorial1_pythia6(Int_t nEvents = 10, TString mcEngine = "TGeant3")
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     run->SetPythiaDecayer(pythia6Config);          // Define Pythia6 as decayer
 
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_pythia8.C b/examples/simulation/Tutorial1/macros/run_tutorial1_pythia8.C
@@ -64,6 +64,7 @@ void run_tutorial1_pythia8(Int_t nEvents = 10, TString mcEngine = "TGeant3")
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     run->SetPythiaDecayer(pythia8Config);          // Define Pythia8 as decayer
 
diff --git a/examples/simulation/Tutorial1/macros/run_tutorial1_urqmd.C b/examples/simulation/Tutorial1/macros/run_tutorial1_urqmd.C
@@ -46,6 +46,7 @@ void run_tutorial1_urqmd(Int_t nEvents = 2, TString mcEngine = "TGeant3")
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
     // ------------------------------------------------------------------------
diff --git a/examples/simulation/Tutorial1/src/CMakeLists.txt b/examples/simulation/Tutorial1/src/CMakeLists.txt
@@ -35,6 +35,7 @@ target_link_libraries(${target} PUBLIC
   FairRoot::FairTools # FairLogger
   FairRoot::FairFastSim
   FairRoot::ParBase
+  FairRoot::ExGConfig
   FairRoot::ExPassive
   FairRoot::ExMCStack
   FairRoot::Gen
diff --git a/examples/simulation/Tutorial2/macros/run_background.C b/examples/simulation/Tutorial2/macros/run_background.C
@@ -51,6 +51,7 @@ void run_background(Int_t nEvents = 130)
     // -----   Create simulation run   ----------------------------------------
     auto run = std::make_unique<FairRunSim>();
     run->SetName("TGeant4");                       // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
     // ------------------------------------------------------------------------
diff --git a/examples/simulation/Tutorial2/macros/run_signal.C b/examples/simulation/Tutorial2/macros/run_signal.C
@@ -53,6 +53,7 @@ void run_signal(Int_t fileNr, Int_t nEvents = 10)
     // -----   Create simulation run   ----------------------------------------
     auto run = std::make_unique<FairRunSim>();
     run->SetName("TGeant4");                       // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
     // ------------------------------------------------------------------------
diff --git a/examples/simulation/Tutorial2/macros/run_tutorial2.C b/examples/simulation/Tutorial2/macros/run_tutorial2.C
@@ -61,6 +61,7 @@ void run_tutorial2(Int_t nEvents = 10, TString mcEngine = "TGeant4", Bool_t isMT
     // -----   Create simulation run   ----------------------------------------
     auto run = std::make_unique<FairRunSim>();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/simulation/Tutorial4/gconfig/g3Config.yaml b/examples/simulation/Tutorial4/gconfig/g3Config.yaml
@@ -0,0 +1,40 @@
+G3_TRIG:          1
+G3_SWIT:          [4, 100]
+G3_DEBU:          [0, 0, 1]
+G3_RAYL:          1
+G3_STRA:          1
+G3_AUTO:          0        #Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
+G3_ABAN:          0        #Restore 3.16 behaviour for abandoned tracks
+G3_OPTI:          2        #Select optimisation level for GEANT geometry searches (0,1,2)
+G3_ERAN:          5.e-7
+G3_CKOV:          1        #Cerenkov photons
+
+MonteCarlo_Process:
+    PAIR:         1        #pair production
+    COMP:         1        #Compton scattering
+    PHOT:         1        #photo electric effect
+    PFIS:         0        #photofission
+    DRAY:         1        #delta-ray
+    ANNI:         1        #annihilation
+    BREM:         1        #bremsstrahlung
+    HADR:         1        #hadronic process
+    MUNU:         1        #muon nuclear interaction
+    DCAY:         1        #decay
+    LOSS:         1        #energy loss
+    MULS:         0        #multiple scattering
+
+MonteCarlo_Cut:
+    CUTGAM:       1.0E-3   #gammas (GeV)
+    CUTELE:       1.0E-3   #electrons (GeV)
+    CUTNEU:       1.0E-3   #neutral hadrons (GeV)
+    CUTHAD:       1.0E-3   #charged hadrons (GeV)
+    CUTMUO:       1.0E-3   #muons (GeV)
+    BCUTE:        1.0E-3   #electron bremsstrahlung (GeV)
+    BCUTM:        1.0E-3   #muon and hadron bremsstrahlung (GeV)
+    DCUTE:        1.0E-3   #delta-rays by electrons (GeV)
+    DCUTM:        1.0E-3   #delta-rays by muons (GeV)
+    PPCUTM:       1.0E-3   #direct pair production by muons (GeV)
+    TOFMAX:       1.E10    #time of flight cut in seconds
+
+FairStack_StoreSecondaries: true
+FairStack_MinPoints:        0
diff --git a/examples/simulation/Tutorial4/gconfig/g4Config.yaml b/examples/simulation/Tutorial4/gconfig/g4Config.yaml
@@ -0,0 +1,40 @@
+Geant4_UserGeometry:     "geomRoot"
+Geant4_PhysicsList:      "QGSP_BERT_EMV"
+Geant4_SpecialProcess:   "stepLimiter+specialCuts+specialControls"
+Geant4_SpecialStacking:  false
+Geant4_Multithreaded:    false
+
+Geant4_MaxNStep:         10000    #default is 30000
+Geant4_Commands:
+ - /cuts/setHighEdge 100 GeV
+ - /tracking/verbose 0
+
+MonteCarlo_Process:
+    PAIR:         1        #pair production
+    COMP:         1        #Compton scattering
+    PHOT:         1        #photo electric effect
+    PFIS:         0        #photofission
+    DRAY:         1        #delta-ray
+    ANNI:         1        #annihilation
+    BREM:         1        #bremsstrahlung
+    HADR:         1        #hadronic process
+    MUNU:         1        #muon nuclear interaction
+    DCAY:         1        #decay
+    LOSS:         1        #energy loss
+    MULS:         1        #multiple scattering
+
+MonteCarlo_Cut:
+    CUTGAM:       1.0E-3   #gammas (GeV)
+    CUTELE:       1.0E-3   #electrons (GeV)
+    CUTNEU:       1.0E-3   #neutral hadrons (GeV)
+    CUTHAD:       1.0E-3   #charged hadrons (GeV)
+    CUTMUO:       1.0E-3   #muons (GeV)
+    BCUTE:        1.0E-3   #electron bremsstrahlung (GeV)
+    BCUTM:        1.0E-3   #muon and hadron bremsstrahlung (GeV)
+    DCUTE:        1.0E-3   #delta-rays by electrons (GeV)
+    DCUTM:        1.0E-3   #delta-rays by muons (GeV)
+    PPCUTM:       1.0E-3   #direct pair production by muons (GeV)
+    TOFMAX:       1.E10    #time of flight cut in seconds
+
+FairStack_StoreSecondaries: true
+FairStack_MinPoints:        0
diff --git a/examples/simulation/Tutorial4/macros/run_tutorial4.C b/examples/simulation/Tutorial4/macros/run_tutorial4.C
@@ -70,6 +70,7 @@ void run_tutorial4(Int_t nEvents = 10, TString mcEngine = "TGeant3", Bool_t doAl
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/examples/simulation/rutherford/macros/run_rad.C b/examples/simulation/rutherford/macros/run_rad.C
@@ -51,6 +51,7 @@ void run_rad(Int_t nEvents = 100, TString mcEngine = "TGeant4")
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    run->SetSimulationConfig(new FairVMCConfig());
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
     // ------------------------------------------------------------------------
diff --git a/examples/simulation/rutherford/macros/run_rutherford.C b/examples/simulation/rutherford/macros/run_rutherford.C
@@ -53,6 +53,7 @@ void run_rutherford(Int_t nEvents = 10, TString mcEngine = "TGeant4", Bool_t isM
     // -----   Create simulation run   ----------------------------------------
     FairRunSim* run = new FairRunSim();
     run->SetName(mcEngine);                        // Transport engine
+    //    run->SetSimulationConfig(new FairVMCConfig());
     run->SetIsMT(isMT);                            // Multi-threading mode (Geant4 only)
     run->SetSink(new FairRootFileSink(outFile));   // Output file
     FairRuntimeDb* rtdb = run->GetRuntimeDb();
diff --git a/fairtools/MCConfigurator/CMakeLists.txt b/fairtools/MCConfigurator/CMakeLists.txt
@@ -29,6 +29,7 @@ target_link_libraries(${target} PUBLIC
 
   geant321
   geant4vmc
+  Pythia6
   yaml-cpp
 
   ROOT::Core
diff --git a/fairtools/MCConfigurator/FairYamlVMCConfig.cxx b/fairtools/MCConfigurator/FairYamlVMCConfig.cxx
@@ -69,6 +69,7 @@ void FairYamlVMCConfig::SetupPostInit(const char* mcEngine)
     if (!((strcmp(mcEngine, "TGeant4") == 0))) {
         LOG(fatal) << "FairYamlVMCConfig::SetupPostInit() only valid for TGeant4.";
     }
+    LOG(info) << "Loading Geant4 PostInit Config." << endl;
 
     TString work = getenv("VMCWORKDIR");
     TString work_config = work + "/gconfig/";
