pileupCalc Namespace Reference

Functions
def	fillPileupHistogram (lumiInfo, calcOption, hist, minbXsec, Nbins, run, ls)
def	MyErf (input)
def	parseInputFile (inputfilename)

Variables
	action
	choices
	default
	description
	dest
	float
	help
	histFile
	inpf
	inputfilecontent
	inputPileupRange
	inputRange
	int
	LSPUlist
	lumiInfo
	nbins
	options
	output
	parser
	pileupHist
	req_group
	required
	type
	upper
	VERSION

Function Documentation

fillPileupHistogram()

def pileupCalc.fillPileupHistogram	(		lumiInfo,
			calcOption,
			hist,
			minbXsec,
			Nbins,
			run,
			ls
	)

lumiinfo:[intlumi per LS, mean interactions ]

intlumi is the deadtime corrected average integrated lumi per lumisection

Definition at line 57 of file pileupCalc.py.

def fillPileupHistogram (lumiInfo, calcOption, hist, minbXsec, Nbins, 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
 
    #coeff = 0
 
    #if RMSInt > 0:
    #    coeff = 1.0/RMSInt/sqrt(6.283185)
 
    #expon = 2.0*RMSInt*RMSInt
 
    Sqrt2 = sqrt(2)
 
    
 
    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("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:
            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("Run %d, LS %d: significant probability density outside of your histogram" % (run, ls))
            print("Consider using a higher value of --maxPileupBin")
 
 
    return hist
 
 
 

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

MyErf()

def pileupCalc.MyErf

(

input

)

Definition at line 25 of file pileupCalc.py.

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
 
 

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

Referenced by fillPileupHistogram().

parseInputFile()

def pileupCalc.parseInputFile

(

inputfilename

)

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

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

action

pileupCalc.action

Definition at line 173 of file pileupCalc.py.

choices

pileupCalc.choices

Definition at line 180 of file pileupCalc.py.

default

pileupCalc.default

Definition at line 184 of file pileupCalc.py.

description

pileupCalc.description

Definition at line 169 of file pileupCalc.py.

dest

pileupCalc.dest

Definition at line 174 of file pileupCalc.py.

float

pileupCalc.float

Definition at line 184 of file pileupCalc.py.

help

pileupCalc.help

Definition at line 173 of file pileupCalc.py.

histFile

pileupCalc.histFile

Definition at line 251 of file pileupCalc.py.

inpf

pileupCalc.inpf

Definition at line 216 of file pileupCalc.py.

inputfilecontent

pileupCalc.inputfilecontent

Definition at line 217 of file pileupCalc.py.

inputPileupRange

pileupCalc.inputPileupRange

Definition at line 222 of file pileupCalc.py.

inputRange

pileupCalc.inputRange

Definition at line 218 of file pileupCalc.py.

Referenced by SiStripNoiseNormalizedWithApvGainBuilder.analyze(), SiStripApvGainBuilderFromTag.analyze(), SiStripGain.fillNewGain(), cms::ReconstructerFP420.produce(), cms::TrackerizerFP420.produce(), cms::DigitizerFP420.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().

int

pileupCalc.int

Definition at line 187 of file pileupCalc.py.

LSPUlist

pileupCalc.LSPUlist

Definition at line 232 of file pileupCalc.py.

lumiInfo

pileupCalc.lumiInfo

Definition at line 237 of file pileupCalc.py.

Referenced by operator<<().

nbins

pileupCalc.nbins

Definition at line 213 of file pileupCalc.py.

options

pileupCalc.options

Definition at line 195 of file pileupCalc.py.

output

pileupCalc.output

Definition at line 196 of file pileupCalc.py.

parser

pileupCalc.parser

autotoc_md311

## Main Program

autotoc_md313

Definition at line 169 of file pileupCalc.py.

pileupHist

pileupCalc.pileupHist

Definition at line 209 of file pileupCalc.py.

req_group

pileupCalc.req_group

Definition at line 172 of file pileupCalc.py.

required

pileupCalc.required

Definition at line 174 of file pileupCalc.py.

type

pileupCalc.type

Definition at line 184 of file pileupCalc.py.

upper

pileupCalc.upper

Definition at line 214 of file pileupCalc.py.

Referenced by ntupleDataFormat.TrackingParticleMatchInfo.__getattr__(), ntupleDataFormat.TrackMatchInfo.__getattr__(), TrackValidation_cff._algoToSelector(), trackingPlots._summaryBinRename(), data_sources.json_list.as_table(), MomentumDependentPedeLabeler.buildMomentumDependencyMap(), VectorHitBuilderAlgorithm.buildVectorHit(), clip(), tmtt::HTbase.convertCoordRangeToBinRange(), ConstrainedTreeBuilder.covarianceMatrix(), ConstrainedTreeBuilderT.covarianceMatrix(), InterpolatedPulse< 1500U >.derivative(), TableParser.extractPages(), MuonShowerInformationFiller.fillHitsByStation(), cond::persistency::RUN_INFO::Table.getInclusiveRunRange(), cond::persistency::RUN_INFO::Table.getInclusiveTimeRange(), SSDigitizerAlgorithm.getSignalScale(), Grid1D.Grid1D(), MagVolumeOutsideValidity.MagVolumeOutsideValidity(), MVAVariableManager< reco::GsfElectron >.MVAVariableManager(), HipPyOptionParser.parseBoolString(), calibratedEgammas_cff.prefixName(), LHCInfoImpl::LumiSectionFilter.process(), runTauIdMVA.TauIDEmbedder.processDeepProducer(), InterpolatedPulse< 1500U >.secondDerivative(), VectorHit.sharesClusters(), jetTools.AddJetCollection.toolCode(), jetTools.UpdateJetCollection.toolCode(), and MFGrid3D.valueInTesla().

VERSION

pileupCalc.VERSION

Definition at line 4 of file pileupCalc.py.