Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2024-09-27 07:02:37

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