pileupCalc Namespace Reference

Functions
def	fillPileupHistogram
def	parseInputFile
Variables
list	CalculationModeChoices = ['truth', 'observed']
string	default = 'PileupCalc.root'
string	help = 'output root file'
tuple	histFile = ROOT.TFile.Open(output, 'recreate')
tuple	inpf = open(options.inputfile, 'r')
tuple	inputfilecontent = inpf.read()
tuple	inputPileupRange = parseInputFile(options.inputLumiJSON)
tuple	inputRange = selectionParser.selectionParser(inputfilecontent)
list	LSPUlist = inputPileupRange[run]
list	lumiInfo = LSPUlist[LSnumber]
	nbins = options.numPileupBins
list	output = args[0]
tuple	parser
	######################## ## ## ################## ## ## ## ## Main Program ## ## ## ## ################## ## ## ######################## ##
tuple	pileupHist
	type = float,
	upper = options.maxPileupBin
string	VERSION = '1.00'

---

Function Documentation

def pileupCalc::fillPileupHistogram	(		lumiInfo,
			calcOption,
			hist,
			minbXsec,
			Nbins
	)

lumiinfo:[intlumi per LS, mean interactions ]

intlumi is the deadtime corrected average integraged lumi per lumisection

Definition at line 22 of file pileupCalc.py.

                                                                     :
    '''
    lumiinfo:[intlumi per LS, mean interactions ]

    intlumi is the deadtime corrected average integraged lumi per lumisection
    '''

    LSintLumi = lumiInfo[0]
    RMSInt = lumiInfo[1]*minbXsec
    AveNumInt = lumiInfo[2]*minbXsec

    coeff = 0
    if RMSInt > 0:
        coeff = 1.0/RMSInt/sqrt(6.283185)
    expon = 2.0*RMSInt*RMSInt

    ##Nbins = hist.GetXaxis().GetNbins()

    ProbFromRMS = []

    # First, re-constitute lumi distribution for this LS from RMS:
    if RMSInt > 0:
        for obs in range (Nbins):
            val = hist.GetBinCenter(obs+1)
            prob = coeff*exp(-1.0*(val-AveNumInt)*(val-AveNumInt)/expon)
            ProbFromRMS.append(prob)
    else:
        obs = hist.FindBin(AveNumInt)
        for bin in range (Nbins):
            ProbFromRMS.append(0.0)
        ProbFromRMS[obs] = 1.0
        if AveNumInt < 1.0E-5:
           ProbFromRMS[obs] = 0.  # just ignore zero values
        
    if calcOption == 'true':  # Just put distribution into histogram
        if RMSInt > 0:
            totalProb = 0
            for obs in range (Nbins):
                prob = ProbFromRMS[obs]
                val = hist.GetBinCenter(obs+1)
                # print obs, val, RMSInt,coeff,expon,prob
                totalProb += prob
                hist.Fill (val, prob * LSintLumi)
            if totalProb < 1:
                hist.Fill (val, (1 - totalProb) * LSintLumi)
        else:
            hist.Fill(AveNumInt,LSintLumi)
    else: # have to convolute with a poisson distribution to get observed Nint
        totalProb = 0
        Peak = 0
        BinWidth = hist.GetBinWidth(1)
        for obs in range (Nbins):
            Peak = hist.GetBinCenter(obs+1)
            RMSWeight = ProbFromRMS[obs]
            for bin in range (Nbins):
                val = hist.GetBinCenter(bin+1)-0.5*BinWidth
                prob = ROOT.TMath.Poisson (val, Peak)
                totalProb += prob
                hist.Fill (val, prob * LSintLumi * RMSWeight)
        if totalProb < 1:
            hist.Fill (Peak, (1 - totalProb) * LSintLumi)

    return hist

def pileupCalc::parseInputFile

(

inputfilename

)

output ({run:[ls:[inlumi, meanint]]})

Definition at line 10 of file pileupCalc.py.

                                 :
    '''
    output ({run:[ls:[inlumi, meanint]]})
    '''
    selectf=open(inputfilename,'r')
    inputfilecontent=selectf.read()
    p=pileupParser.pileupParser(inputfilecontent)                            
    
#    p=inputFilesetParser.inputFilesetParser(inputfilename)
    runlsbyfile=p.runsandls()
    return runlsbyfile

---

Variable Documentation

list pileupCalc::CalculationModeChoices = ['truth', 'observed']

Definition at line 102 of file pileupCalc.py.

string pileupCalc::default = 'PileupCalc.root'

Definition at line 113 of file pileupCalc.py.

string pileupCalc::help = 'output root file'

Definition at line 114 of file pileupCalc.py.

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

Definition at line 208 of file pileupCalc.py.

tuple pileupCalc::inpf = open(options.inputfile, 'r')

Definition at line 170 of file pileupCalc.py.

tuple pileupCalc::inputfilecontent = inpf.read()

