ReadRPCRun2DataFile Namespace Reference

Classes
class	Decodes
class	Mioct
class	MuonGeometry
class	ROI
class	Sector
class	TopoCell

Functions
def	rpcMioctAndSectorInfo (sectorID)
def	read2015Geometry (fn)
def	read2015RPCCodeMapping (geometry)
def	feetRegionMapping (sectorId, roiId)
def	read2016RPCGeomData (fn, oldMapping)
def	updateGeometryTGC (newGeometry, oldGeometry)
def	sectorAsXML_2 (sector, depth, stats)
def	roiAsXML_2 (roi, depth)
def	sectorAsXML (xxx_todo_changeme, depth, stats)
def	roiAsXML (roi, depth)
def	main (args)

Variables
	parser
	description
	formatter_class
	dest
	default
	type
	help
	args

Function Documentation

feetRegionMapping()

def ReadRPCRun2DataFile.feetRegionMapping	(		sectorId,
			roiId
	)

provides the new mapping for the feet and elevator region
The numbers can be found in doc/roi_map_R2.pdf

Definition at line 330 of file ReadRPCRun2DataFile.py.

def feetRegionMapping(sectorId, roiId):
    """ provides the new mapping for the feet and elevator region
    The numbers can be found in doc/roi_map_R2.pdf
    """
 
    etacode = -1
    phicode = -1
    
    if sectorId in [23,24,55,56]:
        # elevator region
 
        # eta
        etaMap = { 23 : [ [0,1,2,3,5,7], [4,6,8,9,10,11], range(12,20) + [21,23], [20,22] + range(24,28) ],
                   24 : [ [0,1,2,3,4,6], [5,7,8,9,10,11], range(12,20) + [20,22], [21,23] + range(24,28) ] }
        etaMap[55] = etaMap[24]
        etaMap[56] = etaMap[23]
                
        for ec, tt in enumerate(etaMap[sectorId]):
            if roiId in tt:
                etacode = ec
                break
 
        # phi
        if sectorId in [23,55]:
            if roiId in [2,3,6,7,10,11,18,19,22,23,27]:
                phicode = 2
            if roiId in [0,1,4,5,8,9,12,13,14,15,16,17,20,21,25]:
                phicode = 3
        elif sectorId in [24,56]:
            if roiId in [0,1,4,5,8,9,12,13,14,15,16,17,20,21,24]:
                phicode = 4
            if roiId in [2,3,6,7,10,11,18,19,22,23,26]:
                phicode = 5
 
        if sectorId == 23:
            if roiId == 27: phicode = 2
            elif roiId == 25: phicode = 3
        elif sectorId == 24:
            if roiId == 24: phicode = 4
            elif roiId == 26: phicode = 5
        elif sectorId == 55:
            if roiId == 24: phicode = 3
            elif roiId == 26: phicode = 2
        elif sectorId == 56:
            if roiId == 27: phicode = 5
            elif roiId == 25: phicode = 4
 
 
 
    elif [21,22,25,26,53,54,57,58]:
        # feet region
 
        # eta
        etaMap = { 21 : [ range(0,8), range(8,16) + [17,19], [16,18] + range(20,32) ],
                   22 : [ range(0,8), range(8,16) + [16,18], [17,19] + range(20,32) ] }
        etaMap[25] = etaMap[21]
        etaMap[26] = etaMap[22]
        etaMap[53] = etaMap[22]
        etaMap[54] = etaMap[21]
        etaMap[57] = etaMap[53]
        etaMap[58] = etaMap[54]
                
        for ec, tt in enumerate(etaMap[sectorId]):
            if roiId in tt:
                etacode = ec
                break
 
        # phi
        if sectorId in [21,53]:
            if roiId in [2,3,6,7,10,11,14,15,18,19,22,23,26,27,30,31]:
                phicode = 6
            elif roiId in [0,1,4,5,8,9,12,13,16,17,20,21,24,25,28,29]:
                phicode = 7
        elif sectorId in [22,54]:
            if roiId in [0,1,4,5,8,9,12,13,16,17,20,21,24,25,28,29]:
                phicode = 0
            elif roiId in [2,3,6,7,10,11,14,15,18,19,22,23,26,27,30,31]:
                phicode = 1
        elif sectorId in [25,57]:
            if roiId in [2,3,6,7,10,11,14,15,18,19,22,23,26,27,30,31]:
                phicode = 6
            elif roiId in [0,1,4,5,8,9,12,13,16,17,20,21,24,25,28,29]:
                phicode = 7
        elif sectorId in [26,58]:
            if roiId in [0,1,4,5,8,9,12,13,16,17,20,21,24,25,28,29]:
                phicode = 0
            elif roiId in [2,3,6,7,10,11,14,15,18,19,22,23,26,27,30,31]:
                phicode = 1
 
    if etacode == -1:
        raise RuntimeError("No etacode for SL %i and ROI %i" % (sectorId, roiId))
    if phicode == -1:
        raise RuntimeError("No phicode for SL %i and ROI %i" % (sectorId, roiId))
 
    return (etacode,phicode)
 
 
 
 

