Jpp
JRadiation.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to plot radiation cross sections, shower energy, range and b(E)
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 if [ -z $JPP_DIR ]; then
15  echo "Variable JPP_DIR undefined."
16  exit
17 fi
18 
19 source $JPP_DIR/setenv.sh $JPP_DIR
20 
21 if ( do_usage $* ); then
22  usage "$script"
23 fi
24 
25 set_variable DEBUG 2
26 set_variable WORKDIR ./
27 
28 
29 if (( 1 )); then
30 
31  JRadiation \
32  -o $WORKDIR/lambda.root \
33  -O length \
34  -d $DEBUG
35 
36  for TYPE in eerad Brems gnrad; do
37 
38  JPlot1D \
39  -f "$WORKDIR/lambda.root:${TYPE}" \
40  -y "1e-1 1e6" \
41  -XXY \
42  -\^ "\lambda [m]" \
43  -> "E [GeV]" \
44  -L TR -T "interaction length" \
45  -o lamdba-$TYPE.gif
46  done
47 
48 fi
49 
50 
51 if (( 1 )); then
52 
53  JRadiation \
54  -o $WORKDIR/radiation.root \
55  -O energy \
56  -n 40000 \
57  -d $DEBUG
58 
59  for TYPE in eerad Brems gnrad; do
60 
61  JPlot1D \
62  -f "$WORKDIR/radiation.root:$TYPE" \
63  -y "0.1 1e9" \
64  -XXY \
65  -\^ "<E_{s}> [GeV]" \
66  -> "E [GeV]" \
67  -L TL -T "shower energy" \
68  -o Es-$TYPE.gif
69  done
70 
71 fi
72 
73 
74 if (( 1 )); then
75 
76  JRadiation \
77  -o $WORKDIR/radiation.root \
78  -O range \
79  -n 1000 \
80  -d $DEBUG
81 
82  JPlot1D \
83  -f $WORKDIR/radiation.root:R\.\* \
84  -y "0 40" \
85  -XX \
86  -> "E [GeV]" \
87  -\^ "R [km]" \
88  -L TL -T "muon range" \
89  -o range.gif
90 
91 fi
92 
93 
94 if (( 1 )); then
95 
96  JRadiation \
97  -o $WORKDIR/radiation.root \
98  -O eloss \
99  -n 4000000 \
100  -d $DEBUG
101 
102  JPlot1D \
103  -f $WORKDIR/radiation.root:hb \
104  -y "0 0.5" \
105  -XX \
106  -> "E [GeV]" \
107  -\^ "b [km^{-1}]" \
108  -T "Energy loss" \
109  -o b.gif
110 
111 fi
Generated by   doxygen 1.8.16