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

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

Variables
	parser
	dest
	default
	type
	help
	args = parser.parse_args()

Function Documentation

◆ feetRegionMapping()

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()

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()

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()

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()

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()

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()

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()

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()

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()

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()

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 = parser.parse_args()

Definition at line 563 of file ReadRPCRun2DataFile.py.

◆ default

ReadRPCRun2DataFile.default

Definition at line 557 of file ReadRPCRun2DataFile.py.

◆ dest

ReadRPCRun2DataFile.dest

Definition at line 557 of file ReadRPCRun2DataFile.py.

◆ help

ReadRPCRun2DataFile.help

Definition at line 558 of file ReadRPCRun2DataFile.py.

◆ parser

ReadRPCRun2DataFile.parser

Initial value:

1= argparse.ArgumentParser( description=__doc__,

2 formatter_class = argparse.RawTextHelpFormatter)

Definition at line 554 of file ReadRPCRun2DataFile.py.

◆ type

ReadRPCRun2DataFile.type

Definition at line 557 of file ReadRPCRun2DataFile.py.

Classes

Functions

Variables

Function Documentation

◆ feetRegionMapping()

◆ main()

◆ read2015Geometry()

◆ read2015RPCCodeMapping()

◆ read2016RPCGeomData()

◆ roiAsXML()

◆ roiAsXML_2()

◆ rpcMioctAndSectorInfo()

◆ sectorAsXML()

◆ sectorAsXML_2()

◆ updateGeometryTGC()

Variable Documentation

◆ args

◆ default

◆ dest

◆ help

◆ parser

◆ type