18{
19
20
21
22
23
26
27
28 TGeant4* geant4
29 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
30
31 cout << "Geant4 has been created." << endl;
32
33
34
35 geant4->ProcessGeantMacro("g4config.in");
36
37
38 geant4->ProcessGeantCommand("/tracking/verbose 1");
39
40
41 geant4->ProcessGeantCommand("/mcPhysics/emModel/setEmModel PAI");
42 geant4->ProcessGeantCommand("/mcPhysics/emModel/setRegions liquidArgon");
43 geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles all");
44
45
46 geant4->ProcessGeantCommand("/mcPhysics/emModel/setEmModel SpecialUrbanMsc");
47 geant4->ProcessGeantCommand("/mcPhysics/emModel/setRegions liquidArgon");
48 geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles e- e+");
49
50
51 geant4->ProcessGeantCommand("/mcVerbose/biasingConfigurationManager 3");
52 geant4->ProcessGeantCommand("/mcPhysics/biasing/setModel inclxx");
53 geant4->ProcessGeantCommand("/mcPhysics/biasing/setRegions Lead");
54 geant4->ProcessGeantCommand("/mcPhysics/biasing/setParticles proton neutron pi+ pi-");
55}