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(