Jpp - the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
hydrophone-XY:run.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 if [ -z $JPP_DIR ]; then
5  echo "Variable JPP_DIR undefined."
6  exit
7 fi
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 zmodload zsh/mathfunc
12 
13 set_variable: DEBUG ACOUSTICS_DEBUG 2
14 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
15 set_variable: QUALITY_FACTOR ACOUSTICS_QUALITYFACTOR 0
16 
17 if do_usage $*; then
18  usage "$script <string identifier> <detector file> (input file)+ (toashort file)+"\
19  "\nNote that the input files and toashort files should be one-to-one related."
20 fi
21 
22 if (( ${#} < 4 || ${#}%2 != 0 )); then
23  fatal "Wrong number of arguments."
24 fi
25 
26 let N="(${#} - 2) / 2"
27 
30 set_array INPUT_FILES $argv[3,2+$N]
31 set_array TOASHORT_FILES $argv[-$N,-1]
32 
33 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
34 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
35 
36 source JAcoustics.sh $DETECTOR_ID
37 
38 CHECK_EXIT_CODE
39 
40 typeset -A TRIPODS
41 
42 get_tripods $WORKDIR/tripod.txt TRIPODS
43 
44 if (( ${#TRIPODS} < 2 )); then
45  fatal "Not enough tripods " ${#TRIPODS}
46 fi
47 
48 getModule -a $DETECTOR -L "$STRINGS[1] 0" >& /dev/null
49 
50 if (( $? != 0 )); then
51 
52  JCookie.sh
53 
54  JDetectorDB \
55  -a $DETECTOR \
56  -o $DETECTOR \
57  -W --!
58 fi
59 
60 if [[ -z "${STRINGS[(r)$STRING]}" ]]; then
61  fatal "Invalid string $STRING; "\
62  "possible values: ${STRINGS[*]}"
63 fi
64 
65 set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
66 set_variable HYDROPHONE_TXT $WORKDIR/hydrophonexy_${STRING}.txt
67 set_array QUANTILES 0.333 0.666
68 set_variable TMPDIR $WORKDIR/.$$
69 
70 mkdir -p $TMPDIR
71 
72 cp hydrophone.txt $TMPDIR/.
73 
74 for (( X = -2.0; $X <= +2.0; X += 0.1 )); do
75  for (( Y = -2.0; $Y <= +2.0; Y += 0.1 )); do
76 
77  JEditHydrophone -f $TMPDIR/hydrophone.txt -S "$STRING set $X $Y 0.0"
78 
79  $JPP_DIR/examples/JAcoustics/JHydrophone \
80  -a $DETECTOR \
81  -f "$INPUT_FILES[*]" \
82  -i "$TOASHORT_FILES[*]" \
83  -o $TMPDIR/hydrophone.root \
84  -H $TMPDIR/hydrophone.txt \
85  -T $WORKDIR/tripod.txt \
86  -V $WORKDIR/sound_velocity.txt \
87  -W $WORKDIR/waveform.txt \
88  -Q $QUALITY_FACTOR \
89  -d $DEBUG --!
90 
91  XMEAN=(0.0 0.0)
92  STDEV=(0.0 0.0)
93  TOTAL=0.0
94 
95  for TRIPOD in ${(@k)TRIPODS}; do
96 
97  JPrintQuantiles \
98  -f "$TMPDIR/hydrophone.root:${MODULE}\.${TRIPOD}" \
99  -Q "$QUANTILES[*]" \
100  -d 0 | read -A RESULT
101 
102  VALUE=$((0.5e3 * ($RESULT[-1] + $RESULT[1]))) # average [ms]
103  TOTAL=$(($TOTAL + 1.0))
104 
105  if (( $TOTAL == 1.0 )); then
106  XMEAN[1]=$VALUE
107  STDEV[1]=0.0
108  else
109  XMEAN[2]=$(($XMEAN[1] + ($VALUE - $XMEAN[1]) / $TOTAL))
110  STDEV[2]=$(($STDEV[1] + ($VALUE - $XMEAN[1]) * ($VALUE - $XMEAN[2])))
111  XMEAN[1]=$XMEAN[2]
112  STDEV[1]=$STDEV[2]
113  fi
114  done
115 
116  printf "%7.3f %7.3f %9.5f\n" $X $Y $((sqrt($TOTAL * $STDEV[1] / ($TOTAL * ($TOTAL - 1.0))))) >> $HYDROPHONE_TXT
117  done
118 done
119 
120 rm -rf $TMPDIR
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Q(UTCMax_s-UTCMin_s)-livetime_s
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
bool read(Vec &v, std::istream &is)
Read a Vec(tor) from a stream.
Definition: io_ascii.hh:141
static const double H
Planck constant [eV s].
exit
Definition: JPizza.sh:36
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
then JPizza f
Definition: JPizza.sh:46
then echo
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then break fi done getCenter read X Y Z let X
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
* usage
set_array INPUT_FILES
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
do set_variable SIGMA_NS set_variable OUTLIERS set_variable OUTPUT_FILE matrix[${ALPHA_DEG}\deg\] root $JPP JMatrixNZ a $DETECTOR f $INPUT_FILE o $OUTPUT_FILE S
Definition: JMatrixNZ.sh:56
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then JCalibrateToT a
Definition: JTuneHV.sh:116
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
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
then display $WORKDIR
Definition: plot-Domino.sh:127
then cp
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable HYDROPHONE_TXT $WORKDIR hydrophonexy_
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
set_variable DETECTOR
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
then usage $script[input file[working directory[option]]] nWhere option can be N
Definition: JMuonPostfit.sh:35
script
Definition: JAcoustics.sh:2
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
esac done
Definition: JAddHDE.sh:21
then QUANTILES