billabong-signal: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
 
X=(-7500.0 +7500.0)                 # generation space
Y=(-7500.0 +7500.0)
Z=(    0.0 +2500.0)
 
let "NUMBER_OF_EVENTS  =  10000"
 
if do_usage $*; then
    usage "$script <detector> <energy [GeV]>"
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR  $argv[1]
set_variable  E_GEV     $argv[2]
 
let "VOLUME =  ($X[2] - $X[1]) * ($Y[2] - $Y[1]) * ($Z[2] - $Z[1]) * 1.0e-9"
    
set_variable  DIR  ${E_GEV}_GeV
 
if [[ ! -d $DIR ]]; then
    mkdir -p $DIR
fi
 
set_variable  VOLUME_TXT      volume.txt
set_variable  OUTPUT_FILE     $DIR/signal.root
 
echo "Writing data in directory $DIR.\n"\
     "Writing results to $VOLUME_TXT."
 
$JPP_DIR/examples/JAcoustics/JSignalToAWriter \
    -a $DETECTOR                               \
    -E $E_GEV                                  \
    -o $OUTPUT_FILE                            \
    -R 1000                                    \
    -W $WAVEFORM                               \
    -T 100.0                                   \
    -n $NUMBER_OF_EVENTS                       \
    -X "$X[*]"                                 \
    -Y "$Y[*]"                                 \
    -Z "$Z[*]"                                 \
    -d $DEBUG --!                              
 
set_variable  INPUT_FILE      $OUTPUT_FILE
set_variable  OUTPUT_FILE     $DIR/signal_trigger.root
 
$JPP_DIR/examples/JAcoustics/billabong:run.sh  $DETECTOR  $INPUT_FILE  $OUTPUT_FILE
 
let "NUMBER_OF_TRIGGERS = $(JPrintTree -f $OUTPUT_FILE -@ "type = JACOUSTICS::JEvent" -k number_of_entries)"
 
V=($(($VOLUME * $NUMBER_OF_TRIGGERS / $NUMBER_OF_EVENTS))  $(($VOLUME * sqrt($NUMBER_OF_TRIGGERS) / $NUMBER_OF_EVENTS)))
 
flock $VOLUME_TXT  printf "%5.2f  %6.2f  %5.2f\n"  $((log10($E_GEV)))  $V[1]  $V[2]  >>  $VOLUME_TXT