In This Package:

__init__.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
usage example:

  nuwa.py -A -n 2 -o output.root -m "ElecTutorial.GenKinematics -i 3500 -n DayaBayAD1_ACU_A_Center_LED -z 400.0"

  The LED names are:
  [Site][Det]_[Location]_LED

  Examples:
  DayaBayAD1_ACU_A_Center_LED
  DayaBayAD1_ACU_B_GdlsEdge_LED
  DayaBayAD1_ACU_C_GammaCatcher_LED
  DayaBayAD1_AD_WallUpper_LED
  DayaBayAD1_AD_WallCenter_LED
  DayaBayAD1_AD_WallLower_LED
  '''

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

import os, math

def configure(argv = []):
    """Configure this module with LED position"""

    import GaudiKernel.SystemOfUnits as units
    import sys, getopt
    from time import gmtime, mktime, strftime, strptime, timezone
    opts,args = getopt.getopt(argv,
                              "p:w:z:n:i:f:m:s:",["no-geometry"])
    wallTime = 0
    axis = 0
    xpos = 0.0
    ypos = 0.0
    zpos = 0.0
    intensity = 35000 # photons per LED flash
    intensityMode = "Gaus" # use Gaussian distribution to smear LED intensity
    intensitySpread = 350 # sigma/FWHM (for Gaus/Lorentz distribution, respectively) of intensity smearing in photons
    ledName = "DayaBayAD1_ACU_A_Center_LED"
    ledFrequency = 500.0 * units.hertz
    pmtDataPath = None
    placeGeometry = True
    for opt,arg in opts:
        if opt == "-p":
            pmtDataPath = arg
        if opt == "-i":
            intensity = int(arg)
            print "======================================================"
            print "Diffuser Ball intensity  = ",intensity," photons per flash"
            print "======================================================"
        if opt == "-m":
            intensityMode = arg
            print "======================================================"
            print "Diffuser Ball intensity is spread using a ",intensityMode," distribution"
            print "======================================================"
        if opt == "-s":
            intensitySpread = int(arg)
            print "======================================================"
            print "Diffuser Ball intensity is spread with sigma/FWHM = ",intensitySpread," photons"
            print "======================================================"
        if opt == "-f":
            ledFrequency = float(arg)
            print "======================================================"
            print "Diffuser Ball frequency  = ",ledFrequency," hertz"
            print "======================================================"
        if opt == "-z":
            zpos = float(arg)
            print "======================================================"
            print "Diffuser Ball Z position = ", zpos, " cm"
            print "======================================================"
            zpos *= units.cm
        if opt == "-n":
            ledName = arg
            print "======================================================"
            print "LED Source = ", ledName
            print "======================================================"
        if opt == "--no-geometry":
            placeGeometry = False
            print "======================================================"
            print "Photons only; ACU diffuser ball geometry will not be added."
            print "======================================================"

    detectorMap = {"DayaBayAD1":"/dd/Structure/AD/db-oil1",
                   "DayaBayAD2":"/dd/Structure/AD/db-oil2",
                   "LingAoAD1":"/dd/Structure/AD/la-oil1",
                   "LingAoAD2":"/dd/Structure/AD/la-oil2",
                   "FarAD1":"/dd/Structure/AD/far-oil1",
                   "FarAD2":"/dd/Structure/AD/far-oil2",
                   "FarAD3":"/dd/Structure/AD/far-oil3",
                   "FarAD4":"/dd/Structure/AD/far-oil4"
                   }

    # Get LED Source information
    from GaudiPython import gbl
    CalibSourceId = gbl.DayaBay.CalibSourceId
    ledId = CalibSourceId( ledName )
    xpos = ledId.nominalX()
    ypos = ledId.nominalY()
    volume = None
    if detectorMap.has_key( ledId.detName() ):
        # Set the correct coordinate system volume
        volume = detectorMap[ ledId.detName() ]
    else:
        print "ERROR: Unknown detector: ", ledId.detName()
        sys.exit(1)

    if ledId.isWallMounted():
        zpos = ledId.nominalZ()
        placeGeometry = False

    import GenTools
    from GenTools.Helpers import DiffuserBall
    ball = DiffuserBall( name=ledName,
                         useGeometry = placeGeometry )
    ball.setVolume( volume )
    ball.setPosition( [xpos, ypos, zpos] )
    ball.ball.PhotonsPerEvent = intensity
    ball.ball.PhotonsPerEventMode = intensityMode
    ball.ball.PhotonsPerEventSpread = intensitySpread
    ball.timerator.LifeTime = 1.0 / ledFrequency
    gtc = GenTools.Configure()
    gtc.generator.TimeStamp = int(wallTime)
    gtc.register(ball)

    # Set run info
    from RunDataSvc.RunDataSvcConf import RunDataSvc
    runDataSvc = RunDataSvc()
    runDataSvc.SimRunType = "Calibration"
    runDataSvc.SimCalibSources = [ ledName ]
    runDataSvc.SimLedFrequency = { ledName : ledFrequency }
    runDataSvc.SimCalibZPosition = { ledName : zpos }    

    import DetSim
    detsim = DetSim.Configure(physlist = DetSim.physics_list_basic)

    #import ElecSim
    #elecsim = ElecSim.Configure()
    #if pmtDataPath != None:
    #    # change PMT properties
    #    elecsim.dataSvc.setPmtSimData( pmtDataPath )

    # Don't save the information for every generated optical photon;
    # the data volume is too large.
    
    from DybAlg.DybAlgConf import DybGenPrune
    genPrune = DybGenPrune()
    from Gaudi.Configuration import ApplicationMgr
    appMgr = ApplicationMgr()
    appMgr.TopAlg += [genPrune]

def run(app):
    pass

---