Jpp - the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JMuonMCEvt.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 test JGandalf - JEnergy - JStart on Monte Carlo true muons.
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 $JPP_DATA
23 set_variable INPUT_FILE $WORKDIR/trigger_efficiency+background.root
24 set_variable DETECTOR $WORKDIR/km3net_reference.detx
25 set_variable PDF $JPP_DATA/J%p.dat
26 set_variable ENERGY_CORRECTION energy_correction_arca.txt
27 
28 if do_usage $*; then
29  usage "$script [detector file [input file [working directory]]]"
30 fi
31 
32 case $# in
33  3) set_variable WORKDIR $3;&
34  2) set_variable INPUT_FILE $2;&
35  1) set_variable DETECTOR $1;;
36 esac
37 
38 set_variable OUTPUT_FILE $WORKDIR/%-mc.root
39 
40 if (( 1 )); then
41 
42  for APP in event gandalf start energy; do
43  rm -f ${OUTPUT_FILE/\%/${APP}}
44  done
45 
46  JMuonMCEvt \
47  -f $INPUT_FILE \
48  -o ${OUTPUT_FILE/\%/event} \
49  -d $DEBUG --!
50 
51  PARAMETERS=(
52  TTS_ns 2
53  E_GeV 1.0e3
54  TMin_ns -50.0
55  TMax_ns +450.0
56  R_Hz 6.0e3
57  roadWidth_m 175.0
58  numberOfPrefits 0)
59 
60  eval JMuonGandalf \
61  -a $DETECTOR \
62  -f ${OUTPUT_FILE/\%/event} \
63  -o ${OUTPUT_FILE/\%/gandalf} \
64  -P $PDF \
65  -@ `make_equation PARAMETERS` \
66  -d $DEBUG --!
67 
68  PARAMETERS=(
69  roadWidth_m 160.0
70  R_Hz 6.0e3
71  TMin_ns -25.0
72  TMax_ns +75.0
73  Pmin1 1.0e-3
74  Pmin2 1.0e-2
75  numberOfPrefits 1)
76 
77  eval JMuonStart \
78  -a $DETECTOR \
79  -f ${OUTPUT_FILE/\%/gandalf} \
80  -o ${OUTPUT_FILE/\%/start} \
81  -P $PDF \
82  -@ `make_equation PARAMETERS` \
83  -d $DEBUG --!
84 
85  PARAMETERS=(
86  roadWidth_m 200.0
87  R_Hz 6.0e3
88  EMin_log 0.0
89  EMax_log 8.0
90  TMin_ns -50.0
91  TMax_ns +450.0
92  mestimator 1
93  numberOfPrefits 1)
94 
95  eval JMuonEnergy \
96  -a $DETECTOR \
97  -f ${OUTPUT_FILE/\%/start} \
98  -o ${OUTPUT_FILE/\%/energy} \
99  -P $PDF \
100  -E energy_correction_arca.txt \
101  -@ `make_equation PARAMETERS` \
102  -d $DEBUG --!
103 fi
104 
105 if (( 1 )); then
106 
107  JMuonPostfit \
108  -f ${OUTPUT_FILE/\%/energy} \
109  -O LOGE \
110  -A 4 \
111  -o $WORKDIR/postfit\[MC\].root
112 fi
113 
114 if (( 1 )); then
115 
116  JPlot1D \
117  -f $WORKDIR/postfit\[MC\].root:hx \
118  -\> "angle [deg]" \
119  -\^ "number of events [a.u.]" \
120  -XX \
121  -L TR \
122  -o angle-truth.gif
123 fi
124 
125 if (( 1 )); then
126 
127  JPlot1D \
128  -f $WORKDIR/postfit\[MC\].root:e0 \
129  -f $WORKDIR/postfit\[MC\].root:e1 \
130  -f $WORKDIR/postfit\[MC\].root:e2 \
131  -\^ "number of events [a.u.]" \
132  -\> "log10(E_{fit}/GeV)" \
133  -XX
134 fi
135 
136 if (( 1 )); then
137 
138  JPlot2D \
139  -f $WORKDIR/postfit\[MC\].root:ee \
140  -O COLZ \
141  -XX -YY \
142  -\> "E_{true} [GeV]" \
143  -\< "E_{fit} [GeV]"
144 fi
145 
146 if (( 1 )); then
147 
148  JFit \
149  -f $WORKDIR/postfit\[MC\].root:er \
150  -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))" \
151  -@ "p0 = 0.8*GetMaximum" \
152  -@ "p1 = 0.0" \
153  -@ "p2 = 0.4" \
154  -x "-0.8 +0.8"
155 
156  JPlot1D \
157  -f fit.root:er \
158  -\> "log10(E_{fit}/E_{true})" \
159  -o er-truth.gif
160 fi
161 
162 if (( 1 )); then
163 
164  JPlot2D \
165  -f $WORKDIR/postfit\[MC\].root:ha \
166  -OCOLZ \
167  -\> "R^{2} [m^{2}]" \
168  -\< "z [m]"
169 fi
170 
171 if (( 1 )); then
172 
173  JPlot1D \
174  -f $WORKDIR/postfit\[MC\].root:hz \
175  -\> "z [m]" \
176  -\^ "number of events [a.u.]" \
177  -L TR \
178  -o z.gif
179 fi
180 
181 if (( 0 )); then
182 
183  JPlot1D \
184  -f $WORKDIR/postfit\[MC\].root:hx \
185  -f $WORKDIR/postfit\[MC\].root:hx \
186  -\> "angle [deg]" \
187  -\^ "number of events [a.u.]" \
188  -XX \
189  -L TR \
190  -o angle-truth.gif
191 fi
192 
193 if (( 0 )); then
194 
195  JSum1D -f $WORKDIR/postfit\[MC\].root:hx -o $WORKDIR/sum\[MC\].root -N
196  JSum1D -f $WORKDIR/postfit\[gandalf\].root:hx -o $WORKDIR/sum\[gandalf\].root -N
197 
198  JPlot1D \
199  -f $WORKDIR/sum\[gandalf\].root:hx \
200  -f $WORKDIRsum\[MC\].root:hx \
201  -\> "angle [deg]" \
202  -\^ "integral of events" \
203  -XX \
204  -L TR \
205  -o sum-truth.gif
206 
207 fi
PDF
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable PDF $JPP_DATA J p dat set_variable HISTOGRAM if do_usage *then usage $script< detector file >< inputfile >< outputfile > fi case set_variable PDF
Definition: JARCACalibrateMuon.sh:23
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
mestimator
mestimator
Definition: acoustics-fit.sh:290
exit
exit
Definition: JPizza.sh:36
f
then JPizza f
Definition: JPizza.sh:46
echo
then echo
Definition: JDOMDAQDriver.sh:56
APP
do set_variable APP
Definition: JCalibrateHeight.sh:80
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
o
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
F
then awk F
Definition: getMIIseacurrent.sh:82
usage
* usage
Definition: JLigier-local.sh:34
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
energy
then for APP in event gandalf start energy
Definition: JMuonMCEvt.sh:42
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
Definition: module-Z:fit.sh:116
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:127
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR tmp set_variable JDAQ_TIMESLICE JDAQTimeslice if do_usage *then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:23
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
in
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
N
then usage $script[input file[working directory[option]]] nWhere option can be N
Definition: JMuonPostfit.sh:35
version
version
Definition: JCalibratePMT.sh:7
A
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
Definition: JShellParser.csh:15
E
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:35
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
done
esac done
Definition: JAddHDE.sh:21
Generated by   doxygen 1.8.5