# File: timeFit.mac
# Brief: Short macro showing how to run the time likelihood fitter
#
# Date: 19-05-2014
# Contact: W Heintzelman, <billh@hep.upenn.edu>
# Revisions:
#       19-05-2014 I. Coulter <icoulter@hep.upenn.edu> Initial commit.
#
# This is a macro which shows how to run a time likelihood fitter.
# This uses the timeLikelihood method, the simple Timing PDF, a Minuit
# optimiser and the null PMTselector which uses all PMTs. It uses a
# seed set a fixed position.
#

/run/initialize

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

# Choose a fixed position for the seed
/rat/proc fitter
/rat/procset method "fixedSeed"
/rat/procset method.position 0 100 0 200

# Run the position likelihood fitter
/rat/proc fitter
/rat/procset method "timeLikelihood"
/rat/procset optimiser "minuit"
/rat/procset pdf "simpleTiming"
# Probably intended for Cherenkov source PDF but doesn't seem to work?
/rat/procset seed "fixedSeed"
/rat/procset selector "null"

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

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

/rat/run/start 10
exit