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

# check whether script is already running

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

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

#
# Method to print usage of script.
#
function usage()
{
    print -u2 "$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\
        -w download files from git server                             \n\
        -R <maximal distance between muon and PMT>                    \
        "
    exit
}

#
# Parse command line options.
#
while getopts "d:W:PMCcBw: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  \&;;
        w) set_variable WGET        $OPTARG;;
        h) usage; exit;;
        A) set_variable absorptionLengthFactor  $OPTARG;;
        S) set_variable scatteringLengthFactor  $OPTARG;;
        R) set_variable RMAX_M      $OPTARG;;
    esac
done

check_variable ROOTSYS JPP_DIR WORKDIR

mkdir -p $WORKDIR

status "Start $script"


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


if (( $MAKE_PDF )); then

    check_process JMakePDF JMakePDG

    if (( 1 )); then

        if [[ ! -f $WORKDIR/J1p.dat ]]; then
            eval nohup JMakePDF -F1 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J1p.dat -d $DEBUG $BACKGROUND
        fi
        if [[ ! -f $WORKDIR/J2p.dat ]]; then
            eval nohup 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 nohup JMakePDF -F3 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J3p-$R.dat -d $DEBUG $BACKGROUND
                fi
                if [[ ! -f $WORKDIR/J4p.dat ]]; then
                    eval nohup 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 nohup 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 nohup JMakePDF -F6 -A $absorptionLengthFactor -S $scatteringLengthFactor -R $R -o $WORKDIR/J6p-$R.dat -d $DEBUG $BACKGROUND
                fi
            fi
        done

        if [[ ! -f $WORKDIR/J13p.dat ]]; then
            eval nohup 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 nohup 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 nohup JMakePD0 -F23 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J23p.dat -d $DEBUG $BACKGROUND
        fi
        if [[ ! -f $WORKDIR/J24p.dat ]]; then
            eval nohup JMakePD0 -F24 -A $absorptionLengthFactor -S $scatteringLengthFactor -o $WORKDIR/J24p.dat -d $DEBUG $BACKGROUND
        fi
    fi

fi


if (( $MERGE_PDF )); then

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

        FILE_LIST=`ls $WORKDIR/J3p-*.dat`

        JMergePDF -f $FILE_LIST -o $WORKDIR/J3p.dat  -d $DEBUG
    fi
        
    if [[ ! -f $WORKDIR/J4p.dat ]]; then
        
        FILE_LIST=`ls $WORKDIR/J4p-*.dat`
        
        JMergePDF -f $FILE_LIST -o $WORKDIR/J4p.dat  -d $DEBUG
    fi
    
    if [[ ! -f $WORKDIR/J6p.dat ]]; then
        
        FILE_LIST=`ls $WORKDIR/J6p-*.dat`

        JMergePDF -f $FILE_LIST -o $WORKDIR/J6p.dat  -d $DEBUG
    fi

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

        JMergePDG -f $FILE_LIST -o $WORKDIR/J14p.dat -d $DEBUG
    fi
fi


if (( $MAKE_CDF )); then

    check_process JMakeCDF JMakeCDG

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

if (( $CLEAN )); then

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

status "End $script"