40 const std::vector<G4double>& expa,
41 const std::vector<G4double>& expcs,
42 const std::vector<G4double>& experr,
45 : nuclz(z), izotop_chsq(0.), average_ratio(0.), aver_rat_err(0.),
46 aver_lhood(0.), aver_matched(0.), exper_as(expa), exper_cs(expcs),
47 exper_err(experr), checkable(check), nucleable(nucl) {}
56 for (
G4int i = 0; i<simulatedAsSize && !here; i++) {
71 for(
G4int i = 0; i < simulatedAsSize ; i++) {
85 for(
G4int iz = 0; iz < experAsSize; iz++) {
90 return std::pair<G4double, G4double>(cs,
err);
97 for(
G4int iz = 0; iz < simulatedAsSize; iz++) {
102 return std::pair<G4double, G4double>(cs,
err);
108 G4cout <<
"\n ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "
118 std::vector<G4bool> not_used(
simulated_cs.size(),
true);
124 for (
G4int iz = 0; iz < experAsSize; iz++) {
132 for (
G4int i = 0; i<simulatedAsSize && !found; i++) {
136 lhood += std::log10(rat) * std::log10(rat);
147 <<
" ratio " << rat <<
" err " << rat_err << G4endl
151 izotop_chsq += (rat - 1.0) * (rat - 1.0) / rat_err / rat_err;
157 if (found) nmatched--;
159 G4cout <<
" not found exper.: A " << a <<
" exp.cs " <<
exper_cs[iz]
163 G4cout <<
" not found in simulations " << nmatched << G4endl
164 <<
" not found in exper: " << nused <<
G4endl;
167 for(
G4int i = 0; i < simulatedAsSize; i++) {
183 lhood = std::pow(10.0, std::sqrt(lhood/matched));
190 <<
" lhood " << lhood <<
G4endl;
197 G4cout <<
" exper. cs " << exp_cs <<
" err " << exp_cs_err << G4endl
198 <<
" inucl. cs " << inucl_cs <<
" err " << inucl_cs_err << G4endl
199 <<
" ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "