Functions
def	fillPileupHistogram (lumiInfo, calcOption, hist, minbXsec, Nbins)

def	MyErf (input)

def	parseInputFile (inputfilename)

Variables
	action

	args

	CalculationModeChoices

	default

	description

	dest

	help

	histFile

	inpf

	inputfilecontent

	inputPileupRange

	inputRange

	LSPUlist

	lumiInfo

	nbins

	options

	output

	parser
	## Main Program More...

	pileupHist

	type

	upper

	VERSION

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 54 of file pileupCalc.py.

References MyErf(), and mathSSE.sqrt().

 def fillPileupHistogram (lumiInfo, calcOption, hist, minbXsec, Nbins):
     '''
     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
 
     Sqrt2 = sqrt(2)
 
     ##Nbins = hist.GetXaxis().GetNbins()
 
     ProbFromRMS = []
     BinWidth = hist.GetBinWidth(1)
 
     # First, re-constitute lumi distribution for this LS from RMS:
     if RMSInt > 0:
 
         AreaLnew = -10.
         AreaL = 0
 
         for obs in range (Nbins):
             #Old Gaussian normalization; inaccurate for small rms and large bins
             #val = hist.GetBinCenter(obs+1)
             #prob = coeff*exp(-1.0*(val-AveNumInt)*(val-AveNumInt)/expon)
             #ProbFromRMS.append(prob)
             
             left = hist.GetBinLowEdge(obs+1)
             right = left+BinWidth
 
             argR = (AveNumInt-right)/Sqrt2/RMSInt
             AreaR = MyErf(argR)
 
             if AreaLnew<-5.:
                 argL = (AveNumInt-left)/Sqrt2/RMSInt
                 AreaL = MyErf(argL)
             else:
                 AreaL = AreaLnew
                 AreaLnew = AreaR  # save R bin value for L next time
 
             NewProb = (AreaL-AreaR)*0.5
 
             ProbFromRMS.append(NewProb)
 
             #print left, right, argL, argR, AreaL, AreaR, NewProb
 
     else:
         obs = hist.FindBin(AveNumInt)
         for bin in range (Nbins):
             ProbFromRMS.append(0.0)
         if obs<Nbins+1:            
             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 1.0-totalProb > 0.01:
                 print "Significant probability density outside of your histogram"
                 print "Consider using a higher value of --maxPileupBin"
                 print "Mean %f, RMS %f, Integrated probability %f" % (AveNumInt,RMSInt,totalProb)
             #    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 1.0-totalProb > 0.01:
             print "Significant probability density outside of your histogram"
             print "Consider using a higher value of --maxPileupBin"
 
 
     return hist
 
 
 

def pileupCalc.MyErf ( input )

Definition at line 22 of file pileupCalc.py.

References funct.abs(), and JetChargeProducer_cfi.exp.

Referenced by fillPileupHistogram().

 def MyErf(input):
 
     # Abramowitz and Stegun approximations for Erf (equations 7.1.25-28)
     X = abs(input)
 
     p = 0.47047
     b1 = 0.3480242
     b2 = -0.0958798
     b3 = 0.7478556
 
     T = 1.0/(1.0+p*X)
     cErf = 1.0 - (b1*T + b2*T*T + b3*T*T*T)*exp(-1.0*X*X)
     if input<0:
         cErf = -1.0*cErf
 
     # Alternate Erf approximation:
     
     #A1 = 0.278393
     #A2 = 0.230389
     #A3 = 0.000972
     #A4 = 0.078108
 
     #term = 1.0+ A1*X+ A2*X*X+ A3*X*X*X+ A4*X*X*X*X
     #denom = term*term*term*term
 
     #dErf = 1.0 - 1.0/denom
     #if input<0:
     #    dErf = -1.0*dErf
         
     return cErf
 
 

def pileupCalc.parseInputFile ( inputfilename )

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

Definition at line 10 of file pileupCalc.py.

 def parseInputFile(inputfilename):
     '''
     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

