lcover_p.Fill(cval, lval, float(lcoverp_ok) / cover_tot) mtot = max(tot, mtot) if not args.asimov else 1 fcn = BasePlotter.parse_formula( 'amplitude*TMath::Gaus(x, mean, width)', 'amplitude[100, 0, 1000], mean[1, 0, 2], width[0.5, 0, 2]', [0, 2]) if args.conly: failing.Fill(cval, failed) else: failing.Fill(cval, lval, failed) xb = cbias.xaxis.find_bin(cval) if not args.conly: yb = cbias.yaxis.find_bin(lval) if args.usemean: bias = ufloat(csf.GetMean(), 0) elif args.asimov: bias = ufloat(csf, 0) else: csf.Fit(fcn, 'LQN0') bias = ufloat(fcn['mean'].value, fcn['mean'].error) bias -= cval bias /= cval bin = cbias[xb] if args.conly else cbias[xb, yb] bin.value = bias.nominal_value bin.error = bias.std_dev if not args.conly: if args.usemean: bias = ufloat(lsf.GetMean(), 0) elif args.asimov: