# binedges1[i] = 10**(log10(2E-4) + i*logwidth/nbins1)
binedges1[i] = epsilon + i*1E-2/nbins1
# print '%2d: %e' % (i, binedges1[i])
for i in range(nbins2+1):
if i <= 15:
binedges2[i] = epsilon + i*2E-4
elif i <= 19:
binedges2[i] = binedges2[15] + (i - 15)*5E-4
else:
binedges2[i] = binedges2[19] + (i - 19)*1E-3
# uniformbins = RooBinning(nbins1, binedges, 'uniformbins')
# Decay binned
decayfn = decaywdt.asTF(RooArgList(time), RooArgList(), timeargset)
# decayhist = decay.createHistogram('decayhist', time,
# RooFit.Binning(uniformbins),
# RooFit.ConditionalObservables(dtargset))
# decayhist.Sumw2()
# timehist = decayhist.Clone('timehist')
# timehist.Reset("ICESM")
timehist = TH1F('timehist', 'Lifetime distribution', nbins1, binedges1)
timehist.Sumw2()
timehist.SetLineColor(kAzure)
canvas1 = TCanvas('canvas1', 'Lifetime', 800, 600)
# canvas1.Divide(2, 1)
# canvas1.cd(1)
