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

fix 'possible msr2data in global+ mode bug' bug report #16 gitea
Build and Deploy Documentation / build-and-deploy (push) Successful in 30s
Details

Browse Source
This commit is contained in:
suter_a
2026-08-23 09:51:46 +02:00
parent 08559fdb59
commit 3236df1678
1 changed files with 14 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/classes/PMsr2Data.cpp
+14 -14
View File
@@ -1453,14 +1453,22 @@ bool PMsr2Data::PrepareGlobalInputFile(unsigned int tempRun, const std::string &
//std::cout << "Number of run specific parameters: " << fNumSpecParam << std::endl;
if (newRunNumber.str().compare(tempRunNumber.str())) { // first run number does not match the template run number
// in global+ mode, read in the single run msr-file of the corresponding run
if (globalPlus) {
singleRunMsrFile = GetSingleRunMsrFile();
if (singleRunMsrFile == nullptr)
return false;
// in global+ mode, read in the single-run msr-file of the first run and use its
// fitted values/steps as starting values for the run-specific parameters -
// also when the first run coincides with the template run
if (globalPlus) {
singleRunMsrFile = GetSingleRunMsrFile();
if (singleRunMsrFile == nullptr)
return false;
for (unsigned int l(fNumGlobalParam); l < msrParamList->size(); ++l) {
fMsrHandler->SetMsrParamValue(l, singleRunMsrFile->GetMsrParamList()->at(l).fValue);
fMsrHandler->SetMsrParamStep(l, singleRunMsrFile->GetMsrParamList()->at(l).fStep);
}
delete singleRunMsrFile;
singleRunMsrFile = nullptr;
}
if (newRunNumber.str().compare(tempRunNumber.str())) { // first run number does not match the template run number
// substitute the template run-numbers in the parameter names
for (unsigned int l(fNumGlobalParam); l < msrParamList->size(); ++l) {
tempParamName = msrParamList->at(l).fName.Data();
@@ -1468,20 +1476,12 @@ bool PMsr2Data::PrepareGlobalInputFile(unsigned int tempRun, const std::string &
if ( loc != std::string::npos ) {
tempParamName.replace(loc, fRunNumberDigits, newRunNumber.str());
msrParamList->at(l).fName = tempParamName;
if (globalPlus && singleRunMsrFile) {
fMsrHandler->SetMsrParamValue(l, singleRunMsrFile->GetMsrParamList()->at(l).fValue);
fMsrHandler->SetMsrParamStep(l, singleRunMsrFile->GetMsrParamList()->at(l).fStep);
}
} else {
std::cerr << std::endl << ">> msr2data: **ERROR** The indices of the run specific parameters do not match the template run number!";
std::cerr << std::endl << ">> msr2data: **ERROR** This should not happen! Please report a bug!";
std::cerr << std::endl;
}
}
if (singleRunMsrFile) {
delete singleRunMsrFile;
singleRunMsrFile = nullptr;
}
// substitute the template run-numbers in the RUN names
for (unsigned int i(0); i < fNumTempRunBlocks; ++i) {