billabong-random:run.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0:t}
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
zmodload zsh/mathfunc
 
set_variable   DEBUG           2
set_variable   WAVEFORM        1    # waveform
 
if do_usage $*; then
    usage "$script <detector> <rate [Hz]>"
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR  $argv[1]
set_variable  RATE_HZ   $argv[2]
 
set_variable  T_S             $((100 * (5.0/$RATE_HZ)**3))
set_variable  DIR             ${RATE_HZ}_Hz
 
if [[ ! -d $DIR ]]; then
    mkdir -p $DIR
fi
 
set_variable  NOISE_TXT       noise.txt
set_variable  CPU_TXT         cpu.txt
set_variable  OUTPUT_FILE     $DIR/random.root
 
echo "Writing data in directory $DIR.\n"\
     "Writing results to $NOISE_TXT and $CPU_TXT."
 
$JPP_DIR/examples/JAcoustics/JRandomToAWriter \
    -a $DETECTOR                               \
    -o $OUTPUT_FILE                            \
    -B $RATE_HZ                                \
    -R 1000                                    \
    -W $WAVEFORM                               \
    -T $T_S                                    \
    -d $DEBUG --!                              
 
set_variable  INPUT_FILE      $OUTPUT_FILE
set_variable  OUTPUT_FILE     $DIR/random_trigger.root
 
timer_start
    
$JPP_DIR/examples/JAcoustics/billabong:run.sh  $DETECTOR  $INPUT_FILE  $OUTPUT_FILE
 
timer_stop
    
let "NUMBER_OF_TRIGGERS = $(JPrintTree -f $OUTPUT_FILE -@ "type = JACOUSTICS::JEvent" -k number_of_entries)"
let "ELAPSED_TIME       = $(timer_print | awk '/Elapsed time/ { print $3}')"
 
R=($(($NUMBER_OF_TRIGGERS / $T_S))  $((sqrt($NUMBER_OF_TRIGGERS + 1.0) / $T_S)))
U=($(($ELAPSED_TIME       / $T_S))  $((1.0                             / $T_S)))
 
flock $NOISE_TXT  printf "%5.2f  %12.5f  %12.5f\n"  $RATE_HZ  $R[1]  $R[2]  >>  $NOISE_TXT
flock $CPU_TXT    printf "%5.2f  %12.5f  %12.5f\n"  $RATE_HZ  $U[1]  $U[2]  >>  $CPU_TXT