Jpp  18.6.0-rc.1
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 
28 set_variable INPUT_FILE $argv[1,-1]
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 JPlot1D \
71  -f "$WORKDIR/parramatta.root:hx" \
72  -f "$WORKDIR/parramatta.root:hy" \
73  -\^ "#DeltaTx/y [mrad]" \
74  -\> "string" \
75  -T "" \
76  -o $WORKDIR/txy.$FORMAT $BATCH
77 
78 typeset -Z4 STRING
79 
80 for H1 in `JPrintName -f "$WORKDIR/parramatta.root:^\[[0-9][0-9]*\]\.tiltdeviation"`; do
81 
82  set_variable STRING `echo $H1 | sed 's/.*\[\([0-9]*\)\].*/\1/'`
83 
84  JPlot2D \
85  -f "$WORKDIR/parramatta.root:^\[${STRING}\]\.tiltdeviation" \
86  -\< "Ty - <Ty> [mrad]" \
87  -\> "Tx - <Tx> [mrad]" \
88  -O COLZ \
89  -G X \
90  -G Y \
91  -N "X 505" \
92  -N "Y 505" \
93  -s 101110 \
94  -T "[${STRING}]" \
95  -o $WORKDIR/tiltdeviation_${STRING}.$FORMAT $BATCH
96 done
97 
98 if (( 0 )); then
99 
100  set_variable FORMULA "[0] + [1]*x + [2]*x*x"
101 
102  for H1 in `JPrintName -f "$WORKDIR/parramatta.root:stretching\[[0-9][0-9]*\]"`; do
103 
104  set_variable STRING `echo $H1 | sed 's/.*\[\([0-9]*\)\].*/\1/'`
105 
106  H1=${H1/\[/\\\[}
107  H1=${H1/\]/\\\]}
108 
109  JFit \
110  -f "$WORKDIR/parramatta.root:${H1}" \
111  -o ${TMPDIR:-/tmp}/fit.root \
112  -F "$FORMULA" \
113  -@ "p0 = 0.0" \
114  -@ "p1 = 1.0e-5" \
115  -@ "p2 = 1.0e-6" \
116  -d $DEBUG
117 
118  JPlot1D \
119  -f ${TMPDIR:-/tmp}/fit.root:\.\* \
120  -\> "|T| [mrad]" \
121  -\^ "stretching [%]" \
122  -y "-0.25 +0.25" \
123  -N "X 505" \
124  -T "$STRING" \
125  -o $WORKDIR/profile_${STRING}.$FORMAT $BATCH
126  done
127 fi
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:867
$argv
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
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!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
exit
exit
Definition: JPizza.sh:36
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
then
then
Definition: datalogs.sh:34
o
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
Definition: software/JDetector/compare-detector.sh:31
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
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: archive-get-calibration.sh:14
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
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
FORMULA
&set_variable FORMULA
Definition: JEvtReweightMupageParameterScan.sh:60
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonStart.sh:47
N
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Definition: JMuonPostfit.sh:40
in
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
chi2
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
X
no fit printf nominal n $STRING awk v X
Definition: detector-XY:fit.sh:102
Z
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
Definition: JAlignDetector.sh:81
echo
then echo
Definition: JQAQC.sh:92
JGIZMO::TIMESTAMP
static const char *const TIMESTAMP
Time stamp of earliest UTC time.
Definition: JGizmoToolkit.hh:95
$WORKDIR
then getFile sh $JPP_ARCHIVE $DETECTOR_ID $RUNS[1] $KEY $VERSION $WORKDIR
Definition: JAcousticsTriggerProcessor.sh:46
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
done
esac done
Definition: JAddHDE.sh:21
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
g1
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25
Generated by   doxygen 1.8.5