Commit 6bfee4f0 authored by Mike Hance's avatar Mike Hance
Browse files

some updates to be able to run sleptons with intermediate charginos

parent 5c5a8468
Pipeline #3140314 passed with stages
in 15 minutes and 16 seconds
This diff is collapsed.
......@@ -38,6 +38,7 @@ sleptonopts=""
stopopts=""
skip_pythia=false
kfactor=-1
chargino=-1
# some modifications based on run parameters
lumi=1
......@@ -63,7 +64,7 @@ fi
# get command line options
while getopts "E:M:P:p:N:m:x:s:e:c:GDAglaB:b:j:J:S:y:k:d:C:iL:f:F:X:h:I:nvTrK:" opt; do
while getopts "E:M:P:p:N:m:x:s:e:c:GDAglaB:b:j:J:S:y:k:d:C:iL:f:F:X:h:I:nvTrK:O:" opt; do
case "${opt}" in
E) ecms=$OPTARG;;
M) mass=$OPTARG;;
......@@ -108,6 +109,8 @@ while getopts "E:M:P:p:N:m:x:s:e:c:GDAglaB:b:j:J:S:y:k:d:C:iL:f:F:X:h:I:nvTrK:"
deltaeta=0
pythia_card="pythia8_card.dat"
stopopts="-r";;
O) chargino=$OPTARG;;
*) exit;;
esac
done
......@@ -115,7 +118,12 @@ done
# construct the tag.
tag="VBFSUSY_${ecms}_${params}_${mass}_mmjj_${mmjj}_${mmjjmax}${suffix}"
if [[ $mmjj == 0.0 || $mmjj == 0 ]]; then
tag="SUSY_${ecms}_${params}_${mass}_${dM}_${proc}_${suffix}"
if [[ $proc == "isrslep" && $chargino != -1 ]]; then
tag="SUSY_${ecms}_${params}_${mass}_${dM}_${chargino}_${proc}_${suffix}"
sleptonopts="${sleptonopts} -O ${chargino}"
else
tag="SUSY_${ecms}_${params}_${mass}_${dM}_${proc}_${suffix}"
fi
fi
pythia_onoff=""
......
......@@ -28,7 +28,7 @@ masssplitting=20
proc="isrslep"
params="SleptonBino"
while getopts "E:M:S:N:c:d:f:P:p:J:L:F:s:glab:" opt; do
while getopts "E:M:S:N:c:d:f:P:p:J:L:F:s:glab:Z" opt; do
case "${opt}" in
E) ecms=$OPTARG;;
M) mass=$OPTARG;;
......@@ -47,6 +47,7 @@ while getopts "E:M:S:N:c:d:f:P:p:J:L:F:s:glab:" opt; do
f) simpleanalysis=$OPTARG;;
F) likelihood=$OPTARG;;
b) database=$OPTARG;;
Z) withchargino=true;;
*) exit;;
esac
done
......@@ -60,6 +61,12 @@ elif $clobber_ana; then
clobberopts="-a"
fi
newmodelopts=""
if $withchargino; then
charginomass=$(bc <<< "scale=0; ${mass}-${masssplitting}/2")
newmodelopts="-v -O $charginomass"
fi
./run_VBFSUSY_standalone.sh \
-E ${ecms} \
-K ${kfactor} \
......@@ -83,4 +90,5 @@ fi
-n \
-b ${database} \
-I "-2.9" \
${newmodelopts} \
${clobberopts}
......@@ -32,8 +32,9 @@ slepton=false
stop=false
sneutrino=false
skip_pythia=false
chargino=-1
while getopts "E:M:P:p:N:m:x:e:c:GgB:b:S:y:k:sd:j:J:X:I:vTr" opt; do
while getopts "E:M:P:p:N:m:x:e:c:GgB:b:S:y:k:sd:j:J:X:I:vTrO:" opt; do
case "${opt}" in
E) ecms=$OPTARG;;
M) mass=$OPTARG;;
......@@ -59,6 +60,7 @@ while getopts "E:M:P:p:N:m:x:e:c:GgB:b:S:y:k:sd:j:J:X:I:vTr" opt; do
d) seed=$OPTARG;;
I) MGversion=$OPTARG;;
T) skip_pythia=true;;
O) chargino=$OPTARG;;
*) exit;;
esac
done
......@@ -76,6 +78,9 @@ if [[ ${slepton} == true ]]; then
mSNU=$(bc <<< "scale=2; ${mass}-${dM}/2")
massopts="${massopts} -m MSNU ${mSNU}"
fi
if [[ $chargino != -1 ]]; then
massopts="${massopts} -m MC1 ${chargino} -m MN2 ${chargino}"
fi
elif [[ ${stop} == true ]]; then
massopts="-m MSTOP ${mass}"
elif [[ ${params} == WinoBino ]]; then
......@@ -115,7 +120,6 @@ fi
-t ${tag} \
${pythia_onoff} ${clobber_opts}
if [[ $? == 0 || ${clobber_mgpy} == true ]]; then
docker run \
--log-driver=journald \
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment