plot-position.sh

Go to the documentation of this file.

#!/usr/bin/env 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       1                    
set_variable   WORKDIR     ./                   
set_variable   APP         gandalf              
set_variable   GRAPH_TXT   graph-position.txt   
set_variable   GRAPH_ROOT  graph-position.root  
 
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
 
    set_variable T   T=0
    set_variable X   X=\*
    set_variable Y   Y=\*
    set_variable Z   Z=0
 
    for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${T}:${X}:${Y}:${Z}`; 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 $X | grep -vF $Y | grep -vF $Z | sed 's/.*=//'`)
        V=(`JQuality -f "$INPUT_FILES[*]" |& sed -n "s/Total\ quality\ *//p"`)
        
        printf "%5.1f %5.1f %12.1f %12.1f\n" $U $V >> $GRAPH_TXT
 
    done
fi
 
if (( 1 )); then
 
    JGraph2D \
        -f $GRAPH_TXT              \
        -o $GRAPH_ROOT
    
    JFit2D \
        -f ${GRAPH_ROOT}:\.\*      \
        -o fit.root                \
        -F "[0] - [1]*(sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([4]*[4])) - 1.0)"  \
        -@"p0 = GetMaximum"        \
        -@"p1 = GetMaximum * 0.01" \
        -@"p2 =   0.0e0"           \
        -@"p3 =   0.0e0"           \
        -@"p4 =  30.0e0"           \
        -d 3
 
    JPlot2D \
        -f ${GRAPH_ROOT}:\.\*      \
        -f fit.root:user           \
        -> "#Deltax [m]"          \
        -< "#Deltay [m]"          \
        -\^ "Likelihood"           \
        -O SURF                    \
        -T ""
 
fi