JMuonCompass.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:  A_RAD             RECONSTRUCTION_A_RAD                1.0
 
if do_usage $*; then
    usage "$script <detector file> <input file>"\
          "\nExample script to test compass 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`
 
if [[ ! -f detector_A.datx ]]; then
 
    JConvertDetectorFormat \
        -a $DETECTOR                          \
        -o detector_A.datx
 
    if (( 0 )); then
        for STRING in $STRINGS[*]; do
    
            JEditDetector \
                -a detector_A.datx                \
                -s "$STRING randrot $A_RAD"       \
                -o detector_A.datx -q
        done
    fi
    if (( 1 )); then
        for MODULE in $MODULES[*]; do
 
            getLocation -a detector_A.datx -M $MODULE -d 0  |  read  STRING FLOOR
 
            if (( $FLOOR != 0 )); then
    
                JEditDetector \
                    -a detector_A.datx                \
                    -M "$MODULE randrot $A_RAD"       \
                    -o detector_A.datx -q
            fi
        done
    fi
 
    JConvertDetectorFormat \
        -a detector_A.datx                    \
        -o detector_B.datx
fi
 
if [[ ! -f sirene.root ]]; then
    JSirene.sh                  $DETECTOR  $INPUT_FILE    sirene.root
fi
if [[ ! -f trigger_efficieny.root ]]; then
    JTriggerEfficiency.sh       $DETECTOR  detector_A.datx  sirene.root  trigger_efficieny.root  $JPP_DATA/trigger_parameters_arca.txt  $JPP_DATA/PMT_parameters.txt
fi
if [[ ! -f jmuon.root ]]; then
    JARCAMuonReconstruction.sh  detector_A.datx  trigger_efficieny.root  %.root  $JPP_DATA/PMT_parameters.txt
fi
 
JConvertDetectorFormat \
    -a detector_B.datx                    \
    -o detector_A.datx
 
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   1)
 
timer_start
 
eval JMuonCompass \
     -a detector_A.datx                   \
     -f jmuon.root                        \
     -o compass.root                      \
     -P $PDF                              \
     -@ `make_equation PARAMETERS`        \
     -d $DEBUG                            
 
timer_stop
timer_print
 
$JPP_DIR/examples/JReconstruction/JMuonCompass \
    -f compass.root                      \
    -o hobbit.root                       \
    -a detector_A.datx                   \
    -A                                   \
    -d $DEBUG
    
popd