#/**
#*	Example Macro for generating pileup C14 and Te130 2n2b events.
#*	Additionally the pileup classifiers are included in the event loop.
#*
#*				pileup.mac
#*	Date:       2014-01-08
#*	contact:    P. G. Jones, QMUL
#*/

# Don't care about hadrons or muons so quicker not to initialise these processes
/rat/physics_list/OmitMuonicProcesses true
/rat/physics_list/OmitHadronicProcesses true

# BEGIN EVENT LOOP
/run/initialize
/rat/proc frontend
/rat/proc trigger
/rat/proc eventbuilder
/rat/proc calibratePMT

/rat/proc scintFitter
# Now the pileup classifiers
/rat/proc classifier
/rat/procset classifier "preTriggerHits"

/rat/proc classifier
/rat/procset classifier "timingPeaks"

/rat/proc classifier
/rat/procset classifier "meanTime"
/rat/procset seed "scintFitter"

/rat/proc classifier
/rat/procset classifier "earlyTime"
/rat/procset seed "scintFitter"

/rat/proc classifier
/rat/procset classifier "isoRegions"
/rat/procset seed "scintFitter"

/rat/proc count
/rat/procset update 10
/rat/proclast outroot
/rat/procset file "pileup.root"
# END EVENT LOOP

# As the timing is important do both orderings
# 2b2n + C14
/generator/add coincidence decay0:fill
/generator/vtx/set 2beta Te130 0 4
/generator/pos/set 0 0 0

/generator/coincidence/add decay0:fill
/generator/coincidence/vtx/set backg C14
/generator/coincidence/pos/set 0 0 0

/generator/coincidence/timewindow 450 # 450ns is longer than a full event window, thus is conservative
/generator/rate/set 1

# C14 + 2b2n
/generator/add coincidence decay0:fill
/generator/vtx/set backg C14
/generator/pos/set 0 0 0

/generator/coincidence/add decay0:fill
/generator/coincidence/vtx/set 2beta Te130 0 4
/generator/coincidence/pos/set 0 0 0

/generator/coincidence/timewindow 450 # 450ns is longer than a full event window, thus is conservative
/generator/rate/set 1

/rat/run/start 10
exit