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

fix of wrong offset in time domain data before Fourier.

Browse Source
This commit is contained in:
suter_a 2019-04-11 16:46:13 +02:00
parent ac1d52a054
commit 1c5069bc48
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/classes/PFourier.cpp
View File

@ -893,29 +893,29 @@ void PFourier::PrepareFFTwInputData(UInt_t apodizationTag)
Double_t mean = 0.0; Double_t mean = 0.0;
if (fDCCorrected) { if (fDCCorrected) {
for (UInt_t i=start; i<fNoOfData; i++) { for (UInt_t i=start; i<start+fNoOfData; i++) {
mean += fData->GetBinContent(i); mean += fData->GetBinContent(i);
} }
mean /= (Double_t)(fNoOfData-start); mean /= (Double_t)fNoOfData;
} }
// 2nd fill fIn // 2nd fill fIn
if (fUseFFTW) { if (fUseFFTW) {
for (UInt_t i=0; i<fNoOfData-start; i++) { for (UInt_t i=0; i<fNoOfData; i++) {
fIn[i][0] = fData->GetBinContent(i+start) - mean; fIn[i][0] = fData->GetBinContent(i+start) - mean;
fIn[i][1] = 0.0; fIn[i][1] = 0.0;
} }
for (UInt_t i=fNoOfData-start; i<fNoOfBins; i++) { for (UInt_t i=fNoOfData; i<fNoOfBins; i++) {
fIn[i][0] = 0.0; fIn[i][0] = 0.0;
fIn[i][1] = 0.0; fIn[i][1] = 0.0;
} }
} else { } else {
for (UInt_t i=0; i<fNoOfData-start; i++) { for (UInt_t i=0; i<fNoOfData; i++) {
#ifdef HAVE_DKS #ifdef HAVE_DKS
fInDKS[i] = fData->GetBinContent(i+start) - mean; fInDKS[i] = fData->GetBinContent(i+start) - mean;
#endif #endif
} }
for (UInt_t i=fNoOfData-start; i<fNoOfBins; i++) { for (UInt_t i=fNoOfData; i<fNoOfBins; i++) {
#ifdef HAVE_DKS #ifdef HAVE_DKS
fInDKS[i] = 0.0; fInDKS[i] = 0.0;
#endif #endif