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  %.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
    R_Hz              6.0e3
    roadWidth_m       175.0
    numberOfPrefits   0)
 
eval JPerth \
     -a detector.datx                     \
     -f jmuon.root                        \
     -P $PDF                              \
     -@ `make_equation PARAMETERS`        \
     -% "number_of_extra_steps = 25"      \
     -% "T_ns                  = 2.5"     \
     -A                                   \
     -N $THREADS                          \
     -d $DEBUG                            >& perth_A.log
 
JConvertDetectorFormat \
    -a detector.datx                      \
    -o detector_B.datx
 
rm -f jmuon.root jmuon.offline.root
 
JARCAMuonReconstruction.sh  detector.datx  trigger_efficieny.root  %.root  $JPP_DATA/PMT_parameters.txt
 
eval JPerth \
     -a detector.datx                     \
     -f jmuon.root                        \
     -P $PDF                              \
     -@ `make_equation PARAMETERS`        \
     -% "number_of_extra_steps = 25"      \
     -% "T_ns                  = 2.5"     \
     -A                                   \
     -N $THREADS                          \
     -d $DEBUG                            >& perth_B.log
 
popd