benchmarks/b0_tracker/forward_protons.sh

0001 #!/bin/bash
0002
0003 if [[ ! -n  "${DETECTOR}" ]] ; then
0004   export DETECTOR="topside"
0005 fi
0006
0007 if [[ ! -n  "${JUGGLER_N_EVENTS}" ]] ; then
0008   export JUGGLER_N_EVENTS=100
0009 fi
0010
0011 export FILE_NAME_TAG="forward_protons"
0012 export JUGGLER_GEN_FILE="${LOCAL_DATA_PATH}/${FILE_NAME_TAG}.hepmc"
0013 export JUGGLER_SIM_FILE="${LOCAL_DATA_PATH}/sim_${FILE_NAME_TAG}.edm4hep.root"
0014
0015 echo "JUGGLER_N_EVENTS = ${JUGGLER_N_EVENTS}"
0016 echo "DETECTOR = ${DETECTOR}"
0017
0018 # Generate the input events
0019 root -b -q "benchmarks/b0_tracker/scripts/gen_${FILE_NAME_TAG}.cxx(${JUGGLER_N_EVENTS}, \"${FILE_NAME_TAG}.hepmc\")"
0020 if [[ "$?" -ne "0" ]] ; then
0021   echo "ERROR running script: generating input events"
0022   exit 1
0023 fi
0024
0025 # Run geant4 simulations
0026 ddsim --runType batch \
0027       -v WARNING \
0028       --part.minimalKineticEnergy 0.5*GeV  \
0029       --filter.tracker edep0 \
0030       --numberOfEvents ${JUGGLER_N_EVENTS} \
0031       --compactFile ${DETECTOR_PATH}/${DETECTOR_CONFIG}.xml \
0032       --inputFiles ${FILE_NAME_TAG}.hepmc \
0033       --outputFile ${JUGGLER_SIM_FILE}
0034
0035 if [[ "$?" -ne "0" ]] ; then
0036   echo "ERROR running npdet"
0037   exit 1
0038 fi
0039
0040 # Directory for plots
0041 mkdir -p results/far_forward/b0
0042
0043 rootls -t ${JUGGLER_SIM_FILE}
0044 # Plot the input events
0045 root -b -q "benchmarks/b0_tracker/analysis/b0_tracker_hits.cxx+(\"${JUGGLER_SIM_FILE}\")"
0046 if [[ "$?" -ne "0" ]] ; then
0047   echo "ERROR running script: events"
0048   exit 1
0049 fi
0050