estimatePileup2 Namespace Reference

Classes
class	pileupParameters

Functions
def	fillPileupHistogram

Variables
string	action = 'store_true'
list	amodetagChoices = [ "PROTPHYS","IONPHYS" ]
list	beamChoices = ['PROTPHYS','IONPHYS']
list	beamModeChoices = [ "","stable", "quiet", "either"]
	## Main Program More...
list	bxlist = perlsdata[8]
list	bxlumiinfo = []
	convert lumiData to lumiDict format #{run:[[cmslsnum,avg]]} More...
list	cmslsnum = perlsdata[1]
tuple	dbGroup = optparse.OptionGroup(parser, "Database Options")
	debug = options.debugLumi)
string	default = 'frontier://LumiCalc/CMS_LUMI_PROD'
list	deliveredlumi = perlsdata[5]
tuple	events = open(options.csvInput, 'r')
tuple	f = open(options.inputfile, 'r')
tuple	fileparsingResult = csvSelectionParser.csvSelectionParser(options.inputfile)
	finecorrections = None
string	help = 'connect string to lumiDB ,default %default'
tuple	hist
tuple	histFile = ROOT.TFile.Open(output, 'recreate')
list	histList = []
tuple	inputfilecontent = f.read()
tuple	inputGroup = optparse.OptionGroup(parser, "Input Options")
dictionary	inputRange = {int(options.runnumber):None}
tuple	lumiData = lumiCalcAPI.lumiForRange(schema,inputRange,beamstatus=options.beamstatus,withBXInfo=True,bxAlgo=options.algoname,xingMinLum=options.xingMinLum,withBeamIntensity=False,datatag=None,finecorrections=finecorrections)
list	output = args[0]
tuple	parser
tuple	pieces = sepRE.split(line.strip())
tuple	pileupGroup = optparse.OptionGroup(parser, "Pileup Options")
tuple	pileupHist
list	recordedlumi = perlsdata[6]
tuple	rruns = inputRange.keys()
dictionary	runDict = {}
tuple	schema = session.nominalSchema()
tuple	sepRE = re.compile(r'[\s,;:]+')
tuple	session = svc.openSession(isReadOnly=True,cpp2sqltype=[('unsigned int','NUMBER(10)'),('unsigned long long','NUMBER(20)')])
tuple	svc = sessionManager.sessionManager(options.connect,authpath=options.authpath,debugON=options.debug)
string	VERSION = '2.00'
list	xingAlgoChoices = [ "OCC1","OCC2","ET"]
list	xingIdx = [int(myidx) for myidx in pieces[4::2] ]
list	xingVal = [float(myval) for myval in pieces[5::2] ]

Function Documentation

def estimatePileup2.fillPileupHistogram	(		bxlumiinfo,
			pileupHistName,
			maxPileupBin,
			runNumber = 0,
			hist = None,
			debug = False
	)

bxlumiinfo:[[cmslsnum(0),avgdelivered(1),avgrecorded(2),bxlumiarray[3]]]

Given luminfo , deadfrac info 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 19 of file estimatePileup2.py.

                         runNumber=0, hist = None, debug = False):
    '''
    bxlumiinfo:[[cmslsnum(0),avgdelivered(1),avgrecorded(2),bxlumiarray[3]]]

    Given luminfo , deadfrac info 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' % (pileupHistName, runNumber)
        hist = ROOT.TH1D (histname, histname, maxPileupBin + 1,
                          -0.5,maxPileupBin + 0.5)
        upper = maxPileupBin
    p=pileupParameters()
    for perlsinfo in bxlumiinfo:
        cmslsnum=perlsinfo[0]
        avgdelivered=perlsinfo[1]
        avgrecorded=perlsinfo[2]
        bxdata=perlsinfo[3]
        bxidx=bxdata[0]
        bxvaluelist=bxdata[1]
        #calculate livefrac
        livetime=1
        if avgrecorded<0:
            avgrecorded=0
        if avgdelivered:
            livetime=avgrecorded/avgdelivered
        else:
            livetime=0
        for idx,bxvalue in enumerate(bxvaluelist):
            xingIntLumi=bxvalue * p.lumiSectionLen * livetime
            if options.minBiasXsec:
                mean = bxvalue * options.minBiasXsec * p.rotationTime
            else:
                mean = bxvalue * p.minBiasXsec * p.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, bxvalue)
                else:
                    print "mean number of pileup events > 100 for lum %d: m %f l %f" % \
                          (cmslsnum, mean, bxvalue)
            totalProb = 0
            for obs in range (upper):
                prob = ROOT.TMath.Poisson (obs, mean)
                totalProb += prob
                hist.Fill (obs, prob * xingIntLumi)
            if debug:
                xing=bxidx[idx]
                print "ls", lumiSection, "xing", xing, "inst", bxvalue, \
                      "mean", mean, "totalProb", totalProb, 1 - totalProb
                print "  hist mean", hist.GetMean()
            if totalProb < 1:
                hist.Fill (obs, (1 - totalProb) * xingIntLumi)
    return hist