# File: waterfitter.mac
# Brief: This uses the waterFitter
#
# Date: 2014-05-29
# Contact: M Mottram, <m.mottram@qmul.ac.uk>
# Revisions:
#       2014-05-29 : P Jones, added header info.
#
# This macro runs the waterFitter on simulated scintillator data.

/rat/physics_list/OmitMuonicProcesses true
/rat/physics_list/OmitHadronicProcesses true

# Choose an Nhit threshold
#/rat/db/set DAQ_RUN_LEVEL n100_hi_thresh 20.0
#/rat/db/set DAQ_RUN_LEVEL n100_med_thresh 18.0
#/rat/db/set DAQ_RUN_LEVEL n100_lo_thresh 16.0
#/rat/db/set DAQ_RUN_LEVEL n20_thresh 20.0
#/rat/db/set DAQ_RUN_LEVEL n20_lb_thresh 20.0
#/rat/db/set DAQ_RUN_LEVEL eshi_thresh 160.0
#/rat/db/set DAQ_RUN_LEVEL eslo_thresh 120.0
/rat/db/set MC thin_factor 3.5

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

/run/initialize
/rat/proc frontend
/rat/proc trigger
/rat/proc eventbuilder
/rat/proc count
/rat/procset update 10
/rat/proc calibratePMT
/rat/proc waterFitter
/rat/proc outroot
/rat/procset file "waterfitter.root"

# Generate 3 classes of events:
#   monoenergetic center
#   monoenergetic distributed
#   energy range center
/generator/add combo gun:point
/generator/vtx/set e- 0 0 0 6.0
/generator/pos/set 0 0 0

/generator/add combo gun:fillshell
/generator/vtx/set e- 0 0 0 6.0
/generator/pos/set 0 0 0 0 0 0 1.0 6005.0

/generator/add combo gun2:point
/generator/vtx/set e- 0 0 0 0 6.0 10
/generator/pos/set 0 0 0

/rat/run/start 1000