pileupCalc.action

Definition at line 182 of file pileupCalc.py.

pileupCalc.args

Definition at line 210 of file pileupCalc.py.

pileupCalc.CalculationModeChoices

Definition at line 172 of file pileupCalc.py.

pileupCalc.default

Definition at line 183 of file pileupCalc.py.

pileupCalc.description

Definition at line 169 of file pileupCalc.py.

pileupCalc.dest

Definition at line 182 of file pileupCalc.py.

pileupCalc.help

Definition at line 184 of file pileupCalc.py.

pileupCalc.histFile

Definition at line 278 of file pileupCalc.py.

pileupCalc.inpf

Definition at line 240 of file pileupCalc.py.

pileupCalc.inputfilecontent

Definition at line 241 of file pileupCalc.py.

pileupCalc.inputPileupRange

Definition at line 247 of file pileupCalc.py.

pileupCalc.inputRange

Definition at line 242 of file pileupCalc.py.

pileupCalc.LSPUlist

Definition at line 257 of file pileupCalc.py.

pileupCalc.lumiInfo

Definition at line 262 of file pileupCalc.py.

Referenced by edm.BranchTypeToInfoTreeName(), and edm.operator<<().

pileupCalc.nbins

Definition at line 237 of file pileupCalc.py.

Referenced by npstat::HistoND< Numeric, Axis >.accumulateBinsLoop(), FedTimingAlgorithm.analyse(), ApvTimingAlgorithm.analyse(), FastFedCablingAlgorithm.analyse(), OptoScanAlgorithm.analyse(), DaqScopeModeAlgorithm.analyse(), BigEventsDebugger< T >.analyze(), ecaldqm::binning::AxisSpecs.AxisSpecs(), BigEventsDebugger< T >.BigEventsDebugger(), FastFedCablingTask.book(), ApvTimingTask.book(), FedTimingTask.book(), VpspScanTask.book(), OptoScanTask.book(), PedestalsTask.book(), FedCablingTask.book(), ecaldqm::MESetEcal.book(), DTTriggerEfficiencyTask.bookChamberHistos(), DTLocalTriggerEfficiencyTest.bookChambHistos(), DTLocalTriggerSynchTest.bookChambHistos(), DTTriggerEfficiencyTest.bookChambHistos(), DQMMessageLogger.bookHistograms(), CTPPSPixelDQMSource.bookHistograms(), TopDiLeptonDQM.bookHistograms(), L1TDTTF.bookHistograms(), PSMonitor.bookHistograms(), L1TGMT.bookHistograms(), JetMETHLTOfflineSource.bookHistograms(), HLTObjectsMonitor.bookHistograms(), DTLocalTriggerSynchTask.bookHistos(), DTLocalTriggerTask.bookHistos(), dqm::TrackAnalyzer.bookHistosForEfficiencyFromHitPatter(), TrackerOfflineValidation.bookHists(), TriggerDQMBase.bookME(), BPHMonitor.bookME(), PhotonMonitor.bookME(), MuonMonitor.bookME(), RazorMonitor.bookME(), JetMonitor.bookME(), METMonitor.bookME(), HTMonitor.bookME(), NoBPTXMonitor.bookNoBPTX(), TrackerOfflineValidation.bookSummaryHists(), BTagEntry.BTagEntry(), L1TDTTFClient.buildPhiEtaPlotO(), L1TDTTFClient.buildPhiEtaPlotOFC(), PixelLumiDQM.calculateBunchMask(), CentralityBins.CentralityBins(), EcalSelectiveReadoutValidation.cIndex2iTtPhi(), SamplingAlgorithm.correctProfile(), reco::DiscretizedEnergyFlow.DiscretizedEnergyFlow(), FitterFuncs::PulseShapeFunctor.EvalPulse(), HistogramManager.executeExtend(), hcaldqm::ContainerSingle1D.extendAxisRange(), hcaldqm::ContainerSingle2D.extendAxisRange(), hcaldqm::Container1D.extendAxisRange(), PedestalsTask.fill(), FineDelayTask.fill(), PedsFullNoiseTask.fill(), SiStripSummaryCreator.fillHistos(), FastTimerServiceClient.fillPlotsVsLumi(), ThroughputServiceClient.fillSummaryPlots(), ecaldqm::binning.getBinningMEM_(), PVValHelper.getMAD(), TrackerOfflineValidationSummary.getMedian(), TrackerOfflineValidation.getMedian(), heppy::BTagSF.getSFb(), heppy::BTagSF.getSFc(), MuonResidualsFitter.histogramChi2GaussianFit(), MSLayersKeeperX0Averaged.init(), Benchmark.isInRange(), PhiScaleHelper.makeBinnedScale(), DTLocalTriggerEfficiencyTest.makeEfficiencyME(), MuIsoValidation.MakeLogBinsForProfile(), AlignmentMonitorMuonSystemMap1D::MuonSystemMapPlot1D.MuonSystemMapPlot1D(), ecaldqm::binning::AxisSpecs.operator=(), DTSegmentAnalysisTest.performClientDiagnostic(), DistortedMuonProducer.produce(), DistortedPFCandProducer.produce(), ISRWeightProducer.produce(), SamplingAlgorithm.pruneProfile(), CSCCrosstalkGenerator.ratio(), RecoTauPlotDiscriminator.RecoTauPlotDiscriminator(), MuonResidualsFitter.residuals_end(), Comp2RefEqualH.runTest(), NoisyChannel.runTest(), CompareToMedian.runTest(), MonitorElement.ShiftFillLast(), smartGausProfile(), smartProfile(), th1ToFormulaLin(), and JetMETHLTOfflineSource.TriggerPosition().

