In This Package:

SimGe68noGEO.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
Script to simulate Ge68 events at the center without background sources and without the source geometry.

usage example:

  nuwa.py -n 2 -o output.root -m "Calibration.SimGe68wBG -z 0.0 -a A"

'''

from DybPython.DybPythonAlg import DybPythonAlg
from DetHelpers.DetHelpersConf import AutoPositionerTool

import os, math

def configure(argv = []):
    """Configure this module with source position"""
    
    import sys, getopt
    import GaudiKernel.SystemOfUnits as units
    from time import gmtime, mktime, strftime, strptime, timezone
    opts,args = getopt.getopt(argv,
                              "z:a:")
    from DybPython.Control import nuwa
    wallTime = nuwa.opts.time
    runNumber = nuwa.opts.run
    #wallTime = 0
    axis = 'A'
    acuName = 'ACU_A_Center'
    xpos = 0.0
    ypos = 0.0
    zpos = 0.0
    pmtDataPath = None
    for opt,arg in opts:
        if opt == "-z":
            zpos = float(arg) * units.cm
            print "======================================================"
            print "Source Z position = ", zpos / units.cm, " cm"
            print "======================================================"
        if opt == "-a":
            axis = arg
            print "======================================================"
            print "ACU Axis = ", axis
            print "======================================================"
            if axis == 'B':
                xpos = 135.0 * math.cos(112.5 * math.pi / 180.) * units.cm
                ypos = 135.0 * math.sin(112.5 * math.pi / 180.) * units.cm
                acuName = 'ACU_B_GdlsEdge'
            elif axis == 'C':
                xpos = 177.25*math.cos( (112.5 + 180) * math.pi / 180.)*units.cm
                ypos = 177.25*math.sin( (112.5 + 180) * math.pi / 180.)*units.cm
                acuName = 'ACU_C_GammaCatcher'
    
    # Set run info
    from RunDataSvc.RunDataSvcConf import RunDataSvc
    runDataSvc = RunDataSvc()
    runDataSvc.SimRunType = "Calibration"
    sourceName = "DayaBayAD1_"+acuName+"_Germanium_68"
    runDataSvc.SimCalibSources = [ sourceName ]
    runDataSvc.SimCalibZPosition = { sourceName : zpos }    

    # Add pull-mixing simulation
    from Stage import Configure as StageConfigure
    stageCfg = StageConfigure()
    stageCfg.addStages( ['Kinematic','Detector','Electronic','TrigRead',
                         'SingleLoader'] )

    # Add kinematic generators
    from Gnrtr.GnrtrConf import Gnrtr
    from GenTools.Helpers import Gun
    from GenTools.Helpers import HepEVT

    SimTime=nuwa.opts.executions * 2 * 0.020

    # Ge68
    
    seed = runNumber
    Ge68lifetime=0.020
    Ge68events= SimTime // Ge68lifetime
   
    mygenGe68 = HepEVT("Ge68.exe -n %d -seed %d |" % (Ge68events,seed), name = "Ge68")
    mygenGe68.positioner.Volume = "/dd/Structure/AD/db-oil1"
    mygenGe68.transformer.Volume = "/dd/Structure/AD/db-oil1"
    mygenGe68.positioner.Position = [xpos, ypos, zpos]
    mygenGe68.positioner.Strategy = "FullVolume"
    mygenGe68.positioner.Mode = "Fixed"
    mygenGe68.positioner.Spread = 3.0 * units.cm
    mygenGe68.timerator.LifeTime = Ge68lifetime*units.second

    gnrtrGe68 = Gnrtr("gnrtrGe68");
    gnrtrGe68.GenTools = mygenGe68.tools()
    gnrtrGe68.ThisStageName = "Kinematic"
    gnrtrGe68.TimeStamp = int(wallTime)
    stageCfg.KinematicSequence.Members.append(gnrtrGe68)
    
    # Rest of Simulation stages
    import DetSim
    detsim = DetSim.Configure(use_push_algs = False)
    detsim.historian(trackSelection="(pdg == -11)",\
                     vertexSelection="(pdg == -11)")
    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",
        'Positron': "pdg == -11",
        '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],
            ["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],
            ["tEnd_Trk1",    "t" , "%(track1)s and %(lastvtx)s"%params],
            ["xEnd_Trk1",    "x", "%(track1)s and %(lastvtx)s"%params],
            ["yEnd_Trk1",    "y", "%(track1)s and %(lastvtx)s"%params],
            ["zEnd_Trk1",    "z", "%(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",   "vx","%(track1)s and %(start)s"%params],
            ["vy_Trk1",   "vy","%(track1)s and %(start)s"%params],
            ["vz_Trk1",   "vz","%(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]
            ])

    from DetSimProc.DetSimProcConf import DetSimProc
    dsp = DetSimProc()
    dsp.ThisStageName = "Detector"
    dsp.LowerStageName = "Kinematic"
    stageCfg.DetectorSequence.Members.append(dsp)

    import ElecSim
    elecsim = ElecSim.Configure(use_push_algs = False)
    from ElecSimProc.ElecSimProcConf import ElecSimProc
    esp = ElecSimProc()
    esp.ThisStageName = "Electronic"
    esp.LowerStageName = "Detector"
    stageCfg.ElectronicSequence.Members.append(esp)
    
    #import TrigSim
    #trigsim = TrigSim.Configure()
    from TrigReadProc.TrigReadProcConf import TrigReadProc
    tsp = TrigReadProc()
    tsp.ThisStageName = "TrigRead"
    tsp.LowerStageName = "Electronic"
    stageCfg.TrigReadSequence.Members.append(tsp)

    #import ReadoutSim
    #rosim = ReadoutSim.Configure()
    from SingleLoader.SingleLoaderConf import SingleLoader
    sll = SingleLoader()
    sll.ThisStageName = "SingleLoader"
    sll.LowerStageName = "TrigRead"
    stageCfg.SingleLoaderSequence.Members.append(sll)

    from Stage.StageConf import Sim15
    sim15=Sim15()
    sim15.TopStage="SingleLoader"
        
    from Gaudi.Configuration import ApplicationMgr
    theApp = ApplicationMgr()
    theApp.TopAlg.append(sim15)
    
def run(app):
    pass

---