Definition at line 171 of file pileupCalc.py.

tuple pileupCalc::inputPileupRange = parseInputFile(options.inputLumiJSON)

Definition at line 177 of file pileupCalc.py.

tuple pileupCalc::inputRange = selectionParser.selectionParser(inputfilecontent)

Definition at line 172 of file pileupCalc.py.

list pileupCalc::LSPUlist = inputPileupRange[run]

Definition at line 187 of file pileupCalc.py.

list pileupCalc::lumiInfo = LSPUlist[LSnumber]

Definition at line 192 of file pileupCalc.py.

Referenced by edm::BranchTypeToInfoTreeName().

pileupCalc::nbins = options.numPileupBins

Definition at line 167 of file pileupCalc.py.

Referenced by ApvTimingAlgorithm::analyse(), FastFedCablingAlgorithm::analyse(), FedTimingAlgorithm::analyse(), OptoScanAlgorithm::analyse(), DaqScopeModeAlgorithm::analyse(), XmasToDQMSource::analyze(), TopDiLeptonDQM::beginJob(), DQMMessageLogger::beginJob(), L1TDTTF::beginJob(), BTagHLTOfflineSource::beginRun(), JetMETHLTOfflineSource::beginRun(), FedTimingTask::book(), OptoScanTask::book(), PedestalsTask::book(), VpspScanTask::book(), FastFedCablingTask::book(), ApvTimingTask::book(), FedCablingTask::book(), L1TGMT::book_(), DTTriggerEfficiencyTest::bookChambHistos(), DTLocalTriggerSynchTask::bookHistos(), TrackerOfflineValidation::bookHists(), TrackerOfflineValidation::bookSummaryHists(), L1TDTTFClient::buildPhiEtaPlotO(), L1TDTTFClient::buildPhiEtaPlotOFC(), CentralityBins::CentralityBins(), L1TauAnalyzer::convertToIntegratedEff(), SamplingAlgorithm::correctProfile(), reco::DiscretizedEnergyFlow::DiscretizedEnergyFlow(), PedsFullNoiseTask::fill(), FineDelayTask::fill(), PedestalsTask::fill(), SiStripSummaryCreator::fillHistos(), ecaldqm::getAverageFromTProfile(), SiPixelActionExecutor::getData(), TrackerOfflineValidation::getMedian(), TrackerOfflineValidationSummary::getMedian(), HiEvtPlaneFlatCalib::HiEvtPlaneFlatCalib(), MuonResidualsFitter::histogramChi2GaussianFit(), MSLayersKeeperX0Averaged::init(), ExpressionHisto< T >::initialize(), PhiScaleHelper::makeBinnedScale(), DTLocalTriggerEfficiencyTest::makeEfficiencyME(), MuIsoValidation::MakeLogBinsForProfile(), AlignmentMonitorMuonSystemMap1D::MuonSystemMapPlot1D::MuonSystemMapPlot1D(), DTSegmentAnalysisTest::performClientDiagnostic(), HLTTauDQMSummaryPlotter::plotIntegratedEffHisto(), DistortedMuonProducer::produce(), ISRWeightProducer::produce(), DistortedPFCandProducer::produce(), SamplingAlgorithm::pruneProfile(), CSCCrosstalkGenerator::ratio(), MuScleFitBase::readProbabilityDistributionsFromFile(), RecoTauPlotDiscriminator::RecoTauPlotDiscriminator(), CompareToMedian::runTest(), NoisyChannel::runTest(), Comp2RefEqualH::runTest(), HcalDigiMonitor::setup(), ecaldqm::shiftAxis(), MonitorElement::ShiftFillLast(), smartGausProfile(), smartProfile(), BTagHLTOfflineSource::TriggerPosition(), and JetMETHLTOfflineSource::TriggerPosition().

list pileupCalc::output = args[0]

Definition at line 149 of file pileupCalc.py.

tuple pileupCalc::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")

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

Definition at line 98 of file pileupCalc.py.

tuple pileupCalc::pileupHist

Initial value:

ROOT.TH1D(options.pileupHistName,options.pileupHistName,
                      options.numPileupBins,
                      0., options.maxPileupBin)

Definition at line 163 of file pileupCalc.py.

pileupCalc::type = float,

Definition at line 122 of file pileupCalc.py.

pileupCalc::upper = options.maxPileupBin

Definition at line 168 of file pileupCalc.py.

Referenced by LHEAnalyzer::analyze(), MomentumDependentPedeLabeler::buildMomentumDependencyMap(), ConstrainedTreeBuilder::covarianceMatrix(), ConstrainedTreeBuilderT::covarianceMatrix(), WebInterface::DrawGif(), MagVolumeOutsideValidity::MagVolumeOutsideValidity(), and MFGrid3D::valueInTesla().

string pileupCalc::VERSION = '1.00'

Definition at line 2 of file pileupCalc.py.