Jpp  16.0.0-rc.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
acoustics-fit_00000049.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 D_ORCA006 (49) 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 7600 7700 7800 7900 8000 8100 8200
19 set_variable: MUL ACOUSTICS_MUL +0.0e-2
20 set_variable: STAGES ACOUSTICS_STAGES 1 2 3
21 
22 if do_usage $*; then
23  usage "$script [user directory]"\
24  "\nThe user directory should be specified when this script is submitted in a batch queue."
25 fi
26 
28 
29 if (( $# == 1 )); then
31 else
33 fi
34 
35 set_variable TRIPOD_INITIAL tripod_initial.txt
36 set_variable TRIPOD tripod.txt
37 set_variable DETECTOR_INITIAL detector_initial.datx
39 
40 
41 # Run in HOMEDIR
42 
43 mkdir -p $HOMEDIR
44 
45 cd $HOMEDIR
46 
47 if [[ ! -f $DETECTOR_INITIAL ]]; then
48 
49  JDetectorDB \
50  -D $DETECTOR_ID \
51  -r $RUNS[1] \
52  -WW \
53  -o $DETECTOR_INITIAL
54 fi
55 
56 if [[ ! -f $TRIPOD_INITIAL ]]; then
57 
58  cat>$TRIPOD_INITIAL<<EOF
59 # Data are from V. Bertin, private communications.
60 1 256877.5 4743716.7 -2438.4
61 2 256815.5 4743395.0 -2435.5
62 3 257096.2 4743636.0 -2439.5
63 EOF
64 fi
65 
66 # Apply modifications to the detector, if any.
67 
68 JEditDetector \
69  -a $DETECTOR_INITIAL \
70  -@ "z = -2441.5" \
71  -s "-1 mul $MUL" \
72  -o $DETECTOR
73 
74 cp -p $TRIPOD_INITIAL $TRIPOD
75 
76 JAcoustics.sh $DETECTOR_ID
77 
78 JAcousticsEventBuilder.sh $DETECTOR $RUNS[*]
79 
80 
81 # Run in WORKDIR
82 
83 cd $WORKDIR
84 
85 if [ ! $HOMEDIR -ef $WORKDIR ]; then
86 
87  cp -p $HOMEDIR/$TRIPOD $WORKDIR
88  cp -p $HOMEDIR/$DETECTOR $WORKDIR
89 
90  for RANGE in $RUNS[*]; do
91  for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
92  cp -p $HOMEDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root $WORKDIR
93  done
94  done
95 
96  JAcoustics.sh $DETECTOR_ID
97 
98 fi
99 
100 for STAGE in `echo $STAGES`; do
101 
102  $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*event.root > $HOMEDIR/stage-${STAGE}.log
103 
104  if [ ! $HOMEDIR -ef $WORKDIR ]; then
105  cp -vru * $HOMEDIR
106  fi
107 done
108 
set_array RUNS
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
then usage $script[< detector file >< tripod file >][user directory] nThe user directory should be specified when this script is submitted in a batch queue nThe optional detector file and tripod file correspond to a pre calibrated sub detector fi set_variable WORKDIR pwd set_variable HOMEDIR $WORKDIR case set_variable HOMEDIR
data_type r[M+1]
Definition: JPolint.hh:758
then echo
*fatal invalid stage $STAGE(possible stages:1, 2, 3, A, D)"
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null if do_usage *then usage $script< detector identifier > fi case set_variable DETECTOR_ID
Definition: JTDC.sh:11
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:20
then rm i $OUTPUT_FILE fi let RUN
&set_variable TRIPOD
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
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
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
then display $WORKDIR
Definition: plot-Domino.sh:128
then cp
set_variable DETECTOR
do set_variable INPUT_FILE $WORKDIR KM3NeT_
do set_variable DETECTOR_TXT $WORKDIR detector
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
script
Definition: JAcoustics.sh:2
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
then cat $TRIPOD_INITIAL<< EOF1 256877.5 4743716.7-2438.42 256815.5 4743395.0-2435.53 257096.2 4743636.0-2439.5EOFfiJEditDetector-a $DETECTOR_INITIAL-@"z = -2441.5"-s"-1 mul $MUL"-o $DETECTORcp-p $TRIPOD_INITIAL $TRIPODJAcoustics.sh $DETECTOR_IDJAcousticsEventBuilder.sh $DETECTOR $RUNS[*]cd $WORKDIRif[!$HOMEDIR-ef $WORKDIR];then cp-p $HOMEDIR/$TRIPOD $WORKDIR cp-p $HOMEDIR/$DETECTOR $WORKDIR for RANGE in $RUNS[*];do for((RUN=${RANGE%%-*};$RUN<=${RANGE##*-};RUN+=1));do cp-p $HOMEDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ ${(l:8::0::0:) RUN}_event.root $WORKDIR done done JAcoustics.sh $DETECTOR_IDfifor STAGE in`echo $STAGES`;do $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ *event.root > $HOMEDIR stage
then fatal Invalid tripod $TRIPOD
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null