33 #ifdef G4MULTITHREADED
56 G4cerr <<
" ./exampleGB07 [-m macro ] "
57 <<
" [-b biasing {'on','off'}]"
58 <<
"\n or\n ./exampleGB07 [macro.mac]"
65 int main(
int argc,
char** argv)
76 if ( argc == 2 ) macro = argv[1];
79 for (
G4int i=1; i<argc; i=i+2 )
81 if (
G4String(argv[i]) ==
"-m" ) macro = argv[i+1];
82 else if (
G4String(argv[i]) ==
"-b" ) onOffBiasing = argv[i+1];
91 if ( onOffBiasing ==
"" ) onOffBiasing =
"on";
100 #ifdef G4MULTITHREADED
102 G4cout <<
" ********** Run Manager constructed in MT mode ************ " <<
G4endl;
107 G4cout <<
" ********** Run Manager constructed in sequential mode ************ " <<
G4endl;
117 if ( onOffBiasing ==
"on" )
121 std::vector< G4String >
122 piPlusProcessesToBias, piMinusProcessesToBias,
123 protonProcessesToBias, antiProtonProcessesToBias,
124 neutronProcessesToBias, antiNeutronProcessesToBias;
125 piPlusProcessesToBias .push_back(
"pi+Inelastic");
126 piMinusProcessesToBias .push_back(
"pi-Inelastic");
127 protonProcessesToBias .push_back(
"protonInelastic");
128 antiProtonProcessesToBias .push_back(
"anti_protonInelastic");
129 neutronProcessesToBias .push_back(
"neutronInelastic");
130 neutronProcessesToBias .push_back(
"nCapture" );
131 antiNeutronProcessesToBias.push_back(
"anti_neutronInelastic");
132 biasingPhysics->
PhysicsBias(
"pi+", piPlusProcessesToBias);
133 biasingPhysics->
PhysicsBias(
"pi-", piMinusProcessesToBias);
134 biasingPhysics->
PhysicsBias(
"proton", protonProcessesToBias);
135 biasingPhysics->
PhysicsBias(
"anti_proton", antiProtonProcessesToBias);
136 biasingPhysics->
PhysicsBias(
"neutron", neutronProcessesToBias);
137 biasingPhysics->
PhysicsBias(
"anti_neutron", antiNeutronProcessesToBias);
140 std::vector< G4String > gammaProcessesToBias,
141 electronProcessesToBias, positronProcessesToBias;
142 gammaProcessesToBias .push_back(
"conv" );
143 gammaProcessesToBias .push_back(
"photonNuclear" );
144 electronProcessesToBias .push_back(
"electronNuclear" );
145 positronProcessesToBias .push_back(
"annihil" );
146 positronProcessesToBias .push_back(
"positronNuclear" );
147 biasingPhysics->
PhysicsBias(
"gamma", gammaProcessesToBias);
148 biasingPhysics->
PhysicsBias(
"e-", electronProcessesToBias);
149 biasingPhysics->
PhysicsBias(
"e+", positronProcessesToBias);
151 std::vector< G4String > pi0ProcessesToBias;
152 pi0ProcessesToBias .push_back(
"Decay" );
153 biasingPhysics->
PhysicsBias(
"pi0", pi0ProcessesToBias);
156 G4cout <<
" ********************************************************* " <<
G4endl;
157 G4cout <<
" ********** processes are wrapped for biasing ************ " <<
G4endl;
158 G4cout <<
" ********************************************************* " <<
G4endl;
162 G4cout <<
" ************************************************* " <<
G4endl;
163 G4cout <<
" ********** processes are not wrapped ************ " <<
G4endl;
164 G4cout <<
" ************************************************* " <<
G4endl;
183 G4String command =
"/control/execute ";