LightFlavorRatios/swimming_correction/run_cutefficiency_systematics.sh

0001 #!/bin/bash
0002
0003 source /opt/sphenix/core/bin/sphenix_setup.sh -n new
0004
0005 export USER="$(id -u -n)"
0006 export LOGNAME=${USER}
0007 export HOME=/sphenix/u/${LOGNAME}
0008
0009 export SPHENIX=${HOME}/sPHENIX
0010 export MYINSTALL=$SPHENIX/install
0011 export LD_LIBRARY_PATH=$MYINSTALL/lib:$LD_LIBRARY_PATH
0012 export ROOT_INCLUDE_PATH=$MYINSTALL/include:$ROOT_INCLUDE_PATH
0013
0014 mkdir -p /tmp/log_mjpeters/
0015
0016 PROCESS=$1
0017
0018 #LAMBDA_INFILE_BASE=/sphenix/tg/tg01/hf/mjpeters/lambdaKshortMB/lambdaKshort_20260422_DetroitMB_CR_2_mode_pTref_1p4/ppi_reco/outputKFParticle_ppi_reco
0019 #KSHORT_INFILE_BASE=/sphenix/tg/tg01/hf/mjpeters/lambdaKshortMB/lambdaKshort_20260422_DetroitMB_CR_2_mode_pTref_1p4/pipi_reco/outputKFParticle_pipi_reco
0020
0021 NOMINAL_CUTEFFICIENCY_INFILE=/sphenix/tg/tg01/hf/mjpeters/LightFlavorProduction/cutEfficiencyCorrection/merged_cut_efficiency.root
0022
0023 LAMBDA_INFILE_BASE=/sphenix/tg/tg01/hf/mjpeters/LightFlavorProduction/closureTestSample/ppi_reco/outputKFParticle_ppi_reco
0024 KSHORT_INFILE_BASE=/sphenix/tg/tg01/hf/mjpeters/LightFlavorProduction/closureTestSample/pipi_reco/outputKFParticle_pipi_reco
0025
0026 OUTFILE_MIN_BASE=/sphenix/tg/tg01/hf/mjpeters/LightFlavorProduction/cutEfficiencyCorrection/systematics_min/cut_efficiency
0027 OUTFILE_MAX_BASE=/sphenix/tg/tg01/hf/mjpeters/LightFlavorProduction/cutEfficiencyCorrection/systematics_max/cut_efficiency
0028
0029 LAMBDA_PDGID=3122
0030 LAMBDA_DAUGHTERS={-211,2212}
0031 LAMBDABAR_DAUGHTERS={211,-2212}
0032
0033 KSHORT_PDGID=310
0034 KSHORT_DAUGHTERS={211,-211}
0035
0036 for i in {0..9}; do
0037         SEGMENT=`printf "%06d" $(( 10*PROCESS + i ))`
0038
0039         rm ${OUTFILE_MIN_BASE}_${SEGMENT}.root
0040         rm ${OUTFILE_MAX_BASE}_${SEGMENT}.root
0041         rm ${OUTFILE_MIN_BASE}_pos_${SEGMENT}.root
0042         rm ${OUTFILE_MAX_BASE}_pos_${SEGMENT}.root
0043         rm ${OUTFILE_MIN_BASE}_neg_${SEGMENT}.root
0044         rm ${OUTFILE_MAX_BASE}_neg_${SEGMENT}.root
0045
0046         # (lambda + lambdabar)/2Kshort
0047         root -l -q -b CutEfficiency_systematics.C\(\"${NOMINAL_CUTEFFICIENCY_INFILE}\",\"Lambda0\",${LAMBDA_PDGID},${LAMBDA_DAUGHTERS},\"${LAMBDA_INFILE_BASE}_${SEGMENT}.root\",true,StandardCuts::MC_Lambda0_cuts,BinInfo::Lambda0_MC_mass_bins,\"K_S0\",${KSHORT_PDGID},${KSHORT_DAUGHTERS},\"${KSHORT_INFILE_BASE}_${SEGMENT}.root\",false,StandardCuts::MC_K_S0_cuts,BinInfo::K_S0_MC_mass_bins,\"${OUTFILE_MIN_BASE}_${SEGMENT}.root\",\"${OUTFILE_MAX_BASE}_${SEGMENT}.root\"\)
0048         # lambda/Kshort
0049         root -l -q -b CutEfficiency_systematics.C\(\"${NOMINAL_CUTEFFICIENCY_INFILE}\",\"Lambda0\",${LAMBDA_PDGID},${LAMBDA_DAUGHTERS},\"${LAMBDA_INFILE_BASE}_${SEGMENT}.root\",false,StandardCuts::MC_Lambda0_pos_cuts,BinInfo::Lambda0_MC_mass_bins,\"K_S0\",${KSHORT_PDGID},${KSHORT_DAUGHTERS},\"${KSHORT_INFILE_BASE}_${SEGMENT}.root\",false,StandardCuts::MC_K_S0_cuts,BinInfo::K_S0_MC_mass_bins,\"${OUTFILE_MIN_BASE}_pos_${SEGMENT}.root\",\"${OUTFILE_MAX_BASE}_pos_${SEGMENT}.root\"\)
0050         # lambdabar/kshort
0051         root -l -q -b CutEfficiency_systematics.C\(\"${NOMINAL_CUTEFFICIENCY_INFILE}\",\"Lambda0\",-${LAMBDA_PDGID},${LAMBDABAR_DAUGHTERS},\"${LAMBDA_INFILE_BASE}_${SEGMENT}.root\",false,StandardCuts::MC_Lambda0_neg_cuts,BinInfo::Lambda0_MC_mass_bins,\"K_S0\",${KSHORT_PDGID},${KSHORT_DAUGHTERS},\"${KSHORT_INFILE_BASE}_${SEGMENT}.root\",false,StandardCuts::MC_K_S0_cuts,BinInfo::K_S0_MC_mass_bins,\"${OUTFILE_MIN_BASE}_neg_${SEGMENT}.root\",\"${OUTFILE_MAX_BASE}_neg_${SEGMENT}.root\"\)
0052 done