In This Package:

OpStack.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
usage example:

  nuwa.py  -n NEvents -o output.root -m " Detsim.OpStack Whichsite MuonRandomseed"

'''

DATETIME_FORMAT = '%Y-%m-%dT%H:%M:%S'

import os, math

test=False

def configure(argv=[]):
    """ script for fast neutron background study"""

    site=argv[0]
    seed=argv[1]

    volume = "/dd/Structure/Pool/db-ows"
    musicsite = 'DYB'
    if 'far' == site:
        volume = "/dd/Structure/Pool/far-ows"
        musicsite = 'Far'
    if 'dayabay' == site:
        volume = "/dd/Structure/Pool/db-ows"
        musicsite = 'DYB'
    if 'lingao' == site:
        volume = "/dd/Structure/Pool/la-ows"
        musicsite = 'LA'

    import GaudiKernel.SystemOfUnits as units
    
    if test:
        import GenTools
        from GenTools.Helpers import Gun
        mygun = Gun()
        mygun.gun.ParticleName = 'mu-'
        mygun.gun.Momentum = 3.0*units.GeV
        mygun.setVolume(volume)
        mygun.positioner.Position = [300*units.cm, 300*units.cm, 0*units.cm]
        mygun.timerator.LifeTime = 1*units.second
        gtc = GenTools.Configure()
        gtc.register(mygun)

    else:
        MuonDataPath = os.getenv('MuonDataPath')
        if MuonDataPath is None :
            print "DetSim.OpStack Muon data path ($MuonDataPath) is not defined."
            print "Might need to get data file  from http://dayabay.ihep.ac.cn/svn/dybsvn/data/trunk/NewMuonGenerator/data/"
            sys.exit()
        else:
            print "Read muon data from ",MuonDataPath


        source = "Muon.exe -n 1000 -s %s -seed %s -r Yes -music_dir %s |"%(musicsite, seed, MuonDataPath)

        print " musicSite: %s, and RandomSeed: %s"%(musicsite,seed)
        
        from GenTools.Helpers import HepEVT
        hepevt = HepEVT(source)
        hepevt.positioner.Volume = volume
        hepevt.positioner.Mode = "Relative"
        hepevt.positioner.Position = [0,0,0]
        hepevt.timerator.LifeTime = 1*units.second
        hepevt.transformer.Volume = volume
        hepevt.transformer.Offset = [0., 0., (0.042)*units.meter]

        import GenTools
        gt = GenTools.Configure(helper=hepevt)
    

    ###Detsim Config
    
    import DetSim
    detsim = DetSim.Configure(site=site)

    from DetSim.DetSimConf import DsPhysConsOptical
    optical = DsPhysConsOptical("GiGa.GiGaPhysListModular.DsPhysConsOptical")
    optical.CerenPhotonScaleWeight = 3.0
    optical.ScintPhotonScaleWeight = 3.0
    optical.UseScintillation = False


    from DetSim.DetSimConf import DsOpStackAction
    opstack = DsOpStackAction("GiGa.DsOpStackAction")
    detsim.giga.StackingAction=opstack
    opstack.TightCut = False
    opstack.PhotonCut = False
    opstack.MaxPhoton = 1e6

    
### record all the muon and neutron track and vertices
#    ts="(pdg==13 or pdg==-13) or (pdg==2112)" 
#    vs = "(MaterialName == '/dd/Materials/MineralOil' or MaterialName == '/dd/Materials/LiquidScintillator') or (MaterialName == '/dd/Materials/GdDopedLS' or MaterialName == '/dd/Materials/Acrylic')" 

    ts="(pdg!=20022)"  ## record all non-optical photon

    detsim.historian(trackSelection=ts,vertexSelection= ts )
    
#    ismuon = "(pdg==13 or pdg==-13) and "
#    isneutron = "(pdg==2112) and "
#    material = "(MaterialName == '/dd/Materials/%s')"

    params = {
        'ismuon':"(pdg==13 or pdg==-13)",
        'isneutron': "(pdg==2112)",
        'isstop': "(IsStopping>0)",
        'start' :"(start > 0)",
        'track1':"(id==1 and ProcessType==1)",
        'track2':"(id==2 and ProcessType==1)",
        'GD':    "MaterialName == '/dd/Materials/GdDopedLS'",
        'LS':    "MaterialName == '/dd/Materials/LiquidScintillator'",
        'Oil':   "MaterialName == '/dd/Materials/MineralOil'",
        'AV':   "MaterialName == '/dd/Materials/Acrylic'",
        'OW':   "MaterialName == '/dd/Materials/OwsWater'",
        'IW':   "MaterialName == '/dd/Materials/IwsWater'"
        }

    stat = [
        # for track 1
        ["pdgId_Trk1","pdg","%(track1)s and %(start)s"%params],
        ["t_Trk1",    "t" , "%(track1)s and %(start)s"%params],
        ["x_Trk1",    "x", "%(track1)s and %(start)s"%params],
        ["y_Trk1",    "y", "%(track1)s and %(start)s"%params],
        ["z_Trk1",    "z", "%(track1)s and %(start)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_LS_Trk1",  "dx","%(track1)s and %(LS)s"%params],
        ["TrkLength_AV_Trk1", "dx","%(track1)s and %(AV)s"%params],
        ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(Oil)s"%params],
        # for track 2 now
        ["pdgId_Trk2","pdg","%(track2)s and %(start)s"%params],
        ["t_Trk2",    "t" , "%(track2)s and %(start)s"%params],
        ["x_Trk2",    "x", "%(track2)s and %(start)s"%params],
        ["y_Trk2",    "y", "%(track2)s and %(start)s"%params],
        ["z_Trk2",    "z", "%(track2)s and %(start)s"%params],
        ["e_Trk2",    "E",  "%(track2)s and %(start)s"%params],
        ["p_Trk2",    "p",  "%(track2)s and %(start)s"%params],
        ["ke_Trk2",   "KE", "%(track2)s and %(start)s"%params],
        ["vx_Trk2",   "lvx","%(track2)s and %(start)s"%params],
        ["vy_Trk2",   "lvy","%(track2)s and %(start)s"%params],
        ["vz_Trk2",   "lvz","%(track2)s and %(start)s"%params],
        ["TrkLength_GD_Trk2",  "dx","%(track2)s and %(GD)s"%params],
        ["TrkLength_LS_Trk2",  "dx","%(track2)s and %(LS)s"%params],
        ["TrkLength_AV_Trk2", "dx","%(track2)s and %(AV)s"%params],
        ["TrkLength_Oil_Trk2", "dx","%(track2)s and %(Oil)s"%params],
        ["MuonTrkLengthInOws","dx","%(ismuon)s and %(OW)s"%params],
        ["MuonTrkLengthInIws", "dx","%(ismuon)s and %(IW)s"%params],
        ["MuonTrkLengthInLS", "dx", "%(ismuon)s and %(LS)s"%params],
        ["MuonTrkLengthInGdLS","dx","%(ismuon)s and %(GD)s"%params],
        ["MuonTrkLengthInOil","dx", "%(ismuon)s and %(Oil)s"%params],
        ["MuonTrkLengthInAV","dx", "%(ismuon)s and %(AV)s"%params],
        ["MuonStop", "dx", "%(ismuon)s and %(isstop)s"%params],
        ["NeutronTrkLengthInOws","dx","%(isneutron)s and %(OW)s"%params],
        ["NeutronTrkLengthInIws", "dx", "%(isneutron)s and %(IW)s"%params],
        ["NeutronTrkLengthInLS", "dx", "%(isneutron)s and %(LS)s"%params],
        ["NeutronTrkLengthInGdLS","dx", "%(isneutron)s and %(GD)s"%params],
        ["NeutronTrkLengthInOil","dx", "%(isneutron)s and %(Oil)s"%params],
        ["NeutronTrkLengthInAV","dx", "%(isneutron)s and %(AV)s"%params],
        ["NeutronStop", "dx", "%(isneutron)s and %(isstop)s"%params]
        ]

    detsim.unobserver(stats=stat)
    

#    from GiGa.GiGaConf import GiGa
#    giga = GiGa()
#    giga.StackingAction = "DsOpStackAction"


    return
    
def run(app):
    pass


#if '__main__' == __name__:
#    pass

---