LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 4 Activity

error feed fStep rather than fPosError for expected chisq. Tweaked the param rounding for expected chisq.

Browse Source
This commit is contained in:
suter_a
2025-07-04 16:29:33 +02:00
parent d9016d82a8
commit 2e951d92f1
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/classes/PFitter.cpp
View File

@ -2178,7 +2178,7 @@ Bool_t PFitter::ExecuteSave(Bool_t firstSave)
for (UInt_t i=0; i<fParams.size(); i++) { for (UInt_t i=0; i<fParams.size(); i++) {
param.push_back(fParams[i].fValue); param.push_back(fParams[i].fValue);
err.push_back(fParams[i].fPosError); err.push_back(fParams[i].fStep);
} }
// CalcExpectedChiSquare handles both, chisq and mlh // CalcExpectedChiSquare handles both, chisq and mlh
@ -2944,7 +2944,7 @@ PDoubleVector PFitter::ParamRound(const PDoubleVector &par, const PDoubleVector
double dval; double dval;
for (unsigned int i=0; i<par.size(); i++) { for (unsigned int i=0; i<par.size(); i++) {
if (err[i] != 0.0) { if (err[i] != 0.0) {
exp = floor(log10(fabs(err[i]))-2); exp = floor(log10(fabs(err[i]))-1);
dval = round(par[i]*pow(10.0, -exp))/pow(10.0, -exp); dval = round(par[i]*pow(10.0, -exp))/pow(10.0, -exp);
par_r[i] = dval; par_r[i] = dval;
} else { } else {