applyFakeRateWeights="4L",
chargeSumSelections=chargeSumSelections,
central_or_shifts=central_or_shifts,
lep_mva_wp=lep_mva_wp,
jet_cleaning_by_index=jet_cleaning_by_index,
gen_matching_by_index=gen_matching_by_index,
max_files_per_job=files_per_job,
era=era,
use_lumi=True,
lumi=lumi,
check_output_files=check_output_files,
running_method=running_method,
num_parallel_jobs=num_parallel_jobs,
executable_addBackgrounds="addBackgrounds",
executable_addBackgroundJetToTauFakes="addBackgroundLeptonFakes",
histograms_to_fit=get_histograms_to_fit(),
select_rle_output=True,
dry_run=dry_run,
isDebug=debug,
use_nonnominal=use_nonnominal,
hlt_filter=hlt_filter,
use_home=use_home,
keep_logs=keep_logs,
blacklist=blacklist,
submission_cmd=sys.argv,
)
if mode.find("forBDTtraining") != -1:
analysis.set_BDT_training(hadTau_selection_relaxed)
job_statistics = analysis.create()
leptonChargeSelections = leptonChargeSelections,
central_or_shifts = central_or_shifts,
lep_mva_wp = lep_mva_wp,
jet_cleaning_by_index = jet_cleaning_by_index,
gen_matching_by_index = gen_matching_by_index,
max_files_per_job = files_per_job,
era = era,
use_lumi = True,
lumi = lumi,
invert_ZbosonMassVeto = True,
check_output_files = check_output_files,
running_method = running_method,
num_parallel_jobs = num_parallel_jobs,
executable_addBackgrounds = "addBackgrounds",
executable_addBackgroundJetToTauFakes = "addBackgroundLeptonFakes",
histograms_to_fit = get_histograms_to_fit("numJets", "dihiggsVisMass", "HT", "STMET"),
select_rle_output = True,
dry_run = dry_run,
isDebug = debug,
use_nonnominal = use_nonnominal,
hlt_filter = hlt_filter,
use_home = use_home,
keep_logs = keep_logs,
blacklist = blacklist,
submission_cmd = sys.argv,
)
if mode == "forBDTtraining":
analysis.set_BDT_training()
#analysis.set_BDT_training(hadTau_selection_relaxed)
applyFakeRateWeights="2lepton",
central_or_shifts=central_or_shifts,
lep_mva_wp=lep_mva_wp,
jet_cleaning_by_index=jet_cleaning_by_index,
gen_matching_by_index=gen_matching_by_index,
max_files_per_job=files_per_job,
era=era,
use_lumi=True,
lumi=lumi,
check_output_files=check_output_files,
running_method=running_method,
num_parallel_jobs=num_parallel_jobs,
executable_addBackgrounds="addBackgrounds",
executable_addFakes="addBackgroundLeptonFakes",
executable_addFlips="addBackgroundLeptonFlips",
histograms_to_fit=get_histograms_to_fit("dihiggsVisMass"),
select_rle_output=True,
dry_run=dry_run,
isDebug=debug,
use_nonnominal=use_nonnominal,
hlt_filter=hlt_filter,
use_home=use_home,
keep_logs=keep_logs,
blacklist=blacklist,
submission_cmd=sys.argv,
)
if mode == "forBDTtraining":
analysis.set_BDT_training()
job_statistics = analysis.create()
"mT_lep2_WZctrl_2lss",
"dR_ll_WZctrl_2lss",
"max_lep_eta_WZctrl_2lss",
"MVAOutput_allBMs",
]
for histogram_to_fit in histograms_to_fit_list:
histograms_to_fit[histogram_to_fit] = {}
else:
histograms_to_fit = get_histograms_to_fit(
"EventCounter",
"m3l",
"dihiggsVisMass_sel",
"mSFOS2l_closestToZ",
"dr_LeptonIdx3_AK4jNear_Approach2",
"dr_LeptonIdx3_2j_inclusive1j_Approach2",
"dr_los_min",
"dr_los_max",
"numSameFlavor_OS_3l",
"met_LD",
"numElectrons",
"numMuons",
"nJetAK4",
"nJetAK8_wSelectorAK8_Wjj"
)
analysis = analyzeConfig_hh_3l(
configDir = os.path.join("/scratch-persistent", getpass.getuser(), "hhAnalysis", era, version),
localDir = os.path.join("/home", getpass.getuser(), "hhAnalysis", era, version),
outputDir = os.path.join("/hdfs/local", getpass.getuser(), "hhAnalysis", era, version),
executable_analyze = "analyze_hh_3l",
cfgFile_analyze = "analyze_hh_3l_cfg.py",
samples = samples,
if __name__ == '__main__':
logging.info(
"Running the jobs with the following systematic uncertainties enabled: %s" % \
', '.join(central_or_shifts)
)
if not use_preselected:
logging.warning('Running the analysis on fully inclusive samples!')
if sample_filter:
samples = filter_samples(samples, sample_filter)
histograms_to_fit = get_histograms_to_fit(
"leptonPairMass", "dihiggsVisMass", "met_LD", "dR_ll", "dR_l_Wjets_min", "dR_l_leadWjet_min", "dR_l_Wjets_max", "dR_l_leadWjet_max", "dR_2j_fromW1", "mT_lep1_met",
"lep1_conePt", "lep2_conePt", "mindr_lep1_jet", "mindr_lep2_jet", "HT", "mht", "mass_2j_fromW1"
)
analysis = analyzeConfig_hh_2lss_leq1tau(
configDir = os.path.join("/scratch-persistent", getpass.getuser(), "hhAnalysis", era, version),
localDir = os.path.join("/home", getpass.getuser(), "hhAnalysis", era, version),
outputDir = os.path.join("/hdfs/local", getpass.getuser(), "hhAnalysis", era, version),
executable_analyze = "analyze_hh_2lss_leq1tau",
cfgFile_analyze = "analyze_hh_2lss_leq1tau_cfg.py",
samples = samples,
hadTauVeto_selection = hadTau_selection_veto,
applyFakeRateWeights = "2lepton",
central_or_shifts = central_or_shifts,
lep_mva_wp = lep_mva_wp,
jet_cleaning_by_index = jet_cleaning_by_index,
gen_matching_by_index = gen_matching_by_index,