Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JDeltaRays.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to determine dE/dx due to visible delta-rays
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23
24set_variable DEBUG 2
25set_variable WORKDIR ./
26set_variable TMIN_GEV 0.000915499
27
28set_variable+ NUMERICAL DELTARAYS_NUMERICAL -N
29
30set_variable: FORMAT GRAPHICS_FORMAT gif
31set_variable+ BATCH GRAPHICS_BATCH -B
32
33if do_usage $*; then
34 usage "$script [working directory]"
35fi
36
37case $# in
38 1) set_variable WORKDIR $1;;
39esac
40
41
42if (( 1 )); then
43
44 for LEPTON in muon tau positron electron; do
45
46 $JPP_DIR/examples/JPhysics/JDeltaRays \
47 -o $WORKDIR/delta-rays\[${LEPTON}\].root \
48 -L $LEPTON \
49 -O Q \
50 -T "$TMIN_GEV 1e10" \
51 ${NUMERICAL} -d $DEBUG
52
53 JPlot1D \
54 -f $WORKDIR/delta-rays\[${LEPTON}\].root:h0 \
55 -XX \
56 -y "0 2.5" \
57 -> "E [GeV]" \
58 -\^ "dE/dx [MeV g^{-1} cm^{2}]" \
59 -L TL -T "" \
60 -S 0.5 \
61 -O "][" \
62 -o $LEPTON.$FORMAT $BATCH
63 done
64fi
65
66if (( 1 )); then
67
68 JPlot1D \
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 \
73 -XX \
74 -y "0 2.5" \
75 -> "E [GeV]" \
76 -\^ "dE/dx [MeV g^{-1} cm^{2}]" \
77 -L TL -T "" \
78 -S 0.5 \
79 -O "][" \
80 -o delta-rays.$FORMAT $BATCH
81fi
82
83if (( 1 )); then
84
85 JPlot1D \
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 \
90 -XX \
91 -> "E [GeV]" \
92 -\^ "dN/dx [g^{-1} cm^{2}]" \
93 -L BR -T "" \
94 -O "][" \
95 -o delta-rays_n.$FORMAT $BATCH
96fi
97