In This Package:

NeutronSim.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
usage example:

  nuwa.py -n 2 -o output.root -m "MDC09a.runNeutron -k 1.0"

  Neutron is generated uniformly in the AD
  "-k argument specify momentum of neutron in eV"
  
'''

def configure(argv = []):
    """Configure this module with neutron energy"""
    
    import sys, getopt
    opts,args = getopt.getopt(argv,"k:")
    neutronE = 100.0
    for opt,arg in opts:
        if opt == "-k":
            neutronE = float(arg)
            print "======================================================"
            print "Neutron energy = ", neutronE, " eV"
            print "======================================================"

    # Set run info
    from RunDataSvc.RunDataSvcConf import RunDataSvc
    runDataSvc = RunDataSvc()
    runDataSvc.SimRunType = "Physics"

    # Generator
    import GaudiKernel.SystemOfUnits as units
    from GenTools.Helpers import Gun
    from GenTools.GenToolsConf import GtGunGenTool, GtPositionerTool
    gun = Gun(name = "neutronGun",
              volume = "/dd/Structure/AD/db-oil1",
              gun = GtGunGenTool("gun",
                                 ParticlesPerEvent = 1,
                                 ParticleName = "n0",
                                 Momentum = neutronE*units.eV,
                                 MomentumMode = "Fixed",
                                 MomentumSpread = 0.*units.MeV,
                                 DirectionMode = "Uniform",
                                 Direction = [ 1, 0, 0 ],
                                 DirectionSpread = 3),
              positioner = GtPositionerTool("pos",
                                            Strategy = "FullVolume",
                                            Mode = "Uniform",
                                            Spread = 2.5*units.meter,
                                            Position = [0,0,0*units.meter])
              )
    gun.timerator.LifeTime = 0.1*units.second
    import GenTools
    GenTools.Configure().register(gun)
    
    # Detector Simulation
    import DetSim
    detsim = DetSim.Configure()
    detsim.historian(trackSelection="(pdg == 2112)",\
                     vertexSelection="(pdg == 2112)")
    params = {
        'start' :"(start > 0)",
        'track1':"(id==1)",
        'track2':"(id==2)",
        'GD':    "MaterialName == '/dd/Materials/GdDopedLS'",
        'LS':    "MaterialName == '/dd/Materials/LiquidScintillator'",
        'MO':   "MaterialName == '/dd/Materials/MineralOil'",
        'IAV':   "DetectorElementName == 'db-iav1'",
        'OAV':   "DetectorElementName == 'db-oav1'",
        'IWS': "MaterialName == '/dd/Materials/IwsWater'",
        'OWS': "MaterialName == '/dd/Materials/OwsWater'",
        'lastvtx': "IsStopping == 1",
        'firstvtx': "IsStarting == 1",
        'Neutron': "pdg == 2112",
        'NeutronMom': "creator == 2112",
        'Gamma': "pdg == 22",
        'Muon': "(pdg == 13 or pdg == -13)"
        }

    detsim.unobserver(stats=[
            ["EDepInGdLS", "dE", "%(GD)s"%params],
            ["EDepInLS", "dE", "%(LS)s"%params],
            ["EDepInIAV", "dE", "%(IAV)s"%params],
            ["EDepInOAV", "dE", "%(OAV)s"%params],
            ["EDepInOIL", "dE", "%(MO)s"%params],

            ["QEDepInGdLS", "qdE", "%(GD)s"%params],
            ["QEDepInLS", "qdE", "%(LS)s"%params],
            ["QEDepInIAV", "qdE", "%(IAV)s"%params],
            ["QEDepInOAV", "qdE", "%(OAV)s"%params],
            ["QEDepInOIL", "qdE", "%(MO)s"%params],

            ["tQESumGdLS", "qEt", "%(GD)s"%params],
            ["xQESumGdLS", "qEx", "%(GD)s"%params],
            ["yQESumGdLS", "qEy", "%(GD)s"%params],
            ["zQESumGdLS", "qEz", "%(GD)s"%params],

            ["tQESumLS", "qEt", "%(LS)s"%params],
            ["xQESumLS", "qEx", "%(LS)s"%params],
            ["yQESumLS", "qEy", "%(LS)s"%params],
            ["zQESumLS", "qEz", "%(LS)s"%params],

            ["tQESumMO", "qEt", "%(MO)s"%params],
            ["xQESumMO", "qEx", "%(MO)s"%params],
            ["yQESumMO", "qEy", "%(MO)s"%params],
            ["zQESumMO", "qEz", "%(MO)s"%params],

            # track 1
            ["pdgId_Trk1","pdg","%(track1)s and %(start)s"%params],
            ["t_Trk1",    "t" , "%(track1)s and %(start)s"%params],
            ["lx_Trk1",    "lx", "%(track1)s and %(start)s"%params],
            ["ly_Trk1",    "ly", "%(track1)s and %(start)s"%params],
            ["lz_Trk1",    "lz", "%(track1)s and %(start)s"%params],
            ["tEnd_Trk1",    "t" , "%(track1)s and %(lastvtx)s"%params],
            ["lxEnd_Trk1",    "lx", "%(track1)s and %(lastvtx)s"%params],
            ["lyEnd_Trk1",    "ly", "%(track1)s and %(lastvtx)s"%params],
            ["lzEnd_Trk1",    "lz", "%(track1)s and %(lastvtx)s"%params],
            ["e_Trk1",    "E",  "%(track1)s and %(start)s"%params],
            ["p_Trk1",    "p",  "%(track1)s and %(start)s"%params],
            ["ke_Trk1",   "KE", "%(track1)s and %(start)s"%params],
            ["vx_Trk1",   "lvx","%(track1)s and %(start)s"%params],
            ["vy_Trk1",   "lvy","%(track1)s and %(start)s"%params],
            ["vz_Trk1",   "lvz","%(track1)s and %(start)s"%params],
            ["TrkLength_GD_Trk1",  "dx","%(track1)s and %(GD)s"%params],
            ["TrkLength_iAV_Trk1", "dx","%(track1)s and %(IAV)s"%params],
            ["TrkLength_LS_Trk1",  "dx","%(track1)s and %(LS)s"%params],
            ["TrkLength_oAV_Trk1", "dx","%(track1)s and %(OAV)s"%params],
            ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(MO)s"%params],

            ["capTarget", "capTargetZ","%(track1)s and %(lastvtx)s"%params]
            
            ])
            
    import ElecSim
    elecsim = ElecSim.Configure()
    
    import TrigSim
    trigsim = TrigSim.Configure()

    import ReadoutSim
    rosim = ReadoutSim.Configure()

def run(app):
    pass

---