JMakePDF.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
# utility script to produce PDF and CDF tables.
script=${0##*/}


# ------------------------------------------------------------------------------------------
#
#                 Utility script to produce PDF and CDF tables.
#
# ------------------------------------------------------------------------------------------

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

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable   DEBUG                   2          
set_variable   WORKDIR                 $JPP_DATA  
set_variable   MAKE_PDF                0          
set_variable   MERGE_PDF               0          
set_variable   MAKE_CDF                0          
set_variable   CLEAN                   0          
set_variable   BACKGROUND                         
set_variable   WGET                               
set_variable   RMAX_M                  999        
set_variable   absorptionLengthFactor  1.0        
set_variable   scatteringLengthFactor  1.0        
set_variable:  TTS_NS                  PHYSICS_TTS_NS     2.0
set_variable   EXECUTOR                "nohup"

#
# Method to print usage of script.
#
function __usage__()
{
    set_variable  DEBUG  3

    usage "$script\n"\
        "-d <debug level>                                                                      \n"\
        "-W <working directory PDF and CDF tables>                                             \n"\
        "-P create  PDF tables                                                                 \n"\
        "-M merge   PDF tables                                                                 \n"\
        "-C convert PDF tables to CDF tables                                                   \n"\
        "-c clean   PDF tables                                                                 \n"\
        "-B parallel running of jobs in background                                             \n"\
        "-E execute command to be evaluated (standard: nohup, example: echo, using -B as well) \n"\
        "-w download files from GIT server (e.g. \"http://sftp.km3net.de/data/latest/\")       \n"\
        "-R <maximal distance between muon and PMT>                                            \n"\
        "Note that the elongated EM shower PDF is created at option -M.                        \n"\
        "Script to generate the PDF/CDFs and merge them;                                       \n"\
        "standard example: JMakePDF.sh -W \${JPP_DATA} -PMCc                                   \n"\
        "print all necessary programs: JMakePDF.sh -W \${JPP_DATA} -E \"echo\" -PMCB"
}

#
# Parse command line options.
#
while getopts "d:W:PMCcBE:w:hA:S:R:" option; do
    case "$option" in
        d) set_variable DEBUG       $OPTARG;;
        W) set_variable WORKDIR     $OPTARG;;
        P) set_variable MAKE_PDF    1;;
        M) set_variable MERGE_PDF   1;;
        C) set_variable MAKE_CDF    1;;
        c) set_variable CLEAN       1;;
        B) set_variable BACKGROUND  \&;;
        E) set_variable EXECUTOR    $OPTARG;;
        w) set_variable WGET        $OPTARG;;
        h) __usage__;;
        A) set_variable absorptionLengthFactor  $OPTARG;;
        S) set_variable scatteringLengthFactor  $OPTARG;;
        R) set_variable RMAX_M      $OPTARG;;
    esac
done

# check whether script is already running

BUFFER=`get_pid $script`
BUFFER=${BUFFER/${$}/}

if [[ -n "${BUFFER/ */}" ]]; then
    notice "${script} already running <${BUFFER}>"
    exit 0
fi

check_variable ROOTSYS JPP_DIR WORKDIR

mkdir -p $WORKDIR

if [[ ! -d $WORKDIR ]]; then
    fatal "Directory $WORKDIR does not exist."
fi

status "Start $script"


if [[ -n $WGET ]]; then
    cd $WORKDIR
    (rm -f *.dat) >& /dev/null
    wget -A dat,txt,root,evt.gz -r -l 1 -nd -q $WGET 
fi


