JORCAShowerPositionPDF.sh

Go to the documentation of this file.

#!/bin/zsh                                                                                    
#                                                                                              
# \author adomi                                                                                
#                                                                                              

version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#                                                                                
#    Utility script to create a ROOT TH2D PDF for E.M. Shower Position Reconstruction.
#                                                                 
# ------------------------------------------------------------------------------------------ 

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

if ( do_usage $* ); then
    usage "$script [detector file [input file [output file]]]"
fi

SUB=${${RUN}:0:1}

DetectorFile="/afs/in2p3.fr/throng/km3net/detectors/orca_115strings_av23min20mhorizontal_18OMs_alt9mvertical_v1.detx"

set_variable  WORKDIR           $JPP_DATA
set_variable  FINAL_HISTOGRAM   $WORKDIR/PDFShowerPosition.root
set_variable  DETECTOR          $DetectorFile
set_variable  DEBUG             1

# DOWNLOAD FILES FROM IRODS AND EXTRACT THEM

MCFilesIRODS=/in2p3/km3net/mc/atm_neutrino/KM3NeT_ORCA_115_23m_9m/v1.1.1/JTE/elec-CC_3-100GeV_1-50.root.tar.gz 

MCFiles=elec-CC_3-100GeV_1-50.root.tar.gz

iget $MCFilesIRODS $WORKDIR

tar -xvzf $WORKDIR/$MCFiles -C $WORKDIR

FILENAME=JTE_elec-CC_3-100GeV_

for file in $(/bin/ls $WORKDIR/$FILENAME*)
do

nrun=$((nrun + 1))

echo $nrun

set_variable  INPUT_FILE        $file
set_variable  OUTPUT_HISTOGRAM  $WORKDIR/shower_position_histogram_$nrun.root

case $# in
    3) OUTPUT_HISTOGRAM=$3;&
    2) INPUT_FILE=$2;&
    1) DETECTOR=$1;;
esac

print_variable     DETECTOR INPUT_FILE OUTPUT_HISTOGRAM
check_input_file   $DETECTOR $INPUT_FILE 
check_output_file  $OUTPUT_HISTOGRAM

if (( 1 )); then

#   default input values                                                         
    set_variable T_RES_NS            "-100 100"
    set_variable D_M                     "0 80"
    set_variable COS_T                   "-1 1"

    timer_start

    JORCAShowerPositionPDF      \
   -a $DETECTOR            \
        -f $INPUT_FILE          \
        -o $OUTPUT_HISTOGRAM    \
        -T ${T_RES_NS}          \
        -D ${D_M}               \
        -C ${COS_T}             \
        -d ${DEBUG} --!

    timer_stop
    timer_print

fi

done

rm $FINAL_HISTOGRAM

hadd $WORKDIR/$FINAL_HISTOGRAM $WORKDIR/shower_position_histogram_*.root

rm $WORKDIR/shower_position_histogram_*.root

# CONVERT HISTOGRAM INTO A PDF (NORMALISE IT)

set_variable  OPTION        "NXYW"    

$JPP_DIR/software/JGizmo/JConvertToPDF2D     \
    -O $OPTION                               \
    -f $FINAL_HISTOGRAM":hPDF2Dist"          \
    -o $OUTPUT_PDF

rm $WORKDIR/$FILENAME*