macros/wEEC_3D_unfolding_kinematics/run_wEEC_unfolding.sh

0001 #!/bin/bash
0002 # ─────────────────────────────────────────────────────────────────────────────
0003 #  run_wEEC_unfolding.sh  —  runs wEEC_doUnfolding.C
0004 #
0005 #  Usage:  ./run_wEEC_unfolding.sh <DIR> <MODE> [NITER_WEEC]
0006 #    DIR        — output directory
0007 #    MODE       — kFull | kHalf | kData
0008 #    NITER_WEEC — Bayes iterations for wEEC unfolding (default: 4)
0009 #
0010 #  Can be run independently of run_unfolding.sh — does not need jet_pT file.
0011 #  For kData: uses fullClosure response for MC corrections,
0012 #             data_measured-all.root for measured wEEC.
0013 # ─────────────────────────────────────────────────────────────────────────────
0014
0015 export USER="$(id -u -n)"
0016 export LOGNAME=${USER}
0017 export HOME=/sphenix/u/${LOGNAME}/macros/detectors/sPHENIX/
0018 export MYINSTALL=/sphenix/user/bkimelman/sPHENIX/install/
0019
0020 source /opt/sphenix/core/bin/sphenix_setup.sh -n
0021 source /opt/sphenix/core/bin/setup_local.sh $MYINSTALL
0022
0023 DIR=/sphenix/tg/tg01/jets/bkimelman/VandyDSTs_wEEC_3D_unfolding_kinematics_Apr29_2026/
0024 MODE=$1
0025 NITER=${2:-4}
0026
0027 if [ -z "$DIR" ] || [ -z "$MODE" ]; then
0028     echo "Usage: $0 <DIR> <MODE> [NITER_WEEC]"
0029     exit 1
0030 fi
0031
0032 case $MODE in
0033     kFull) LABEL="fullClosure" ;;
0034     kHalf) LABEL="halfClosure" ;;
0035     kData) LABEL="dataClosure" ;;
0036     *)
0037         echo "Unknown mode: $MODE  (expected kFull, kHalf, or kData)"
0038         exit 1
0039         ;;
0040 esac
0041
0042 OUT_FILE="${DIR}/wEEC-${LABEL}.root"
0043
0044 if [ "$MODE" = "kData" ]; then
0045     RESP_FILE="${DIR}/response-all-${LABEL}.root"
0046     MEAS_FILE="${DIR}/data_measured-all.root"
0047     echo "Running wEEC_doUnfolding: mode=${MODE}, resp=${RESP_FILE}, meas=${MEAS_FILE}, out=${OUT_FILE}"
0048     root -b -q "wEEC_doUnfolding.C(\"${RESP_FILE}\",\"${OUT_FILE}\",${NITER},Mode::${MODE},\"${MEAS_FILE}\")"
0049 else
0050     RESP_FILE="${DIR}/response-all-${LABEL}.root"
0051     echo "Running wEEC_doUnfolding: mode=${MODE}, resp=${RESP_FILE}, out=${OUT_FILE}"
0052     root -b -q "wEEC_doUnfolding.C(\"${RESP_FILE}\",\"${OUT_FILE}\",${NITER},Mode::${MODE})"
0053 fi
0054
0055 echo "Done: ${OUT_FILE}"