I’m trying to morph a set of templates I have based on histograms containing pT fractions. I’m currently passing these histograms into RooMomentMorph as RooDataHists, but the resulting morphed template has some jagged regions https://byates.web.cern.ch/byates/Top2016/2016/test/mtop/tree/sfs/fix/meson/morph/TemplateMorph_172v5_d0_rB_755-1055.png. When I try to produce an NLL the scan is not smooth and often get “function value is NAN”. It seems like RooIntegralMorph or even PiecewiseInterpolation are better suited for this, but they only seem to accept two templates (low end and high end) and I have four. I’ve also head people mention HistFactory, but this seems to be only for templates that are +/- n*sigma.
Does anyone have experience with smoothing morphing binned histograms?
Thanks.