Hi,
I think you have continued in another post, No matching function for call to 'bindFunction'.
You should probably keep the same one.
Trying your code, the main problem is that your pdf function you want to convolute with a Gaussian is too sharp, it is almost as a delta function and you cannot perform a numerical convolution with such a peak. I think in this case you can probably assume it is a delta function and then your final pdf will be a Gaussian centred at the peak location.
Lorenzo