# File: waterFitterSmear.mac
# Brief: Short macro showing how to run the waterFitter with 
#      reconstruction smearing turned on
#
# Date: 11-01-2017
# Contact: Morgan Askins, <maskins@ucdavis.edu>
# Revisions:
#       11-01-2017 M. Askins Initial commit.
#
# This macro runs the waterFitter with a smearing function applied
# to the position, direction, and energy reconstruction.


# Load a water phase geometry
/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 count
/rat/procset update 10
/rat/proc calibratePMT

# Run the waterFitter
/rat/proc waterFitter
# Energy
/rat/procset smear.energy "on"
/rat/procset smear.function "TMath::Gaus(x, [0], 1)"
/rat/procset smear.low -100.0
/rat/procset smear.high 100.0
# Position
/rat/procset smear.position "on"
/rat/procset smear.function "TMath::Gaus(x, [0], 1000)"
/rat/procset smear.low -10000.0
/rat/procset smear.high 10000.0
# Direction
/rat/procset smear.direction "on"
/rat/procset smear.function "TMath::Gaus(x, [0], 0.15)"
/rat/procset smear.low -100.0
/rat/procset smear.high 100.0

/rat/proclast outroot
/rat/procset file "waterFitSmeared.root"
# END EVENTLOOP

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

/rat/run/start 10
exit