JPerth.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable:  DEBUG             RECONSTRUCTION_DEBUG                2
set_variable:  WORKDIR           RECONSTRUCTION_WORKDIR              ./
set_variable:  PDF               RECONSTRUCTION_PDF                  $JPP_DATA/J%p.dat
set_variable:  THREADS           RECONSTRUCTION_THREADS              5
set_variable:  T_NS              RECONSTRUCTION_T_NS                 10.0
 
if do_usage $*; then
    usage "$script <detector file> <input file>"\
              "\nExample script to test inter-string calibration."\
              "\nThe input file corresponds to an input file for the detector simulation (e.g. JSirene)."
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR        $argv[1]
set_variable INPUT_FILE      $argv[2]
 
mkdir -p $WORKDIR
pushd    $WORKDIR
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
JConvertDetectorFormat \
    -a $DETECTOR                          \
    -o detector.datx
 
for STRING in $STRINGS[*]; do
    
    JEditDetector \
            -a detector.datx                  \
            -s "$STRING randadd $T_NS"        \
            -o detector.datx
done
 
JConvertDetectorFormat \
    -a detector.datx                      \
    -o detector_A.datx
 
JSirene.sh                  $DETECTOR  $INPUT_FILE    sirene.root
JTriggerEfficiency.sh       $DETECTOR  detector.datx  sirene.root  trigger_efficieny.root  $JPP_DATA/trigger_parameters_arca.txt  $JPP_DATA/PMT_parameters.txt
 
JARCAMuonReconstruction.sh  detector.datx  trigger_efficieny.root  jmuon.root  $JPP_DATA/PMT_parameters.txt
 
typeset -A PARAMETERS
 
PARAMETERS=(
    TTS_ns            2
    E_GeV             1.0e3
    TMin_ns            -50.0
    TMax_ns           +450.0
    ZMin_m            -1.0e4
    ZMax_m            +1.0e4
    R_Hz              6.0e3
    roadWidth_m       175.0
    numberOfPrefits   0)
 
eval JPerth \
     -a detector.datx                       \
     -f jmuon.root                          \
     -F $PDF                                \
     -@ `make_equation PARAMETERS`          \
     -% \"number_of_iterations  = 10000\"   \
     -% \"number_of_extra_steps = 5\"       \
     -% \"T_ns                  = 1.0\"     \
     -S -1                                  \
     -A                                     \
     -N $THREADS                            \
     -d $DEBUG                              >& perth_A.log
 
JConvertDetectorFormat \
    -a detector.datx                        \
    -o detector_B.datx
 
rm -f jmuon.root
 
JARCAMuonReconstruction.sh  detector.datx  trigger_efficieny.root  jmuon.root  $JPP_DATA/PMT_parameters.txt
 
eval JPerth \
     -a detector.datx                       \
     -f jmuon.root                          \
     -F $PDF                                \
     -@ `make_equation PARAMETERS`          \
     -% \"number_of_iterations  = 10000\"   \
     -% \"number_of_extra_steps = 5\"       \
     -% \"T_ns                  = 1.0\"     \
     -S -1                                  \
     -A                                     \
     -N $THREADS                            \
     -d $DEBUG                              >& perth_B.log
 
popd