# File: socFitter.mac
# Brief: Short macro showing how to run the soc fitter
#
# Date: 19-05-2014
# Contact: P G Jones, <p.g.jones@qmul.ac.uk>
# Revisions:
#       19-05-2014 I. Coulter <icoulter@hep.upenn.edu> Initial commit.
#       19-05-2015 J. Maneira <maneira@lip.pt> Updates to generator
#                  and DB interface, turn on extasy trigger and conditional
#                  processing for triggered events only.
#
# This is a macro which shows how to run the socFitter.
# Details?
#


/rat/physics_list/OmitMuonicProcesses true
/rat/physics_list/OmitHadronicProcesses true

# To use simulated external asynch trigger
/rat/db/set DAQ_RUN_LEVEL trigger_mask 32768
/rat/db/set DAQ_RUN_LEVEL trigger_enable 35967

/rat/db/set LASERBALL_SIMULATION intensity 10000
# Position is not an integer -> make sure there's a "."
/rat/db/set LASERBALL_SIMULATION position [0.0, 0.0, 0.0]
/rat/db/set LASERBALL_SIMULATION dye_name "BIS-MSB"
/rat/db/set LASERBALL_SIMULATION lb_id "OCT03"

/run/initialize

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

# Require soc data
/rat/proc/if trigTypeSelector
/rat/procset trigType "EXTASY"
    /rat/proc count
    /rat/procset update 1
    /rat/proc socdata
/rat/proc/endif

# First need to run the peak finder
/rat/proc socpeakfind

# Run socFitter
/rat/proc socfitter
/rat/procset method "socPositionTimeChiSquared"
/rat/procset optimiser "minuit"
/rat/procset selector "null"
# Changing fitter name
/rat/procset name "socfit"

/rat/proc outsoc
/rat/procset file "soc.root"

# END EVENT LOOP

# Choose the laserball generator from fixed point, uniform in time
/generator/add combo vertexlaserball:poslaserball:poisson
# Do not use the generator/pos/set or /generator/vtx/set commands.
# All variables should be set as ratdb variables, so that other
# routines, in addition to the generator, can access them

/generator/rate/set 100

/rat/run/start 10
exit