for i in xrange(graph.GetN()) :
ggH_list.append(float(graph.GetX()[i]))
bbH_list.append(float(graph.GetY()[i]))
ggH_list = sorted(set(ggH_list))
bbH_list = sorted(set(bbH_list))
ggH_bins=len(ggH_list)
bbH_bins=len(bbH_list)
print "r_ggH_list: ", ggH_list, "Total number: ", ggH_bins
print "r_bbH_list: ", bbH_list, "Total number: ", bbH_bins
#Create canvas and TH2D for each component
plot.ModTDRStyle(width=600, l=0.12)
#Slightly thicker frame to ensure contour edges dont overlay the axis
ROOT.gStyle.SetFrameLineWidth(3)
c1=ROOT.TCanvas()
axis = plot.makeVarBinHist2D('hist2d', ggH_list, bbH_list)
axis.GetXaxis().SetTitle("#sigma#font[42]{(gg#phi)}#upoint#font[52]{B}#font[42]{(#phi#rightarrow#tau#tau)} [pb]")
axis.GetYaxis().SetTitle("#sigma#font[42]{(bb#phi)}#upoint#font[52]{B}#font[42]{(#phi#rightarrow#tau#tau)} [pb]")
#Create two pads, one is just for the Legend
pad1 = ROOT.TPad("pad1","pad1",0,0.5,1,1)
pad1.SetFillStyle(4000)
pad1.Draw()
pad2 = ROOT.TPad("pad2","pad2",0,0,1,1)
pad2.SetFillStyle(4000)
pad2.Draw()
pads=[pad1,pad2]
pads[1].cd()
bfit = plot.bestFit(tree, "r_ggH", "r_bbH", "")
#Note the binning of the TH2D for the interpolation should ~ match the initial input grid
#Could in future implement variable binning here
for i in xrange(graph.GetN()):
ggH_list.append(float(graph.GetX()[i]))
bbH_list.append(float(graph.GetY()[i]))
ggH_list = sorted(set(ggH_list))
bbH_list = sorted(set(bbH_list))
ggH_bins = len(ggH_list)
bbH_bins = len(bbH_list)
print "r_ggH_list: ", ggH_list, "Total number: ", ggH_bins
print "r_bbH_list: ", bbH_list, "Total number: ", bbH_bins
#Create canvas and TH2D for each component
plot.ModTDRStyle(width=600, l=0.12)
#Slightly thicker frame to ensure contour edges dont overlay the axis
ROOT.gStyle.SetFrameLineWidth(3)
c1 = ROOT.TCanvas()
axis = plot.makeVarBinHist2D('hist2d', ggH_list, bbH_list)
axis.GetXaxis().SetTitle(
"#sigma#font[42]{(gg#phi)}#upoint#font[52]{B}#font[42]{(#phi#rightarrow#tau#tau)} [pb]"
)
axis.GetYaxis().SetTitle(
"#sigma#font[42]{(bb#phi)}#upoint#font[52]{B}#font[42]{(#phi#rightarrow#tau#tau)} [pb]"
)
#Create two pads, one is just for the Legend
pad1 = ROOT.TPad("pad1", "pad1", 0, 0.5, 1, 1)
pad1.SetFillStyle(4000)
pad1.Draw()
pad2 = ROOT.TPad("pad2", "pad2", 0, 0, 1, 1)
pad2.SetFillStyle(4000)
pad2.Draw()
pads = [pad1, pad2]
pads[1].cd()
tanb_list.append(float(graph_exp.GetY()[i]))
tanb_list = sorted(set(tanb_list))
mA_list = sorted(set(mA_list))
tanb_bins = len(tanb_list)
mA_bins = len(mA_list)
if int(args.verbosity) > 0:
print "mA_list: ", mA_list, "Total number: ", mA_bins
print "tanb_list: ", tanb_list, "Total number: ", tanb_bins
#Create canvas and TH2D for each component
plot.ModTDRStyle(width=600, l=0.12)
#Slightly thicker frame to ensure contour edges dont overlay the axis
ROOT.gStyle.SetFrameLineWidth(2)
c1 = ROOT.TCanvas()
#axis = plot.makeHist2D('hist2d', mA_bins, tanb_bins, graph_exp)
axis = plot.makeVarBinHist2D('hist2d', mA_list, tanb_list)
axis.GetYaxis().SetTitle("tan#beta")
axis.GetXaxis().SetTitle("m_{A} (GeV)")
#Create two pads, one is just for the Legend
pad1 = ROOT.TPad("pad1", "pad1", 0, 0.82, 1, 1)
pad1.SetFillStyle(4000)
pad1.Draw()
pad2 = ROOT.TPad("pad2", "pad2", 0, 0, 1, 0.82)
pad2.SetFillStyle(4000)
pad2.Draw()
pads = [pad1, pad2]
pads[1].cd()
#Note the binning of the TH2D for the interpolation should ~ match the initial input grid
#Could in future implement variable binning here
#h_exp = plot.makeVarBinHist2D("h_exp", mA_list, tanb_list)
tanb_list.append(float(graph_exp.GetY()[i]))
tanb_list = sorted(set(tanb_list))
mA_list = sorted(set(mA_list))
tanb_bins=len(tanb_list)
mA_bins=len(mA_list)
if int(args.verbosity) > 0 :
print "mA_list: ", mA_list, "Total number: ", mA_bins
print "tanb_list: ", tanb_list, "Total number: ", tanb_bins
#Create canvas and TH2D for each component
plot.ModTDRStyle(width=600, l=0.12)
#Slightly thicker frame to ensure contour edges dont overlay the axis
ROOT.gStyle.SetFrameLineWidth(2)
c1=ROOT.TCanvas()
#axis = plot.makeHist2D('hist2d', mA_bins, tanb_bins, graph_exp)
axis = plot.makeVarBinHist2D('hist2d', mA_list, tanb_list)
axis.GetYaxis().SetTitle("tan#beta")
axis.GetXaxis().SetTitle("m_{A} (GeV)")
#Create two pads, one is just for the Legend
pad1 = ROOT.TPad("pad1","pad1",0,0.82,1,1)
pad1.SetFillStyle(4000)
pad1.Draw()
pad2 = ROOT.TPad("pad2","pad2",0,0,1,0.82)
pad2.SetFillStyle(4000)
pad2.Draw()
pads=[pad1,pad2]
pads[1].cd()
#Note the binning of the TH2D for the interpolation should ~ match the initial input grid
#Could in future implement variable binning here
#h_exp = plot.makeVarBinHist2D("h_exp", mA_list, tanb_list)
