Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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 set_variable DEBUG 2
22 set_variable WORKDIR /tmp/
23 
24 if ( do_usage $* ); then
25  usage "$script [working directory]"
26 fi
27 
28 case $# in
29  1) set_variable WORKDIR $1;;
30 esac
31 
32 if (( 1 )); then
33 
34  JRadiation \
35  -o $WORKDIR/lambda.root \
36  -O length \
37  -d $DEBUG
38 
39  for TYPE in eerad Brems gnrad; do
40 
41  JPlot1D \
42  -f "$WORKDIR/lambda.root:${TYPE}" \
43  -y "1e-1 1e6" \
44  -XXY \
45  -\^ "\lambda [m]" \
46  -\> "E [GeV]" \
47  -L TR -T "interaction length" \
48  -o lamdba-$TYPE.gif
49  done
50 
51 fi
52 
53 
54 if (( 1 )); then
55 
56  JRadiation \
57  -o $WORKDIR/radiation.root \
58  -O energy \
59  -n 40000 \
60  -d $DEBUG
61 
62  for TYPE in eerad Brems gnrad; do
63 
64  JPlot1D \
65  -f "$WORKDIR/radiation.root:$TYPE" \
66  -y "0.1 1e9" \
67  -XXY \
68  -\^ "<E_{s}> [GeV]" \
69  -\> "E [GeV]" \
70  -L TL -T "shower energy" \
71  -o Es-$TYPE.gif
72  done
73 
74 fi
75 
76 
77 if (( 1 )); then
78 
79  JRadiation \
80  -o $WORKDIR/radiation.root \
81  -O range \
82  -n 1000 \
83  -d $DEBUG
84 
85  JPlot1D \
86  -f $WORKDIR/radiation.root:R\.\* \
87  -y "0 40" \
88  -XX \
89  -\> "E [GeV]" \
90  -\^ "R [km]" \
91  -L TL -T "muon range" \
92  -o range.gif
93 
94 fi
95 
96 
97 if (( 1 )); then
98 
99  JRadiation \
100  -o $WORKDIR/radiation.root \
101  -O eloss \
102  -n 4000000 \
103  -d $DEBUG
104 
105  JPlot1D \
106  -f $WORKDIR/radiation.root:hb \
107  -y "0 0.5" \
108  -XX \
109  -\> "E [GeV]" \
110  -\^ "b [km^{-1}]" \
111  -T "Energy loss" \
112  -o b.gif
113 
114 fi
JSYSTEM::set_variable
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
Definition: JSystemToolkit.hh:404
in
esac print_variable DETECTOR INPUT_FILE OUTPUT_FILE CDF for TYPE in
Definition: JSirene.sh:45
exit
exit
Definition: JPizza.sh:36
gnrad
then JRadiation o $WORKDIR lambda root O length d $DEBUG for TYPE in eerad Brems gnrad
Definition: JRadiation.sh:39
f
then JPizza f
Definition: JPizza.sh:46
echo
then echo
Definition: JDOMDAQDriver.sh:56
tmp
do montage tile geometry tmp
Definition: plot-PMTParameters2D.sh:116
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
d
then print_variable DETECTOR INPUT_FILE INTERMEDIATE_FILE check_input_file $DETECTOR $INPUT_FILE check_output_file $INTERMEDIATE_FILE $OUTPUT_FILE JMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JPath.sh:52
usage
* usage
Definition: JLigier-local.sh:33
R
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
energy
then for APP in event gandalf start energy
Definition: JMuonMCEvt.sh:42
n
alias put_queue eval echo n
Definition: qlib.csh:19
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:142
version
version
Definition: JCalibratePMT.sh:7
WORKDIR
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
done
esac done
Definition: JAddHDE.sh:21
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5