Jpp  18.3.0-209-g56ce19a
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JParramatta.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 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13 set_variable: ZMAX_M ACOUSTICS_ZMAX_M 0.0
14 set_variable: TIMESTAMP GRAPHICS_TIMESTAMP utc
15 set_variable: FORMAT GRAPHICS_FORMAT gif
16 set_variable: SIZE GRAPHICS_SIZE 0.5
17 set_variable+ BATCH GRAPHICS_BATCH -B
18 
19 if do_usage $*; then
20  usage "$script (input file)+"\
21  "\nInput files correspond to the output of JKatoomba[.sh]."
22 fi
23 
24 if (( $# < 1 )); then
25  fatal "Wrong number of arguments."
26 fi
27 
29 
30 JParramatta \
31  -f "$INPUT_FILE[*]" \
32  -o $WORKDIR/parramatta.root \
33  -Z $ZMAX_M \
34  -d $DEBUG --!
35 
36 echo "Creating graphics output in directory $WORKDIR"
37 
38 JPlot1D \
39  -w 1200x600 \
40  -f $WORKDIR/parramatta.root:g1 \
41  -\^ "#chi^{2}/NDF" \
42  -N "X 505" \
43  -t "$TIMESTAMP" \
44  -S $SIZE \
45  -T "" \
46  -o $WORKDIR/chi2-time.$FORMAT $BATCH
47 
48 JPlot1D \
49  -w 1200x600 \
50  -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].amplitude" \
51  -\^ "amplitude [mrad]" \
52  -N "X 505" \
53  -t "$TIMESTAMP" \
54  -S $SIZE \
55  -T "" \
56  -L TR \
57  -o $WORKDIR/amplitude.$FORMAT $BATCH
58 
59 JPlot1D \
60  -w 1200x600 \
61  -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].stretching" \
62  -\^ "stretching [%]" \
63  -N "X 505" \
64  -t "$TIMESTAMP" \
65  -S $SIZE \
66  -T "" \
67  -L BR \
68  -o $WORKDIR/stretching.$FORMAT $BATCH
69 
70 typeset -Z4 STRING
71 
72 for H1 in `JPrintName -f "$WORKDIR/parramatta.root:stretching\[[0-9][0-9]*\]"`; do
73 
74  set_variable STRING `echo $H1 | sed 's/.*\[\([0-9]*\)\].*/\1/'`
75 
76  JPlot2D \
77  -f "$WORKDIR/parramatta.root:^\[${STRING}\]\.tiltdeviation" \
78  -\< "Ty - <Ty> [mrad]" \
79  -\> "Tx - <Tx> [mrad]" \
80  -O COLZ \
81  -N "X 505" \
82  -N "Y 505" \
83  -s 101110 \
84  -T "[${STRING}]" \
85  -o $WORKDIR/tiltdeviation_${STRING}.$FORMAT $BATCH
86 
87  JPrintResult \
88  -f "$WORKDIR/parramatta.root:^\[${STRING}\]\.tiltdeviation" \
89  -F "GetMean(1)" \
90  -F "GetMean(2)" \
91  -F "GetStdDev(1)" \
92  -F "GetStdDev(2)" | read -A RESULT
93 
94  printf "%04d %6.3f %6.3f %6.3f %6.3f\n" $STRING $RESULT[*]
95 done
96 
97 if (( 0 )); then
98 
99  set_variable FORMULA "[0] + [1]*x + [2]*x*x"
100 
101  for H1 in `JPrintName -f "$WORKDIR/parramatta.root:stretching\[[0-9][0-9]*\]"`; do
102 
103  set_variable STRING `echo $H1 | sed 's/.*\[\([0-9]*\)\].*/\1/'`
104 
105  H1=${H1/\[/\\\[}
106  H1=${H1/\]/\\\]}
107 
108  JFit \
109  -f "$WORKDIR/parramatta.root:${H1}" \
110  -o ${TMPDIR:-/tmp}/fit.root \
111  -F "$FORMULA" \
112  -@ "p0 = 0.0" \
113  -@ "p1 = 1.0e-5" \
114  -@ "p2 = 1.0e-6" \
115  -d $DEBUG
116 
117  JPlot1D \
118  -f ${TMPDIR:-/tmp}/fit.root:\.\* \
119  -\> "|T| [mrad]" \
120  -\^ "stretching [%]" \
121  -y "-0.25 +0.25" \
122  -N "X 505" \
123  -T "$STRING" \
124  -o $WORKDIR/profile_${STRING}.$FORMAT $BATCH
125  done
126 fi
data_type w[N+1][M+1]
Definition: JPolint.hh:867
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
Definition: JDataMonitor.sh:24
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
exit
Definition: JPizza.sh:36
then
Definition: datalogs.sh:34
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
then fatal Invalid string $STRING
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
do set_variable OUTPUT_DIRECTORY $WORKDIR T
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
then awk F
* usage
set_variable INPUT_FILE
&set_variable FORMULA
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Definition: JMuonPostfit.sh:40
*fatal Invalid number of arguments esac set_variable DETID getDetector D $DETID O string set_variable FILENAME $WORKDIR
Definition: datalogs.sh:27
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
then fatal The output file must have the wildcard in the e g root fi 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
then if[[!-f $DETECTOR]] then JDetector sh $DETECTOR fi cat $WORKDIR trigger_parameters txt<< EOFtrigger3DMuon.enabled=1;trigger3DMuon.numberOfHits=5;trigger3DMuon.gridAngle_deg=1;ctMin=0.0;TMaxLocal_ns=15.0;EOF set_variable TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY INPUT_FILES=() for((i=1;$i<=$NUMBER_OF_RUNS;++i));do JSirene.sh $DETECTOR $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz $WORKDIR/sirene_ ${i}.root JTriggerEfficiency.sh $DETECTOR $DETECTOR $WORKDIR/sirene_ ${i}.root $WORKDIR/trigger_efficiency_ ${i}.root $WORKDIR/trigger_parameters.txt $JPP_DATA/PMT_parameters.txt INPUT_FILES+=($WORKDIR/trigger_efficiency_ ${i}.root) done for ANGLE_DEG in $ANGLES_DEG[*];do set_variable SIGMA_NS 3.0 set_variable OUTLIERS 3 set_variable OUTPUT_FILE $WORKDIR/matrix\[${ANGLE_DEG}\deg\].root $JPP_DIR/examples/JReconstruction-f"$INPUT_FILES[*]"-o $OUTPUT_FILE-S ${SIGMA_NS}-A ${ANGLE_DEG}-O ${OUTLIERS}-d ${DEBUG}--!fiif[[$OPTION=="plot"]];then if((0));then for H1 in h0 h1;do JPlot1D-f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"-y"1 2e3"-Y-L TR-T""-\^"number of events [a.u.]"-> o chi2
Definition: JMatrixNZ.sh:106
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
then echo
Definition: JQAQC.sh:96
static const char *const TIMESTAMP
Time stamp of earliest UTC time.
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
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