Jpp  15.0.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JNarrabri.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 DYNAMICS_DEBUG 2
12 set_variable: WORKDIR DYNAMICS_WORKDIR ./
13 set_variable: TMAX_S DYNAMICS_TMAX_S 1000.0
14 set_variable: YMIN DYNAMICS_YMIN 0.00
15 set_variable: YMAX DYNAMICS_YMAX 0.02
16 set_variable DIR $JPP_DIR/examples/JDynamics/
17 set_variable: TIMESTAMP GRAPHICS_TIMESTAMP utc
18 set_variable: FORMAT GRAPHICS_FORMAT gif
19 set_variable+ BATCH GRAPHICS_BATCH -B
20 
21 if do_usage $*; then
22  usage "$script <detector file> (input file)+"\
23  "\nInput files correspond to the output of JKatoomba[.sh]."
24 fi
25 
26 if (( $# < 2 )); then
27  fatal "Wrong number of arguments."
28 fi
29 
30 set_variable DETECTOR $argv[1]
31 set_variable INPUT_FILE $argv[2,-1]
32 
33 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
34 
35 
36 $DIR/JNarrabri \
37  -a $DETECTOR \
38  -f "$INPUT_FILE[*]" \
39  -o $WORKDIR/narrabri.root \
40  -T $TMAX_S \
41  -d $DEBUG --!
42 
43 
44 typeset -Z 4 STRING
45 
46 for STRING in $STRINGS[*]; do
47 
48  JPrintResult \
49  -f "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt" \
50  -F "GetStdDev(1)" \
51  -F "GetStdDev(2)" \
52  -F "GetEntries - Integral" \
53  -F "GetEntries" | read -A RESULT
54 
55  printf "%04d %4.2f %4.2f %5.0f / %8.0f\n" $STRING $RESULT[*]
56 
57 done
58 
59 echo "Creating graphics output in directory $WORKDIR"
60 
61 JPlot1D \
62  -f $WORKDIR/narrabri.root:time \
63  -\> "#Deltat [s]" \
64  -\^ "number of events [a.u.]" \
65  -XX \
66  -L TL \
67  -T "" \
68  -o $WORKDIR/time.$FORMAT $BATCH
69 
70 JPlot1D \
71  -w 1200x600 \
72  -f $WORKDIR/narrabri.root:G\.\*orientation \
73  -\^ "orientation [rad]" \
74  -y "-3.20 +3.20" \
75  -N "X 505" \
76  -t "$TIMESTAMP" \
77  -L BR \
78  -T "" \
79  -o $WORKDIR/orientation.$FORMAT $BATCH
80 
81 JPlot1D \
82  -w 1200x600 \
83  -f $WORKDIR/narrabri.root:G\.\*amplitude \
84  -\^ "amplitude [rad]" \
85  -N "X 505" \
86  -t "$TIMESTAMP" \
87  -L TR \
88  -T "" \
89  -o $WORKDIR/amplitude.$FORMAT $BATCH
90 
91 JPlot1D \
92  -w 1200x600 \
93  -f $WORKDIR/narrabri.root:H\.\*orientation \
94  -\^ "orientation [rad]" \
95  -y "-3.20 +3.20" \
96  -N "X 505" \
97  -t "$TIMESTAMP" \
98  -L BR \
99  -T "" \
100  -o $WORKDIR/orientation-L.$FORMAT $BATCH
101 
102 JPlot1D \
103  -w 1200x600 \
104  -f $WORKDIR/narrabri.root:H\.\*amplitude \
105  -\^ "amplitude [rad]" \
106  -N "X 505" \
107  -t "$TIMESTAMP" \
108  -L TR \
109  -T "" \
110  -o $WORKDIR/amplitude-L.$FORMAT $BATCH
111 
112 for STRING in $STRINGS[*]; do
113 
114  JPlot1D \
115  -w 1200x600 \
116  -f "$WORKDIR/narrabri.root:G\[${STRING}\]\.orientation" \
117  -f "$WORKDIR/narrabri.root:H\[${STRING}\]\.orientation" \
118  -\^ "orientation [rad]" \
119  -y "-3.20 +3.20" \
120  -N "X 505" \
121  -t "$TIMESTAMP" \
122  -T "[${STRING}]" \
123  -o $WORKDIR/orientation_${STRING}.$FORMAT $BATCH
124 
125  JPlot1D \
126  -w 1200x600 \
127  -f "$WORKDIR/narrabri.root:G\[${STRING}\]\.amplitude" \
128  -f "$WORKDIR/narrabri.root:H\[${STRING}\]\.amplitude" \
129  -\^ "amplitude [rad]" \
130  -y "$YMIN $YMAX" \
131  -N "X 505" \
132  -t "$TIMESTAMP" \
133  -T "[${STRING}]" \
134  -o $WORKDIR/amplitude_${STRING}.$FORMAT $BATCH
135 
136  JPlot2D \
137  -f "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt" \
138  -\> "dTx/dt [mrad/10 min]" \
139  -\< "dTy/dt [mrad/10 min]" \
140  -O COLZ \
141  -N "X 505" \
142  -N "Y 505" \
143  -s 101010 \
144  -T "[${STRING}]" \
145  -o $WORKDIR/tilt_${STRING}.$FORMAT $BATCH
146 done
$argv
then usage $script< detector file >< detectorfile > nIf the range of floors is the first detector file is aligned to the second before the comparison nIn only modules within the given range of floors are used fi case set_variable RANGE $argv[3]
Definition: JCompareDetector.sh:31
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:741
JCOMPAREHISTOGRAMS::read
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
Definition: JTestSummary.hh:222
STRING
clean eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY set_variable STRING
Definition: examples/JTrigger/JTriggerEfficiencyRunByRun.sh:64
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
N
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
Definition: JShowerPostfit.sh:95
JPHYSICS::H
static const double H
Planck constant [eV s].
Definition: JPhysics/JConstants.hh:28
exit
exit
Definition: JPizza.sh:36
DIR
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:54
YMIN
fatal Invalid number of arguments esac let YMIN
Definition: testToT.sh:39
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
$STRING
then fatal Invalid string $STRING
Definition: detector-XY:run.sh:35
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
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
F
then awk F
Definition: getMIIseacurrent.sh:86
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
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
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
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
JGIZMO::TIMESTAMP
static const char *const TIMESTAMP
Time stamp of earliest UTC time.
Definition: JGizmoToolkit.hh:85
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
Generated by   doxygen 1.8.5