Jpp  16.0.0-rc.2
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 set_variable: FORMAT GRAPHICS_FORMAT gif
28 set_variable+ BATCH GRAPHICS_BATCH -B
29 
30 if do_usage $*; then
31  usage "$script [detector file [input file [working directory]]]"
32 fi
33 
34 case $# in
35  3) set_variable WORKDIR $3;&
36  2) set_variable INPUT_FILE $2;&
37  1) set_variable DETECTOR $1;;
38 esac
39 
40 set_variable OUTPUT_FILE $WORKDIR/%-mc.root
41 
42 if (( 1 )); then
43 
44  for APP in event gandalf start energy; do
45  rm -f ${OUTPUT_FILE/\%/${APP}}
46  done
47 
48  JMuonMCEvt \
49  -f $INPUT_FILE \
50  -o ${OUTPUT_FILE/\%/event} \
51  -d $DEBUG --!
52 
53  PARAMETERS=(
54  TTS_ns 2
55  E_GeV 1.0e3
56  TMin_ns -50.0
57  TMax_ns +450.0
58  R_Hz 6.0e3
59  roadWidth_m 175.0
60  numberOfPrefits 0)
61 
62  eval JMuonGandalf \
63  -a $DETECTOR \
64  -f ${OUTPUT_FILE/\%/event} \
65  -o ${OUTPUT_FILE/\%/gandalf} \
66  -P $PDF \
67  -@ `make_equation PARAMETERS` \
68  -d $DEBUG --!
69 
70  PARAMETERS=(
71  roadWidth_m 160.0
72  R_Hz 6.0e3
73  TMin_ns -25.0
74  TMax_ns +75.0
75  Pmin1 1.0e-3
76  Pmin2 1.0e-2
77  numberOfPrefits 1)
78 
79  eval JMuonStart \
80  -a $DETECTOR \
81  -f ${OUTPUT_FILE/\%/gandalf} \
82  -o ${OUTPUT_FILE/\%/start} \
83  -P $PDF \
84  -@ `make_equation PARAMETERS` \
85  -d $DEBUG --!
86 
87  PARAMETERS=(
88  roadWidth_m 200.0
89  R_Hz 6.0e3
90  EMin_log 0.0
91  EMax_log 8.0
92  TMin_ns -50.0
93  TMax_ns +450.0
94  mestimator 1
95  numberOfPrefits 1)
96 
97  eval JMuonEnergy \
98  -a $DETECTOR \
99  -f ${OUTPUT_FILE/\%/start} \
100  -o ${OUTPUT_FILE/\%/energy} \
101  -P $PDF \
102  -E energy_correction_arca.txt \
103  -@ `make_equation PARAMETERS` \
104  -d $DEBUG --!
105 fi
106 
107 if (( 1 )); then
108 
109  JMuonPostfit \
110  -f ${OUTPUT_FILE/\%/energy} \
111  -O LOGE \
112  -A 4 \
113  -o $WORKDIR/postfit\[MC\].root
114 fi
115 
116 if (( 1 )); then
117 
118  JPlot1D \
119  -f $WORKDIR/postfit\[MC\].root:hx \
120  -\> "angle [deg]" \
121  -\^ "number of events [a.u.]" \
122  -XX \
123  -L TR \
124  -o angle-truth.$FORMAT $BATCH
125 fi
126 
127 if (( 1 )); then
128 
129  JPlot1D \
130  -f $WORKDIR/postfit\[MC\].root:e0 \
131  -f $WORKDIR/postfit\[MC\].root:e1 \
132  -f $WORKDIR/postfit\[MC\].root:e2 \
133  -\^ "number of events [a.u.]" \
134  -\> "log10(E_{fit}/GeV)" \
135  -XX \
136  -o e0.$FORMAT $BATCH
137 fi
138 
139 if (( 1 )); then
140 
141  JPlot2D \
142  -f $WORKDIR/postfit\[MC\].root:ee \
143  -O COLZ \
144  -XX -YY \
145  -\> "E_{true} [GeV]" \
146  -\< "E_{fit} [GeV]" \
147  -o ee.$FORMAT $BATCH
148 fi
149 
150 if (( 1 )); then
151 
152  JFit \
153  -f $WORKDIR/postfit\[MC\].root:er \
154  -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))" \
155  -@ "p0 = 0.8*GetMaximum" \
156  -@ "p1 = 0.0" \
157  -@ "p2 = 0.4" \
158  -x "-0.8 +0.8"
159 
160  JPlot1D \
161  -f fit.root:er \
162  -\> "log10(E_{fit}/E_{true})" \
163  -o er-truth.$FORMAT $BATCH
164 fi
165 
166 if (( 1 )); then
167 
168  JPlot2D \
169  -f $WORKDIR/postfit\[MC\].root:ha \
170  -OCOLZ \
171  -\> "R^{2} [m^{2}]" \
172  -\< "z [m]" \
173  -o ha.$FORMAT $BATCH
174 fi
175 
176 if (( 1 )); then
177 
178  JPlot1D \
179  -f $WORKDIR/postfit\[MC\].root:hz \
180  -\> "z [m]" \
181  -\^ "number of events [a.u.]" \
182  -L TR \
183  -o z.$FORMAT $BATCH
184 fi
185 
186 if (( 0 )); then
187 
188  JPlot1D \
189  -f $WORKDIR/postfit\[MC\].root:hx \
190  -f $WORKDIR/postfit\[MC\].root:hx \
191  -\> "angle [deg]" \
192  -\^ "number of events [a.u.]" \
193  -XX \
194  -L TR \
195  -o angle-truth.$FORMAT $BATCH
196 fi
197 
198 if (( 0 )); then
199 
200  JSum1D -f $WORKDIR/postfit\[MC\].root:hx -o $WORKDIR/sum\[MC\].root -N
201  JSum1D -f $WORKDIR/postfit\[gandalf\].root:hx -o $WORKDIR/sum\[gandalf\].root -N
202 
203  JPlot1D \
204  -f $WORKDIR/sum\[gandalf\].root:hx \
205  -f $WORKDIRsum\[MC\].root:hx \
206  -\> "angle [deg]" \
207  -\^ "integral of events" \
208  -XX \
209  -L TR \
210  -o sum-truth.$FORMAT $BATCH
211 
212 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
E
then usage E
Definition: JMuonPostfit.sh:35
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
N
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
Definition: JShowerPostfit.sh:95
mestimator
mestimator
Definition: acoustics-fit.sh:494
exit
exit
Definition: JPizza.sh:36
YY
set_variable LOGY YY
Definition: JMuonPostfit.sh:73
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
done
set_variable NUMBER_OF_ITERATIONS set_variable EPSILON cat acoustics_fit_parameters txt<< EOF $CONFIGURATION[*]Nmin=3;sigma_s=100.0e-6;stdev=10.0;mestimator=0;fixStrings=0;EOF for STRING in $STRINGS[*];do#fit stretching and(z) position of given string set_variable DETECTOR_TMP ${TMPDIR:-/tmp}/detector_A.datx JEditDetector-a $DETECTOR-o $DETECTOR_TMP-r $STRING JEditDetector-a $DETECTOR-o $DETECTOR-k $STRING for MUL in 0.005 0.001;do DX_M=0.2 for((N=0;$N< $NUMBER_OF_ITERATIONS;++N));do CHI2[3]=$CHI2[1] fitPositionOfString $STRING Z $DX_M fitStretchingOfString $STRING $MUL if(($CHI2[3]-$CHI2[1]< $EPSILON));then break fi done if(($N >=$NUMBER_OF_ITERATIONS));then printf"warning: reached maximum number of iterations %d - converenge %7.3f\n"$N $(($CHI2[3]-$CHI2[1])) fi done JMergeDetector-a $DETECTOR-a $DETECTOR_TMP-o $DETECTOR rm-f $DETECTOR_TMP JConvertDetectorFormat-a $DETECTOR-o $DETECTOR-r-d 0 > &dev null done
Definition: acoustics-fit.sh:563
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
XX
set_variable LOG XX
Definition: JMuonPostfit.sh:71
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
energy
then for APP in event gandalf start energy
Definition: JMuonMCEvt.sh:44
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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:118
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:128
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:42
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
version
version
Definition: JCalibratePMT.sh:7
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
Generated by   doxygen 1.8.5