# File: QEnergy.mac
# Brief: Short macro showing how to run the QEnergy.mac fitter
#
# Date: 06-09-2019
# Contact: Yang Zhang, <yzhang23@ualberta.ca>
#
# This is a macro which shows how to run the QEnergy fitter.

# Load a partial fill phase geometry
/rat/db/set DETECTOR geo_file "geo/snoplus_partial.geo"
/rat/db/set GEO[inner_av] split_z 0.0

/run/initialize

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

# Requires a position seed
# Use the partial fitter.
/rat/proc partialFitter

/rat/proc fitter
/rat/procset method "QpartialEnergy"
/rat/procset seed "partialFitter"

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

#/generator/add combo gun:fill:poisson
/generator/add combo gun:point:poisson
/generator/vtx/set e- 0 0 0 1.0
/generator/pos/set 0 0 2000
#/generator/pos/set inner_av_top
/generator/rate/set 1

/rat/run/start 10
exit