MCScenario_CRAFT1_22X Namespace Reference

Classes
class	Alignable
class	CSCChamber
class	DTChamber
class	Operation
class	Position

Functions
def	cmsRun (fileName)
def	CSCsorter (a, b)
def	define_scenario ()
def	DTsorter (a, b)
def	make_scenario_sqlite ()
	called once at the end of this script More...
def	sorter (a, b)
def	write_check_cfg (fileName, dbFileName, xmlFileName)
def	write_conversion_cfg (fileName, xmlFileName, dbFileName)
	that's it! everything this uses is defined below More...
def	write_xml (scenario, fileName)

Variables
dictionary	CSCpreferred_order = {"endcap":1, "station":2, "ring":3, "chamber":4, "layer":5}
dictionary	DTpreferred_order = {"wheel":1, "station":2, "sector":3, "superlayer":4, "layer":5}
	writing a scenario in XML ############################################################## More...
dictionary	preferred_order = {"x":1, "y":2, "z":3, "phix":4, "phiy":5, "phiz":6}

Function Documentation

cmsRun()

def MCScenario_CRAFT1_22X.cmsRun

(

fileName

)

Definition at line 65 of file MCScenario_CRAFT1_22X.py.

def cmsRun(fileName):
    os.system("cmsRun %(fileName)s" % vars())
 

Referenced by make_scenario_sqlite().

CSCsorter()

def MCScenario_CRAFT1_22X.CSCsorter	(		a,
			b
	)

Definition at line 74 of file MCScenario_CRAFT1_22X.py.

def CSCsorter(a, b): return cmp(CSCpreferred_order[a], CSCpreferred_order[b])
 
# an instance of this class corresponds to one <DTChamber ... /> or <CSCStation ... />, etc.

define_scenario()

def MCScenario_CRAFT1_22X.define_scenario

(

)

Definition at line 147 of file MCScenario_CRAFT1_22X.py.

