In This Package:

__init__.py

Go to the documentation of this file.

#!/usr/bin/env python

__all__ = ['SimHistsExampleConf','Configure', 'Tutorial']

import GaudiKernel.SystemOfUnits as units


class Configure:
    '''
    Do default configuration for SimHistsExample package
    '''

    def __init__(self,do_detsim = True, volume = '/dd/Structure/AD/far-oil1', particle = 'e+'):

        print "Configuring GenTools..."
        import GenTools
        from GenTools.Helpers import Gun
        gtc = GenTools.Configure()
        gunner = Gun()
        gtc.register(gunner)

        # Non-default values for the gen tools:

        gunner.timerator.LifeTime = 1*units.microsecond

        gunner.positioner.Strategy = "FullVolume"
        gunner.positioner.Volume = volume
        gunner.positioner.Mode = "Smeared"
        gunner.positioner.Spread = 1*units.meter
        gunner.positioner.Position = [0,0,0]

        gunner.gun.ParticlesPerEvent = 1
        gunner.gun.ParticleName = particle
        gunner.gun.Momentum = 4*units.MeV
        gunner.gun.MomentumMode = "Smeared"
        gunner.gun.MomentumSpread = 1*units.MeV
        gunner.gun.DirectionMode = "Uniform"
        gunner.gun.DirectionSpread = 180*units.degree
        gunner.gun.PolarizeMode = "Random"

        # Add the dumper to see log output of kinematics
        dump = GenTools.Dumper()
        print "... GenTools done."

        print 'Configure the THistSvc...'
        from GaudiSvc.GaudiSvcConf import THistSvc
        histsvc = THistSvc()
        histsvc.Output = [
            "file1 DATAFILE='simhists.root' OPT='RECREATE' TYP='ROOT' "
            ]
        print '... THistSvc config done.'

        # Get the app so we can add to TopAlg
        from Gaudi.Configuration import ApplicationMgr
        theApp = ApplicationMgr()

        print 'Configure the GenHists alg...'
        from SimHistsExample.SimHistsExampleConf import GenHists
        gh = GenHists()
        theApp.TopAlg.append(gh)
        gh.Location = "/Event/Gen/GenHeader"
        gh.Volume = volume
        gh.FilePath = "/file1/gen"
        gh.MaxEnergy = 5
        gh.EnergyUnits = units.MeV
        print '... GenHists alg configured.'

        if not do_detsim:
            print "Not doing DetSim configuration, only GenTools level."
            return

        print "Configuring DetSim with only basic physics list..."
        import DetSim
        from DetSim.Default import physics_list_basic
        detsim = DetSim.Configure(physlist = physics_list_basic)
        print '... done.'

        print 'Configure the SimHists alg...'
        from SimHistsExample.SimHistsExampleConf import SimHists
        sh = SimHists()
        theApp.TopAlg.append(sh)
        sh.Location = "/Event/Sim/SimHeader"
        sh.FilePath = "/file1/sim"
        print '... SimHists alg configured.'


        return
    pass

---