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

fixed a problem showing up at SL4

Browse Source
This commit is contained in:
nemu 2008-03-03 10:24:14 +00:00
parent 12c1aef4cd
commit 0d03962415
1 changed files with 17 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
src/tests/skewedGaussianTest/fakeData.cpp
View File

@ -103,10 +103,10 @@ int main(int argc, char *argv[])
cout << endl << ">> read p(B) file ..."; cout << endl << ">> read p(B) file ...";
// read pB input file and fill B and pB // read pB input file and fill B and pB
ifstream f; ifstream fpB;
f.open(pBFileName.Data(), iostream::in); fpB.open(pBFileName.Data(), iostream::in);
if (!f.is_open()) { if (!fpB.is_open()) {
cout << endl << "**ERROR**: Sorry, couldn't open input file (p(B) file): " << pBFileName.Data(); cout << endl << "**ERROR**: Sorry, couldn't open input file (p(B) file): " << pBFileName.Data();
cout << endl << " Will quit"; cout << endl << " Will quit";
cout << endl; cout << endl;
@ -116,9 +116,9 @@ int main(int argc, char *argv[])
char str[256]; char str[256];
Double_t bb, ppb; Double_t bb, ppb;
int status, lineNo=0; int status, lineNo=0;
while (!f.eof()) { while (!fpB.eof()) {
// read a line // read a line
f.getline(str, sizeof(str)); fpB.getline(str, sizeof(str));
// ignore comments or empty lines // ignore comments or empty lines
if ((str[0] == '#') || (strlen(str)==0)) if ((str[0] == '#') || (strlen(str)==0))
@ -130,7 +130,7 @@ int main(int argc, char *argv[])
cout << endl << "**ERROR**: Problems while reading the input file (line no " << lineNo << "), will quit."; cout << endl << "**ERROR**: Problems while reading the input file (line no " << lineNo << "), will quit.";
cout << endl << " status = " << status; cout << endl << " status = " << status;
cout << endl; cout << endl;
f.close(); fpB.close();
return -1; return -1;
} }
@ -141,7 +141,7 @@ int main(int argc, char *argv[])
lineNo++; lineNo++;
} }
f.close(); fpB.close();
// normalize p(B) // normalize p(B)
Double_t sum = 0.0; Double_t sum = 0.0;
@ -163,9 +163,12 @@ int main(int argc, char *argv[])
cout << endl << ">> read parameter input file ..."; cout << endl << ">> read parameter input file ...";
// read parameter input file
ifstream fparam;
// open parameter input file and extract the parameters // open parameter input file and extract the parameters
f.open(paramInputFileName.Data(), iostream::in); fparam.open(paramInputFileName.Data(), iostream::in);
if (!f.is_open()) { if (!fparam.is_open()) {
cout << endl << "**ERROR**: Sorry, couldn't open parameter input file: " << paramInputFileName.Data(); cout << endl << "**ERROR**: Sorry, couldn't open parameter input file: " << paramInputFileName.Data();
cout << endl << " Will quit"; cout << endl << " Will quit";
cout << endl; cout << endl;
@ -179,9 +182,9 @@ int main(int argc, char *argv[])
Int_t ival; Int_t ival;
lineNo = 0; lineNo = 0;
while (!f.eof()) { while (!fparam.eof()) {
// read a line // read a line
f.getline(str, sizeof(str)); fparam.getline(str, sizeof(str));
lineNo++; lineNo++;
// ignore comments or empty lines // ignore comments or empty lines
@ -261,7 +264,7 @@ int main(int argc, char *argv[])
} }
} }
f.close(); fparam.close();
cout << endl << ">> number of pB's = " << pB.size(); cout << endl << ">> number of pB's = " << pB.size();
cout << endl << ">> number of channels = " << noOfChannels; cout << endl << ">> number of channels = " << noOfChannels;
@ -444,9 +447,9 @@ int main(int argc, char *argv[])
runHeader->SetTimeResolution(fval); runHeader->SetTimeResolution(fval);
runHeader->SetNChannels(noOfChannels); runHeader->SetNChannels(noOfChannels);
runHeader->SetNHist(histoData.size()); runHeader->SetNHist(histoData.size());
Int_t *t0array = new Int_t[histoData.size()]; Double_t *t0array = new Double_t[histoData.size()];
for (UInt_t i=0; i<histoData.size(); i++) for (UInt_t i=0; i<histoData.size(); i++)
t0array[i] = t0[i]; t0array[i] = (Double_t)t0[i];
runHeader->SetTimeZero(t0array); runHeader->SetTimeZero(t0array);
if (t0array) if (t0array)
delete t0array; delete t0array;