pmtSimData.py

Go to the documentation of this file.

#!/usr/bin/env python
#
# Generate a random table of pmt properties
# for use in the electronics simulation

import random

import makeFeeCableMap

def makePmtProperties():
    "Return a list of PMTs and their properties"
    properties = []
    pmtIds = makeFeeCableMap.getAllPmtIds()
    #print "nDetectors = ",len(pmtIds)
    detectors = pmtIds.keys()
    detectors.sort()
    for detector in detectors:
        #print "Detector",detector.detName(),"nPmts = ",len(pmtIds[detector])
        for pmtId in pmtIds[detector]:
            properties.append( getProperties(pmtId) )
    return properties

def getProperties(pmtId):
    "Generate simulation properties for this PMT"
    gain           = -1
    sigmaGain      = -1
    efficiency     = -1
    timeSpread     = -1
    timeOffset     = -1
    prePulseProb   = -1
    afterPulseProb = -1
    darkRate       = -1

    if pmtId.detName().find("AD") > 0 or pmtId.detName().find("WS") > 0:
        # Add PMT variation in AD
        while gain <= 0:
            gain = random.gauss(1, 0.06)
        
        while sigmaGain <= 0:
            sigmaGain = random.gauss(0.30, 0.02) 

        while timeOffset <= 0:
            timeOffset = random.gauss(40.0, 3.0)

        while timeSpread <= 0:
            timeSpread = random.gauss(2.0, 0.2)

        while efficiency <= 0 or efficiency > 1.0:
            efficiency = random.gauss(0.90, 0.05)
    
        while prePulseProb <= 0 or prePulseProb >=1.0:
            prePulseProb = random.gauss(0.004, 0.0004)

        while afterPulseProb <= 0 or afterPulseProb >=1.0:
            afterPulseProb = random.gauss(0.016, 0.0016)
        
        while darkRate <= 0:
            darkRate = random.gauss(3000., 1000.)

        description = makeFeeCableMap.pmtDescription(pmtId)
        format = "%-2d %30s %-4.3f %-4.3f %-4.2f %-4.2f %-4.3f %-4.5f %-4.4f %-4.1f"
        line = format  % (pmtId.fullPackedData(),
                          description,
                          gain,
                          sigmaGain,
                          timeOffset,
                          timeSpread,
                          efficiency,
                          prePulseProb,
                          afterPulseProb,
                          darkRate)
        #print line
        return line
    else:
        return ""

def printPmtProperties():
    "Print table of PMT properties"
    print "# Simulation input data for pmt properties"
    print "# PmtId description gain sigmaGain tOffset tSpread effic prePuls afterPuls darkRate"
    lines = makePmtProperties()
    for line in lines:
        print line

if __name__ == "__main__":
    printPmtProperties()