c.SaveAs(region + "/mmtt_pt_histo_" + region + ".root") if args.nofit and args.nofitkk: ## Phi fitting kcanvas = TCanvas("kcanvas", "kcanvas", 1400, 800) pullpad = TPad("pullpad", "pullpad", 0.0, 0.05, 1.0, 0.33) plotpad = TPad("histopad", "histopad", 0.0, 0.35, 1.0, 1.0) plotpad.SetFillStyle(4004) pullpad.SetFillStyle(4004) plotpad.Draw() pullpad.Draw() traKFitData = theData.Clone("fitTrakData") if binnedfit: tt_mass.setBins(30) traKFitData = theData.binnedClone("binnedTrakData") phimean = 1.019 gammavalue = 0.0012 fitphimin = 1.00 fitphimax = 1.04 kkSigma = RooRealVar("#sigma", "#sigma", 0.0013) kkGamma = RooRealVar("#Gamma", "#Gamma", gammavalue, 0.001, 0.015) kkMean = RooRealVar("m_{kk}", "m_{kk}", phimean, phimean - 0.005, phimean + 0.005)
hist_mm_trig.Draw("SAME") gPad.BuildLegend(0.65, 0.5, 0.85, 0.65) c.SaveAs("mm_mass_histos.png") c.SaveAs("mm_mass_histos.root") if args.nofit: sys.exit() #### #### Done plotting #### #### Fitting ### Trak trak # traKFitData = trakData.Clone("fitTrakData") if binnedfit: tt_mass.setBins(100) traKFitData = trakData.binnedClone("binnedTrakData") phimean = 1.019 gammavalue = 0.01 fitphimin = 1.0 fitphimax = 1.04 if args.nofitkk and args.noreco: ## RECO