Assumes PMs are a convolution of normal and exponential. So we
observe X+Y where X is background and Y is signal. bg.adjust
returns E[Y|X+Y, Y>0] as our background corrected PM.
bg.parameters provides ad hoc estimates of the parameters of the
normal and exponential distributions.