IBD.py

Go to the documentation of this file.

#!/usr/bin/env python
'''
  Configure background generators 

  Examples for how to use.
  For Fifteen package usage:

    from SourceRate.IBD import dyb_ibd_lso
    ibd_lso =dyb_ibd_lso()

    from Gnrtr.GnrtrConf import Gnrtr
    gnrtrIBDlso = Gnrtr("gnrtrIBDlso");
    gnrtrIBDlso.GenTools = ibd_lso.tools()
    gnrtrIBDlso.ThisStageName = "Kinematic"
    gnrtrIBDlso.TimeStamp = int(wallTime)
    stageCfg.KinematicSequence.Members.append(gnrtrIBDlso)

  For Single usage:

    from SourceRate.IBD import dyb_ibd_lso
    ibd_lso=dyb_ibd_lso()

    import GenTools
    gtc = GenTools.Configure(genname="IBD_lso")
    gtc.generator.TimeStamp = int(wallTime)
    gtc.register(ibd_lso)

'''
from GaudiKernel import SystemOfUnits as units

class EvtGenerator:
    '''
    set up the kinematic generator in the given volume with chosen options
    using the InvBetaDecay generator.
    '''

    def __init__(self,
                 name="IBD",
                 seed = 99999,
                 finalstate = 'FullIbd',
                 volume='/dd/Structure/AD/db-gds1',
                 material = 'GdDopedLS',
                 strategy = 'FullVolume',
                 mode = 'Uniform',
                 lifetime = 9*units.second):

        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)
        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.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
        
        self.ibd = decay

    def tools(self):
        return self.ibd.tools()
    
    pass

    
#The following rates are based on daya bay site
class dyb_ibd_oil:
    ''' IBD in oil '''

    def __init__(self,
                 seed=99999):    

        ibd = EvtGenerator(name     = 'IBD_oil',
                           seed     = seed,
                           volume   = '/dd/Structure/AD/db-oil1',
                           strategy = 'AvoidDaughters',
                           #strategy = 'Material',
                           material = 'MineralOil',
                           mode     = 'Uniform',
                           lifetime = 36.3*units.second #daya bay site
                          )
        
        self.ibd = ibd

    def tools(self):
        return self.ibd.tools()
    
    pass

class dyb_ibd_gds:
    ''' IBD in gds '''

    def __init__(self,
                 seed=99999):    

        ibd = EvtGenerator(name     = 'IBD_gds',
                           seed     = seed,
                           volume   = '/dd/Structure/AD/db-oil1',
                           strategy = 'Material',
                           material = 'GdDopedLS',
                           mode     = 'Uniform',
                           lifetime = 78.4*units.second, #daya bay site
                           )
        
        self.ibd = ibd

    def tools(self):
        return self.ibd.tools()
    
    pass

class dyb_ibd_lso:
    ''' IBD in lso '''

    def __init__(self,
                 seed=99999):    

        ibd = EvtGenerator(name     = 'IBD_lso',
                           seed     = seed,
                           volume   = '/dd/Structure/AD/db-oil1',
                           strategy = 'Material',
                           material = 'LiquidScintillator',
                           mode     = 'Uniform',
                           lifetime = 73.0*units.second #daya bay site
                           )
        
        self.ibd = ibd

    def tools(self):
        return self.ibd.tools()
    
    pass

class dyb_ibd_acrylic:
    ''' IBD in acrylic '''

    def __init__(self,
                 seed=99999):    

        ibd = EvtGenerator(name     = 'IBD_acrylic',
                           seed     = seed,
                           volume   = '/dd/Structure/AD/db-oil1',
                           strategy = 'Material',
                           material = 'Acrylic',
                           mode     = 'Uniform',
                           lifetime = 901.6*units.second #daya bay site
                           )
        
        self.ibd = ibd

    def tools(self):
        return self.ibd.tools()
    
    pass