JFitVolume.sh

Go to the documentation of this file.

#!/bin/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               


if do_usage $*; then
    usage "$script [detector file [trigger file [input file]]]"
fi

case $# in
    3) INPUT_FILE=$3;&
    2) TRIGGER_FILE=$2;&
    1) 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.gif