Classes
class	EquidistantBinning

Functions
def	fillPileupHistogram (lumiInfo, calcOption, binning, hContents, minbXsec, run, ls)

def	MyErf (xInput)

def	parseInputFile (inputfilename)

def	poisson (x, par)

Variables
	action

	binning

	choices

	default

	description

	dest

	float

	hContents

	help

	histFile

	inpf

	inputfilecontent

	inputPileupRange

	inputRange

	int

	LSPUlist

	lumiInfo

	options

	output

	parser
	## Main Program More...

	pileupHist

	req_group

	required

	Sqrt2

	type

	VERSION

Function Documentation

◆ fillPileupHistogram()

def pileupCalc.fillPileupHistogram	(	lumiInfo,
		calcOption,
		binning,
		hContents,
		minbXsec,
		run,
		ls
	)

lumiinfo:[intlumi per LS, mean interactions ]

intlumi is the deadtime corrected average integrated lumi per lumisection

Definition at line 69 of file pileupCalc.py.

References MyErf(), poisson(), and print().

 def fillPileupHistogram(lumiInfo, calcOption, binning, hContents, minbXsec, run, ls):
     '''
     lumiinfo:[intlumi per LS, mean interactions ]
 
     intlumi is the deadtime corrected average integrated lumi per lumisection
     '''
 
     LSintLumi = lumiInfo[0]
     RMSInt = lumiInfo[1]*minbXsec
     AveNumInt = lumiInfo[2]*minbXsec
 
     # First, re-constitute lumi distribution for this LS from RMS:
     if RMSInt > 0:
         areaAbove = MyErf((AveNumInt-binning.edges)/Sqrt2/RMSInt)
         
         ProbFromRMS = .5*(areaAbove[:-1]-areaAbove[1:])
     else:
         ProbFromRMS = np.zeros(hContents.shape)
         obs = binning.find(AveNumInt)
         if ( obs < binning.num ) and ( AveNumInt >= 1.0E-5 ): # just ignore zero values
             ProbFromRMS[obs] = 1.0
         
     if calcOption == 'true':  # Just put distribution into histogram
         if RMSInt > 0:
             hContents += ProbFromRMS*LSintLumi
             totalProb = np.sum(ProbFromRMS)
 
             if 1.0-totalProb > 0.01:
                 print("Run %d, LS %d: Significant probability density outside of your histogram (mean %.2f," % (run, ls, AveNumInt))
                 print("rms %.2f, integrated probability %.3f). Consider using a higher value of --maxPileupBin." % (RMSInt, totalProb))
         else:
             hContents[obs] += LSintLumi 
     else: # have to convolute with a poisson distribution to get observed Nint
         if not hasattr(binning, "poissConv"): 
             binning.poissConv = poisson(
                 binning.edges[:-1,np.newaxis], 
                 binning.centers[np.newaxis,:]) 
         prob = np.sum(binning.poissConv * ProbFromRMS[np.newaxis,:], axis=1)
         hContents += prob*LSintLumi
         #if ( not np.all(ProbFromRMS == 0) ) and 1.0-np.sum(prob) > 0.01:
         #    print("Run %d, LS %d: significant probability density outside of your histogram, %f" % (run, ls, np.sum(prob)))
         #    print("Consider using a higher value of --maxPileupBin")
 

◆ MyErf()

def pileupCalc.MyErf ( xInput )

Definition at line 24 of file pileupCalc.py.

Referenced by fillPileupHistogram().

 def MyErf(xInput):
     # Abramowitz and Stegun approximations for Erf (equations 7.1.25-28)
     X = np.abs(xInput)
 
     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)*np.exp(-1.0*X*X)
 
     # 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
 
     return np.where(xInput < 0, -cErf, cErf)
 

◆ parseInputFile()

def pileupCalc.parseInputFile ( inputfilename )

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

Definition at line 12 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
 

◆ poisson()

def pileupCalc.poisson	(	x,
		par
	)

