61int main(
int argc,
char** argv)
72 if ( argc == 2 ) macro = argv[1];
75 for ( G4int i=1; i<argc; i=i+2 )
77 if (
G4String(argv[i]) ==
"-m" ) macro = argv[i+1];
78 else if (
G4String(argv[i]) ==
"-b" ) onOffBiasing = argv[i+1];
87 if ( onOffBiasing ==
"" ) onOffBiasing =
"on";
96 auto* runManager = G4RunManagerFactory::CreateRunManager();
97 runManager->SetNumberOfThreads(4);
101 runManager->SetUserInitialization(detector);
106 if ( onOffBiasing ==
"on" )
110 std::vector< G4String > piPlusProcessesToBias, piMinusProcessesToBias,
111 protonProcessesToBias, antiProtonProcessesToBias, neutronProcessesToBias,
112 antiNeutronProcessesToBias;
113 piPlusProcessesToBias .push_back(
"pi+Inelastic");
114 piMinusProcessesToBias .push_back(
"pi-Inelastic");
115 protonProcessesToBias .push_back(
"protonInelastic");
116 antiProtonProcessesToBias .push_back(
"anti_protonInelastic");
117 neutronProcessesToBias .push_back(
"neutronInelastic");
118 neutronProcessesToBias .push_back(
"nCapture" );
119 antiNeutronProcessesToBias.push_back(
"anti_neutronInelastic");
120 biasingPhysics->PhysicsBias(
"pi+", piPlusProcessesToBias);
121 biasingPhysics->PhysicsBias(
"pi-", piMinusProcessesToBias);
122 biasingPhysics->PhysicsBias(
"proton", protonProcessesToBias);
123 biasingPhysics->PhysicsBias(
"anti_proton", antiProtonProcessesToBias);
124 biasingPhysics->PhysicsBias(
"neutron", neutronProcessesToBias);
125 biasingPhysics->PhysicsBias(
"anti_neutron", antiNeutronProcessesToBias);
128 std::vector< G4String > gammaProcessesToBias, electronProcessesToBias,
129 positronProcessesToBias;
130 gammaProcessesToBias .push_back(
"conv" );
131 gammaProcessesToBias .push_back(
"photonNuclear" );
132 electronProcessesToBias .push_back(
"electronNuclear" );
133 positronProcessesToBias .push_back(
"annihil" );
134 positronProcessesToBias .push_back(
"positronNuclear" );
135 biasingPhysics->PhysicsBias(
"gamma", gammaProcessesToBias);
136 biasingPhysics->PhysicsBias(
"e-", electronProcessesToBias);
137 biasingPhysics->PhysicsBias(
"e+", positronProcessesToBias);
139 std::vector< G4String > pi0ProcessesToBias;
140 pi0ProcessesToBias .push_back(
"Decay" );
141 biasingPhysics->PhysicsBias(
"pi0", pi0ProcessesToBias);
143 physicsList->RegisterPhysics(biasingPhysics);
144 G4cout <<
" ********************************************************* " << G4endl;
145 G4cout <<
" ********** processes are wrapped for biasing ************ " << G4endl;
146 G4cout <<
" ********************************************************* " << G4endl;
150 G4cout <<
" ************************************************* " << G4endl;
151 G4cout <<
" ********** processes are not wrapped ************ " << G4endl;
152 G4cout <<
" ************************************************* " << G4endl;
154 runManager->SetUserInitialization(physicsList);
159 runManager->Initialize();
164 visManager->Initialize();
167 G4UImanager* UImanager = G4UImanager::GetUIpointer();
171 G4String command =
"/control/execute ";
172 UImanager->ApplyCommand(command+macro);
176 UImanager->ApplyCommand(
"/control/execute vis.mac");