# File: isotropy.mac
# Brief: Short macro showing how to run the isotropy classifier
#
# Date: 19-11-2014
# Contact: E. W. Beier, <geneb@hep.upenn.edu>
#
# This is a macro which shows how to run the Isotropy classifier.
# This returns six classifications, thetaij, snobeta14, beta1, beta2,
# beta3 and beta4.
#

/run/initialize

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

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

# Now run the isotropy classifier
# Choose a different time interval from that in CLASSIFIER_ISOTROPY.ratdb
/rat/proc classifier
/rat/procset classifier "isotropy"
/rat/procset seed "quad"
/rat/procset selector "timeResidualCut"
/rat/procset selector.lowLimit  -9.
/rat/procset selector.highLimit +1.5

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