main()

def ReadRPCRun2DataFile.main

(

args

)

Definition at line 528 of file ReadRPCRun2DataFile.py.

def main(args):
 
    # read 2015 geometry xml file
    muonGeometry2015 = read2015Geometry( args.infile2015 )
 
    muonGeometry2015.printStats()
 
    # maps roi number to etacode and phicode
    rpcCodeMapping = read2015RPCCodeMapping( muonGeometry2015 )
 
 
    muonGeometry2016 = read2016RPCGeomData( args.infile, rpcCodeMapping )
    muonGeometry2016.printStats()
 
    
    updateGeometryTGC(muonGeometry2016, muonGeometry2015)
    muonGeometry2016.printStats()
 
    # write out the RPC info and the TGC info
    muonGeometry2016.writeXML("TestMioctGeometry2016.xml")
 
 
 
 

read2015Geometry()

def ReadRPCRun2DataFile.read2015Geometry

(

fn

)

Definition at line 299 of file ReadRPCRun2DataFile.py.

def read2015Geometry( fn ):
    xmlgeometry = MioctGeometryXMLReader(fn)
    miocts = []
    for MIOCT in xmlgeometry.getMIOCTs():
        sectors = []
        for sector in MIOCT.Sectors:
            #if sn.startswith('B'): continue
            rois = []
            for roiElem in sector.ROIs:
                rois += [ ROI.fromROIelement(roiElem) ]
            sectors += [ Sector(name=sector['name'], connector=int(sector['connector']), rois=rois) ]
        miocts += [ Mioct( id=int(MIOCT['id']), slot=int(MIOCT['slot']), sectors=sectors) ]
        decodes = Decodes()
        for tc in MIOCT.Decode.TopoCells:
            decodes.addTopoCell( TopoCell.fromTopoCellElement(tc) )
        miocts[-1].decodes=decodes
    return MuonGeometry("full geometry 2015", miocts=miocts)
 
 
 

read2015RPCCodeMapping()

def ReadRPCRun2DataFile.read2015RPCCodeMapping

(

geometry

)

Definition at line 319 of file ReadRPCRun2DataFile.py.

def read2015RPCCodeMapping( geometry ):
    mapping = {}  
    for mioct in geometry.miocts:
        for sector in mioct.sectors:
            if sector.isBarrel():
                mapping[sector.id] = dict( [ (roi.roiid, (roi.etacode, roi.phicode)) for roi in sector.rois ] )
    return mapping
 
 
 
 

read2016RPCGeomData()

def ReadRPCRun2DataFile.read2016RPCGeomData	(		fn,
			oldMapping
	)

Definition at line 429 of file ReadRPCRun2DataFile.py.

def read2016RPCGeomData(fn, oldMapping):
 
    rpcGeometry2016 = MuonGeometry("RPC geometry 2016")
 
    print("read2016RPCGeomData: Reading %s", fn)
    f = open(fn,"read")
    for line in f:
        if line.lstrip().startswith("#"): # remove comments
            continue
        ls = line.split()
        side     = int(ls[0])
        sector   = int(ls[1]) 
        roiid    = int(ls[2]) 
        etamin   = float(ls[3])
        etamax   = float(ls[4])
        phimin   = float(ls[5])
        phimax   = float(ls[6])
        if max(phimin, phimax) < 0:
            phimin += 2 * PI
            phimax += 2 * PI
        eta = (etamin + etamax) / 2
        phi = (phimin + phimax) / 2
 
        sectorId = sector + 32 * side
        if sectorId in [21,22,23,24,25,26,53,54,55,56,57,58]:
            (etacode, phicode) = feetRegionMapping(sectorId, roiid)
        else:
            (etacode, phicode) = oldMapping[sectorId][roiid]
 
        roi = ROI( eta=eta, phi=phi, phimin=phimin, phimax=phimax, etamin=etamin, etamax=etamax, etacode=etacode, phicode=phicode, roiid=roiid )
 
        # now find MIOCT and Sector where to add the ROI to
        mioctId, mioctSlot, sectorConnector, sectorName = rpcMioctAndSectorInfo(sectorId)
        mioct = rpcGeometry2016.getMioct(mioctId)
        if not mioct:
            mioct = rpcGeometry2016.addMioct(Mioct(id=mioctId, slot=mioctSlot))
 
        sector = mioct.getSector(sectorConnector)
        if not sector:
            sector = mioct.addSector( Sector(name=sectorName, connector=sectorConnector) )
            
        sector.addROI( roi )
 
    return rpcGeometry2016
 
 

roiAsXML()

def ReadRPCRun2DataFile.roiAsXML	(		roi,
			depth
	)

Definition at line 514 of file ReadRPCRun2DataFile.py.

