Jpp  16.0.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JCompareDetector.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Utility script to compare two detectors and to create graphics output.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
19 
20 set_variable WORKDIR ${TMPDIR:-/tmp}
21 set_variable DEBUG 1
22 set_variable: TIME PRECISION_TIME 0.01 # [ns]
23 set_variable: ANGLE PRECISION_ANGLE 0.01 # [deg]
24 set_variable: POSITION PRECISION_POSITION 0.01 # [m]
25 set_variable: ORIENTATION PRECISION_ORIENTATION 0.01 # [rad]
26 set_variable: FORMAT GRAPHICS_FORMAT gif
27 set_variable+ BATCH GRAPHICS_BATCH -B
28 
29 if do_usage $*; then
30  usage "$script <detector file> <detector file> [<lower floor>-<upper floor>]"\
31  "\nIf the range of floors is specified, the first detector file is aligned to the second before the comparison."\
32  "\nIn this, only modules within the given range of floors are used."
33 fi
34 
35 case $# in
36  3) set_variable RANGE $argv[3];&
37  2) set_array DETECTOR $argv[1] $argv[2];;
38  *) fatal "Wrong number of arguments."
39 esac
40 
41 if [[ "$RANGE" == [0-9]*"-"[0-9]* ]]; then
42 
43  JConvertDetectorFormat \
44  -a $DETECTOR[1] \
45  -o ${TMPDIR:-/tmp}/detector.datx \
46  -d $DEBUG
47 
48  echo "" > ${TMPDIR:-/tmp}/tripod.txt
49 
50  JAlignDetector \
51  -a ${TMPDIR:-/tmp}/detector.datx \
52  -b $DETECTOR[2] \
53  -AO \
54  -r "${RANGE%%-*} ${RANGE##*-}" \
55  -T ${TMPDIR:-/tmp}/tripod.txt \
56  -d $DEBUG --!
57 
58  DETECTOR[1]=${TMPDIR:-/tmp}/detector.datx
59 
60  rm -f ${TMPDIR:-/tmp}/tripod.txt
61 fi
62 
63 JCompareDetector \
64  -a $DETECTOR[1] \
65  -b $DETECTOR[2] \
66  -p "time = $TIME" \
67  -p "angle = $ANGLE" \
68  -p "position = $POSITION" \
69  -p "orientation = $ORIENTATION" \
70  -o $WORKDIR/detector.root \
71  -d $DEBUG --!
72 
73 typeset -A HISTOGRAMS
74 
75 HISTOGRAMS[X2]="#Deltax [m]"
76 HISTOGRAMS[Y2]="#Deltay [m]"
77 HISTOGRAMS[Z2]="#Deltaz [m]"
78 #HISTOGRAMS[R2]="angle [deg]"
79 
80 for KEY VALUE in ${(@kv)HISTOGRAMS}; do
81 
82  JPlot2D \
83  -f $WORKDIR/detector.root:${KEY} \
84  -\> "string" \
85  -\< "floor" \
86  -\^ "$VALUE" \
87  -O COLZ \
88  -T "$KEY" \
89  -o ${KEY}.$FORMAT $BATCH
90 
91 done
92 
93 eval `JPrintDetector -a $DETECTOR[1] -O SUMMARY`
94 
95 for (( i = 1; $i <= ${#STRINGS}; ++i )); do
96 
97  set_variable STRING $STRINGS[$i]
98 
99  let XMIN="$i - 1 - 0.5"
100  let XMAX="$i - 1 + 0.5"
101 
102  JProject2D \
103  -f $WORKDIR/detector.root:Z2 \
104  -P Y \
105  -x "$XMIN $XMAX" \
106  -o ${TMPDIR:-/tmp}/py_\[${STRING}\].root \
107  -d 0
108 
109 done
110 
111 JPlot1D \
112  -f${TMPDIR:-/tmp}/py_\[${^STRINGS}\].root:\.\* \
113  -\> "floor" \
114  -\^ "#Deltaz [m]" \
115  -x "0.5 18.5" \
116  -O "HIST ][" \
117  -T "" \
118  -L BR \
119  -o Z1.$FORMAT $BATCH
HISTOGRAMS
do echo n Creating graphics for string $STRING HISTOGRAMS
Definition: plot-canberra.sh:74
$argv
*set_variable DETECTOR_ID $argv[1] set_array RUNS $argv[2,-2]
Definition: JDataQuality.sh:24
STRING
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
Definition: examples/JTrigger/JTriggerEfficiencyRunByRun.sh:64
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
exit
Definition: JPizza.sh:36
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:758
Y
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Definition: detector-XY:fit1d.sh:33
echo
then echo
Definition: JDOMDAQDriver.sh:56
$STRINGS
do JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T $WORKDIR tripod txt V $WORKDIR sound_velocity txt M $WORKDIR mechanics txt H $WORKDIR hydrophone txt E $EMITTER $DISABLE d $DEBUG!done kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR for STRING in $STRINGS[*]
Definition: JCanberra.sh:64
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
TMPDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
Definition: JEnergyCorrection.sh:21
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
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
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
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
XMIN
then let XMIN
Definition: detector-XY:fit.sh:67
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
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
script
script
Definition: JAcoustics.sh:2
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
$i
do alias $i
Definition: JDOMDAQDriver.sh:23
Generated by   doxygen 1.8.5