Li9He8.py

Go to the documentation of this file.

#!/usr/bin/env python
import GaudiKernel.SystemOfUnits as units

def Li9He8(stage = 'null',
           name = 'Li9He8_gds',
           volume = '/dd/Structure/AD/db-gds1',
           strategy = "AvoidDaughters",
           material = "StainlessSteel",
           mode = "Uniform",
           Li9Fraction = 0.9,
           CompleteDecay = True,
           Lifetime = 1*units.second,
           start_time = 0):
    '''
    Add a Li9 generator
    '''
        
    # set up Gnrtr itself
    from Gnrtr.GnrtrConf import Gnrtr
    gnrtr = Gnrtr(name)
    
    gnrtr.GenTools = [ "Li9He8Decayerator/"+name+"Gen",                       
                       "GtPositionerTool/"+name+"Pos",
                       "GtTimeratorTool/"+name+"Tim",
                       "GtTransformTool/"+name+"Tra"]
    
    gnrtr.TimeStamp = start_time
    gnrtr.GenName = name
        
    # set up each tools
    from Li9He8Decay.Li9He8DecayConf import Li9He8Decayerator
    from GenTools.GenToolsConf import GtPositionerTool, GtTransformTool, GtTimeratorTool

    #generator
    gen=Li9He8Decayerator(name+'Gen')
    gen.Li9fraction = Li9Fraction
    gen.CompleteDecay = CompleteDecay
        
    # Set up positioner
    pos=GtPositionerTool(name+'Pos',Volume=volume)
    pos.Strategy = strategy 
    if pos.Strategy == 'VolumeType':
        pos.FillVolumes = [fillvolumes]
    if pos.Strategy == 'Material':
        pos.FillMaterials = [material]
    pos.Mode = mode
    pos.Spread = 25*units.m
    pos.Position = [0,0,0*units.m]

    # time
    tim = GtTimeratorTool(name+"Tim")
    tim.LifeTime = Lifetime

    # transform
    tra=GtTransformTool(name+'Tra',Volume=volume)

    return gnrtr