if (( $MAKE_PDF )); then

    check_process JMakePDF JMakePDG pgrep

    if (( 1 )); then

        if [[ ! -f $WORKDIR/J1p.dat ]]; then
            eval $EXECUTOR JMakePDF -F1 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J1p.dat -d $DEBUG $BACKGROUND
        fi
        if [[ ! -f $WORKDIR/J2p.dat ]]; then
            eval $EXECUTOR JMakePDF -F2 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J2p.dat -d $DEBUG $BACKGROUND
        fi

        for R in 0.1 0.3 0.5 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 \
            11.0 12.0 13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 \
            22.0 24.0 26.0 28.0 30.0 32.0 34.0 36.0 38.0 40.0 42.0 44.0 46.0 48.0 50.0 \
            55.0 60.0 65.0 70.0 75.0 80.0 85.0 90.0 95.0 100.0 110.0 120.0 130.0 140.0 150.0 \
            170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0; do
            
            if (( $R <= $RMAX_M )); then
                    
                if [[ ! -f $WORKDIR/J3p.dat ]]; then
                    eval $EXECUTOR JMakePDF -F3 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J3p-$R.dat -d $DEBUG $BACKGROUND
                fi
                if [[ ! -f $WORKDIR/J4p.dat ]]; then
                    eval $EXECUTOR JMakePDF -F4 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J4p-$R.dat -d $DEBUG $BACKGROUND
                fi
            fi
        done
        
        if [[ ! -f $WORKDIR/J5p.dat ]]; then
            eval $EXECUTOR JMakePDF -F5 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J5p.dat -d $DEBUG $BACKGROUND
        fi
        
        for R in 0.1 0.3 0.5 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 \
                     11.0 12.0 13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 \
                     22.0 24.0 26.0 28.0 30.0 32.0 34.0 36.0 38.0 40.0 42.0 44.0 46.0 48.0 50.0 \
                     55.0 60.0 65.0 70.0 75.0 80.0 85.0 90.0 95.0 100.0 110.0 120.0 130.0 140.0 150.0 \
                     170.0 190.0 210.0 230.0 250.0; do
            
            if (( $R <= $RMAX_M )); then

                if [[ ! -f $WORKDIR/J6p.dat ]]; then
                    eval $EXECUTOR JMakePDF -F6 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J6p-$R.dat -d $DEBUG $BACKGROUND
                fi
            fi
        done

        let "DMAX_M = ${RMAX_M}*1.5"

        for D in 0.1 0.5 1.0 5.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 \
            120.0 150.0 170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0; do

            if (( $D <= $DMAX_M )); then
                if [[ ! -f $WORKDIR/J12p.dat ]]; then
                    eval $EXECUTOR JMakePDG -F12 -A $absorptionLengthFactor -S $scatteringLengthFactor -D $D -o $WORKDIR/J12p-$D.dat -d $DEBUG $BACKGROUND
                fi
            fi
        done

        if [[ ! -f $WORKDIR/J13p.dat ]]; then
            eval $EXECUTOR JMakePDG -F13 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J13p.dat -d $DEBUG $BACKGROUND
        fi
        
        let "DMAX_M = ${RMAX_M}*1.5"

        for D in 0.1 0.5 1.0 5.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 \
            120.0 150.0 170.0 190.0 210.0 230.0 250.0 270.0 290.0 310.0 \
            340.0 370.0 400.0 430.0 460.0 490.0 520.0 550.0 600.0 650.0 700.0 750.0 800.0; do

            if (( $D <= $DMAX_M )); then
                if [[ ! -f $WORKDIR/J14p.dat ]]; then
                    eval $EXECUTOR JMakePDG -F14 -A $absorptionLengthFactor -S $scatteringLengthFactor -D $D -o $WORKDIR/J14p-$D.dat -d $DEBUG $BACKGROUND
                fi
            fi
        done

        if [[ ! -f $WORKDIR/J23p.dat ]]; then
            eval $EXECUTOR JMakePD0 -F23 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J23p.dat -d $DEBUG $BACKGROUND
        fi
        if [[ ! -f $WORKDIR/J24p.dat ]]; then
            eval $EXECUTOR JMakePD0 -F24 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J24p.dat -d $DEBUG $BACKGROUND
        fi

        if [[ "$BACKGROUND" != "" ]]; then
            wait
        fi
    fi
fi


