estimatePileup Namespace Reference

Functions
def	fillPileupHistogram

Variables
string	action = 'store_true'
list	beamModeChoices = ["","stable", "quiet", "either"]
dictionary	csvDict = {}
tuple	dbGroup = optparse.OptionGroup(parser, "Database Options")
list	deadTable = runDTarray[2]
	debug = options.debugLumi,
tuple	events = open(options.csvInput, 'r')
tuple	f = open(options.inputfile, 'r')
tuple	fileparsingResult = csvSelectionParser.csvSelectionParser(options.inputfile)
string	help = 'connect string to lumiDB (optional, default to frontier://LumiCalc/CMS_LUMI_PROD)'
tuple	hist
tuple	histFile = ROOT.TFile.Open(output, 'recreate')
list	histList = []
tuple	inputfilecontent = f.read()
tuple	inputGroup = optparse.OptionGroup(parser, "Input Options")
	inputRange = options.runnumber
	Get input source. More...
string	mode = 'csv'
list	output = args[0]
tuple	parameters = LumiQueryAPI.ParametersObject()
	## Main Program More...
tuple	parser
tuple	pieces = sepRE.split(line.strip())
tuple	pileupGroup = optparse.OptionGroup(parser, "Pileup Options")
tuple	pileupHist
	Save what we need in the parameters object. More...
tuple	recordedData = LumiQueryAPI.recordedLumiForRange(session, parameters, inputRange)
dictionary	runLumiDict = {}
list	runNumber = runDTarray[0]
	pprint (recordedData) More...
tuple	sepRE = re.compile(r'[\s,;:]+')
list	xingInstLumiArray

Function Documentation

def estimatePileup.fillPileupHistogram	(		deadTable,
			parameters,
			runNumber = 0,
			hist = None,
			debug = False,
			mode = 'deadtable'
	)

Given a deadtable and run number, will (create if necessary
and) fill histogram with expected pileup distribution.  If a
histogram is created, it is owned by the user and is his/her
responsibility to clean up the memory.

Definition at line 14 of file estimatePileup.py.

                         mode='deadtable'):
    '''Given a deadtable and run number, will (create if necessary
    and) fill histogram with expected pileup distribution.  If a
    histogram is created, it is owned by the user and is his/her
    responsibility to clean up the memory.'''
    if hist:
        maxBin = hist.GetNbinsX()
        upper = int( hist.GetBinLowEdge(maxBin) + \
                     hist.GetBinWidth(maxBin) + 0.25 )
    else:
        histname = '%s_%s' % (parameters.pileupHistName, runNumber)
        hist = ROOT.TH1D (histname, histname, parameters.maxPileupBin + 1,
                          -0.5, parameters.maxPileupBin + 0.5)
        upper  = parameters.maxPileupBin
    for lumiSection, deadArray in sorted (deadTable.iteritems()):
        if mode == 'csv':
            numerator     = float (deadArray[1])
            denominator   = float (deadArray[0])
            instLumiArray =        deadArray[2]
            livetime = 1
            if numerator < 0:
                numerator = 0
            if denominator:
                livetime = numerator / denominator
        else:
            # we got everything from lumiDB
            if len(deadArray) <= parameters.xingIndex:
                # for some reason the xing instantaneous luminosity
                # information isn't there.  Print a warning and then skip
                # it:
                if parameters.noWarnings:
                    continue
                if runNumber:
                    print "No Xing Instantaneous luminosity information for run %d, lumi section %d" \
                          % (runNumber, lumiSection)
                else:
                    print "No Xing Instantaneous luminosity information for lumi section %d" \
                          % lumiSection
                continue
            numerator   = float (deadArray[0])
            denominator = float (deadArray[2] * deadArray[4])
            xingInstLumiArray = deadArray[parameters.xingIndex]
            # here we only want the instantaneous luminosities and don't
            # care which crosings they fall in.  So we only want the odd 
            instLumiArray = [(xingInstLumiArray[index], xingInstLumiArray[index + 1]) \
                             for index in xrange( 0, len (xingInstLumiArray), 2 ) ]
            livetime = 1
            if numerator < 0:
                numerator = 0
            if denominator:
                livetime = 1 - numerator / denominator
        # totalInstLumi = reduce(lambda x, y: x+y, instLumiArray) # not needed
        for xing, xingInstLumi in instLumiArray:
            xingIntLumi = xingInstLumi * parameters.lumiSectionLen * livetime
            mean = xingInstLumi * parameters.minBiasXsec * \
                   parameters.rotationTime
            if mean > 100:
                if runNumber:
                    print "mean number of pileup events > 100 for run %d, lum %d : m %f l %f" % \
                          (runNumber, lumiSection, mean, xingInstLumi)
                else:
                    print "mean number of pileup events > 100 for lum %d: m %f l %f" % \
                          (lumiSection, mean, xingInstLumi)

            for obs in range (upper):
                prob = ROOT.TMath.Poisson (obs, mean)
                totalProb += prob
                hist.Fill (obs, prob * xingIntLumi)
            if debug:
                print "ls", lumiSection, "xing", xing, "inst", xingInstLumi, \
                      "mean", mean, "totalProb", totalProb, 1 - totalProb
                print "  hist mean", hist.GetMean()
            if totalProb < 1:
                hist.Fill (obs, (1 - totalProb) * xingIntLumi)
    return hist