def define_scenario():
    # this will be a list of operations to write to an XML file
    scenario = []
 
    # Uncertainty in DT chamber positions comes in two parts:
    #    1. positions within sectors
    #    2. positions of the sector-groups
 
    # Aligned chambers (wheels -1, 0, +1 except sectors 1 and 7)
    # uncertainty within sectors:
    #    x: 0.08 cm (from segment-matching)  phix: 0.0007 rad (from MC)
    #    y: 0.10 cm (from MC)                phiy: 0.0007 rad (from segment-matching)
    #    z: 0.10 cm (from MC)                phiz: 0.0003 rad (from MC)
    # uncertainty of sector-groups (depends on choice of pT cut, not well understood):
    #    x: 0.05 cm
 
    # Unaligned chambers uncertainty within sectors:
    #    x: 0.08 cm (same as above)          phix: 0.0016 rad
    #    y: 0.24 cm                          phiy: 0.0021 rad
    #    z: 0.42 cm with a -0.35 cm bias     phiz: 0.0010 rad
    # uncertainty of sector-groups:
    #    x: 0.65 cm
    # These come from actual alignments measured in the aligned
    # chambers (we assume that the unaligned chambers have
    # misalignments on the same scale)
 
    # Also, superlayer z uncertainty is 0.054 cm
 
    # Before starting, let's build a list of chambers
    DTchambers = []
    for wheel in -2, -1, 0, 1, 2:
        for station in 1, 2, 3, 4:
            if station == 4: nsectors = 14
            else: nsectors = 12
            for sector in range(1, nsectors+1):
                DTchambers.append(DTChamber(wheel = wheel, station = station, sector = sector))
 
    # the superlayers
    for dtchamber in DTchambers:
        for superlayer in 1, 2, 3:
            if superlayer == 2 and dtchamber.station == 4: continue
 
            alignable = Alignable("DTSuperLayer", wheel = dtchamber.wheel, station = dtchamber.station, sector = dtchamber.sector, superlayer = superlayer)
            position = Position(x = 0, y = 0, z = random.gauss(0, 0.054), phix = 0, phiy = 0, phiz = 0)
            scenario.append(Operation(alignable, position))
 
    sector_errx = {}
 
    # sector-groups for aligned chambers:
    for wheel in -1, 0, 1:
        for sector in 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14:
            sector_errx[wheel, sector] = random.gauss(0., 0.05)
 
    # sector-groups for unaligned chambers:
    for wheel in -1, 0, 1:
        for sector in 1, 7:
            sector_errx[wheel, sector] = random.gauss(0., 0.65)
    for wheel in -2, 2:
        for sector in 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14:
            sector_errx[wheel, sector] = random.gauss(0., 0.65)
 
    for dtchamber in DTchambers:
        # within sectors for aligned chambers:
        if dtchamber.wheel in (-1, 0, 1) and dtchamber.sector in (2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14):
            errx = random.gauss(0, 0.08)
            erry = random.gauss(0, 0.10)
            errz = random.gauss(0, 0.10)
            errphix = random.gauss(0, 0.0007)
            errphiy = random.gauss(0, 0.0007)
            errphiz = random.gauss(0, 0.0003)
 
        # within sectors for unaligned chambers:
        else:
            errx = random.gauss(0, 0.08)
            erry = random.gauss(0, 0.24)
            errz = random.gauss(-0.35, 0.42)
            errphix = random.gauss(0, 0.0016)
            errphiy = random.gauss(0, 0.0021)
            errphiz = random.gauss(0, 0.0010)
 
        errx += sector_errx[dtchamber.wheel, dtchamber.sector]
 
        # now turn this into an operation
        alignable = Alignable("DTChamber", wheel = dtchamber.wheel, station = dtchamber.station, sector = dtchamber.sector)
        position = Position(x = errx, y = erry, z = errz, phix = errphix, phiy = errphiy, phiz = errphiz)
        scenario.append(Operation(alignable, position))
 
    # Uncertainty in CSC chamber positions comes in 5 parts:
    #    1. 0.0092 cm layer x misalignments observed with beam-halo tracks
    #    2. isotropic photogrammetry uncertainty of 0.03 cm (x, y, z) and 0.00015 rad in phiz
    #    3. 0.0023 rad phiy misalignment observed with beam-halo tracks
    #    4. 0.1438 cm z and 0.00057 rad phix uncertainty between rings from SLM (from comparison in 0T data with PG)
    #    5. 0.05 cm (x, y, z) disk misalignments and 0.0001 rad rotation around beamline
 
    # Before starting, let's build a list of chambers
    CSCchambers = []
    for endcap in 1, 2:
        for station, ring in (1, 1), (1, 2), (1, 3), (1, 4), (2, 1), (2, 2), (3, 1), (3, 2), (4, 1):
            if station > 1 and ring == 1:
                nchambers = 18
            else:
                nchambers = 36
 
            for chamber in range(1, nchambers+1):
                CSCchambers.append(CSCChamber(endcap = endcap, station = station, ring = ring, chamber = chamber))
            
    # First, the layer uncertainties: x only for simplicity, observed 0.0092 cm in overlaps alignment test
    for chamber in CSCchambers:
        for layer in 1, 2, 3, 4, 5, 6:
            alignable = Alignable("CSCLayer", endcap = chamber.endcap, station = chamber.station, ring = chamber.ring, chamber = chamber.chamber, layer = layer)
            position = Position(x = random.gauss(0, 0.0092), y = 0, z = 0, phix = 0, phiy = 0, phiz = 0)
            scenario.append(Operation(alignable, position))
 
    # Next, the ring errors from DCOPS (derived from comparison with photogrammetry)
    CSCrings = []
    for endcap in 1, 2:
        for station, ring in (1, 1), (1, 2), (1, 3), (1, 4), (2, 1), (2, 2), (3, 1), (3, 2), (4, 1):
            CSCrings.append(CSCChamber(endcap = endcap, station = station, ring = ring, z = random.gauss(0, 0.1438), phix = random.gauss(0, 0.00057)))
 
    # Next, the chamber errors
    for chamber in CSCchambers:
        errx = random.gauss(0, 0.03)
        erry = random.gauss(0, 0.03)
        errz = random.gauss(0, 0.03)
        errphix = random.gauss(0, 0.00057)
        errphiy = random.gauss(0, 0.0023)
        errphiz = random.gauss(0, 0.00015)
        
        for ring in CSCrings:
            if ring.endcap == chamber.endcap and ring.station == chamber.station and ring.ring == chamber.ring:
                errz += ring.z
                errphix += ring.phix
                break
 
        alignable = Alignable("CSCChamber", endcap = chamber.endcap, station = chamber.station, ring = chamber.ring, chamber = chamber.chamber)
        position = Position(x = errx, y = erry, z = errz, phix = errphix, phiy = errphiy, phiz = errphiz)
        scenario.append(Operation(alignable, position))
 
    # Finally, the disk errors
    for endcap in 1, 2:
        for station in 1, 2, 3, 4:
            alignable = Alignable("CSCStation", endcap = endcap, station = station)
            position = Position(x = random.gauss(0, 0.05), y = random.gauss(0, 0.05), z = random.gauss(0, 0.05), phix = 0., phiy = 0., phiz = random.gauss(0, 0.0001))
            scenario.append(Operation(alignable, position))
 
    return scenario
 
# run it all!
make_scenario_sqlite()

References make_scenario_sqlite(), and FastTimerService_cff.range.

Referenced by make_scenario_sqlite().

DTsorter()

def MCScenario_CRAFT1_22X.DTsorter	(		a,
			b
	)

Definition at line 73 of file MCScenario_CRAFT1_22X.py.

def DTsorter(a, b): return cmp(DTpreferred_order[a], DTpreferred_order[b])

make_scenario_sqlite()