Definition at line 49 of file pileupCalc.py.

Referenced by fillPileupHistogram().

 def poisson(x, par):
     
     return np.where(x == 0., np.exp(-par),
            np.exp(x*np.log(par)-loggamma(x+1)-par))
 

Variable Documentation

◆ action

pileupCalc.action

Definition at line 128 of file pileupCalc.py.

◆ binning

pileupCalc.binning

Definition at line 163 of file pileupCalc.py.

Referenced by V0Monitor.bookHisto1D(), TriggerDQMBase.bookME(), TagProbeFitTreeAnalyzer.calculateEfficiency(), JME::JetParameters.createVector(), and HLTMuonMatchAndPlot.fillEdges().

◆ choices

pileupCalc.choices

Definition at line 135 of file pileupCalc.py.

◆ default

pileupCalc.default

Definition at line 139 of file pileupCalc.py.

◆ description

pileupCalc.description

Definition at line 124 of file pileupCalc.py.

◆ dest

pileupCalc.dest

Definition at line 129 of file pileupCalc.py.

◆ float

pileupCalc.float

Definition at line 139 of file pileupCalc.py.

◆ hContents

pileupCalc.hContents

Definition at line 164 of file pileupCalc.py.

◆ help

pileupCalc.help

Definition at line 128 of file pileupCalc.py.

◆ histFile

pileupCalc.histFile

Definition at line 209 of file pileupCalc.py.

◆ inpf

pileupCalc.inpf

Definition at line 166 of file pileupCalc.py.

◆ inputfilecontent

pileupCalc.inputfilecontent

Definition at line 167 of file pileupCalc.py.

◆ inputPileupRange

pileupCalc.inputPileupRange

Definition at line 172 of file pileupCalc.py.

◆ inputRange

pileupCalc.inputRange

Definition at line 168 of file pileupCalc.py.

Referenced by SiStripNoiseNormalizedWithApvGainBuilder.analyze(), SiStripApvGainBuilderFromTag.analyze(), SiStripGain.fillNewGain(), PixelDigiCollection.put(), and CTPPSPixelDigiCollection.put().

◆ int

pileupCalc.int

Definition at line 142 of file pileupCalc.py.

◆ LSPUlist

pileupCalc.LSPUlist

Definition at line 182 of file pileupCalc.py.

◆ lumiInfo

pileupCalc.lumiInfo

Definition at line 187 of file pileupCalc.py.

Referenced by operator<<().

◆ options

pileupCalc.options

Definition at line 150 of file pileupCalc.py.

◆ output

pileupCalc.output

Definition at line 151 of file pileupCalc.py.

◆ parser

pileupCalc.parser

## Main Program

Definition at line 124 of file pileupCalc.py.

◆ pileupHist

pileupCalc.pileupHist

Definition at line 204 of file pileupCalc.py.

◆ req_group

pileupCalc.req_group

Definition at line 127 of file pileupCalc.py.

◆ required

pileupCalc.required

Definition at line 129 of file pileupCalc.py.

◆ Sqrt2

pileupCalc.Sqrt2

Definition at line 67 of file pileupCalc.py.

◆ type

pileupCalc.type

Definition at line 139 of file pileupCalc.py.

◆ VERSION

pileupCalc.VERSION

Definition at line 4 of file pileupCalc.py.

Classes

Functions

Variables

## Main Program

Function Documentation

◆ fillPileupHistogram()

◆ MyErf()

◆ parseInputFile()

◆ poisson()

Variable Documentation

◆ action

◆ binning

◆ choices

◆ default

◆ description

◆ dest

◆ float

◆ hContents

◆ help

◆ histFile

◆ inpf

◆ inputfilecontent

◆ inputPileupRange

◆ inputRange

◆ int

◆ LSPUlist

◆ lumiInfo

◆ options

◆ output

◆ parser

## Main Program

◆ pileupHist

◆ req_group

◆ required

◆ Sqrt2

◆ type

◆ VERSION