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

make sure that exported MusrRoot files have proper stats in the histo objects.

Browse Source
This commit is contained in:
suter_a 2023-08-31 08:22:37 +02:00
parent 0b3fc43123
commit 4d1c7c10c8
1 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/classes/PRunDataHandler.cpp
View File

@ -4647,9 +4647,13 @@ Bool_t PRunDataHandler::WriteMusrRootFile(TString fln)
str.Form("hDecay%03d", dataSet->GetHistoNo()); str.Form("hDecay%03d", dataSet->GetHistoNo());
length = dataSet->GetData()->size(); length = dataSet->GetData()->size();
histo = new TH1F(str.Data(), str.Data(), length+1, -0.5, static_cast<Double_t>(length)+0.5); histo = new TH1F(str.Data(), str.Data(), length+1, -0.5, static_cast<Double_t>(length)+0.5);
Int_t sum=0, entries=0;
for (UInt_t j=0; j<length; j++) { for (UInt_t j=0; j<length; j++) {
histo->SetBinContent(j+1, dataSet->GetData()->at(j)); entries = dataSet->GetData()->at(j);
histo->SetBinContent(j+1, entries);
sum += entries;
} }
histo->SetEntries(sum);
histos.push_back(histo); histos.push_back(histo);
} }
} else { // rebin > 1 } else { // rebin > 1
@ -4666,14 +4670,18 @@ Bool_t PRunDataHandler::WriteMusrRootFile(TString fln)
length = dataSet->GetData()->size(); length = dataSet->GetData()->size();
histo = new TH1F(str.Data(), str.Data(), static_cast<Int_t>(length/fAny2ManyInfo->rebin)+1, -0.5, static_cast<Double_t>(static_cast<Int_t>(length/fAny2ManyInfo->rebin))+0.5); histo = new TH1F(str.Data(), str.Data(), static_cast<Int_t>(length/fAny2ManyInfo->rebin)+1, -0.5, static_cast<Double_t>(static_cast<Int_t>(length/fAny2ManyInfo->rebin))+0.5);
dataCount = 0; dataCount = 0;
Int_t sum=0, entries=0;
for (UInt_t j=0; j<length; j++) { for (UInt_t j=0; j<length; j++) {
if ((j > 0) && (j % fAny2ManyInfo->rebin == 0)) { if ((j > 0) && (j % fAny2ManyInfo->rebin == 0)) {
dataCount++; dataCount++;
histo->SetBinContent(dataCount, dataRebin); histo->SetBinContent(dataCount, dataRebin);
dataRebin = 0; dataRebin = 0;
} }
dataRebin += static_cast<UInt_t>(dataSet->GetData()->at(j)); entries = dataSet->GetData()->at(j);
sum += entries;
dataRebin += static_cast<UInt_t>(entries);
} }
histo->SetEntries(sum);
histos.push_back(histo); histos.push_back(histo);
} }
} }