plot-orientation.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-orientation-%.txt   
set_variable   GRAPH_ROOT  graph-orientation-%.root  
set_variable   ID          ID=2                      
set_variable   R           R=\*                      

if do_usage $*; then
    usage "$script  <string identifier> [working directory]"
fi

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

set_variable  GRAPH_TXT           ${GRAPH_TXT/\%/${ID}}
set_variable  GRAPH_ROOT          ${GRAPH_ROOT/\%/${ID}}

if ( ! reuse_file $GRAPH_TXT ); then

    for OUTPUT_DIRECTORY in `eval echo $WORKDIR/${ID}:${R}`; do

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

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

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

        printf "%5.1f %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]*cos(x-[2])"  \
        -@"p0 = GetMaximum"        \
        -@"p1 = GetMaximum * 0.01" \
        -@"p2 =   0.0e0"           \
        -d3

    JPlot1D \
        -f fit.root:\.\*           \
        -x "-3.14 +3.14"           \
        -\> "#phi [rad]"           \
        -\^ "Likelihood"           \
        -T "String ${ID}"

fi