BiasScaleFactor = False)

    res_model = MultiHistEfficiencyModel(Name = "RMHE", Original = sig_t.pdf(), Observable = t,
                                         ConditionalCategories = True, UseSingleBinConstraint = False,
                                         ResolutionModel = tres.model(), **spec)
    pdf = Single_Exponent_Time(Name = 'pdf', time = t, resolutionModel = res_model)
    pdf = pdf.pdf()

    pdf.Print('v')        
    data = pdf.generate([t, hlt1_excl_biased_dec, hlt2_unbiased, hlt2_biased], 30000)

## Fit
print 'fitting data'
## from profiler import profiler_start, profiler_stop
## profiler_start("acceptance.log")
result = pdf.fitTo(data, **fitOpts)
## profiler_stop()

from ROOT import kDashed, kRed, kGreen, kBlue, kBlack
from ROOT import TCanvas, RooBinning
canvas = {}
print 'plotting'

states_signal = set([(state, label) for d in valid_definition for state, label in d])
def sort_combination(combination):
    valid_def = valid_definition[:]
    valid_def.reverse()
    level_left = 0
    n = 0
    c = set(combination)
    for level in valid_def: