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