JFitVolume.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author lquinn
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                Utility script to test JVolume1D and JFitEfficiency
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
 
set_variable   DEBUG         1                                             
set_variable   WORKDIR       ./                                            
set_variable   TRIGGER_FILE  $JPP_DATA/trigger_efficiency+background.root  
set_variable   INPUT_FILE    $JPP_DATA/start.root                          
set_variable   DETECTOR      $JPP_DATA/km3net_reference.detx               
set_variable:  FORMAT        GRAPHICS_FORMAT        gif
set_variable+  BATCH         GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [detector file [trigger file [input file]]]"
fi
 
case $# in
    3) set_variable  INPUT_FILE    $3;&
    2) set_variable  TRIGGER_FILE  $2;&
    1) set_variable  DETECTOR      $1;;
esac
 
if (( 1 )); then
 
  JPrintTree -f $TRIGGER_FILE
 
  echo "Check that the Monte Carlo and DAQ trees have different numbers of events. Else run JTriggerEfficiency without the -O option.\n"
 
fi
 
if (( 1 )); then 
 
  print_variable TRIGGER_FILE INPUT_FILE DETECTOR
 
  $JPP_DIR/examples/JTrigger/JVolume1D \
      -f $TRIGGER_FILE          \
      -o "volume[trigger].root" \
      -X                        \
      -d $DEBUG
 
  $JPP_DIR/examples/JFit/JFitEfficiency \
      -f $INPUT_FILE          \
      -o efficiency.root      \
      -a $DETECTOR            \
      -X -C                   \
      -W "2 0"                \
      -Q 0.5                  \
      -d $DEBUG
fi
 
JOpera1D  \
  -f "volume[trigger].root:hV"  \
  -f efficiency.root:he         \
  -o "volume[reco].root"        \
  -u "Multiply"                 \
  -d $DEBUG
 
 
JPlot1D \
    -f "volume[reco].root:Multiply"  \
    -f "volume[trigger].root:hV"     \
    -XX                              \
    -> "E [GeV]"                    \
    -\^ "V [km^{3}]"                 \
    -LL                              \
    -o volume.$FORMAT                $BATCH