Jpp  debug
the software that should make you happy
JPlotPDF_L1.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 tabulated PDFs.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable NORTH "0.00 0.00"
25 set_variable EAST "1.57 0.00"
26 set_variable SOUTH "3.14 0.00"
27 set_variable WEST "1.57 3.14"
28 
29 set_variable WORKDIR ${TMPDIR:-/tmp}
30 set_variable E 1e3
31 set_variable R 50
32 set_variable T_NS 10
33 set_variable: YMAX PHYSICS_YMAX 2e-0
34 set_variable: YMIN PHYSICS_YMIN 1e-5
35 set_variable: HISTOGRAM PHYSICS_HISTOGRAM "1040 -20 +500"
36 
37 set_variable PDF $JPP_DATA/J%p.dat
38 
39 set_variable: FORMAT GRAPHICS_FORMAT gif
40 set_variable+ BATCH GRAPHICS_BATCH -B
41 
42 if do_usage $*; then
43  usage "$script [energy [distance]]"
44 fi
45 
46 case $# in
47  2) set_variable R $2;&
48  1) set_variable E $1;;
49 esac
50 
51 
52 if (( 1 )); then
53 
54  for dir in NORTH EAST SOUTH WEST; do
55 
56  echo "Generating $dir."
57 
58  eval D=\$$dir
59 
60  $JPP_DIR/examples/JPhysics/JPlotPDF_L1 \
61  -f $PDF \
62  -D "${D}" -R $R -E $E \
63  -H "$HISTOGRAM" \
64  -T $T_NS \
65  -o $WORKDIR/F\[${dir}\].root
66 
67  done
68 
69 fi
70 
71 
72 if (( 1 )); then
73 
74  JPlot1D \
75  -f $WORKDIR/F\[NORTH\].root:h0 \
76  -f $WORKDIR/F\[EAST\].root:h0 \
77  -f $WORKDIR/F\[SOUTH\].root:h0 \
78  -f $WORKDIR/F\[WEST\].root:h0 \
79  -y "$YMIN $YMAX" -Y \
80  -> "#Deltat [ns]" \
81  -\^ "dP/dt [npe/ns]" \
82  -L TR \
83  -T "E = ${E} GeV#; R = ${R} m" \
84  -o PDF.$FORMAT $BATCH
85 
86 fi
Generated by doxygen 1.9.1