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 26 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 53 of file studyJets.py.

References getJetProperties(), edm.print(), and str.

def printJetProperties(etaRange):

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

    for size in jetEtaSizes:

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

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



print("\n")
print("==============================================================================================================================================================================================================================")
print("Eta-dependence of jet sizes & seed centrality")
print("First number  = jet diameter in eta, normalised to 0.8 (offline)")
print("Second number = jet size on forward side of jet seed / jet size on central side of jet seed")
print("Ideally, best jet choice is where both numbers are closest to unity!")
print("9xN cen = larger area on side of jet further from beam pipe")
print("9xN fwd  = larger area on side of jet closer to beam pipe")
print("==============================================================================================================================================================================================================================")