Warning
UNDER CONSTRUCTION
The skimming is performed by the hhskim.py script.
Run the skims on the grid (after setting up the panda client and your VOMS proxy with the phys-higgs production role):
./skim --yall mc11_hadhad mc12_hadhad data11_hadhad data12_hadhad embed11_hadhad embed12_hadhadAfter the skims are finished and downloaded, update the paths in higgstautau/datasets_config.yml and update the datasets database:
./dsdb --reset hhThen launch the batch jobs that create all the analysis ntuples (nominal and systematics) with:
./run_all_hh.shAutomatically organize ROOT and log files with:
init-ntupThe above also merges the output from the subjobs for embedding and data.
Add the Higgs pt weights:
make higgs-ptThis creates copies of all signal ROOT files ntuples/running/hhskim/hhskim*tautauhh*.root at ntuples/running/hhskim/weighted.hhskim*tautauhh*.root.
Make a backup of the original signal files:
mkdir ntuples/running/hhskim/higgs_unweighted
mv ntuples/running/hhskim/hhskim*tautauhh*.root ntuples/running/hhskim/higgs_unweightedThen remove weighted. from the weighted signal files:
rename weighted. "" ntuples/running/hhskim/weighted.hhskim*tautauhh*.rootThen move the hhskim running directory to production (replace XX with the skim version number):
mkdir ntuples/prod_vXX
mv ntuples/running/hhskim ntuples/prod_vXXUpdate the production path in the Makefile (HHNTUP) and mva/__init__.py (NTUPLE_PATH).
And finally create the merged hhskim.root and hhskim.h5:
make ntupGenerate the cache of all background normalizations:
make normsCreate all validation plots with:
make plotsRun the batch jobs that train the BDTs at each mass point with:
make trainGo get some coffee.
Create all the BDT validation plots with:
make mva-plotsRun the batch jobs to determine the optimal binning for each mass point in each category and year:
make binningGo get some coffee.
Run the batch jobs that create the workspaces with:
make mva-workspaces
make cuts-workspacesGo get some coffee.
Move to the directory containing all the MVA workspaces to combine:
cd workspaces/hh_nos_nonisol_ebz_mvaCombine workspaces across years with:
for mass in $(seq 100 5 150); do
combine hh_11_vbf_$mass hh_12_vbf_$mass --name hh_vbf_$mass;
combine hh_11_boosted_$mass hh_12_boosted_$mass --name hh_boosted_$mass;
combine hh_11_combination_$mass hh_12_combination_$mass --name hh_combination_$mass;
doneCheck your email.
Move to the directory containing all the CBA workspaces to combine:
cd workspaces/hh_nos_nonisol_ebz_cutsCreate VBF and boosted combinations for each year, and a combination across years:
for mass in $(seq 100 5 150); do
combine hh_11_cuts_boosted_loose_$mass hh_11_cuts_boosted_tight_$mass --name hh_11_cuts_boosted_$mass;
combine hh_12_cuts_boosted_loose_$mass hh_12_cuts_boosted_tight_$mass --name hh_12_cuts_boosted_$mass;
combine hh_11_cuts_vbf_lowdr_$mass hh_11_cuts_vbf_highdr_$mass --name hh_11_cuts_vbf_$mass;
combine hh_12_cuts_vbf_lowdr_$mass hh_12_cuts_vbf_highdr_loose_$mass hh_12_cuts_vbf_highdr_tight_$mass --name hh_12_cuts_vbf_$mass;
combine hh_11_cuts_boosted_$mass hh_12_cuts_boosted_$mass --name hh_cuts_boosted_$mass;
combine hh_11_cuts_vbf_$mass hh_12_cuts_vbf_$mass --name hh_cuts_vbf_$mass;
combine hh_11_combination_$mass hh_12_combination_$mass --name hh_combination_$mass;
doneApply all of the HSG4 workspace fixes with:
cd workspaces
fix-workspace --quiet --symmetrize --prune-shapes --chi2-thresh 0.9 hh_nos_nonisol_ebz_mva
fix-workspace --quiet --symmetrize --prune-shapes --chi2-thresh 0.9 hh_nos_nonisol_ebz_cutsConstruct the profile of every nuisance parameter (NP):
# submit a batch job for each NP. If --submit is omitted simply print the command.
multinp scans_fit --file path_to_measurement_file.root --submit
# merge all the output in a single file and compute the nominal NLL for normalisation
multinp merge --file path_to_measurement_file.root --jobs -1
# Clean the directory from the individual pickle files (keep only the master)
multinp clean --file path_to_measurement_file.rootUpdate the paths in plot-nuis and plot the profiles with:
plot-nuisCompute the pull of each nuisance parameter with:
multinp pulls --file path_to_measurement_file.root --jobs -1Update the path in plot-ranking and plot the ranking with:
plot-rankingCompute the expected significance (bkg. only hypothesis) with:
# Walk trough the directory and subdirectory and look for workspaces
multisig path_to_directory_containing_workspacesCompute the postfit histograms and errors with:
# --fit_var bdt_score/mmc_mass
plot-postfit path_to_measurement_file.root --fit-var bdt_score --force-fit --jobs -1
# If the fit has already been performed
plot-postfit path_to_measurement_file.root --fit-var bdt_score