62{
63
64
65 if ( argc > 5 ) {
67 return 1;
68 }
69
72 if ( argc == 2 ) macro = argv[1];
73 else
74 {
75 for ( G4int i=1; i<argc; i=i+2 )
76 {
77 if (
G4String(argv[i]) ==
"-m" ) macro = argv[i+1];
78 else if (
G4String(argv[i]) ==
"-b" ) onOffBiasing = argv[i+1];
79 else
80 {
82 return 1;
83 }
84 }
85 }
86
87 if ( onOffBiasing == "" ) onOffBiasing = "on";
88
89
91 if ( macro == "" ) {
93 }
94
95
96 auto* runManager = G4RunManagerFactory::CreateRunManager();
97 runManager->SetNumberOfThreads(4);
98
99
101 runManager->SetUserInitialization(detector);
102
104
106 if ( onOffBiasing == "on" )
107 {
108
109
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);
126
127
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);
138
139 std::vector< G4String > pi0ProcessesToBias;
140 pi0ProcessesToBias .push_back( "Decay" );
141 biasingPhysics->PhysicsBias("pi0", pi0ProcessesToBias);
142
143 physicsList->RegisterPhysics(biasingPhysics);
144 G4cout << " ********************************************************* " << G4endl;
145 G4cout << " ********** processes are wrapped for biasing ************ " << G4endl;
146 G4cout << " ********************************************************* " << G4endl;
147 }
148 else
149 {
150 G4cout << " ************************************************* " << G4endl;
151 G4cout << " ********** processes are not wrapped ************ " << G4endl;
152 G4cout << " ************************************************* " << G4endl;
153 }
154 runManager->SetUserInitialization(physicsList);
155
157
158
159 runManager->Initialize();
160
161
163
164 visManager->Initialize();
165
166
167 G4UImanager* UImanager = G4UImanager::GetUIpointer();
168
169 if ( macro != "" )
170 {
171 G4String command =
"/control/execute ";
172 UImanager->ApplyCommand(command+macro);
173 }
174 else
175 {
176 UImanager->ApplyCommand("/control/execute vis.mac");
177
178
179 ui->SessionStart();
180 delete ui;
181 }
182
183 delete visManager;
184 delete runManager;
185}
Action initialization class.