I’m having trouble getting RooFFTConvPdf to work with a RooProdPdf, when trying to convolve the RooProdPdf with a RooGaussian. The RooProdPdf is the product of two RooGenericPdf’s (a box function and a “spectrum”).
As the attached macro can demonstrate, the convolution seems to work fine when I use the RooNumConvPdf class (except for a dependence on whether the Gaussian is the first or second pdf argument). However, when I try to use RooFFTConvPdf, the code returns a large number of errors, and produces no final result. The FFT convolution works fine for the Gaussian with each of the box and spectrum functions, just not the product of the two.
Would someone be able to point out what I’m doing wrong? or if this is a bug?
prodConv.C (1.5 KB)