109 G4cout <<
"\n ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "
110 <<
"\n **** izotop Z **** " << nuclz <<
G4endl;
119 std::vector<G4bool> not_used(simulated_cs.size(),
true);
120 G4int nmatched = exper_as.size();
121 G4int nused = simulated_cs.size();
123 G4int experAsSize = exper_as.size();
128 exp_cs += exper_cs[
iz];
129 exp_cs_err += exper_err[
iz];
132 G4int simulatedAsSize = simulated_as.size();
133 for (
G4int i = 0; i<simulatedAsSize && !found; i++) {
134 if (std::fabs(simulated_as[i] - a) < small) {
135 G4double rat = simulated_cs[i] / exper_cs[
iz];
137 lhood += std::log10(rat) * std::log10(rat);
140 = std::sqrt(simulated_errors[i]*simulated_errors[i] +
141 exper_err[
iz]*exper_err[
iz] * rat*rat) / exper_cs[
iz];
142 average_ratio += rat;
143 aver_rat_err += rat_err;
145 G4cout <<
" A " << a <<
" exp.cs " << exper_cs[
iz] <<
" err "
146 << exper_err[
iz] << G4endl <<
" sim. cs " << simulated_cs[i]
147 <<
" err " << simulated_errors[i] << G4endl
148 <<
" ratio " << rat <<
" err " << rat_err << G4endl
149 <<
" simulated production rate " << simulated_prob[i] <<
G4endl;
152 izotop_chsq += (rat - 1.0) * (rat - 1.0) / rat_err / rat_err;
158 if (found) nmatched--;
160 G4cout <<
" not found exper.: A " << a <<
" exp.cs " << exper_cs[
iz]
161 <<
" err " << exper_err[
iz] <<
G4endl;
164 G4cout <<
" not found in simulations " << nmatched << G4endl
165 <<
" not found in exper: " << nused <<
G4endl;
167 G4int simulatedAsSize = simulated_as.size();
168 for(
G4int i = 0; i < simulatedAsSize; i++) {
169 inucl_cs += simulated_cs[i];
170 inucl_cs_err += simulated_errors[i];
173 G4cout <<
" extra simul.: A " << simulated_as[i] <<
" sim. cs "
174 << simulated_cs[i] <<
" err " << simulated_errors[i] <<
G4endl;
176 G4cout <<
" simulated production rate " << simulated_prob[i] <<
G4endl;
179 G4int matched = exper_as.size() - nmatched;
183 aver_matched = matched;
184 lhood = std::pow(10.0, std::sqrt(lhood/matched));
186 G4cout <<
" matched " << matched <<
" CHSQ " << std::sqrt(izotop_chsq) / matched
188 <<
" raw chsq " << izotop_chsq << G4endl
189 <<
" average ratio " << average_ratio / matched
190 <<
" err " << aver_rat_err / matched << G4endl
191 <<
" lhood " << lhood <<
G4endl;
198 G4cout <<
" exper. cs " << exp_cs <<
" err " << exp_cs_err << G4endl
199 <<
" inucl. cs " << inucl_cs <<
" err " << inucl_cs_err << G4endl
200 <<
" ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "
G4GLOB_DLL std::ostream G4cout