# File: partialfitter.mac
# Brief: This uses the partialFitter.
#
# Date: 2022-03-22
# Contact: W. Parker, <william.parker@physics.ox.ac.uk>
# Revisions:
#
# This macro runs the partialFitter on simulated partial fill data.
# It's based on the scintFitter rat test macro

# Load a partial fill phase geometry
/rat/db/set DETECTOR geo_file "geo/snoplus_partial.geo"
/rat/db/set GEO[inner_av] material_top "labppo_0p5_scintillator" # Measured 0.53 g/L PPO
/rat/db/set GEO[inner_av] split_z -5950

/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 2.5

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

# Generate 3 classes of events:
#   monoenergetic at (0,0,1)m
#   monoenergetic distributed in inner_av_top
#   energy range at (0,0,1)m
/generator/add combo gun:point
/generator/vtx/set e- 0 0 0 1.0
/generator/pos/set 0 0 1000

/generator/add combo gun:fill
/generator/vtx/set e- 0 0 0 1.0
/generator/pos/set 0 0 1000

/generator/add combo gun2:point
/generator/vtx/set e- 0 0 0 0 1.0 3.5
/generator/pos/set 0 0 1000

/rat/run/start 1000