JPerth.sh

Go to the documentation of this file.

#!/bin/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`        \
     -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`        \
     -A                                   \
     -N $THREADS                          \
     -d $DEBUG                            >& perth_B.log

popd