Beispiel #1
0
st_means = {'2011' : 0.350, '2012' : 0.349}
st_means[year] = sig_sdata.mean(st._target_())

if options.simultaneous:
    split_pars = [[]]
    if options.wpv:
        for par in time_pdf.Parameters():
            if not par.getAttribute('Yield'):
                continue
            if par.GetName().find(components[-1].GetName()) != -1:
                continue
            split_pars[0].append(par)
    else:
        split_pars[0] = [par for par in time_pdf.Parameters() if par.getAttribute('Yield')]
    split_pars[0] += sig_tres.splitVars()

    if options.wpv and options.wpv_type == 'Gauss':
        split_pars[0].extend([sig_wpv.getYield()])
        ## split_pars[0].extend([wpv_sigma])
    elif options.wpv and signal_MC and options.wpv_type == 'Rest':
        split_pars[0].append(rest.getYield())
        ## split_pars[0].append(rest_tres._left_rlifeSF)
        ## split_pars[0].append(rest_tres._frac_left)
    
    ## Calculate the mean in each bin for the splitting observable
    split_obs = split_util.observables()[0]
    split_obs_mean = st_means[year]
    split_cat = split_cats[0][0]
    if hasattr(split_cat, '_target_'):
        split_cat = split_cat._target_()