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

fixed start != t0 phase issue

Browse Source
This commit is contained in:
suter_a 2018-10-18 13:54:05 +02:00
parent b9e0812699
commit 06ae0a12ea
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/classes/PFourier.cpp
View File

@ -312,7 +312,7 @@ PFourier::PFourier(TH1F *data, Int_t unitTag, Double_t startTime, Double_t endTi
// if endTime == 0 set it to the last time slot
if (fEndTime == 0.0) {
Int_t last = fData->GetNbinsX();
Int_t last = fData->GetNbinsX()-1;
fEndTime = fData->GetBinCenter(last);
}
@ -901,21 +901,21 @@ void PFourier::PrepareFFTwInputData(UInt_t apodizationTag)
// 2nd fill fIn
if (fUseFFTW) {
for (UInt_t i=0; i<fNoOfData-t0bin; i++) {
fIn[i][0] = fData->GetBinContent(i+t0bin) - mean;
for (UInt_t i=0; i<fNoOfData-start; i++) {
fIn[i][0] = fData->GetBinContent(i+start) - mean;
fIn[i][1] = 0.0;
}
for (UInt_t i=fNoOfData-t0bin; i<fNoOfBins; i++) {
for (UInt_t i=fNoOfData-start; i<fNoOfBins; i++) {
fIn[i][0] = 0.0;
fIn[i][1] = 0.0;
}
} else {
for (UInt_t i=0; i<fNoOfData-t0bin; i++) {
for (UInt_t i=0; i<fNoOfData-start; i++) {
#ifdef HAVE_DKS
fInDKS[i] = fData->GetBinContent(i+t0bin) - mean;
fInDKS[i] = fData->GetBinContent(i+start) - mean;
#endif
}
for (UInt_t i=fNoOfData-t0bin; i<fNoOfBins; i++) {
for (UInt_t i=fNoOfData-start; i<fNoOfBins; i++) {
#ifdef HAVE_DKS
fInDKS[i] = 0.0;
#endif