The reason that predicted distributions (esp. of c-cbar) never agreed with generated was that the physics functions used in the prediction took VPDL and σ(vpdl) as inputs. The physics function for c-cbar used a single gaussian for the VPDL distribution, which is not correct. The single gaussian applies to the Lxy distribution.
To make comparisons with binned chi2 fits, need to retain the capability to use VPDL, σ(vpdl) as the event variables. So change meaning of physics functions and resolution variables to take this into account.