# File: qpdt.mac
# Brief: Short macro showing how to run the QPDT classifier
#
# Date: 19-05-2014
# Contact: E Marzec, <marzece@sas.upenn.edu>
# Revisions:
#       19-05-2014 I. Coulter <icoulter@hep.upenn.edu> Initial commit.
#
# This is a macro which shows how to run the QPDT classifier. This
# compares the observed PMT charges for hits that are in an early time
# window to a DiPolya+Exponential distribution to obtain the hit's
# charge probability. This is then used to calculate the event's QPDT
# which is returned along with the number of hits within the early
# time window.
#

# Time window defined in CLASSIFIER_QPDT.ratdb.
# This can be over written by the user as:
/rat/db/set CLASSIFIER_QPDT t1 -75.0d
/rat/db/set CLASSIFIER_QPDT t2 -20.0d

/run/initialize

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

# QPDT classifier requires a position+time seed.
# Use the quadFitter.
/rat/proc fitter
/rat/procset method "quad"

# Now run the QPDT classifier
/rat/proc classifier
/rat/procset classifier "QPDT"
/rat/procset seed "quad"

/rat/proclast outroot
/rat/procset file "qpdt.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