len(param_points))
# Interpret coupling argument list
names = param_point[::2]
values = map(float, param_point[1::2])
modification_dict = {c: v for c, v in zip(names, values)}
# Let's not leave the user in the dark
logger.info("Model: %s", args.model)
logger.info("Process: %s", args.process)
logger.info("Couplings: %s",
", ".join(["%s=%5.4f" % c for c in modification_dict.items()]))
# make process
p = Process(process=args.process, nEvents=args.nEvents, config=config)
# Make grid pack
if args.makeGridpack:
gridpack = p.makeGridpack(modified_couplings=modification_dict,
overwrite=args.overwrite)
# calculate x-sec
if args.calcXSec:
xsec = p.xsec(modified_couplings=modification_dict,
overwrite=args.overwrite)
logger.info("xsec: %s ", repr(xsec))
if not args.keepWorkspace:
config.cleanup()
"color": ROOT.kGreen + 2
}
}
processes = ["ttZ", "ttH", "ttW"]
SM_xsec = {}
modified_couplings = {args.coupling: 0.0}
for proc in processes:
logger.info("Checking SM x-sec:")
config = Configuration(model_name=model_name,
modified_couplings=modified_couplings)
p = Process(process=proc, nEvents=50000, config=config)
SM_xsec[proc] = p.xsec()
logger.info("SM x-sec for %s is %s", proc, SM_xsec[proc])
if SM_xsec[proc].val == 0.: SM_xsec[proc] = u_float(1)
del config
hists = []
fits = []
m = 0
if args.scale:
scale = lambdaSqInv[args.coupling]
else:
scale = 1
for proc in processes:
cans = []
pads = []
SM_xsec = {}
latex1 = ROOT.TLatex()
latex1.SetNDC()
latex1.SetTextSize(0.04)
latex1.SetTextAlign(11)
logger.info("Model: %s", model_name)
for proc in processes[:1]:
logger.info("Starting with process %s", proc)
config = Configuration(model_name=model_name)
p = Process(process=proc, nEvents=50000, config=config)
SM_xsec[proc] = p.xsec(modified_couplings={})
logger.info("SM x-sec: %s", SM_xsec[proc])
for comb in combinations:
logger.info("Making 2D plots for %s", comb)
fixedPoints = list(nonZeroCouplings)
subDir = "{}_vs_{}".format(comb[0], comb[1])
if interpolate:
subDir += "_interpolated"
for c in comb:
fixedPoints.remove(c)
for v in points[fixedPoints[0]]:
for w in points[fixedPoints[1]]:
fileName = "{}={:.2f}_{}={:.2f}".format(