# File: multiPDFMethodPartial.mac
# Brief: Short macro showing how to run the multiPDF position
#        fitting method in partial fill
#
# Date: 06-12-2020
# Contact: W. Parker, <william.parker@physics.ox.ac.uk>
# Revisions:
#
# This is a macro which shows how to run the multiPDF position fitting
# method. Currently, this uses different PDFs depending on the
# pmt height and seeded event position (here from the quad fitter)). 
# They are split into 4 bands of pmt height and two event positions.
#

# Load a partial fill phase geometry
/rat/db/set DETECTOR geo_file "geo/snoplus_partial.geo"

#ppo concentration and fill level the current pdfs apply to
/rat/db/set GEO[inner_av] material "labppo_0p5_scintillator"
/rat/db/set GEO[inner_av] material_top "labppo_0p5_scintillator"
/rat/db/set GEO[inner_av] split_z 747.48

#Set number of position dependent and non-position dependent pdfs
/rat/db/set ET1D_Partial_MultiPDF[labppo_0p5_scintillator]  num_posdep_pdfs 8
/rat/db/set ET1D_Partial_MultiPDF[labppo_0p5_scintillator]  num_nonposdep_pdfs 4

/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
/rat/proc fitter
/rat/procset method "quad"

/rat/proc fitter
/rat/procset method "multipdf-scintwater"
/rat/procset seed "quad"
/rat/procset name "multiPDFFit"

/rat/proc outroot
/rat/procset file "multiPDFPartial.root"
# END EVENTLOOP

/generator/add combo gun:fill

/generator/vtx/set e- 0 0 0 5.0
/generator/pos/set 0 0 3000
/generator/rate/set 1

/rat/run/start 10
exit