# File: unseededAlpha.mac
# Brief: Short macro showing how to run the unseededAlpha classifier
#
# Date: 19-05-2014
# Contact: K Majumdar, <Krishanu.Majumdar@physics.ox.ac.uk>
# Revisions:
#       19-05-2014 I. Coulter <icoulter@hep.upenn.edu> Initial commit.
#
# This is a macro which shows how to run the unseededAlpha classifier.
# Looks at the distribution of PMT trigger times, as a longer tail on
# this distribution indicates an alpha. Will return the ratio of PMT
# hits in the tail to the total number of PMT hits in the event.
#

# The time of the tail and total event window are defined in
# CLASSIFIER_ALPHA_UNSEEDED.ratdb.
# This can be overwritten by the user using the commands:
/rat/db/set CLASSIFIER_ALPHA_UNSEEDED tail_fraction 0.45d
/rat/db/set CLASSIFIER_ALPHA_UNSEEDED pmt_Min 250.0d
/rat/db/set CLASSIFIER_ALPHA_UNSEEDED pmt_Max 500.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

# Run the Unseeded Alpha classifier
/rat/proc classifier
/rat/procset classifier "UnseededAlpha"

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