In This Package:

DefaultFilter.py

Go to the documentation of this file.

#!/usr/bin/env python
#
# Default first-pass filter for reconstructed data.
#
# Usage:
#   nuwa.py -o reconData.root DybAlg.DefaultFilter 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.ClearStore = True  # First filter should clear storage
    reconFilter.StorePath = "/Event/Rec/RecHeader"
    app.addAlgorithm(reconFilter)

    # Save readout data for singles events in all detectors with a prescale.
    singlesReadoutFilter = FilterAlg("SinglesReadoutFilter")
    # Secondary filters shouldn't re-clear storage
    singlesReadoutFilter.ClearStore = False 
    singlesReadoutFilter.StorePath = "/Event/Readout/ReadoutHeader"
    singlesReadoutFilter.Prescale = 0.01
    app.addAlgorithm(singlesReadoutFilter)
    
    # Save AD coincidence readouts based on coincidence time window
    coincFilter = CoincidenceFilterAlg("AdCoincidenceFilter")
    #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 = 2.0 * units.millisecond
    # Identify coincidences in same AD only
    coincFilter.AdOnly = True
    # Save other detector readouts occuring between prompt/delayed readouts
    coincFilter.SaveIntermediate = True
    # Save readouts occuring in other detectors in small time window
    # before prompt AD readout
    coincFilter.SaveNeighbors = True
    coincFilter.NeighborWindow = 1.0 * units.microsecond
    app.addAlgorithm(coincFilter)
    pass

---