studyJets Namespace Reference

Functions
def	getJetProperties (jetSeed, etaFwd, etaCen)
	Python script for studying eta dependence of jet sizes and seed centrality # Author: Aaron Bundock (aaron. More...
def	printJetProperties (etaRange)

Function Documentation

def studyJets.getJetProperties	(		jetSeed,
			etaFwd,
			etaCen
	)

Python script for studying eta dependence of jet sizes and seed centrality # Author: Aaron Bundock (aaron.

*nospamthankyamaam*bundock.ch) # Date: 15/01/18 # #

Definition at line 25 of file studyJets.py.

References funct.abs().

Referenced by printJetProperties().

def getJetProperties(jetSeed,etaFwd,etaCen):

    jetSize = 0
    etaFwdSize = 0
    etaCenSize = 0

    for ring in xrange(jetSeed-etaCen,jetSeed+etaFwd+1):

        if ring < 1:
            ring = abs(ring-1)

        if ring >= len(towerEtaWidths):
            break

        if ring < jetSeed:
            etaFwdSize += towerEtaWidths[ring]
        if ring > jetSeed:
            etaCenSize += towerEtaWidths[ring]

    jetSize = (etaFwdSize + etaCenSize + towerEtaWidths[jetSeed])/0.8
    seedCent = etaCenSize/etaFwdSize

    jetProps = [jetSize, seedCent]

    return jetProps

def studyJets.printJetProperties

(

etaRange

)

Definition at line 52 of file studyJets.py.

References getJetProperties(), reco.if(), edm.print(), and harvestTrackValidationPlots.str.

def printJetProperties(etaRange):

    print "Size  \  eta\t",
    for seedEta in etaRange:
        if(seedEta<29):
            print str(seedEta)+"\t\t",
        else:
            print str(seedEta+1)+"\t\t",
    print

    for size in jetEtaSizes:

        if size%2:
            print "   9x"+str(size)+("\t"),
            for seedEta in etaRange:
                print("\t"),
                etaFwdCen = (size-1)/2
                jetProps = getJetProperties(seedEta,etaFwdCen,etaFwdCen)
                print("%.2f / %.2f" %(jetProps[0],jetProps[1])),
            print

        else:
            print "   9x"+str(size)+" cen",
            for seedEta in etaRange:
                print("\t"),
                etaFwd = size/2-1
                etaCen = size/2
                jetProps = getJetProperties(seedEta, etaFwd, etaCen)
                print("%.2f / %.2f" %(jetProps[0],jetProps[1])),
            print
            print "   9x"+str(size)+" fwd",
            for seedEta in etaRange:
                print("\t"),
                etaFwd = size/2
                etaCen = size/2-1
                jetProps = getJetProperties(seedEta, etaFwd, etaCen)
                print("%.2f / %.2f" %(jetProps[0],jetProps[1])),
            print