TestRPCsim.py

Go to the documentation of this file.

#!/usr/bin/env python

import GaudiKernel.SystemOfUnits as units
from DybPython.Constants import *

class ConfigureTest:

    '''
    simple usage example:

    nuwa.py --history=off -n 10 -o output.root TestRPCsim

    For a list of the options, run

    nuwa.py -m "TestRPCElecSim --help"

    Note: You need to turn off simheader conversion to run this module:

    export NUWA_PRUNESIMHITS=1 (sh)
    setenv NUWA_PRUNESIMHITS 1 (csh)
    '''

    def __init__(self,argv):

        self.parse_args(argv)

        return

    def parse_args(self,argv):

        import optparse

        parser = optparse.OptionParser(usage=self.__doc__)

        parser.add_option("-D","--detector",
                          help="List of detectors to be simulated. [default: %default]",
                          default='DayaBayRPC',
                          type="string")
        parser.add_option("-c","--single_Channel",
                          help="Only simulate 1 channel per detector. [default: %default]",
                          action='store_true',
                          default=False)
        parser.add_option("-t","--delay-Time",
                          help="Delay time of in ns. [default: %default]",
                          default=500,
                          type="float")
        parser.add_option("-w","--weight",
                          help="Weight of the hits. [default: %default]",
                          default=1,
                          type="float")
        parser.add_option("-s","--time-Spread",
                          help="Sigma of Gaussian hit time distribution. [default: %default]",
                          default=0,
                          type="float")

        (options,args) = parser.parse_args(args=argv)
        self.opts = options
        self.args = args


        print "======================================================"
        print "Simulate detector ",self.opts.detector
        print "------------------------------------------------------"
        print "Only simulate 1 channel  = ",self.opts.single_Channel
        print "------------------------------------------------------"
        print "Delay time = ",self.opts.delay_Time, " ns"
        print "------------------------------------------------------"
        print "Hit weight = ",self.opts.weight
        print "------------------------------------------------------"
        print "Gaussian hit time spread = ",self.opts.time_Spread, " ns"
        print "======================================================"

        return

    def configure(self):

        from ElecSim.ElecSimConf import EsTestAlg
        testAlg = EsTestAlg()
        testAlg.Detectors = [self.opts.detector]
        testAlg.DelayTime = self.opts.delay_Time
        testAlg.Weight = self.opts.weight
        testAlg.TimeSmear = self.opts.time_Spread
        testAlg.SingleChannel = self.opts.single_Channel

        from Gaudi.Configuration import ApplicationMgr
        app = ApplicationMgr()
        app.TopAlg.append(testAlg)
        app.HistogramPersistency = 'ROOT'

        import ElecSim
        elecsim = ElecSim.Configure()

        from ElecSim.ElecSimConf import EsIdealFecTool
        fect = EsIdealFecTool()
        #fect.IdealNoise = 1*units.kilohertz

        import TrigSim
        trigsim = TrigSim.Configure()

        import ReadoutSim
        rosim = ReadoutSim.Configure()

        return

def configure(argv=[]):

    cft = ConfigureTest(argv)
    cft.configure()

    return