Jpp - the software that should make you happy
 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 
23 
24 if ( do_usage $* ); then
25  usage "$script [working directory]"
26 fi
27 
28 case $# in
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
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
exit
Definition: JPizza.sh:36
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:79
then JRadiation o $WORKDIR lambda root O length d $DEBUG for TYPE in eerad Brems gnrad
Definition: JRadiation.sh:39
then echo
do set_variable OUTPUT_DIRECTORY $WORKDIR T
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
* usage
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi mv $WORKDIR/fit.root $MODULE_ROOT typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
Definition: module-Z:fit.sh:84
then for APP in event gandalf start energy
Definition: JMuonMCEvt.sh:42
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
alias put_queue eval echo n
Definition: qlib.csh:19
then display $WORKDIR
Definition: plot-Domino.sh:127
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:38
version
Definition: JCalibratePMT.sh:7
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62