def MCScenario_CRAFT1_22X.make_scenario_sqlite

(

)

called once at the end of this script

Definition at line 34 of file MCScenario_CRAFT1_22X.py.

def make_scenario_sqlite():
    scenario = define_scenario()
    write_xml(scenario, "MCScenario_CRAFT1_22X.xml")
    write_conversion_cfg("convert_cfg.py", "MCScenario_CRAFT1_22X.xml", "MCScenario_CRAFT1_22X.db")
    cmsRun("convert_cfg.py")
    write_check_cfg("check_cfg.py", "MCScenario_CRAFT1_22X.db", "MCScenario_CRAFT1_22X_CHECKME.xml")
    cmsRun("check_cfg.py")

References cmsRun(), define_scenario(), write_check_cfg(), write_conversion_cfg(), and write_xml().

Referenced by define_scenario().

sorter()

def MCScenario_CRAFT1_22X.sorter	(		a,
			b
	)

Definition at line 97 of file MCScenario_CRAFT1_22X.py.

def sorter(a, b): return cmp(preferred_order[a], preferred_order[b])
 
# an instance of this class corresponds to one <setposition ... />

Referenced by PFMuonAlgo.getTrackWithSmallestError(), PFHFRecHitCreator.importRecHits(), L1TMuonBarrelKalmanRegionModule.process(), L2TauJetsMerger.produce(), EcalZeroSuppressionProducer.produce(), SimpleForwardNavigableLayer.setInwardLinks(), and SimpleBarrelNavigableLayer.setInwardLinks().

write_check_cfg()

def MCScenario_CRAFT1_22X.write_check_cfg	(		fileName,
			dbFileName,
			xmlFileName
	)

Definition at line 51 of file MCScenario_CRAFT1_22X.py.

def write_check_cfg(fileName, dbFileName, xmlFileName):
    outfile = file(fileName, "w")
    outfile.write("""
from Alignment.MuonAlignment.convertSQLitetoXML_cfg import *
process.PoolDBESSource.connect = "sqlite_file:%(dbFileName)s"
process.MuonGeometryDBConverter.outputXML.fileName = "%(xmlFileName)s"
process.MuonGeometryDBConverter.outputXML.relativeto = "ideal"
process.MuonGeometryDBConverter.outputXML.suppressDTChambers = False
process.MuonGeometryDBConverter.outputXML.suppressDTSuperLayers = False
process.MuonGeometryDBConverter.outputXML.suppressDTLayers = True
process.MuonGeometryDBConverter.outputXML.suppressCSCChambers = False
process.MuonGeometryDBConverter.outputXML.suppressCSCLayers = False
""" % vars())
 

References geometryDiff.file.

Referenced by make_scenario_sqlite().

write_conversion_cfg()

def MCScenario_CRAFT1_22X.write_conversion_cfg	(		fileName,
			xmlFileName,
			dbFileName
	)

that's it! everything this uses is defined below

Definition at line 43 of file MCScenario_CRAFT1_22X.py.

def write_conversion_cfg(fileName, xmlFileName, dbFileName):
    outfile = file(fileName, "w")
    outfile.write("""
from Alignment.MuonAlignment.convertXMLtoSQLite_cfg import *
process.MuonGeometryDBConverter.fileName = "%(xmlFileName)s"
process.PoolDBOutputService.connect = "sqlite_file:%(dbFileName)s"
""" % vars())
 

References geometryDiff.file.

Referenced by make_scenario_sqlite().

write_xml()

def MCScenario_CRAFT1_22X.write_xml	(		scenario,
			fileName
	)

Definition at line 125 of file MCScenario_CRAFT1_22X.py.

def write_xml(scenario, fileName):
    # a scenario is an ordered list of Operations
    XMLlist = ["<MuonAlignment>\n"]
    for operation in scenario:
        XMLlist.append(operation.writeXML())
    XMLlist.append("</MuonAlignment>\n")
    XMLstring = "".join(XMLlist)
        
    outfile = file(fileName, "w")
    outfile.write(XMLstring)
 

References geometryDiff.file, and join().

Referenced by make_scenario_sqlite().

Variable Documentation

CSCpreferred_order

dictionary MCScenario_CRAFT1_22X.CSCpreferred_order = {"endcap":1, "station":2, "ring":3, "chamber":4, "layer":5}

Definition at line 72 of file MCScenario_CRAFT1_22X.py.

DTpreferred_order

dictionary MCScenario_CRAFT1_22X.DTpreferred_order = {"wheel":1, "station":2, "sector":3, "superlayer":4, "layer":5}

writing a scenario in XML ##############################################################

Definition at line 71 of file MCScenario_CRAFT1_22X.py.

preferred_order

dictionary MCScenario_CRAFT1_22X.preferred_order = {"x":1, "y":2, "z":3, "phix":4, "phiy":5, "phiz":6}

Definition at line 96 of file MCScenario_CRAFT1_22X.py.