Jpp  17.0.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
pre-calibration_00000100.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 # --------------------------------------------------------------------------------------------
5 #
6 # Steering script for the global-fit-of-global-fits procedure of the D0ORCA010 (100) detector.
7 #
8 # --------------------------------------------------------------------------------------------
9 
10 if [ -z $JPP_DIR ]; then
11  echo "Variable JPP_DIR undefined."
12  exit
13 fi
14 
15 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
16 
18 set_array RUNS 11342 11345 11347 11374 11376
19 
20 if do_usage $*; then
21  usage "$script <detector file> <tripod file> [user directory]"\
22  "\nThe user directory should be specified when this script is submitted in a batch queue."\
23  "\nThe detector file and tripod file correspond to a pre-calibrated D_ORCA006 detector."
24 fi
25 
28 
29 case $# in
30  3) set_variable HOMEDIR $3;&
33  *) fatal "Wrong number of arguments."
34 esac
35 
36 set_variable TRIPOD_INITIAL tripod_initial.txt
37 set_variable DETECTOR_INITIAL detector_initial.datx
38 
39 
40 # Run in HOMEDIR
41 
42 mkdir -p $HOMEDIR
43 
44 cd $HOMEDIR
45 
46 if [[ ! -f $DETECTOR_INITIAL ]]; then
47 
48  JDetectorDB \
49  -D $DETECTOR_ID \
50  -r $RUNS[1] \
51  -WW \
52  -o $DETECTOR_INITIAL
53 fi
54 
55 if [[ ! -f $TRIPOD_INITIAL ]]; then
56 
57  cat>$TRIPOD_INITIAL<<EOF
58 # Very preliminary acoustic emitter positions from survey
59 2 +256815.500 +4743395.000 -2436.677
60 3 +257096.200 +4743636.000 -2439.354
61 4 +256919.238 +4743354.178 -2439.354
62 EOF
63 fi
64 
65 # Apply modifications.
66 
67 JEditDetector \
68  -a $DETECTOR_INITIAL \
69  -o detector.datx
70 
71 eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`; STRINGS_INITIAL=($STRINGS[*])
72 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
73 
74 JEditDetector -a detector.datx -r "$STRINGS[*]" -o detector.datx
75 JMergeDetector -a detector.datx -a $DETECTOR -o detector.datx
76 JEditDetector -a detector.datx -k "$STRINGS_INITIAL[*]" -o detector.datx
77 
78 FIXED_STRINGS=($STRINGS[*])
79 
80 cp -p $TRIPOD_INITIAL tripod.txt
81 
82 source JAcousticsToolkit.sh
83 
84 typeset -A TRIPODS_INITIAL
85 typeset -A TRIPODS
86 
87 get_tripods $TRIPOD_INITIAL TRIPODS_INITIAL
88 get_tripods $TRIPOD TRIPODS
89 
90 JEditTripod -f tripod.txt -D "${(k)TRIPODS}"
91 JMergeTripod -f tripod.txt -f $TRIPOD -o tripod.txt
92 JEditTripod -f tripod.txt -K "${(k)TRIPODS_INITIAL}"
93 
94 FIXED_TRIPODS=(${(k)TRIPODS[*]})
95 
98 
99 JAcoustics.sh $DETECTOR_ID
100 
101 cat>acoustics_trigger_parameters.txt<<EOF
102 Q = 0.0;
103 TMax_s = 0.015;
104 numberOfHits = 150;
105 EOF
106 
107 JAcousticsEventBuilder.sh $DETECTOR $RUNS[*]
108 
109 INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
110 
111 
112 # Run in WORKDIR
113 
114 cd $WORKDIR
115 
116 if [ ! $HOMEDIR -ef $WORKDIR ]; then
117  cp -p $HOMEDIR/$DETECTOR $WORKDIR
118  cp -p $HOMEDIR/$TRIPOD $WORKDIR
119  cp -p $HOMEDIR/${^INPUT_FILES} $WORKDIR
120  cp -p $HOMEDIR/{acoustics_fit_parameters,acoustics_trigger_parameters,disable,hydrophone,mechanics,sound_velocity,tripod,waveform}.txt \ $WORKDIR
121 fi
122 
123 source $JPP_DIR/examples/JAcoustics/acoustics-fit-toolkit.sh
124 
126 
128 
129 fixStrings $FIXED_STRINGS[*]
130 fixTripods $FIXED_TRIPODS[*]
131 
132 stage_0 >& stage-0.log; timer_lap
133 stage_1B >& stage-1B.log; timer_lap
134 stage_2 >& stage-2a.log; timer_lap
135 
137 
138 stage_2 >& stage-2b.log; timer_lap
139 stage_3 >& stage-3.log; timer_lap
140 
141 if [ ! $HOMEDIR -ef $WORKDIR ]; then
142  cp -vru * $HOMEDIR
143 fi
144 
145 clean
146 
&set_variable TRIPOD
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Q(UTCMax_s-UTCMin_s)-livetime_s
static const uint32_t K[64]
Definition: crypt.cc:77
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null if do_usage *then usage $script< DEETECTOR_ID >< inputdetectorfile >< outputdetectorfile >< firstrun >< lastrun >< namesuffix >< outputdirectory >< datadirectory > fi DETECTOR_ID
Definition: SN_timecalib.sh:15
then JLigier sh continue fi cat
Definition: JDAQDriver.sh:51
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
Definition: JPizza.sh:36
function stage_0()
data_type r[M+1]
Definition: JPolint.hh:758
then usage source $script nThis script contains auxiliary functions for the pre calibration of specific detectors nIt should be sourced by the detector specific scripts fi source JAcousticsToolkit sh set_variable WORKDIR pwd set_variable TMPDIR $WORKDIR typeset A TRIPODS typeset a CHI2 set_variable PRECISION set_variable EPSILON set_variable RADIUS_M function initialise()
clean for KEY in pcal ccal tripod transmitter hydrophone
function fixTripods()
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null if do_usage *then usage source $script nToolkit for acoustics scripts fi function get_tripods()
then echo
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:68
function stage_1B()
then timer_start(JRandomTimesliceWriter-a ${DETECTOR}-o ${PIPE}-n ${NUMBER_OF_SLICES}-P ${PMT_FILE}-d ${DEBUG}-S ${SEED}-B"$BACKGROUND_HZ"-N"$RECYCLING"--!&&for((i=0;$i!=10;++i));do;echo 'EOF' > $PIPE;usleep 100000;done &&rm-f $PIPE)&JTimesliceReprocessor-a $
then fatal Detector file should be binary format fi eval JPrintDetector a $DETECTOR O IDENTIFIER typeset a INPUT_FILES for RANGE in $RUNS[*]
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
case $OPTION in clean clean
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 JConvertDetectorFormat a $DETECTOR[1] o
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:77
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DETECTOR_ID set_array RUNS set_variable MUL if do_usage *then usage $script< detector file >< tripodfile > nThe user directory should be specified when this script is submitted in a batch queue nThe detector file and tripod file correspond to a pre calibrated D0ARCA006 detector fi set_variable WORKDIR pwd set_variable HOMEDIR $WORKDIR case set_variable HOMEDIR
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
*fatal Wrong number of arguments esac eval JPrintDetector a $DETECTOR O IDENTIFIER RUNS
then error Missing CDF files error Run
Definition: JSirene.sh:58
then display $WORKDIR
Definition: plot-Domino.sh:128
function stage_2()
then cp
set_variable DETECTOR
do set_variable DETECTOR_TXT $WORKDIR detector
function fixStrings()
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:46
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
function stage_3()
script
Definition: JAcoustics.sh:2
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
then fatal Invalid tripod $TRIPOD
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
do JPrintDAQHeader f $INPUT_FILE d read DETECTOR_ID RUN FRAME_INDEX TIMESLICE_START set_variable OUTPUT_FILE $WORKDIR KM3NeT_
then usage $script< archive >< detectorfile >< run > nGet tripod