Radioact.py

Go to the documentation of this file.

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

def Radioact(stage = 'null',
             name = 'K40_gds',
             volume = '/dd/Structure/AD/db-gds1',
             nuclide = "K40",
             abundance = 3.01e17,
             strategy = "AvoidDaughters",
             material = "StainlessSteel",
             fillvolumes = "lvPmtHemiVacuum",
             start_time = 0):
    '''
    Configure GenDecay generator
    Originally created by Qing.
    Modified from MDC09b/runIBD15 by Zhe  Dec. 13, 2010
    '''
        
    # set up Gnrtr itself
    from Gnrtr.GnrtrConf import Gnrtr
    gnrtr = Gnrtr(name)
    
    gnrtr.GenTools = [ "GtDecayerator/"+name+"Gen",
                       "GtPositionerTool/"+name+"Pos",
                       "GtTransformTool/"+name+"Tra"]
    
    gnrtr.TimeStamp = start_time
    gnrtr.GenName = name
        
    # set up each tools
    from GenDecay.GenDecayConf import GtDecayerator
    from GenTools.GenToolsConf import GtPositionerTool, GtTransformTool

    #generator
    gen=GtDecayerator(name+'Gen')
    gen.ParentNuclide = nuclide
    gen.ParentAbundance = abundance
    gen.SecularEquilibrium = True
    gen.CorrelationTime = 1*units.second
    
    # 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 = "Uniform"
    pos.Spread = 25*units.m
    pos.Position = [0,0,0*units.m]

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

    return gnrtr