def makeNumPro(method):
c = Canvas('c')
c.cd()
c.SetTicks(1,1)
c.SetLogy()
h_16 = getHist(method,['2016'],'h_num_pro')
h_17 = getHist(method,['2017'],'h_num_pro')
h_18 = getHist(method,['2018'],'h_num_pro')
h_18.SetTitle(''), h_18.GetXaxis().SetTitle('Number Of Reconstructed Protons'), h_18.GetYaxis().SetTitle('Events')
h_18.SetMarkerStyle(20), h_17.SetMarkerStyle(20), h_16.SetMarkerStyle(20)
h_18.SetMarkerColor(ROOT.darkLime), h_17.SetMarkerColor(ROOT.turquois), h_16.SetMarkerColor(ROOT.orangeSoda)
h_18.SetMarkerSize(1.5), h_17.SetMarkerSize(1.5), h_16.SetMarkerSize(1.5)
h_18.SetLineColor(ROOT.darkLime), h_17.SetLineColor(ROOT.turquois), h_16.SetLineColor(ROOT.orangeSoda)
h_18.SetLineWidth(2), h_17.SetLineWidth(2), h_16.SetLineWidth(2)
h_18.SetMinimum(10000.0)
h_18.Draw('p')
h_17.Draw('p same')
h_16.Draw('p same')
c.SetGrid(0,1)
legend = TLegend(0.6,0.7,0.8,0.78)
legend.SetTextSize(0.03)
legend.SetLineColor( 0 )
legend.SetFillColor( 0 )
legend.AddEntry(h_18,"2018",'lp')
legend.AddEntry(h_17,"2017",'lp')
legend.AddEntry(h_16,"2016",'lp')
legend.Draw()
pLabel = prelimLabel('top',True,h_18.GetMaximum())
pLabel.Draw()
c.SaveAs('plots/combined/h_num_pro_%s.png' % (selection))
def makeXiAcceptance(pot,log):
c = Canvas('c')
c.cd()
if log: c.SetLogy()
c.SetTicks(1,1)
c.SetGrid(1,1)
h = ROOT.TH1F('h', '', 100, 0, 0.2)
for year in years:
for pf in protonFiles:
if year == pf[0]: pfile = pf[1]
h.Add( pfile.Get('plots/h_pro_xi_%s' % pot) )
h.Rebin(2)
#h.Scale(1.0/100.0)
for pp in proton_pot:
if pot == pp[0]: v_pot = pp
h.GetXaxis().SetTitle('#xi '+v_pot[1])
h.GetYaxis().SetTitle('Events')
h.SetLineColor(ROOT.kBlack)
h.SetMarkerStyle(24)
minx, miny, maxx, maxy = 0.0001, h.GetMinimum(), 0.2, h.GetMaximum()*10 if log else h.GetMaximum()*1.2
c.DrawFrame(minx, miny, maxx, maxy)
h.Draw('p e2')
h.SetMaximum(maxy)
pLabel, lLabel = prelimLabel('top',log,maxy), lumiLabel(log,years)
lLabel.SetTextSize(0.034)
pLabel.Draw(), lLabel.Draw()
c.SaveAs('plots/%s/h_xi_%s_%s.pdf' % (s_years,side,s_years))
def makeDetType():
c = Canvas('c')
c.cd()
c.SetTicks(1,1)
h_type = ROOT.TH1F('h_type', '', 2, 3, 5)
for year in years:
for pf in protonFiles:
if year == pf[0]: pfile = pf[1]
h_type.Add( pfile.Get('plots/h_detType') )
h_type.Scale(1.0/1000.0)
h_type.GetXaxis().SetBinLabel(1,'Si Strip'), h_type.GetXaxis().SetBinLabel(2,'Pixel')
h_type.GetYaxis().SetTitle('Events/0.001')
h_type.Draw('p')
h_type.SetLineColor(ROOT.kBlack)
h_type.SetLineWidth(2)
h_type.SetMarkerStyle(20)
h_type.SetMinimum(0.0), h_type.SetMaximum( h_type.GetMaximum() * 1.5 )
h_type.SetMarkerSize(0.9)
c.SetGrid(0,1)
pLabel, sLabel, lLabel = prelimLabel('left',False,h_type.GetMaximum()), selectionLabel('singleRP',False,False,h_type.GetMaximum()), lumiLabel(False,years)
pLabel.Draw(), sLabel.Draw(), lLabel.Draw()
c.SaveAs('plots/%s/h_detType_%s_%s.png' % (s_years,selection,s_years))
def makeProtonSide(log):
c = Canvas('c')
c.cd()
c.SetLeftMargin(0.15)
if log: c.SetLogy()
c.SetTicks(1,1)
h_single = ROOT.TH1F('h_single', '', 4, 0, 4)
h_multi = ROOT.TH1F('h_multi', '', 4, 0, 4)
for year in years:
for pf in protonFiles:
if year == pf[0]: pf_single, pf_multi = pf[1], pf[2]
h_single.Add( pf_single.Get('plots/h_proton_side') )
h_multi.Add( pf_multi.Get('plots/h_proton_side') )
denom = h_single.GetEntries()
h_multi.Scale(1.0/denom)
h_multi.GetYaxis().SetTitle('Fraction Of Events')
h_multi.GetXaxis().SetBinLabel(1,'No protons'), h_multi.GetXaxis().SetBinLabel(2,'sector45 only'), h_multi.GetXaxis().SetBinLabel(3,'sector56 only'), h_multi.GetXaxis().SetBinLabel(4,'Both')
h_multi.Draw('p')
h_single.Scale(1.0/denom)
h_single.Draw('p same')
h_single.SetMarkerColor(62), h_multi.SetMarkerColor(207)
h_single.SetLineColor(62), h_multi.SetLineColor(207)
h_single.SetLineWidth(2), h_multi.SetLineWidth(2)
h_single.SetMarkerStyle(20), h_multi.SetMarkerStyle(20)
h_multi.SetMinimum(0.0), h_multi.SetMaximum( h_multi.GetMaximum()*1.5 )
h_single.SetMarkerSize(0.9), h_multi.SetMarkerSize(0.9)
c.SetGrid(0,1)
legend = TLegend(0.6,0.7,0.8,0.78)
legend.SetTextSize(0.03)
legend.SetLineColor( 0 )
legend.SetFillColor( 0 )
legend.AddEntry(h_single,"singleRP",'lp')
legend.AddEntry(h_multi,"multiRP",'lp')
legend.Draw()
pLabel, lLabel = prelimLabel('top',log,h_single.GetMaximum()), lumiLabel(False,years)
pLabel.SetMargin(0.49)
pLabel.Draw(), lLabel.Draw()
c.SaveAs('plots/%s/h_proton_side_%s_%s.pdf' % (s_years,selection,s_years))
def makeXiAcceptanceSide(side,log):
c = Canvas('c')
c.cd()
if log: c.SetLogy()
c.SetTicks(1,1)
c.SetGrid(1,1)
h = ROOT.TH1F('h', '', 100, 0, 0.3)
for year in years:
for pf in protonFiles:
if year == pf[0]: pfile = pf[2]
h.Add( pfile.Get('plots/h_pro_xi%s' % side) )
#h.Rebin(2)
h.GetXaxis().SetTitle('#xi ^{%s}' % ('+' if side == 'p' else '-'))
h.GetYaxis().SetTitle('Events')
h.SetLineColor(ROOT.kBlack)
h.SetMarkerStyle(24)
minx, miny, maxx, maxy = 0.0001, h.GetMinimum(), 0.2, h.GetMaximum()*10 if log else h.GetMaximum()*1.2
c.DrawFrame(minx, miny, maxx, maxy)
h.Draw('p e2')
h.SetMaximum(maxy)
pLabel, lLabel = prelimLabel('top',log,maxy), lumiLabel(log,years)
lLabel.SetTextSize(0.034)
pLabel.Draw(), lLabel.Draw()
c.SaveAs('plots/%s/h_xi%s_%s.pdf' % (s_years,side,s_years))