36 : verboseLevel(0), eventNumber(0.0), averageMultiplicity(0.0),
37 averageProtonNumber(0.0), averageNeutronNumber(0.0),
38 averagePionNumber(0.0), averageNucleonKinEnergy(0.0),
39 averageProtonKinEnergy(0.0), averageNeutronKinEnergy(0.0),
40 averagePionKinEnergy(0.0), averageExitationEnergy(0.0),
41 averageOutgoingNuclei(0.0), fissy_prob(0.0), averagePionPl(0.0),
42 averagePionMin(0.0), averagePion0(0.0), averageA(0.0), averageZ(0.0),
43 inel_csec(0.0), withNuclei(
false) {
105 if (nucleus.size() > 0) {
126 const std::vector<G4InuclElementaryParticle>& particles =
130 for (
G4int i = 0; i <
G4int(particles.size()); i++) {
137 if (particles[i].type() == 1) {
150 }
else if (particles[i].
pion()) {
155 if (particles[i].type() == 3) {
159 }
else if (particles[i].type() == 5) {
163 }
else if (particles[i].type() == 7) {
174 const std::vector<G4InuclElementaryParticle>& particles =
178 for (
G4int i = 0; i <
G4int(particles.size()); i++) {
183 if (particles[i].type() == 1) {
192 }
else if (particles[i].
pion()) {
203 G4cout <<
" >>> G4Analyser::printResultsSimple" <<
G4endl;
221 <<
" average Excitation Energy " <<
256 <<
" average Excitation Energy " <<
268 G4cout <<
" >>> G4Analyser::handleWatcherStatistics" <<
G4endl;
292 std::pair<G4double, G4double> rat_err =
ana_watchers[iw].getAverageRatio();
293 averat += rat_err.first;
294 ave_err += rat_err.second;
296 std::pair<G4double, G4double> cs_err =
ana_watchers[iw].getExpCs();
297 tot_exper += cs_err.first;
298 tot_exper_err += cs_err.second;
299 std::pair<G4double, G4double> inucl_cs_err =
ana_watchers[iw].getInuclCs();
300 tot_inucl += inucl_cs_err.first;
301 tot_inucl_err += inucl_cs_err.second;
304 if (iz_checked > 0.0) {
306 checked += iz_checked;
312 gl_chsq = std::sqrt(gl_chsq) / checked;
315 fgr = std::pow(10.0, std::sqrt(fgr / checked));
319 G4cout <<
" total exper c.s. " << tot_exper <<
" err " << tot_exper_err <<
320 " tot inucl c.s. " << tot_inucl <<
" err " << tot_inucl_err <<
G4endl;
321 G4cout <<
" checked total " << checked <<
" lhood " << fgr << G4endl
322 <<
" average ratio " << averat <<
" err " << ave_err << G4endl
323 <<
" global chsq " << gl_chsq <<
G4endl;
330 G4cout <<
" >>> G4Analyser::printResultsNtuple" <<
G4endl;