Exemple #1
0
                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: