Functions
def	CalcPileup (deadTable, parameters, mode='deadtable')

Variables
	args

	aveLumi

	csvDict

	default

	delivered

	description

	dest

	events

	GapDict

	help

	InGap

	inputGroup

	LastDelivered

	lastindex

	lastLumiS

	LastValidLumi

	lumi

	LumiArray

	LumiString

	mode

	OldRun

	options

	output

	outputfile

	OUTPUTLINE
	Let's start the fun. More...

	parameters
	## Main Program More...

	parser

	peakratio

	pieces

	pileup

	pileupGroup

	record

	recorded

	run

	runLumiDict

	sepRE

	trunc

	type

	xingInstLumiArray

Function Documentation

def estimatePileup_makeJSON.CalcPileup	(	deadTable,
		parameters,
		mode = `'deadtable'`
	)

Given a deadtable, will calculate parameters of pileup distribution. Return formatted
string with LumiSection, LS integrated lumi, RMS of bunch to bunch lumi and pileup.

Definition at line 15 of file estimatePileup_makeJSON.py.

References edm.print(), and mathSSE.sqrt().

 def CalcPileup (deadTable, parameters, mode='deadtable'):
     '''Given a deadtable, will calculate parameters of pileup distribution. Return formatted
     string with LumiSection, LS integrated lumi, RMS of bunch to bunch lumi and pileup.'''
 
     LumiString = ""
     LumiArray = []
 
     for lumiSection, deadArray in sorted (six.iteritems(deadTable)):
         numerator = 0
         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:
             print("no csv input! Doh!")
             return
         # totalInstLumi = reduce(lambda x, y: x+y, instLumiArray) # not needed
         if lumiSection > 0:
              TotalLumi = 0 
              TotalInt = 0
              TotalInt2 = 0
              TotalWeight = 0
              TotalWeight2 = 0
              FilledXings = 0
              for xing, xingInstLumi in instLumiArray:
                 xingIntLumi = xingInstLumi * parameters.lumiSectionLen * livetime
                 mean = xingInstLumi * 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))
 #                print "mean number of pileup events for lum %d: m %f idx %d l %f" % (lumiSection, mean, xing, xingIntLumi)
 
                 if xingInstLumi > 0.1:
 
                     TotalLumi = TotalLumi+xingIntLumi
                     TotalInt+= mean*xingIntLumi
                     FilledXings = FilledXings+1
                     #print "xing inst lumi %f %f %d" % (xingIntLumi,TotalLumi,FilledXings)
 
              #compute weighted mean, then loop again to get weighted RMS       
              MeanInt = 0
              if TotalLumi >0:
                  MeanInt = TotalInt/TotalLumi
              for xing, xingInstLumi in instLumiArray:
                  if xingInstLumi > 0.1:
                      xingIntLumi = xingInstLumi * parameters.lumiSectionLen * livetime
                      mean = xingInstLumi * parameters.rotationTime
                      TotalInt2+= xingIntLumi*(mean-MeanInt)*(mean-MeanInt)
                      TotalWeight+= xingIntLumi
                      TotalWeight2+= xingIntLumi*xingIntLumi
 
 
 
         if ((lumiSection > 0)):
             #print " LS, Total lumi, filled xings %d, %f, %d" %(lumiSection,TotalLumi,FilledXings)
             if FilledXings > 0:
                 AveLumi = TotalLumi/FilledXings
             else:
                 AveLumi = 0
             RMSLumi = 0
             Denom = TotalWeight*TotalWeight-TotalWeight2
             if TotalLumi > 0 and Denom > 0:
                 RMSLumi = sqrt(TotalWeight/(TotalWeight*TotalWeight-TotalWeight2)*TotalInt2)
             LumiString = "[%d,%2.4e,%2.4e,%2.4e]," % (lumiSection, TotalLumi, RMSLumi, MeanInt)
             LumiArray.append(lumiSection)
             LumiArray.append(TotalLumi)  # should really weight by total lumi in LS
             LumiArray.append(RMSLumi)
             LumiArray.append(MeanInt)
             
     return LumiArray
     
 
 