In This Package:

FilterExample.py

Go to the documentation of this file.

#!/usr/bin/env python
#
# Example of filter the recon data from a file containing raw and
# reconstructed data.
#
# Usage:
#   nuwa.py -o reconData.root DybAlg.FilterExample fullData.root
#

# Load GaudiPython
from DybAlg.FilterAlgs import FilterAlg, CoincidenceFilterAlg
import GaudiKernel.SystemOfUnits as units

def configure():
    return

def run(app):
    '''
    Add filtering algorithm to the job
    '''
    # Store all Recon data
    reconFilter = FilterAlg("ReconFilter")
    reconFilter.OutputLevel = 1
    reconFilter.ClearStore = True  # First filter should clear storage
    reconFilter.StorePath = "/Event/Rec/RecHeader"
    app.addAlgorithm(reconFilter)

    # Save simulation and readout data with a prescale.
    singlesReadoutFilter = FilterAlg("SinglesReadoutFilter")
    singlesReadoutFilter.OutputLevel = 1
    # Secondary filters shouldn't re-clear storage
    singlesReadoutFilter.ClearStore = False 
    singlesReadoutFilter.StorePath = "/Event/Readout/ReadoutHeader"
    singlesReadoutFilter.Prescale = 0.01
    app.addAlgorithm(singlesReadoutFilter)
    
    # Save readouts if both within coincidence time window
    coincFilter = CoincidenceFilterAlg("CoincidenceFilter")
    coincFilter.OutputLevel = 1
    #coincFilter.StorePaths = ["/Event/Readout/ReadoutHeader"]
    coincFilter.StorePaths = ["/Event/Readout/ReadoutHeader",
                              "/Event/SimReadout/SimReadoutHeader",
                              "/Event/SimTrig/SimTrigHeader",
                              "/Event/Elec/ElecHeader",
                              "/Event/Sim/SimHeader",
                              "/Event/Gen/GenHeader"]
    coincFilter.CoincidenceWindow = 15.0 * units.millisecond
    # Tag coincidences in same AD only
    # coincFilter.AdOnly = True
    # Tag coincidences in any detector, as long as prompt/delayed in same det.
    # coincFilter.SameDetectorOnly = True
    # Tag coincidences in a list of detectors
    coincFilter.DetectorNames = ["DayaBayAD1", "DayaBayAD2"]
    # Save other detector readouts between prompt/delayed readouts
    coincFilter.SaveIntermediate = True
    # Save readouts occuring in other detectors in small time window
    # before prompt.
    coincFilter.SaveNeighbors = True
    coincFilter.NeighborWindow = 1.0 * units.microsecond
    app.addAlgorithm(coincFilter)
    pass

---