JPDFSupportkit.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#              Utility script to plot PDF configuration.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
 
set_variable   DEBUG       2      
set_variable   WORKDIR     ${TMPDIR:-/tmp}/  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [working directory]"
fi
 
case $# in
    1) set_variable  WORKDIR  $1;;
esac
 
ABSORPTION_LENGTHS=("nemo          1.0 1.0e6"
                    "SmithAndBaker 1.0 1.0e6"
                    "SmithAndBaker 1.0 610.0")
 
for ABSORPTION_LENGTH in $ABSORPTION_LENGTHS[*]; do
 
    $JPP_DIR/examples/JPhysics/JPDFSupportkit \
        -@ "absorptionLength = $ABSORPTION_LENGTH"           \
        -o $WORKDIR/h\["${ABSORPTION_LENGTH}"\].root         \
        -d $DEBUG
done
 
JPlot1D \
    -f$WORKDIR/h\[${^ABSORPTION_LENGTHS}\].root:h1           \
    -y "0.0 100.0"                                           \
    -> "#lambda [nm]"                                       \
    -\^ "#lambda_{abs} [m]"                                  \
    -G XY                                                    \
    -L TL                                                    \
    -T ""                                                    \
    -o absorption-length.$FORMAT                             $BATCH
 
rm -f $WORKDIR/h\[${^ABSORPTION_LENGTHS}\].root
 
 
SCATTERING_LENGTHS=("default  1.0"
                    "default  0.8"
                    "default  1.2")
 
for SCATTERING_LENGTH in $SCATTERING_LENGTHS[*]; do
 
    $JPP_DIR/examples/JPhysics/JPDFSupportkit \
        -@ "scatteringLength = $SCATTERING_LENGTH"           \
        -o $WORKDIR/h\["${SCATTERING_LENGTH}"\].root         \
        -d $DEBUG
done
 
JPlot1D \
    -f$WORKDIR/h\[${^SCATTERING_LENGTHS}\].root:h2           \
    -y "0.0 160.0"                                           \
    -> "#lambda [nm]"                                       \
    -\^ "#lambda_{s} [m]"                                    \
    -G XY                                                    \
    -L TL                                                    \
    -T ""                                                    \
    -o scattering-length.$FORMAT                             $BATCH
 
rm -f $WORKDIR/h\[${^SCATTERING_LENGTHS}\].root
 
 
QES=("default  1.0"
     "default  0.9"
     "default  1.1")
 
for QE in $QES[*]; do
 
    $JPP_DIR/examples/JPhysics/JPDFSupportkit \
        -@ "QE = $QE"                                        \
        -o $WORKDIR/h\["${QE}"\].root                        \
        -d $DEBUG
done
 
JPlot1D \
    -f$WORKDIR/h\[${^QES}\].root:h3                          \
    -y "0.0 0.40"                                            \
    -> "#lambda [nm]"                                       \
    -\^ "QE"                                                 \
    -G XY                                                    \
    -L TR                                                    \
    -T ""                                                    \
    -o QE.$FORMAT                                            $BATCH
 
rm -f $WORKDIR/h\[${^QES}\].root
 
 
ANGULAR_ACCEPTANCES=("default")
 
for ANGULAR_ACCEPTANCE in $ANGULAR_ACCEPTANCES[*]; do
 
    $JPP_DIR/examples/JPhysics/JPDFSupportkit \
        -@ "angularAcceptance = $ANGULAR_ACCEPTANCE" \
        -o $WORKDIR/h\["${ANGULAR_ACCEPTANCE}"\].root        \
        -d $DEBUG
done
 
JPlot1D \
    -f$WORKDIR/h\[${^ANGULAR_ACCEPTANCES}\].root:h4          \
    -y "0.0 1.2"                                             \
    -> "cos(#theta_{s})"                                    \
    -\^ "P"                                                  \
    -G XY                                                    \
    -L TL                                                    \
    -T ""                                                    \
    -o angular-acceptance.$FORMAT                            $BATCH
 
rm -f $WORKDIR/h\[${^ANGULAR_ACCEPTANCES}\].root
 
 
SCATTERING_PROBABILITYS=("default   1.0"
                         "f4        1.0"
                         "p00075    1.0"
                         "petzhold  1.0")
 
for SCATTERING_PROBABILITY in $SCATTERING_PROBABILITYS[*]; do
 
    $JPP_DIR/examples/JPhysics/JPDFSupportkit \
        -@ "scatteringProbability = $SCATTERING_PROBABILITY" \
        -o $WORKDIR/h\["${SCATTERING_PROBABILITY}"\].root    \
        -d $DEBUG
done
 
JPlot1D \
    -f$WORKDIR/h\[${^SCATTERING_PROBABILITYS}\].root:h5      \
    -y "1.0e-3 1.0e2"                                        \
    -Y                                                       \
    -> "cos(#theta_{s})"                                    \
    -\^ "P"                                                  \
    -G XY                                                    \
    -L TL                                                    \
    -T ""                                                    \
    -o scattering-probability.$FORMAT                        $BATCH
 
rm -f $WORKDIR/h\[${^SCATTERING_PROBABILITYS}\].root