def roiAsXML(roi, depth):
    roi['eta'] = (float(roi['etamin']) + float(roi['etamax'])) / 2
    roi['phi'] =  (float(roi['phimin']) + float(roi['phimax'])) / 2
    s = ' ' * depth + '<ROI eta="%f" phi="%f" etacode="0x%s" phicode="0x%s" etamin="%f" etamax="%f" phimin="%f" phimax="%f" roiid="%i"/>\n' % (roi["eta"], roi["phi"], roi["etacode"], roi["phicode"], roi["etamin"], roi["etamax"], roi["phimin"], roi["phimax"], roi["roiid"] )
    return s
 
 
 
 
 
 
 
 
 

roiAsXML_2()

def ReadRPCRun2DataFile.roiAsXML_2	(		roi,
			depth
	)

Definition at line 497 of file ReadRPCRun2DataFile.py.

def roiAsXML_2(roi, depth):
    attr = ["eta", "phi", "etacode", "phicode", "etamin", "etamax", "phimin", "phimax", "roiid"]
    s = ' ' * depth + "<ROI %s/>\n" % (" ".join(['%s="%s"' % (a, roi[a]) for a in attr]) )
    return s
 
 

rpcMioctAndSectorInfo()

def ReadRPCRun2DataFile.rpcMioctAndSectorInfo

(

sectorID

)

Definition at line 285 of file ReadRPCRun2DataFile.py.

def rpcMioctAndSectorInfo(sectorID):
    sectorName = "B%02i" % sectorID
    sectorConnector = (sectorID+2) % 4
 
    mioctID = ( (sectorID+2) % 32 ) / 4
    mioctSlot = mioctID + 4
    if sectorID >= 32:
        mioctID += 8
        mioctSlot = mioctID + 6
 
    return mioctID, mioctSlot, sectorConnector, sectorName
 
 
 

sectorAsXML()

def ReadRPCRun2DataFile.sectorAsXML	(		xxx_todo_changeme,
			depth,
			stats
	)

Definition at line 503 of file ReadRPCRun2DataFile.py.

def sectorAsXML(xxx_todo_changeme, depth, stats):
    (connector, name, rois) = xxx_todo_changeme
    s = ' '*depth + '<Sector connector="%s" name="%s">\n' % (connector, name)
    s += ' '*depth + '    <!-- contains %i ROIs -->\n' % len(rois)
    s += ' '*depth + '    <!-- mapping from ROI to coding scheme -->\n'
    for roi in rois:
        s += roiAsXML(roi, depth + 4 )
    s += ' '*depth + "</Sector>\n"
    stats['rois'] += len(rois)
    return s
 

sectorAsXML_2()

def ReadRPCRun2DataFile.sectorAsXML_2	(		sector,
			depth,
			stats
	)

Definition at line 486 of file ReadRPCRun2DataFile.py.

def sectorAsXML_2(sector, depth, stats):
    attr = ["connector", "name"]
    s = ' '*depth + "<%s %s>\n" % (sector.tag, " ".join(['%s="%s"' % (a, sector[a]) for a in attr]) )
    s += ' '*depth + '    <!-- contains %i ROIs -->\n' % len(sector.ROIs)
    s += ' '*depth + '    <!-- mapping from ROI to coding scheme -->\n'
    for roi in sorted(sector.ROIs, key=lambda roi: int(roi['roiid'])):
        s += roiAsXML_2(roi, depth + 4 )
    s += ' '*depth + "</%s>\n" % sector.tag
    stats['rois'] += len(sector.ROIs)
    return s
 

updateGeometryTGC()

def ReadRPCRun2DataFile.updateGeometryTGC	(		newGeometry,
			oldGeometry
	)

Definition at line 475 of file ReadRPCRun2DataFile.py.

def updateGeometryTGC(newGeometry, oldGeometry):
    newGeometry.name = "updated full geometry for 2016"
    for mioct in newGeometry.miocts:
        oldMioct = oldGeometry.getMioct(mioct.id)
        for sector in oldMioct.sectors:
            if sector.isBarrel(): continue
            mioct.addSector(sector)
        mioct.fillTopoCells()
 
 
 

Variable Documentation

args

ReadRPCRun2DataFile.args

Definition at line 563 of file ReadRPCRun2DataFile.py.

default

ReadRPCRun2DataFile.default

Definition at line 557 of file ReadRPCRun2DataFile.py.

description

ReadRPCRun2DataFile.description

Definition at line 554 of file ReadRPCRun2DataFile.py.

dest

ReadRPCRun2DataFile.dest

Definition at line 557 of file ReadRPCRun2DataFile.py.

formatter_class

ReadRPCRun2DataFile.formatter_class

Definition at line 555 of file ReadRPCRun2DataFile.py.

help

ReadRPCRun2DataFile.help

Definition at line 558 of file ReadRPCRun2DataFile.py.

parser

ReadRPCRun2DataFile.parser

Definition at line 554 of file ReadRPCRun2DataFile.py.

type

ReadRPCRun2DataFile.type

Definition at line 557 of file ReadRPCRun2DataFile.py.