In This Package:

Default.py

Go to the documentation of this file.

#!/usr/bin/env python

# Available physics lists
physics_list_basic = [
    "DsPhysConsGeneral", 
    "DsPhysConsOptical",
    "DsPhysConsEM",
    ]
physics_list_nuclear = [
    "DsPhysConsElectroNu",
    "DsPhysConsHadron",
    "DsPhysConsIon",
    ]


class Configure:
    '''
    Do default DetSim configuration.
    '''

    # Available geometry broken up by site
    giga_far_items = [
        "/dd/Structure/Sites/far-rock",
        "/dd/Geometry/AdDetails/AdSurfacesAll",
        "/dd/Geometry/AdDetails/AdSurfacesFar",
        "/dd/Geometry/PoolDetails/FarPoolSurfaces",
        "/dd/Geometry/PoolDetails/PoolSurfacesAll",
        ]
    giga_dayabay_items = [
        "/dd/Structure/Sites/db-rock",
        "/dd/Geometry/AdDetails/AdSurfacesAll",
        "/dd/Geometry/AdDetails/AdSurfacesNear",
        "/dd/Geometry/PoolDetails/NearPoolSurfaces",
        "/dd/Geometry/PoolDetails/PoolSurfacesAll",
        ]
    giga_lingao_items = [
        "/dd/Structure/Sites/la-rock",
        "/dd/Geometry/AdDetails/AdSurfacesAll",
        "/dd/Geometry/AdDetails/AdSurfacesNear",
        "/dd/Geometry/PoolDetails/NearPoolSurfaces",
        "/dd/Geometry/PoolDetails/PoolSurfacesAll",
        ]

    def __init__(self,site="far,dayabay,lingao", 
                 physlist = physics_list_basic+physics_list_nuclear,
                 use_push_algs = True,
                 use_sim_subseq=False ):
        '''
        Configure DetSim.  

        "site" can be "far", "dayabay" or "lingao".  Default is all three

        "physlist" specifies the physics lists.  You can use the
        predefined lists in DetSim.configure.: physics_list_basic and
        physics_list_nuclear.  Default is to use both.
        
        After creating this object you may want to call historian() or
        unobserver() to add to their configuration.

        '''

        from GaussTools.GaussToolsConf import GiGaPhysListModular
        import GaudiKernel.SystemOfUnits as units

        # Note: we must name this with "GiGa." as it is assumed later when the
        # properties are looked up.  Really the action of giving it to GiGa
        # should take care of this.  More bugs in Configurables
        physics_list = GiGaPhysListModular("GiGa.GiGaPhysListModular")
        physics_list.CutForElectron = 100*units.micrometer
        physics_list.CutForPositron = 100*units.micrometer
        physics_list.CutForGamma = 1*units.millimeter
        physics_list.PhysicsConstructors = physlist
        self.physics_list = physics_list

        from GiGa.GiGaConf import GiGa
        giga = GiGa()
        giga.PhysicsList = physics_list

        # Start empty step action sequence to hold historian/unobserver
        from GaussTools.GaussToolsConf import GiGaStepActionSequence
        sa = GiGaStepActionSequence('GiGa.GiGaStepActionSequence')
        giga.SteppingAction = sa

        self.giga = giga

        # Tell GiGa the size of the world.
        # Set default world material to be vacuum to speed propagation of
        # particles in regions of little interest.
        from GiGaCnv.GiGaCnvConf import GiGaGeo
        giga_geom = GiGaGeo()
        giga_geom.XsizeOfWorldVolume = 2.4*units.kilometer
        giga_geom.YsizeOfWorldVolume = 2.4*units.kilometer
        giga_geom.ZsizeOfWorldVolume = 2.4*units.kilometer
        giga_geom.WorldMaterial = "/dd/Materials/Vacuum"
        self.gigageo = giga_geom

        # Set up for telling GiGa what geometry to use, but don't
        # actually set that.
        from GaussTools.GaussToolsConf import GiGaInputStream
        giga_items = GiGaInputStream()
        giga_items.ExecuteOnce = True
        giga_items.ConversionSvcName = "GiGaGeo"
        giga_items.DataProviderSvcName = "DetectorDataSvc"
        giga_items.StreamItems = [ ]
        site = site.lower()
        if "far" in site: 
            giga_items.StreamItems += self.giga_far_items
        if "dayabay" in site: 
            giga_items.StreamItems += self.giga_dayabay_items
        if "lingao" in site: 
            giga_items.StreamItems += self.giga_lingao_items
        self.giga_items = giga_items

        # Make sequencer alg to run all this stuff as subalgs
        from GaudiAlg.GaudiAlgConf import GaudiSequencer
        giga_sequence = GaudiSequencer()
        giga_sequence.Members = [ self.giga_items ]
        self.giga_sequence=giga_sequence
        if use_push_algs:
            # DetSim's algs
            from DetSim.DetSimConf import DsPushKine, DsPullEvent
            self.detsim_push_kine = DsPushKine()
            self.detsim_pull_event = DsPullEvent()
            giga_sequence.Members += [self.detsim_push_kine,
                                      self.detsim_pull_event]
            pass
        
        if not use_sim_subseq:
            from Gaudi.Configuration import ApplicationMgr
            theApp = ApplicationMgr()
            theApp.TopAlg.append(giga_sequence)
        
        return

    def historian(self,trackSelection="",vertexSelection="",useFastMuEnergyCut=False):
        ''' Add the Historian to the simulation and configure its
        default track and vertex selection, default is empty.  Its
        configurable is returned in order to allow further, direct
        modification.
        '''
        from Historian.HistorianConf import HistorianStepAction
        # Must name stepping action with anticipation of its final ownership.
        # This works around some inconsistency with configurable tools
        hist = HistorianStepAction('GiGa.GiGaStepActionSequence.HistorianStepAction')
        hist.TrackSelection = trackSelection
        hist.VertexSelection = vertexSelection
        hist.UseFastMuEnergyCut = useFastMuEnergyCut

        self.giga.SteppingAction.Members.append(hist)
        return hist

    def unobserver(self,stats = []):
        ''' Add the Unobserver to the simulation and configure it with
        a set of stats (default is empty).  Its configurable is
        returned in order to allow further, direct modification.
        '''
        from Historian.HistorianConf import UnObserverStepAction
        # Must name stepping action with anticipation of its final ownership.
        # This works around some inconsistency with configurable tools
        unobs = UnObserverStepAction('GiGa.GiGaStepActionSequence.UnObserverStepAction')
        unobs.Stats = stats
        self.giga.SteppingAction.Members.append(unobs)
        return unobs
        
    
def configure(argv=None):
    try:
        style = argv[0]
    except IndexError:
        Configure()
        return

    if style == 'basic':
        Configure(physlist = physics_list_basic)
        return

    Configure()
    return

---