pileupCalc.options

Definition at line 210 of file pileupCalc.py.

pileupCalc.output

Definition at line 219 of file pileupCalc.py.

pileupCalc.parser

## Main Program

Definition at line 168 of file pileupCalc.py.

pileupCalc.pileupHist

Definition at line 233 of file pileupCalc.py.

pileupCalc.type

Definition at line 192 of file pileupCalc.py.

pileupCalc.upper

Definition at line 238 of file pileupCalc.py.

Referenced by TrackValidation_cff._algoToSelector(), trackingPlots._summaryBinRename(), data_sources.json_list.as_table(), MomentumDependentPedeLabeler.buildMomentumDependencyMap(), ConstrainedTreeBuilder.covarianceMatrix(), ConstrainedTreeBuilderT.covarianceMatrix(), HcalDbASCIIIO.createObject< HcalElectronicsMap >(), InterpolatedPulse< 1500U >.derivative(), Vispa.Plugins.ConfigEditor.ConfigEditorTabController.ConfigEditorTabController.dumpPython(), cmsPerfSuiteHarvest.exportIgProfReport(), TableParser.extractPages(), MuonShowerInformationFiller.fillHitsByStation(), GEMTrackMatch.getClosestChambers(), cond::persistency::RUN_INFO::Table.getInclusiveRunRange(), cond::persistency::RUN_INFO::Table.getInclusiveTimeRange(), FileNamesHelper.getJobID_fromIgProfLogName(), MagVolumeOutsideValidity.MagVolumeOutsideValidity(), HipPyOptionParser.parseBoolString(), InterpolatedPulse< 1500U >.secondDerivative(), jetTools.AddJetCollection.toolCode(), jetTools.UpdateJetCollection.toolCode(), MFGrid3D.valueInTesla(), cond::persistency::IRunInfoTable.~IRunInfoTable(), and cond::persistency::RUN_INFO::Table.~Table().

pileupCalc.VERSION

Definition at line 2 of file pileupCalc.py.

Functions

Variables

## Main Program

Function Documentation

Variable Documentation

## Main Program