# File: radius_cut.mac
# Brief: Macro showing how to cut on reconstructed radius.
#
# Date: 2015-11-20
# Contact: J Walker <john.walker@liverpool.ac.uk>
# Revisions:
#
# This macro demonstrates how to cut events with waterFitter position
# reconstruction outside the inner AV.

/rat/db/set DETECTOR geo_file "geo/snoplus_water.geo"

/run/initialize

# BEGIN EVENT LOOP
/rat/proc frontend
/rat/proc trigger
/rat/proc eventbuilder
/rat/proc calibratePMT
/rat/proc count
/rat/procset update 10

/rat/proc waterFitter

/rat/proc/if radiusCut
/rat/procset radius 6005.0
/rat/procset fitName "waterFitter"
/rat/procset mode "greater"
    /rat/proc outroot
    /rat/procset file "radius_cut.root"
/rat/proc/endif
# END EVENT LOOP

/generator/add combo gun:fill
/generator/vtx/set e- 0 0 0 6.0
/generator/pos/set 0 0 0
/generator/rate/set 1

/rat/run/start 10

exit