# 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)