leg_2d_exclusion_v1_v2.AddEntry(graph_exclusion_exp_v1, v1_label + " Exp", "L")
leg_2d_exclusion_v1_v2.AddEntry(graph_exclusion_exp_v2, v2_label + " Exp", "L")
if drawObs:
leg_2d_exclusion_v1_v2.AddEntry(graph_exclusion_obs_v1, v1_label + " Obs",
"L")
leg_2d_exclusion_v1_v2.AddEntry(graph_exclusion_obs_v2, v2_label + " Obs",
"L")
leg_2d_exclusion_v1_v2.Draw()
drawCMS2(myC2D, 13, lumi)
f1_lambda = TF1("f1", "(x+6.00)/1.454", 72.902, 416.78)
A1_lambda = TGaxis(100.0, 0.02, 600.0, 0.02, "f1", 1010, "NI")
A1_lambda.SetLabelFont(42)
A1_lambda.SetLabelSize(0.035)
A1_lambda.SetTextFont(42)
A1_lambda.SetTextSize(1.2)
A1_lambda.SetTitle("#Lambda [TeV]")
A1_lambda.SetTitleSize(0.04)
A1_lambda.SetTitleOffset(0.9)
A1_lambda.Draw()
myC2D.SaveAs(outputDir + "/limits" + "/limit_exclusion_region_2D_" +
year_to_plot + "_" + v1_tag + "_vs_" + v2_tag + ".pdf")
myC2D.SaveAs(outputDir + "/limits" + "/limit_exclusion_region_2D_" +
year_to_plot + "_" + v1_tag + "_vs_" + v2_tag + ".png")
myC2D.SaveAs(outputDir + "/limits" + "/limit_exclusion_region_2D_" +
year_to_plot + "_" + v1_tag + "_vs_" + v2_tag + ".C")
######CMS 8TeV, two photons
mass_cms_8TeV_2g = np.array([198., 227., 256., 256., 227., 198.])
ctau_cms_8TeV_2g = 0.01 * np.array([0.4, 2, 9, 9, 25., 50.])
graph_exclusion_cms_8TeV_2g = TGraph(6, mass_cms_8TeV_2g, ctau_cms_8TeV_2g)
graph_exclusion_cms_8TeV_2g.SetFillColorAlpha(kGray, 0.6)
graph_exclusion_cms_8TeV_2g.SetLineColorAlpha(kGray, 0.6)
####legend
#Lambda axis
f1_lambda = TF1("f1", "(x+6.00)/1.454", 72.902, 416.78)
A1_lambda = TGaxis(100.0, 0.0125, 600.0, 0.0125, "f1", 1010, "NI")
A1_lambda.SetLabelFont(42)
A1_lambda.SetLabelSize(0.035)
A1_lambda.SetTextFont(42)
A1_lambda.SetTextSize(1.3)
A1_lambda.SetTitle("#Lambda (TeV)")
A1_lambda.SetTitleSize(0.04)
A1_lambda.SetTitleOffset(0.9)
### only
graph_exclusion_exp_ = TGraph(
len(lambda_point_boundary_exp_),
np.array(1.454 * lambda_point_boundary_exp_ - 6.0),
0.01 * np.array(ctau_points))
graph_exclusion_exp_p1sig_ = TGraph(
len(lambda_point_boundary_exp_p1sig_),
np.array(1.454 * lambda_point_boundary_exp_p1sig_ - 6.0),
0.01 * np.array(ctau_points))
graph_exclusion_exp_m1sig_ = TGraph(
len(lambda_point_boundary_exp_m1sig_),
