(muclear = 0.028, sigmaclear = 0.005) are defined as variables at the top of the file.
Instead of throwing random to clear erf, use erf value as weight.
Plot dvv and dvvc with the binning used for the fit.
Option to read (muclear, sigmaclear) from argv; overlay dvv with binning; compute and print chi2 of dvvc to dvv.
Plot the relative difference between dvvc and dvv.