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

some additional parameter checks

Browse Source
This commit is contained in:
suter_a 2012-09-12 08:25:27 +00:00
parent da0d0e6ed9
commit 9095f82c33
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/classes/PTheory.cpp
View File

@ -1397,6 +1397,12 @@ Double_t PTheory::DynamicGaussKTLF(register Double_t t, const PDoubleVector& par
if ((val[0] == 0.0) && (val[1] == 0.0) && (val[2] == 0.0))
return 1.0;
// make sure that damping and hopping are positive definite
if (val[1] < 0.0)
val[1] = -val[1];
if (val[2] < 0.0)
val[2] = -val[2];
// check that Delta != 0, if not (i.e. stupid parameter) return 1, which is the correct limit
if (fabs(val[1]) < 1.0e-6) {
return 1.0;
@ -1628,6 +1634,13 @@ Double_t PTheory::DynamicLorentzKTLF(register Double_t t, const PDoubleVector& p
if ((val[0] == 0.0) && (val[1] == 0.0) && (val[2] == 0.0))
return 1.0;
// make sure that damping and hopping are positive definite
if (val[1] < 0.0)
val[1] = -val[1];
if (val[2] < 0.0)
val[2] = -val[2];
Double_t tt;
if (fParamNo.size() == 3) // no tshift
tt = t;