Co60.py

Go to the documentation of this file.

#!/usr/bin/env python

def Co60(volume):

    import GaudiPython as gm
    app = gm.AppMgr()            
    #msg = self.app.service("MessageSvc")
    #msg.OutputLevel = 1

    #processor
    gen = app.algorithm("Co60")
    #gen.OutputLevel = 2
    gen.GenTools = [ "GtHepEvtGenTool/Co60Gen", "GtPositionerTool/Co60Pos", "GtTimeratorTool/Co60Tim" ]
    gen.GenName = "Bang Bang"
    gen.Location = "/Event/Gen/HepMCEvents" # this is default anyways

    #
    from GaudiKernel import SystemOfUnits as units

    # Set up Co60Gen
    seed=1234567
    nevents=1000
    hepevt_source="Co60.exe -seed %(seed)s -n %(nevents)s|"
    # Deal with executable
    if hepevt_source[-1] == "|":
        exe = hepevt_source.split(' ')[0]
        if exe[0] != '/':   # Try to find full path
            import os, os.path
            path = os.getenv('PATH')
            for p in path:
                if (os.path.isfile(path+"/"+exe)):
                    exe = path+"/"+exe
                    break
                continue
            pass
        source = exe + ' ' + ' '.join(hepevt_source.split(' ')[1:])
        if "%" in source:   # Fill in any placemarks
            source = source%{'nevents':str(nevents),
                             'seed':str(seed)}
            pass
        pass
    
    
    app.property('ToolSvc.Co60Gen').OutputLevel = 2
    app.property('ToolSvc.Co60Gen').HepEvtDataSource = source

    # Set up positioner
    app.property('ToolSvc.Co60Pos').OutputLevel = 2
    app.property('ToolSvc.Co60Pos').Volume = volume
    app.property('ToolSvc.Co60Pos').Strategy = "FullVolume" # also "AvoidDaughters" and "Surface"
    app.property('ToolSvc.Co60Pos').Mode = "Uniform"
    app.property('ToolSvc.Co60Pos').Spread = 10*units.cm
    app.property('ToolSvc.Co60Pos').Position = [0,0,0]
    
    # Set up timerator
    app.property('ToolSvc.Co60Tim').LifeTime = int(5*units.ns)
    
    
    return