estimatePileup Namespace Reference

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

Variables
	action
	args
	basename
	beammode
	beamModeChoices
	choices
	csvDict
	dbGroup
	deadTable
	debug
	default
	delivered
	description
	dest
	events
	extension
	f
	fileparsingResult
	help
	hist
	histFile
	histList
	inputfilecontent
	inputGroup
	inputRange
	Get input source. More...
	lumi
	lumiversion
	lumiXing
	maxPileupBin
	minBiasXsec
	mode
	noWarnings
	options
	output
	parameters
	## Main Program More...
	parser
	pieces
	pileupGroup
	pileupHist
	Let's start the fun. More...
	pileupHistName
	recorded
	recordedData
	run
	runLumiDict
	runNumber
	pprint (recordedData) More...
	sepRE
	session
	svc
	type
	verbose
	Save what we need in the parameters object. More...
	xingInstLumiArray
	xingMinLum

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 16 of file estimatePileup.py.

References createfilelist.int, and edm.print().

                         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 (six.iteritems(deadTable)):
        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
    

Variable Documentation

estimatePileup.action

Definition at line 110 of file estimatePileup.py.

estimatePileup.args

Definition at line 151 of file estimatePileup.py.

estimatePileup.basename

Definition at line 244 of file estimatePileup.py.

Referenced by utils_v2.auth_download_file(), fetchall_from_DQM_v2.auth_wget(), compare_using_files_v2.RootFileComparison.compare(), ValidationMatrix_v2.ReleaseComparison.compare(), ExternalLHEAsciiDumper.endRun(), directories2html.get_comparisons(), relmon_rootfiles_spy.get_page(), directories2html.get_title_section(), WatcherStreamFileReader.getInputFile(), utils_v2.init_database(), EcalSimRawData.iTtPhi2cIndex(), directories2html.make_summary_table(), and utils.wget().

estimatePileup.beammode

Definition at line 171 of file estimatePileup.py.

estimatePileup.beamModeChoices

Definition at line 104 of file estimatePileup.py.

estimatePileup.choices

Definition at line 123 of file estimatePileup.py.

estimatePileup.csvDict

Definition at line 196 of file estimatePileup.py.

estimatePileup.dbGroup

Definition at line 107 of file estimatePileup.py.

estimatePileup.deadTable

Definition at line 260 of file estimatePileup.py.

estimatePileup.debug

Definition at line 219 of file estimatePileup.py.

estimatePileup.default

Definition at line 123 of file estimatePileup.py.

estimatePileup.delivered

Definition at line 206 of file estimatePileup.py.

estimatePileup.description

Definition at line 106 of file estimatePileup.py.

estimatePileup.dest

Definition at line 110 of file estimatePileup.py.

estimatePileup.events

Definition at line 195 of file estimatePileup.py.

estimatePileup.extension

Definition at line 244 of file estimatePileup.py.

estimatePileup.f

Definition at line 249 of file estimatePileup.py.

estimatePileup.fileparsingResult

Definition at line 247 of file estimatePileup.py.

estimatePileup.help

Definition at line 111 of file estimatePileup.py.

estimatePileup.hist

Definition at line 217 of file estimatePileup.py.

estimatePileup.histFile

Definition at line 229 of file estimatePileup.py.

estimatePileup.histList

Definition at line 188 of file estimatePileup.py.

estimatePileup.inputfilecontent

Definition at line 250 of file estimatePileup.py.

estimatePileup.inputGroup

Definition at line 108 of file estimatePileup.py.

estimatePileup.inputRange

Get input source.

Definition at line 242 of file estimatePileup.py.

Referenced by SiStripNoiseNormalizedWithApvGainBuilder.analyze(), SiStripApvGainBuilderFromTag.analyze(), SiStripGain.fillNewGain(), cms::ReconstructerFP420.produce(), cms::DigitizerFP420.produce(), cms::TrackerizerFP420.produce(), cms::ClusterizerFP420.produce(), CTPPSPixelDigiCollection.put(), PixelDigiCollection.put(), DigiCollectionFP420.put(), TrackCollectionFP420.put(), ClusterCollectionFP420.put(), RecoCollectionFP420.put(), RecoCollectionFP420.putclear(), TrackCollectionFP420.putclear(), ClusterCollectionFP420.putclear(), DigiCollectionFP420.putclear(), FP420TrackMain.run(), and FP420ClusterMain.run().

estimatePileup.lumi

Definition at line 205 of file estimatePileup.py.

estimatePileup.lumiversion

Definition at line 169 of file estimatePileup.py.

Referenced by lumi::LumiDummy2DB.retrieveData(), LumiSummary.setLumiVersion(), lumi::Lumi2DB.writeAllLumiData(), lumi::Lumi2DB.writeBeamIntensityOnly(), and LumiSummary.~LumiSummary().

estimatePileup.lumiXing

Definition at line 168 of file estimatePileup.py.

estimatePileup.maxPileupBin

Definition at line 175 of file estimatePileup.py.

estimatePileup.minBiasXsec

Definition at line 173 of file estimatePileup.py.

estimatePileup.mode

Definition at line 220 of file estimatePileup.py.

estimatePileup.noWarnings

Definition at line 167 of file estimatePileup.py.

estimatePileup.options

Definition at line 151 of file estimatePileup.py.

estimatePileup.output

Definition at line 159 of file estimatePileup.py.

estimatePileup.parameters

## Main Program

Definition at line 103 of file estimatePileup.py.

estimatePileup.parser

Definition at line 105 of file estimatePileup.py.

estimatePileup.pieces

Definition at line 198 of file estimatePileup.py.

estimatePileup.pileupGroup

Definition at line 109 of file estimatePileup.py.

estimatePileup.pileupHist

Let's start the fun.

Definition at line 185 of file estimatePileup.py.

estimatePileup.pileupHistName

Definition at line 174 of file estimatePileup.py.

estimatePileup.recorded

Definition at line 206 of file estimatePileup.py.

estimatePileup.recordedData

Definition at line 256 of file estimatePileup.py.

estimatePileup.run

Definition at line 205 of file estimatePileup.py.

estimatePileup.runLumiDict

Definition at line 194 of file estimatePileup.py.

estimatePileup.runNumber

pprint (recordedData)

Definition at line 218 of file estimatePileup.py.

estimatePileup.sepRE

Definition at line 193 of file estimatePileup.py.

estimatePileup.session

Definition at line 177 of file estimatePileup.py.

estimatePileup.svc

Definition at line 177 of file estimatePileup.py.

estimatePileup.type

Definition at line 125 of file estimatePileup.py.

estimatePileup.verbose

Save what we need in the parameters object.

Definition at line 166 of file estimatePileup.py.

estimatePileup.xingInstLumiArray

Definition at line 207 of file estimatePileup.py.

estimatePileup.xingMinLum

Definition at line 172 of file estimatePileup.py.