if (( $MERGE_PDF )); then

    check_process JMergePDF pgrep
    
    if [[ ! -f $WORKDIR/J3p.dat ]]; then

        FILE_LIST=(`echo $WORKDIR/J3p-*.dat`)

        eval $EXECUTOR JMergePDF -f${^FILE_LIST} -o $WORKDIR/J3p.dat  -d $DEBUG $BACKGROUND
    fi
        
    if [[ ! -f $WORKDIR/J4p.dat ]]; then
        
        FILE_LIST=(`echo $WORKDIR/J4p-*.dat`)
        
        eval $EXECUTOR JMergePDF -f${^FILE_LIST} -o $WORKDIR/J4p.dat  -d $DEBUG $BACKGROUND
    fi
    
    if [[ ! -f $WORKDIR/J6p.dat ]]; then
        
        FILE_LIST=(`echo $WORKDIR/J6p-*.dat`)

        eval $EXECUTOR JMergePDF -f${^FILE_LIST} -o $WORKDIR/J6p.dat  -d $DEBUG $BACKGROUND
    fi

    if [[ ! -f $WORKDIR/J12p.dat ]]; then
    
        FILE_LIST=(`echo $WORKDIR/J12p-*.dat`)

        eval $EXECUTOR JMergePDG -f${^FILE_LIST} -o $WORKDIR/J12p.dat -d $DEBUG $BACKGROUND
    fi

    if [[ ! -f $WORKDIR/J14p.dat ]]; then
    
        FILE_LIST=(`echo $WORKDIR/J14p-*.dat`)

        eval $EXECUTOR JMergePDG -f${^FILE_LIST} -o $WORKDIR/J14p.dat -d $DEBUG  $BACKGROUND
    fi
    
    if [[ "$BACKGROUND" != "" ]]; then
        wait
    fi

    if [[ ! -f $WORKDIR/J113p.dat ]]; then
    
        eval $EXECUTOR JMakePDE -F $WORKDIR/J%p.dat -o $WORKDIR/J113p.dat -T $TTS_NS -d $DEBUG $BACKGROUND
    fi
    
    if [[ "$BACKGROUND" != "" ]]; then
        wait
    fi
fi


if (( $MAKE_CDF )); then

    check_process JMakeCDF JMakeCDG pgrep

    if [[ ! -f $WORKDIR/I1p.dat ]]; then
        eval $EXECUTOR JMakeCDF -f $WORKDIR/J1p.dat  -e 1e-5 -o $WORKDIR/I1p.dat  -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I2p.dat ]]; then
        eval $EXECUTOR JMakeCDF -f $WORKDIR/J2p.dat  -e 1e-4 -o $WORKDIR/I2p.dat  -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I5p.dat ]]; then
        eval $EXECUTOR JMakeCDF -f $WORKDIR/J5p.dat  -e 1e-5 -o $WORKDIR/I5p.dat  -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I6p.dat ]]; then
        eval $EXECUTOR JMakeCDF -f $WORKDIR/J6p.dat  -e 1e-5 -o $WORKDIR/I6p.dat  -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I12p.dat ]]; then
        eval $EXECUTOR JMakeCDG -f $WORKDIR/J12p.dat -e 1e-5 -o $WORKDIR/I12p.dat -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I13p.dat ]]; then
        eval $EXECUTOR JMakeCDG -f $WORKDIR/J13p.dat -e 1e-5 -o $WORKDIR/I13p.dat -d $DEBUG $BACKGROUND
    fi
    if [[ ! -f $WORKDIR/I14p.dat ]]; then
        eval $EXECUTOR JMakeCDG -f $WORKDIR/J14p.dat -e 1e-5 -o $WORKDIR/I14p.dat -d $DEBUG $BACKGROUND
    fi

    if [[ "$BACKGROUND" != "" ]]; then
        wait
    fi
fi

if (( $CLEAN )); then
    if ( rm -f $WORKDIR/J*p-*.dat ); then; fi    
fi

status "End $script"