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.
string	mode = 'csv'
list	output = args[0]
tuple	parameters = LumiQueryAPI.ParametersObject()
	######################## ## ## ################## ## ## ## ## Main Program ## ## ## ## ################## ## ## ######################## ##
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.
tuple	recordedData = LumiQueryAPI.recordedLumiForRange(session, parameters, inputRange)
dictionary	runLumiDict = {}
list	runNumber = runDTarray[0]
	pprint (recordedData)
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 12 of file estimatePileup.py.

                                          :
    '''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

string estimatePileup::action = 'store_true'

Definition at line 140 of file estimatePileup.py.

list estimatePileup::beamModeChoices = ["","stable", "quiet", "either"]

Definition at line 102 of file estimatePileup.py.

dictionary estimatePileup::csvDict = {}

Definition at line 194 of file estimatePileup.py.

tuple estimatePileup::dbGroup = optparse.OptionGroup(parser, "Database Options")

Definition at line 105 of file estimatePileup.py.

list estimatePileup::deadTable = runDTarray[2]

Definition at line 259 of file estimatePileup.py.

estimatePileup::debug = options.debugLumi,

Definition at line 224 of file estimatePileup.py.

tuple estimatePileup::events = open(options.csvInput, 'r')

Definition at line 193 of file estimatePileup.py.

tuple estimatePileup::f = open(options.inputfile, 'r')

Definition at line 248 of file estimatePileup.py.

tuple estimatePileup::fileparsingResult = csvSelectionParser.csvSelectionParser(options.inputfile)

Definition at line 246 of file estimatePileup.py.

string estimatePileup::help = 'connect string to lumiDB (optional, default to frontier://LumiCalc/CMS_LUMI_PROD)'

Definition at line 109 of file estimatePileup.py.

tuple estimatePileup::hist

Initial value:

fillPileupHistogram(lumiDict, parameters,
                                            runNumber = runNumber,
                                            debug = options.debugLumi,
                                            mode='csv')

Definition at line 215 of file estimatePileup.py.

Referenced by PrimaryVertexAnalyzer4PU::add(), PrimaryVertexAnalyzer4PU::beginJob(), ImpactParameterCalibration::beginJob(), CalibAnalyzer< Ref, Rec, Alg >::beginJob(), APVCyclePhaseMonitor::beginRun(), PhysicalPartsTree::beginRun(), Benchmark::book1D(), Benchmark::book2D(), TouchableToHistory::buildAll(), fwlite::Scanner< Collection >::draw(), examples::TrackAnalysisAlgorithm::draw(), fwlite::Scanner< Collection >::draw2D(), fwlite::Scanner< Collection >::drawProf(), L1TEventInfoClient::dumpContentMonitorElements(), PrimaryVertexAnalyzer::endJob(), EcalPedHists::endJob(), EcalCosmicsHists::endJob(), PrimaryVertexAnalyzer4PU::endJob(), EcalURecHitHists::endJob(), DQMHistNormalizer::endRun(), ExpressionHisto< T >::fill(), DigiCollectionProfiler< T >::fill(), edm::Principal::fillPrincipal(), EgHLTOfflineSummaryClient::getEgHLTSumHist_(), TouchableToHistory::getNavStory(), fwlite::Scanner< Collection >::getSameH1(), fwlite::Scanner< Collection >::getSameH2(), fwlite::Scanner< Collection >::getSameProf(), FWEveViewManager::highlightAdded(), HistoAnalyzer< C >::HistoAnalyzer(), MuonResidualsFitter::histogramChi2GaussianFit(), EcalDisplaysByEvent::init2DEcalHist(), EcalDisplaysByEvent::init3DEcalHist(), ImpactParameterCalibration::initFromFirstES(), EcalURecHitHists::initHists(), EcalCosmicsHists::initHists(), ExpressionHisto< T >::initialize(), main(), LA_Filler_Fitter::make_and_fit_symmchi2(), SiStripSpyDisplayModule::MakeDigiHist_(), ValidateGeometry::makeHistogram(), SiStripSpyDisplayModule::MakeProcessedRawDigiHist_(), SiStripSpyDisplayModule::MakeRawDigiHist_(), DDExpandedView::mergedSpecificsV(), DQMGenericClient::normalizeToEntries(), DQMRivetClient::normalizeToIntegral(), DQMRivetClient::normalizeToLumi(), CSCAlignmentCorrections::plot(), MuonResidualsFitter::plotsimple(), MuonResidualsFitter::plotweighted(), EcalPedHists::readEBdigis(), EcalPedHists::readEEdigis(), DQMRivetClient::scaleByFactor(), Book::size(), DDExpandedView::specificsV(), tfwliteselectortest::ThingsTSelector::terminate(), __class__< T >::terminate(), tfwliteselectortest::ThingsTSelector2::terminate(), sistrip::EnsembleCalibrationLA::write_ensembles_plots(), sistrip::MeasureLA::write_report_plots(), and sistrip::EnsembleCalibrationLA::write_samples_plots().

tuple estimatePileup::histFile = ROOT.TFile.Open(output, 'recreate')

Definition at line 227 of file estimatePileup.py.

list estimatePileup::histList = []

Definition at line 186 of file estimatePileup.py.

tuple estimatePileup::inputfilecontent = f.read()

Definition at line 249 of file estimatePileup.py.

tuple estimatePileup::inputGroup = optparse.OptionGroup(parser, "Input Options")

Definition at line 106 of file estimatePileup.py.

tuple estimatePileup::inputRange = options.runnumber

Get input source.

Definition at line 241 of file estimatePileup.py.

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

string estimatePileup::mode = 'csv'

Definition at line 225 of file estimatePileup.py.

list estimatePileup::output = args[0]

Definition at line 157 of file estimatePileup.py.

tuple estimatePileup::parameters = LumiQueryAPI.ParametersObject()

######################## ## ## ################## ## ## ## ## Main Program ## ## ## ## ################## ## ## ######################## ##

Definition at line 101 of file estimatePileup.py.

tuple estimatePileup::parser

Initial value:

optparse.OptionParser("Usage: %prog [--options] output.root",
                                    description = "Script to estimate pileup distribution using xing instantaneous luminosity information and minimum bias cross section.  Output is TH1D stored in root file")

Definition at line 103 of file estimatePileup.py.

tuple estimatePileup::pieces = sepRE.split(line.strip())

Definition at line 196 of file estimatePileup.py.

tuple estimatePileup::pileupGroup = optparse.OptionGroup(parser, "Pileup Options")

Definition at line 107 of file estimatePileup.py.

tuple estimatePileup::pileupHist

Initial value:

ROOT.TH1D(parameters.pileupHistName, parameters.pileupHistName,
                      1001,
                      0.0, 25.)

Save what we need in the parameters object.

Let's start the fun

Definition at line 183 of file estimatePileup.py.

tuple estimatePileup::recordedData = LumiQueryAPI.recordedLumiForRange(session, parameters, inputRange)

Definition at line 255 of file estimatePileup.py.

dictionary estimatePileup::runLumiDict = {}

Definition at line 192 of file estimatePileup.py.

list estimatePileup::runNumber = runDTarray[0]

pprint (recordedData)

Definition at line 258 of file estimatePileup.py.

tuple estimatePileup::sepRE = re.compile(r'[\s,;:]+')

Definition at line 191 of file estimatePileup.py.

list estimatePileup::xingInstLumiArray

Initial value:

[( int(orbit), float(lum) ) \
                                     for orbit, lum in zip( pieces[9::2],
                                                            pieces[10::2] ) ]

Definition at line 205 of file estimatePileup.py.