theData.plotOn(mmttFrame) mmttFrame.Draw() c.SaveAs("mmtt_mass" + region + ".png") #### X lxy Data print("TrakTrakMuMu data plotting . . .") print("All : " + str(theData.numEntries())) lxyFrame = lxysig.frame() theData.plotOn(lxyFrame) lxyFrame.Draw() c.SaveAs("lxy" + region + ".png") ##Plotting hist_mmtt = (theData.createHistogram( mmtt_mass, mmtt_mass)).ProjectionX("hist_mmtt_mass") hist_mm = (theData.createHistogram(mm_mass, mmtt_mass)).ProjectionX("hist_mm_mass") hist_tt = (theData.createHistogram(tt_mass, mmtt_mass)).ProjectionX("hist_tt_mass") hist_tt_pt = (theData.createHistogram(mm_pt, mmtt_mass)).ProjectionX("hist_mm_pt") hist_mm_pt = (theData.createHistogram(tt_pt, mmtt_mass)).ProjectionX("hist_tt_pt") hist_mmtt_pt = (theData.createHistogram( mmtt_pt, mmtt_mass)).ProjectionX("hist_mmtt_pt") ##Masses hists hist_mmtt.SetTitle("M(#mu#muKK)" + region + " ;M(#mu#muKK) [GeV]; no. entries")
theData.plotOn(mmttFrame) mmttFrame.Draw() c.SaveAs(region + "/mmtt_mass_" + region + ".png") #### X lxy Data print("TrakTrakMuMu data plotting . . .") print("All : " + str(theData.numEntries())) lxyFrame = lxysig.frame() theData.plotOn(lxyFrame) lxyFrame.Draw() c.SaveAs(region + "/lxy_" + region + ".png") ##Plotting hist_mmtt = (theData.createHistogram(mmtt_mass, mmtt_mass, 100, 100)).ProjectionX("hist_mmtt_mass") hist_mm = (theData.createHistogram(mm_mass, mmtt_mass, 100, 100)).ProjectionX("hist_mm_mass") hist_tt = (theData.createHistogram(tt_mass, mmtt_mass, 400, 100)).ProjectionX("hist_tt_mass") hist_mm_pt = (theData.createHistogram(mm_pt, mmtt_mass, int(ptmax) * 2, 100)).ProjectionX("hist_mm_pt") hist_tt_pt = (theData.createHistogram(tt_pt, mmtt_mass, int(ptmax) * 2, 100)).ProjectionX("hist_tt_pt") hist_mmtt_pt = (theData.createHistogram(mmtt_pt, mmtt_mass, int(ptmax) * 2, 100)).ProjectionX("hist_mmtt_pt") ##Masses hists
c.SaveAs("mmtt_mass_trigger" + region + ".png") #### X lxy Data print("TrakTrakMuMu data plotting . . .") print("All : " + str(mmttData.numEntries())) lxyFrame = lxysig.frame() mmttData.plotOn(lxyFrame) lxyFrame.Draw() c.SaveAs("lxy" + region + ".png") mmttcuts = [4.0, 4.2, 4.3, 4.4, 4.6, 4.8, 5.0] trigcands = [] patcands = [] hist_mmtt = (mmttData.createHistogram( mmtt_mass, mmtt_mass_trig)).ProjectionX("hist_mmtt") hist_mmtt_trig = (mmttData.createHistogram( mmtt_mass, mmtt_mass_trig)).ProjectionY("hist_mmtt_trig") mmtt_all = hist_mmtt.Integral(0, hist_mmtt.GetNbinsX()) mmtt_trig_all = hist_mmtt_trig.Integral(0, hist_mmtt_trig.GetNbinsX()) print("####TRIGGER") for i in range(len(mmttcuts)): integral = hist_mmtt_trig.Integral(hist_mmtt_trig.FindBin(mmttcuts[i]), hist_mmtt_trig.GetNbinsX()) print("#No. cands. > " + str(mmttcuts[i]) + " : " + str(integral) + " fraction : " + str(integral / mmtt_trig_all)) trigcands.append(integral) print(trigcands) print("####PAT")