Jpp  18.5.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
plot-acoustics-fit.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 
15 set_variable: PRECISION ACOUSTICS_PRECISION 0.1
16 set_variable: FORMAT GRAPHICS_FORMAT gif
17 set_variable+ BATCH GRAPHICS_BATCH -B
18 
19 if do_usage $*; then
20  usage "$script <directory>"\
21  "\nInput directory corresponds to location of stage files of a pre-calibration procedure."
22 fi
23 
24 if (( $# == 1 )); then
25  set_variable DIR $1
26 else
27  fatal "Wrong number of arguments."
28 fi
29 
30 for INPUT_FILE in `ls -rt $DIR/stage-*.log`; do
31 
32  grep "\(chi2\[.\]\|gradient\)" $INPUT_FILE > $WORKDIR/stage.txt
33 
34  awk 'BEGIN { YMAX = 0.0; YMIN = 1.0e10 } END { printf "%9.5f %9.5f", YMIN, YMAX} /chi2/ { if ($3 < 1.0e10) { if ($3 < YMIN) { YMIN = $3}; if ($3 > YMAX) { YMAX = $3 }}}' $WORKDIR/stage.txt | read YMIN YMAX
35 
36  let "Y = $YMAX - $YMIN"
37  let "YMAX = floor(($YMAX + 0.1 * ($Y)) / $PRECISION + 1.0) * $PRECISION"
38  let "YMIN = floor(($YMIN - 0.1 * ($Y)) / $PRECISION - 0.0) * $PRECISION"
39 
40  if (( $YMIN < 0.0 )); then
41  let "YMIN = 0.0"
42  fi
43 
44  OUTPUT_FILES=()
45 
46  awk '/chi2\[.\]/ { printf "%5d %9.5f\n", NR, $3 }' $WORKDIR/stage.txt > $WORKDIR/g1.txt
47 
48  JGraph \
49  -f $WORKDIR/g1.txt \
50  -o $WORKDIR/g1.root \
51  -T $TITLE
52 
53  OUTPUT_FILES+=$WORKDIR/g1.root
54 
55  awk \
56  -v YMIN=$YMIN \
57  -v YMAX=$YMAX \
58  '/gradient/ { printf "%5d %9.5f %5d %9.5f\n", NR, YMIN, NR, YMAX }' $WORKDIR/stage.txt > $WORKDIR/L1.txt
59 
60  JLine \
61  -p $WORKDIR/L1.txt \
62  -o $WORKDIR/L1.root \
63  -@ "width = 1" \
64  -@ "style = 2"
65 
66  OUTPUT_FILES+=$WORKDIR/L1.root
67 
68  JPlot1D \
69  -w 1200x600 \
70  -f${^OUTPUT_FILES}":.*" \
71  -\> "step" \
72  -\^ "#chi^{2}/NDF" \
73  -y "$YMIN $YMAX" \
74  -N "X 505" \
75  -G Y \
76  -S 0.5 \
77  -T "${${INPUT_FILE:t}:r}" \
78  -o ${${INPUT_FILE:t}:r}.$FORMAT $BATCH
79 
80  rm -f $WORKDIR/stage.txt $WORKDIR/g1.root $WORKDIR/L1.root
81 done
data_type w[N+1][M+1]
Definition: JPolint.hh:867
$WORKDIR stage
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
usr share Modules init zsh export TMPDIR
exit
Definition: JPizza.sh:36
data_type r[M+1]
Definition: JPolint.hh:868
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
Definition: datalogs.sh:31
fatal Invalid number of arguments esac let YMIN
Definition: testToT.sh:39
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DETECTOR_ID set_array RUNS if do_usage *then usage $script< source directory > nThe source directory corresponds to a pre calibrated D0ARCA009 detector fi case set_variable DIR
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
then rm
Definition: sftpput.zsh:30
do set_variable STRING_TXT awk
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then usage set_variable JSON_FILE $argv[-1] JEditDetector a $DETECTOR P set HIGH_RATE_VETO_DISABLE o $WORKDIR detector datx d set_variable DETECTOR $WORKDIR detector datx OUTPUT_FILES
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
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
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
set_variable INPUT_FILE
Auxiliary data structure for line.
Definition: JLine.cc:18
do grep(chi2\[.\]\|gradient\)" $INPUT_FILE > $WORKDIR/stage.txt awk 'BEGIN
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Definition: JMuonPostfit.sh:40
fi set_variable FILENAME $WORKDIR
Definition: datalogs.sh:34
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:48
data_type v[N+1][M+1]
Definition: JPolint.hh:866
then echo
Definition: JQAQC.sh:96
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25