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

References createfilelist.int.

                         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
    

Variable Documentation

estimatePileup.action

Definition at line 108 of file estimatePileup.py.

estimatePileup.args

Definition at line 149 of file estimatePileup.py.

estimatePileup.basename

Definition at line 242 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 169 of file estimatePileup.py.

estimatePileup.beamModeChoices

Definition at line 102 of file estimatePileup.py.

estimatePileup.choices

Definition at line 121 of file estimatePileup.py.

estimatePileup.csvDict

Definition at line 194 of file estimatePileup.py.

estimatePileup.dbGroup

Definition at line 105 of file estimatePileup.py.

estimatePileup.deadTable

Definition at line 258 of file estimatePileup.py.

estimatePileup.debug

Definition at line 217 of file estimatePileup.py.

estimatePileup.default

Definition at line 121 of file estimatePileup.py.

estimatePileup.delivered

Definition at line 204 of file estimatePileup.py.

estimatePileup.description

Definition at line 104 of file estimatePileup.py.

estimatePileup.dest

Definition at line 108 of file estimatePileup.py.

estimatePileup.events

Definition at line 193 of file estimatePileup.py.

estimatePileup.extension

Definition at line 242 of file estimatePileup.py.

estimatePileup.f

Definition at line 247 of file estimatePileup.py.

estimatePileup.fileparsingResult

Definition at line 245 of file estimatePileup.py.

estimatePileup.help

Definition at line 109 of file estimatePileup.py.

estimatePileup.hist

Definition at line 215 of file estimatePileup.py.

estimatePileup.histFile

Definition at line 227 of file estimatePileup.py.

estimatePileup.histList

Definition at line 186 of file estimatePileup.py.

estimatePileup.inputfilecontent

Definition at line 248 of file estimatePileup.py.

estimatePileup.inputGroup

Definition at line 106 of file estimatePileup.py.

estimatePileup.inputRange

Get input source.

Definition at line 240 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 203 of file estimatePileup.py.

estimatePileup.lumiversion

Definition at line 167 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 166 of file estimatePileup.py.

estimatePileup.maxPileupBin

Definition at line 173 of file estimatePileup.py.

estimatePileup.minBiasXsec

Definition at line 171 of file estimatePileup.py.

estimatePileup.mode

Definition at line 218 of file estimatePileup.py.

estimatePileup.noWarnings

Definition at line 165 of file estimatePileup.py.

estimatePileup.options

Definition at line 149 of file estimatePileup.py.

estimatePileup.output

Definition at line 157 of file estimatePileup.py.

estimatePileup.parameters

## Main Program

Definition at line 101 of file estimatePileup.py.

estimatePileup.parser

Definition at line 103 of file estimatePileup.py.

estimatePileup.pieces

Definition at line 196 of file estimatePileup.py.

estimatePileup.pileupGroup

Definition at line 107 of file estimatePileup.py.

estimatePileup.pileupHist

Let's start the fun.

Definition at line 183 of file estimatePileup.py.

estimatePileup.pileupHistName

Definition at line 172 of file estimatePileup.py.

estimatePileup.recorded

Definition at line 204 of file estimatePileup.py.

estimatePileup.recordedData

Definition at line 254 of file estimatePileup.py.

estimatePileup.run

Definition at line 203 of file estimatePileup.py.

estimatePileup.runLumiDict

Definition at line 192 of file estimatePileup.py.

estimatePileup.runNumber

pprint (recordedData)

Definition at line 216 of file estimatePileup.py.

estimatePileup.sepRE

Definition at line 191 of file estimatePileup.py.

estimatePileup.session

Definition at line 175 of file estimatePileup.py.

estimatePileup.svc

Definition at line 175 of file estimatePileup.py.

estimatePileup.type

Definition at line 123 of file estimatePileup.py.

estimatePileup.verbose

Save what we need in the parameters object.

Definition at line 164 of file estimatePileup.py.

estimatePileup.xingInstLumiArray

Definition at line 205 of file estimatePileup.py.

estimatePileup.xingMinLum

Definition at line 170 of file estimatePileup.py.