benchmarks/barrel_hcal/run_hcal_barrel_particles.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 if [[ ! -n  "${E_START}" ]] ; then
0012   export E_START=5.0
0013 fi
0014
0015 if [[ ! -n  "${E_END}" ]] ; then
0016   export E_END=5.0
0017 fi
0018
0019 export PARTICLE=$1
0020 if [[ ! -n  "${PARTICLE}" ]] ; then
0021   export PARTICLE="electron"
0022 fi
0023
0024 if [[ ! -n  "${JUGGLER_FILE_NAME_TAG}" ]] ; then
0025   export JUGGLER_FILE_NAME_TAG="hcal_barrel_${PARTICLE}"
0026 fi
0027
0028 export JUGGLER_GEN_FILE="data/${JUGGLER_FILE_NAME_TAG}.hepmc"
0029
0030 export JUGGLER_SIM_FILE="sim_output/sim_${JUGGLER_FILE_NAME_TAG}.edm4hep.root"
0031 export JUGGLER_REC_FILE="rec_${JUGGLER_FILE_NAME_TAG}.root"
0032
0033 echo "JUGGLER_N_EVENTS = ${JUGGLER_N_EVENTS}"
0034 echo "DETECTOR = ${DETECTOR}"
0035
0036 # Generate the input events
0037 root -b -q "benchmarks/barrel_hcal/scripts/hcal_barrel_particles_gen.cxx+(${JUGGLER_N_EVENTS}, ${E_START}, ${E_END}, \"${PARTICLE}\")"
0038 if [[ "$?" -ne "0" ]] ; then
0039   echo "ERROR running script: generating input events"
0040   exit 1
0041 fi
0042 # Plot the input events
0043 root -b -q "benchmarks/barrel_hcal/scripts/hcal_barrel_particles_reader.cxx+(\"${PARTICLE}\")"
0044 if [[ "$?" -ne "0" ]] ; then
0045   echo "ERROR running script: plotting input events"
0046   exit 1
0047 fi
0048
0049 ls -ltRhL
0050
0051 ddsim --runType batch \
0052       -v WARNING \
0053       --part.minimalKineticEnergy 0.5*GeV  \
0054       --filter.tracker edep0 \
0055       --numberOfEvents ${JUGGLER_N_EVENTS} \
0056       --compactFile ${DETECTOR_PATH}/${DETECTOR_CONFIG}.xml \
0057       --inputFiles ${JUGGLER_GEN_FILE} \
0058       --outputFile ${JUGGLER_SIM_FILE}
0059
0060 if [[ "$?" -ne "0" ]] ; then
0061   echo "ERROR running npdet"
0062   exit 1
0063 fi
0064
0065 # Directory for plots
0066 mkdir -p results