IBD.py

Go to the documentation of this file.

#!/usr/bin/env python

from GaudiKernel import SystemOfUnits as units

def EvtGenerator(name='IBD_gds',
                 seed = 42,
                 finalstate = 'FullIbd',
                 volume='/dd/Structure/AD/db-gds1',
                 material = 'GdDopedLS',
                 strategy = 'FullVolume',
                 mode = 'Uniform',
                 lifetime = 9*units.second,
                 wallTime=0):
    '''
    set up the kinematic generator in the given volume with chosen options
    using the InvBetaDecay generator.

    Originally created by Qing.
    Modified from MDC09b/runIBD15 by Zhe  Dec. 13, 2010
    '''
    
    print name+"'s strategy is: ", strategy
    print name+"'s seed is: ignored"
    print name+"'s volume is: ", volume
    print name+"'s Material is: ", material
        
    # set up each tools
    from GaudiKernel import SystemOfUnits as units
    from InvBetaDecay.Helpers import Decay
    
    decay=Decay(name=name)
    if finalstate == 'PositronOnly':
        decay.decay.PositronOnly = 1
    elif finalstate == 'NeutronOnly':
        decay.decay.NeutronOnly = 1
    decay.positioner.Volume = volume
    decay.positioner.Strategy = strategy
    if decay.positioner.Strategy == 'Material':
        #            decay.positioner.Volume='/dd/Structure/AD/db-ade1'
        decay.positioner.FillMaterials = [material]
    decay.positioner.Mode = mode
    decay.positioner.Position = [0,0,0]
    decay.positioner.Spread = 10*units.m
    decay.timerator.LifeTime = lifetime
    decay.transformer.Volume = volume
    
    # set up Gnrtr itself
    from Gnrtr.GnrtrConf import Gnrtr
    gnrtr = Gnrtr(name);
    gnrtr.TimeStamp = wallTime
    gnrtr.GenName = name
    
    gnrtr.GenTools = decay.tools()

    return gnrtr