# File: BiPoLikelihood.mac
# Brief: Short macro showing how to run the BiPoLikelihoodDiff 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.
#       11-05-2015 M. Mottram <m.mottram@qmul.ac.uk> Modified syntax to invoke classifier.
#
# This is macro showing the commands need to run the likelihood-based
# BiPo classifier. It uses the log-likelihood difference from comparing
# time residual distributions to classify events as either pileup
# 212/214BiPo coincidence events or 130Te 0vBB. The classifier returns
# the difference between the log-likelihood values of the 130Te and BiPo
# events distributions.
#

# The range of the time residual distribution is defined in the
# CLASSIFIER_BIPO_LIKELIHOODDIFF_xxx.ratdb files.
# Also, the trigger window times and the meanPoAlphaNhits.
# These can be overwritten by the user using the commands:
/rat/db/set CLASSIFIER_BIPO_LIKELIHOODDIFF[te_0p3_labppo_scintillator_bisMSB_Dec2013_212BiPo] meanPoAlphaNhits 200.300195868
/rat/db/set CLASSIFIER_BIPO_LIKELIHOODDIFF[te_0p3_labppo_scintillator_bisMSB_Dec2013_212BiPo] triggerWindowStartTime -180.0
/rat/db/set CLASSIFIER_BIPO_LIKELIHOODDIFF[te_0p3_labppo_scintillator_bisMSB_Dec2013_212BiPo] triggerWindowEndTime 220.0

/run/initialize

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

# BiPoLikelihoodClassifier requires a position+time seed.
/rat/proc fitter
/rat/procset method "quad"

# Now run the classifier.
# It requires an optimiser, grid search works best, and a PDF
# as found in one of the CLASSIFIER_BIPO_LIKELIHOODDIFF_xxx.ratdb files.

# Run the classifier in default mode
# This will run the classifier using PDFs for the inner AV material and 214BiPo
/rat/proc classifier
/rat/procset classifier "BiPoLikelihoodDiff"
/rat/procset optimiser "grid-55"
/rat/procset seed "quad"

# Run the classifier again, using a named material and BiPo212
# Also set a name (otherwise, both classifiers have the same name and only one is stored)
/rat/proc classifier
/rat/procset classifier "BiPoLikelihoodDiff-te_0p3_labppo_scintillator_bisMSB_Dec2013-212"
# Note that if specifying the material, the 212/214 index MUST be set
# It is possible to specify just the 212/214 index and use the default material with:
# /rat/procset classifier "BiPoLikelihoodDiff-212"
/rat/procset optimiser "grid-55"
/rat/procset seed "quad"
/rat/procset name "BiPoLikelihoodDiffTe"

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