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

check for parameter names adopted to minuit2 limitation

Browse Source
This commit is contained in:
nemu 2008-05-13 05:56:13 +00:00
parent 1ec5e6e6c0
commit e3ab2a21ab
1 changed files with 10 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/classes/PMsrHandler.cpp
View File

@ -220,6 +220,7 @@ int PMsrHandler::ReadMsrFile()
unsigned int parX, parY; unsigned int parX, parY;
if (!CheckUniquenessOfParamNames(parX, parY)) { if (!CheckUniquenessOfParamNames(parX, parY)) {
cout << endl << "PMsrHandler::ReadMsrFile: **SEVERE ERROR** parameter name " << fParam[parX].fName.Data() << " is identical for parameter no " << fParam[parX].fNo << " and " << fParam[parY].fNo << "!"; cout << endl << "PMsrHandler::ReadMsrFile: **SEVERE ERROR** parameter name " << fParam[parX].fName.Data() << " is identical for parameter no " << fParam[parX].fNo << " and " << fParam[parY].fNo << "!";
cout << endl << "It might be that only the 10 first characters are identical. Due to a current MINUIT2 limitation, this means they are \"identical\", sorry.";
cout << endl << "Needs to be fixed first!"; cout << endl << "Needs to be fixed first!";
error = PMUSR_MSR_SYNTAX_ERROR; error = PMUSR_MSR_SYNTAX_ERROR;
} }
@ -2265,10 +2266,19 @@ void PMsrHandler::FillParameterInUse(PMsrLines &theory, PMsrLines &funcs, PMsrLi
bool PMsrHandler::CheckUniquenessOfParamNames(unsigned int &parX, unsigned int &parY) bool PMsrHandler::CheckUniquenessOfParamNames(unsigned int &parX, unsigned int &parY)
{ {
bool unique = true; bool unique = true;
char name1[11], name2[11];
for (unsigned int i=0; i<fParam.size()-1; i++) { for (unsigned int i=0; i<fParam.size()-1; i++) {
for (unsigned int j=i+1; j<fParam.size(); j++) { for (unsigned int j=i+1; j<fParam.size(); j++) {
// due to the current minuit2 limitation of 10 characters for the parameter name
memset(name1, 0, sizeof(name1));
strncpy(name1, fParam[i].fName.Data(), 10);
memset(name2, 0, sizeof(name2));
strncpy(name2, fParam[j].fName.Data(), 10);
/* as soon as the minuit2 parameter name limit is lifted, the if should become this:
if (fParam[i].fName.CompareTo(fParam[j].fName) == 0) { // equal if (fParam[i].fName.CompareTo(fParam[j].fName) == 0) { // equal
*/
if (!strcmp(name1, name2)) { // equal
unique = false; unique = false;
parX = i; parX = i;
parY = j; parY = j;