plot-height.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
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             ./
set_variable  APP                 gandalf
set_variable  GRAPH_TXT           graph-height.txt
set_variable  GRAPH_ROOT          graph-height.root
set_variable  T                   T=+20.4
set_variable  M                   M=\*

if ( do_usage $* ); then
    usage "$script [working directory]"
fi

case $# in
    1) set_variable WORKDIR $1;;
    0) ;;
    *) fatal "Wrong number of arguments."
esac

if ( ! reuse_file $GRAPH_TXT ); then

    for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${T}:${M}`; do

   INPUT_FILES=($OUTPUT_DIRECTORY/*_${APP}\.root)

   echo "Processing:"
   echo "$INPUT_FILES[*]" | tr ' ' '\n'

   U=`echo ${OUTPUT_DIRECTORY##.*/} | tr ':' '\n' | grep -vF $T | grep -vF $M | sed 's/.*=//'`
   V=(`JQuality -f "$INPUT_FILES[*]" |& sed -n "s/Total\ quality\ *//p"`)

   printf "%5.3f %12.1f %12.1f\n" $U $V >> $GRAPH_TXT

    done
fi

if (( 1 )); then

    JGraph \
   -f $GRAPH_TXT              \
   -o $GRAPH_ROOT
    
    JFit \
   -f ${GRAPH_ROOT}:\.\*      \
   -o fit.root                \
   -F "[0] - [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3])) - 1.0)"  \
   -@"p0 = GetMaximum"        \
   -@"p1 = GetMaximum * 0.01" \
   -@"p2 =   1.0"             \
   -@"p3 =   1.0e-3"          \
   -d3

    JPlot1D \
   -f fit.root:\.\*           \
   -x "0.90 1.10"             \
   -> "height scaling"       \
   -\^ "Likelihood"           \
   -T ""

fi