Beispiel #1
0
    else:
        time_pdf = SimultaneousPdf(time_pdf.GetName() + '_simul'
                                   , MasterPdf       = time_pdf
                                   , SplitCategories = split_cats
                                   , SplitParameters = split_pars)
        st_mean = lambda st: 0.04921 + st * 1.220
        st_sigma = lambda st: 0.012 + st * 0.165

        pars = time_pdf.getParameters(RooArgSet())
        for m, s in zip(means, split_cat):
            for n, f in (('timeResSFMean', st_mean), ('timeResSFSigma', st_sigma)):
               p = pars.find(n + '_' + s.GetName())
               p.setVal(f(m - split_obs_mean))

if (not options.simultaneous and options.mu_param) or 'sigmat' in options.mu_param:
    placeholder = sig_tres.muPlaceHolder()
    placeholder.setVal(st_means[year])

if options.sf_param:
    placeholder = sig_tres.sfPlaceHolder()
    placeholder.setVal(st_means[year])

if options.plot:
    cats = RooArgSet(st)
    if options.simultaneous:
        cats.add(time_pdf.indexCat())
        st_data = fit_data_full.reduce(cats)
    else:
        st_data = fit_data.reduce(cats)
    
if options.reuse_result and options.cache: