Variables | |
| tuple | log = logging.getLogger(__name__) |
| tuple | parser = argparse.ArgumentParser(description="Plot significance from combine output.") |
| list | channels = ["et", "mt", "tt", "em"] |
| tuple | args = parser.parse_args() |
| tuple | fig = plt.figure() |
| tuple | ax = fig.add_subplot(111) |
| list | plots = [] |
| mass = args.mass | |
| tuple | x_range = ([float(a) if "." in a else int(a) for a in args.x_range.split(",")]) |
| tuple | p1s = np.linspace(*x_range) |
| tuple | p1_width = (max(p1s) - min(p1s)) |
| tuple | y_range = ([float(a) if "." in a else int(a) for a in args.y_range.split(",")]) |
| tuple | p2s = np.linspace(*y_range) |
| tuple | p2_width = (max(p2s) - min(p2s)) |
| list | matrix = [[0 for x in range(len(p1s))] for y in range(len(p2s))] |
| tuple | in_raw = os.path.join(args.input_dir, str(p1)+"_"+str(p2), args.folder, str(mass), "higgsCombineTest.ProfileLikelihood.mH125.*.root") |
| tuple | in_file = glob.glob( in_raw) |
| tuple | current_chain = ROOT.TChain("chain_%i_%f"%(p1, p2)) |
| y = line.limit | |
| tuple | npmatrix = np.matrix(matrix) |
| tuple | min_value = np.min(npmatrix[np.nonzero(npmatrix)]) |
| tuple | im = plt.imshow(matrix, origin = "lower", cmap=plt.cm.Reds, interpolation="None", extent=[p1s[0]-p1_width/2, p1s[-1]+p1_width/2,p2s[0]-p2_width/2, p2s[-1]+p2_width/2], clim = (min_value,npmatrix.max())) |
| dictionary | labels = { "mjj" : "$m_{jj}$", "jdeta" : "$\Delta \eta_{jj}$", "hpt" : "$p_T^{H}$", "pt2" : "$p_T^2$"} |
| tuple significance_2d.args = parser.parse_args() |
| tuple significance_2d.ax = fig.add_subplot(111) |
| list significance_2d.channels = ["et", "mt", "tt", "em"] |
| tuple significance_2d.current_chain = ROOT.TChain("chain_%i_%f"%(p1, p2)) |
| tuple significance_2d.fig = plt.figure() |
| tuple significance_2d.im = plt.imshow(matrix, origin = "lower", cmap=plt.cm.Reds, interpolation="None", extent=[p1s[0]-p1_width/2, p1s[-1]+p1_width/2,p2s[0]-p2_width/2, p2s[-1]+p2_width/2], clim = (min_value,npmatrix.max())) |
| tuple significance_2d.in_file = glob.glob( in_raw) |
| tuple significance_2d.in_raw = os.path.join(args.input_dir, str(p1)+"_"+str(p2), args.folder, str(mass), "higgsCombineTest.ProfileLikelihood.mH125.*.root") |
| dictionary significance_2d.labels = { "mjj" : "$m_{jj}$", "jdeta" : "$\Delta \eta_{jj}$", "hpt" : "$p_T^{H}$", "pt2" : "$p_T^2$"} |
| tuple significance_2d.log = logging.getLogger(__name__) |
| significance_2d.mass = args.mass |
| tuple significance_2d.npmatrix = np.matrix(matrix) |
| tuple significance_2d.p1s = np.linspace(*x_range) |
| tuple significance_2d.p2s = np.linspace(*y_range) |
| tuple significance_2d.parser = argparse.ArgumentParser(description="Plot significance from combine output.") |
| list significance_2d.plots = [] |
| tuple significance_2d.x_range = ([float(a) if "." in a else int(a) for a in args.x_range.split(",")]) |
| significance_2d.y = line.limit |
| tuple significance_2d.y_range = ([float(a) if "." in a else int(a) for a in args.y_range.split(",")]) |