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="topside"
0005 fi
0006 
0007 if [[ ! -n  "${JUGGLER_N_EVENTS}" ]] ; then 
0008   export JUGGLER_N_EVENTS=1000
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 JUGGLER_FILE_NAME_TAG="emcal_barrel_uniform_pions"
0020 export JUGGLER_GEN_FILE="${JUGGLER_FILE_NAME_TAG}.hepmc"
0021 
0022 export JUGGLER_SIM_FILE="sim_${JUGGLER_FILE_NAME_TAG}.edm4hep.root"
0023 export JUGGLER_REC_FILE="rec_${JUGGLER_FILE_NAME_TAG}.root"
0024 
0025 echo "JUGGLER_N_EVENTS = ${JUGGLER_N_EVENTS}"
0026 echo "DETECTOR = ${DETECTOR}"
0027 
0028 # Generate the input events
0029 root -b -q "benchmarks/barrel_ecal/scripts/emcal_barrel_pions.cxx(${JUGGLER_N_EVENTS}, ${E_start}, ${E_end}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
0030 if [[ "$?" -ne "0" ]] ; then
0031   echo "ERROR running script: generating input events"
0032   exit 1
0033 fi
0034 # Plot the input events
0035 root -b -q "benchmarks/barrel_ecal/scripts/emcal_barrel_pions_reader.cxx(${E_start}, ${E_end}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
0036 if [[ "$?" -ne "0" ]] ; then
0037   echo "ERROR running script: plotting input events"
0038   exit 1
0039 fi
0040 
0041 # Run geant4 simulations
0042 ddsim --runType batch \
0043       -v WARNING \
0044       --part.minimalKineticEnergy 0.5*GeV  \
0045       --filter.tracker edep0 \
0046       --numberOfEvents ${JUGGLER_N_EVENTS} \
0047       --compactFile ${DETECTOR_PATH}/${DETECTOR_CONFIG}.xml \
0048       --inputFiles ${JUGGLER_FILE_NAME_TAG}.hepmc \
0049       --outputFile sim_output/${JUGGLER_SIM_FILE}
0050 
0051 if [[ "$?" -ne "0" ]] ; then
0052   echo "ERROR running npdet"
0053   exit 1
0054 fi
0055 
0056 # Directory for plots
0057 mkdir -p results
0058 
0059 # Move ROOT output file
0060 #mv ${JUGGLER_REC_FILE} sim_output/
0061