8# ------------------------------------------------------------------------------------------
10# Utility script to determine dE/dx due to visible delta-rays
12# ------------------------------------------------------------------------------------------
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
21source $JPP_DIR/setenv.sh $JPP_DIR
25set_variable WORKDIR ./
26set_variable TMIN_GEV 0.000915499
28set_variable+ NUMERICAL DELTARAYS_NUMERICAL -N
30set_variable: FORMAT GRAPHICS_FORMAT gif
31set_variable+ BATCH GRAPHICS_BATCH -B
34 usage "$script [working directory]"
38 1) set_variable WORKDIR $1;;
44 for LEPTON in muon tau positron electron; do
46 $JPP_DIR/examples/JPhysics/JDeltaRays \
47 -o $WORKDIR/delta-rays\[${LEPTON}\].root \
51 ${NUMERICAL} -d $DEBUG
54 -f $WORKDIR/delta-rays\[${LEPTON}\].root:h0 \
58 -\^ "dE/dx [MeV g^{-1} cm^{2}]" \
62 -o $LEPTON.$FORMAT $BATCH
69 -f $WORKDIR/delta-rays\[muon\].root:h0 \
70 -f $WORKDIR/delta-rays\[tau\].root:h0 \
71 -f $WORKDIR/delta-rays\[positron\].root:h0 \
72 -f $WORKDIR/delta-rays\[electron\].root:h0 \
76 -\^ "dE/dx [MeV g^{-1} cm^{2}]" \
80 -o delta-rays.$FORMAT $BATCH
86 -f $WORKDIR/delta-rays\[muon\].root:h1 \
87 -f $WORKDIR/delta-rays\[tau\].root:h1 \
88 -f $WORKDIR/delta-rays\[positron\].root:h1 \
89 -f $WORKDIR/delta-rays\[electron\].root:h1 \
92 -\^ "dN/dx [g^{-1} cm^{2}]" \
95 -o delta-rays_n.$FORMAT $BATCH