Reference Manual

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

MCScenario_CRAFT1_22X.cmsRun

def cmsRun

Definition: MCScenario_CRAFT1_22X.py:65

def MCScenario_CRAFT1_22X.CSCsorter	(	a,
		b
	)

Definition at line 74 of file MCScenario_CRAFT1_22X.py.

75 def CSCsorter(a, b): return cmp(CSCpreferred_order[a], CSCpreferred_order[b])

# an instance of this class corresponds to one <DTChamber ... /> or <CSCStation ... />, etc.

MCScenario_CRAFT1_22X.CSCsorter

def CSCsorter

Definition: MCScenario_CRAFT1_22X.py:74

def MCScenario_CRAFT1_22X.define_scenario ( )

Definition at line 147 of file MCScenario_CRAFT1_22X.py.

References make_scenario_sqlite(), and sistrip::SpyUtilities.range().

Referenced by make_scenario_sqlite().

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

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

MCScenario_CRAFT1_22X.DTsorter

def DTsorter

Definition: MCScenario_CRAFT1_22X.py:73

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.

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

Referenced by define_scenario().

35 def make_scenario_sqlite():

36 scenario = define_scenario()

37 write_xml(scenario, "MCScenario_CRAFT1_22X.xml")

38 write_conversion_cfg("convert_cfg.py", "MCScenario_CRAFT1_22X.xml", "MCScenario_CRAFT1_22X.db")

39 cmsRun("convert_cfg.py")

40 write_check_cfg("check_cfg.py", "MCScenario_CRAFT1_22X.db", "MCScenario_CRAFT1_22X_CHECKME.xml")

cmsRun("check_cfg.py")

MCScenario_CRAFT1_22X.define_scenario

def define_scenario

Definition: MCScenario_CRAFT1_22X.py:147

MCScenario_CRAFT1_22X.write_check_cfg

def write_check_cfg

Definition: MCScenario_CRAFT1_22X.py:51

MCScenario_CRAFT1_22X.write_xml

def write_xml

Definition: MCScenario_CRAFT1_22X.py:125

MCScenario_CRAFT1_22X.make_scenario_sqlite

def make_scenario_sqlite

called once at the end of this script

Definition: MCScenario_CRAFT1_22X.py:34

MCScenario_CRAFT1_22X.write_conversion_cfg

def write_conversion_cfg

that's it! everything this uses is defined below

Definition: MCScenario_CRAFT1_22X.py:43

MCScenario_CRAFT1_22X.cmsRun

def cmsRun

Definition: MCScenario_CRAFT1_22X.py:65

def MCScenario_CRAFT1_22X.sorter	(	a,
		b
	)

Definition at line 97 of file MCScenario_CRAFT1_22X.py.

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

98 def sorter(a, b): return cmp(preferred_order[a], preferred_order[b])

# an instance of this class corresponds to one <setposition ... />

MCScenario_CRAFT1_22X.sorter

def sorter

Definition: MCScenario_CRAFT1_22X.py:97

def MCScenario_CRAFT1_22X.write_check_cfg	(	fileName,
		dbFileName,
		xmlFileName
	)

Definition at line 51 of file MCScenario_CRAFT1_22X.py.

References mergeVDriftHistosByStation.file.

Referenced by make_scenario_sqlite().

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

mergeVDriftHistosByStation.file

tuple file

Definition: mergeVDriftHistosByStation.py:55

MCScenario_CRAFT1_22X.write_check_cfg

def write_check_cfg

Definition: MCScenario_CRAFT1_22X.py:51

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.

References mergeVDriftHistosByStation.file.

Referenced by make_scenario_sqlite().

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

mergeVDriftHistosByStation.file

tuple file

Definition: mergeVDriftHistosByStation.py:55

MCScenario_CRAFT1_22X.write_conversion_cfg

def write_conversion_cfg

that's it! everything this uses is defined below

Definition: MCScenario_CRAFT1_22X.py:43

def MCScenario_CRAFT1_22X.write_xml	(	scenario,
		fileName
	)

Definition at line 125 of file MCScenario_CRAFT1_22X.py.

References mergeVDriftHistosByStation.file, and join().

Referenced by make_scenario_sqlite().

 
 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)

mergeVDriftHistosByStation.file

tuple file

Definition: mergeVDriftHistosByStation.py:55

MCScenario_CRAFT1_22X.write_xml

def write_xml

Definition: MCScenario_CRAFT1_22X.py:125

join

static std::string join(char **cmd)

Definition: RemoteFile.cc:19

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.

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.

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.

Classes

Functions

Variables

Function Documentation

Variable Documentation

Classes
class	Alignable

class	CSCChamber

class	DTChamber

class	Operation

class	Position

Functions
def	cmsRun

def	CSCsorter

def	define_scenario

def	DTsorter

def	make_scenario_sqlite
	called once at the end of this script More...

def	sorter

def	write_check_cfg

def	write_conversion_cfg
	that's it! everything this uses is defined below More...

def	write_xml

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}