plot_configs = []
for datacard in args.datacards:
cb = ch.CombineHarvester()
cb.SetFlag("workspaces-use-clone", True)
for template in filename_templates:
template_tag = template.split("$")[0]
if template_tag in datacard:
matched_template = os.path.join(datacard[:datacard.index(template_tag)], template).replace("{", "").replace("}", "")
print datacard
print matched_template
cb.QuickParseDatacard(datacard, matched_template)
break
datacards = cpstudiesdatacards.CPStudiesDatacards(cb=cb)
for model in args.models:
model_settings = models.get(model, {})
datacards_workspaces = {}
for freeze_syst_uncs in args.freeze_syst_uncs:
output_dir_base = args.output_dir
if output_dir_base is None:
output_dir_base = os.path.splitext(datacard)[0]
sub_dir_base = os.path.join("projection", model, "statUnc" if freeze_syst_uncs else "totUnc")
output_dir_base = os.path.abspath(os.path.expandvars(os.path.join(output_dir_base, sub_dir_base)))
if args.clear_output_dir and os.path.exists(output_dir_base):
logger.subprocessCall("rm -r " + output_dir_base, shell=True)
cp_mixings_scan = list(numpy.arange(args.cp_mixings[0], args.cp_mixings[-1]+0.001, (args.cp_mixings[-1]-args.cp_mixings[0])/(args.cp_mixing_scan_points-1)))
cp_mixings_combine_range_min = (3*cp_mixings_scan[0]-cp_mixings_scan[1]) / 2.0
cp_mixings_combine_range_max = (3*cp_mixings_scan[-1]-cp_mixings_scan[-2]) / 2.0
# initialisations for plotting
sample_settings = samples.Samples()
binnings_settings = binnings.BinningsDict()
systematics_factory = systematics.SystematicsFactory()
plot_configs = []
output_files = []
merged_output_files = []
hadd_commands = []
datacards = cpstudiesdatacards.CPStudiesDatacards(cp_mixings_str, add_data=True) # add_data=args.add_data)
# initialise datacards
tmp_input_root_filename_template = "input/${ANALYSIS}_${CHANNEL}_${BIN}_${SYSTEMATIC}_${ERA}.root"
input_root_filename_template = "input/${ANALYSIS}_${CHANNEL}_${BIN}_${ERA}.root"
bkg_histogram_name_template = "${BIN}/${PROCESS}"
sig_histogram_name_template = "${BIN}/${PROCESS}${MASS}"
bkg_syst_histogram_name_template = "${BIN}/${PROCESS}_${SYSTEMATIC}"
sig_syst_histogram_name_template = "${BIN}/${PROCESS}${MASS}_${SYSTEMATIC}"
datacard_filename_templates = [
"datacards/individual/${CHANNEL}/${BINID}/${ANALYSIS}_${CHANNEL}_${BINID}_${ERA}.txt",
"datacards/channel/${CHANNEL}/${ANALYSIS}_${CHANNEL}_${ERA}.txt",
"datacards/category/${BINID}/${ANALYSIS}_${BINID}_${ERA}.txt",
"datacards/combined/${ANALYSIS}_${ERA}.txt",
]
output_root_filename_template = "datacards/common/${ANALYSIS}.input_${ERA}.root"