Skip to content

Commit

Permalink
bf
Browse files Browse the repository at this point in the history
  • Loading branch information
rkansal47 committed Mar 7, 2024
1 parent f61a5ba commit a05d3a4
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 20 deletions.
4 changes: 2 additions & 2 deletions src/HHbbVV/postprocessing/bash_scripts/MassPlots.sh
Original file line number Diff line number Diff line change
Expand Up @@ -15,6 +15,7 @@
MAIN_DIR="../../.."
TAG=""
samples="HHbbVV VBFHHbbVV NMSSM_XToYHTo2W2BTo4Q2B_MX-900_MY-80 NMSSM_XToYHTo2W2BTo4Q2B_MX-1200_MY-190 NMSSM_XToYHTo2W2BTo4Q2B_MX-2000_MY-125 NMSSM_XToYHTo2W2BTo4Q2B_MX-3000_MY-250 NMSSM_XToYHTo2W2BTo4Q2B_MX-4000_MY-150"
# samples="HHbbVV VBFHHbbVV NMSSM_XToYHTo2W2BTo4Q2B_MX-900_MY-80"

options=$(getopt -o "" --long "tag:" -- "$@")
eval set -- "$options"
Expand Down Expand Up @@ -45,8 +46,7 @@ if [[ -z $TAG ]]; then
exit 1
fi

# for year in 2016APV 2016 2017 2018
for year in 2016
for year in 2016APV 2016 2017 2018
do
python -u postprocessing.py --control-plots --year $year --resonant \
--data-dir "${MAIN_DIR}/../data/skimmer/24Mar5AllYears" \
Expand Down
19 changes: 1 addition & 18 deletions src/HHbbVV/postprocessing/postprocessing.py
Original file line number Diff line number Diff line change
Expand Up @@ -504,9 +504,6 @@ def main(args):
args.signal_data_dirs[0],
)

# Check for 0 weights - would be an issue for weight shifts
check_weights(events_dict)

print("\nMaking templates")
templates = {}

Expand Down Expand Up @@ -1459,7 +1456,7 @@ def control_plots(
events_dict, shape_var, bb_masks, weight_key=weight_key, selection=selection
)

ylim = np.max([h.values() for h in hists.values()]) if same_ylim else None
ylim = np.max([h.values() for h in hists.values()]) * 1.05 if same_ylim else None

if HEM2d and year == "2018":
hists["HEM2d"] = hists_HEM2d(events_dict, bb_masks, weight_key, selection)
Expand Down Expand Up @@ -1553,20 +1550,6 @@ def plot_bdt_sculpting(
)


def check_weights(events_dict):
# Check for 0 weights - would be an issue for weight shifts
print(
"\nAny 0 weights:",
np.any(
[
np.any(events["weight_nonorm"] == 0)
for key, events in events_dict.items()
if key != data_key
]
),
)


def _get_fill_data(
events: pd.DataFrame, bb_mask: pd.DataFrame, shape_vars: list[ShapeVar], jshift: str = ""
):
Expand Down

0 comments on commit a05d3a4

Please sign in to comment.