test/Visualization/displayMFGeom_cfg.py

0001 #################################################
0002 #
0003 # Please run the script with cmsRun:
0004 #
0005 # cmsRun displayMFGeom_cfg.py
0006 #
0007 #################################################
0008
0009 import FWCore.ParameterSet.Config as cms
0010
0011 process = cms.Process("DISPLAY")
0012
0013 process.load("MagneticField.Engine.volumeBasedMagneticField_160812_cfi")
0014
0015 MFGeom=True;
0016
0017 process.XMLIdealGeometryESSource = cms.ESSource("XMLIdealGeometryESSource",
0018       geomXMLFiles = cms.vstring('Geometry/CMSCommonData/data/normal/cmsextent.xml',
0019         'Geometry/CMSCommonData/data/cms.xml',
0020         'Geometry/CMSCommonData/data/cmsMagneticField.xml',
0021         # for 090322 and older
0022 #        'MagneticField/GeomBuilder/data/MagneticFieldVolumes_1103l.xml',
0023
0024         # for 160812
0025         'MagneticField/GeomBuilder/data/MagneticFieldVolumes_160812_1.xml',
0026         'MagneticField/GeomBuilder/data/MagneticFieldVolumes_160812_2.xml',
0027
0028         'Geometry/CMSCommonData/data/materials.xml'),
0029  #     rootNodeName = cms.string('cmsMagneticField:MAGF')
0030       rootNodeName = cms.string('cms:World')
0031 )
0032
0033
0034
0035 process.maxEvents = cms.untracked.PSet(
0036     input = cms.untracked.int32(1)
0037     )
0038
0039
0040 process.source = cms.Source("EmptySource")
0041
0042 process.EveService = cms.Service("EveService")
0043
0044 ### Extractor of geometry needed to display it in Eve.
0045 ### Required for "DummyEvelyser".
0046 process.add_( cms.ESProducer(
0047         "TGeoMgrFromDdd",
0048         verbose = cms.untracked.bool(False),
0049         level   = cms.untracked.int32(8)
0050 ))
0051
0052 #Z_SECTION=-165  #S3 chimney
0053 #Z_SECTION=165   #S4 chimney
0054 #Z_SECTION=760  #YE1
0055 #Z_SECTION=880  #YE2
0056 #Z_SECTION=979  #YE3
0057 #Z_SECTION=1069 #YE4 edge
0058 #Z_SECTION=1072 #YE4
0059 #Z_SECTION=-1519.4957 #CASTOR
0060
0061 #Z_SECTION=1434 #?
0062 #Z_SECTION=1548 #?
0063
0064
0065
0066
0067 process.dump = cms.EDAnalyzer("DisplayGeom",
0068     nodes = cms.untracked.vstring("cmsMagneticField:MAGF_1"),
0069     level   = cms.untracked.int32(2),
0070     MF_component = cms.untracked.string("B"), #Field map: None, B, AbsBZ, AbsBR, AvsBphi, BR, Bphi
0071     MF_pickable = cms.untracked.bool(False),     #Field map: pickable values
0072
0073 #Field map on ZY plane @ X=0 (Top-right quarter)
0074 #    MF_plane_d0 = cms.untracked.vdouble(0, 0, 0),
0075 #    MF_plane_d1 = cms.untracked.vdouble(0, 0., 2000),
0076 #    MF_plane_d2 = cms.untracked.vdouble(0, 900., 0.),
0077
0078 #Field map on ZY plane @ X=0 (full CMS)
0079     MF_plane_d0 = cms.untracked.vdouble(0, -900, -2400),
0080     MF_plane_d1 = cms.untracked.vdouble(0, -900., 2400),
0081     MF_plane_d2 = cms.untracked.vdouble(0, 900., -2400.),
0082
0083 #Field map on XY plane @ Z=0
0084 #    MF_plane_d0 = cms.untracked.vdouble(-900, -900, 0.),
0085 #    MF_plane_d1 = cms.untracked.vdouble(-900, 900, 0.),
0086 #    MF_plane_d2 = cms.untracked.vdouble(900, -900, 0.),
0087
0088 #Field map on XY plane @ S3chimney
0089 #    MF_plane_d0 = cms.untracked.vdouble(-900, -900, -165.),
0090 #    MF_plane_d1 = cms.untracked.vdouble(-900, 900, -165.),
0091 #    MF_plane_d2 = cms.untracked.vdouble(900, -900, -165.),
0092
0093 #Field map on XY plane @ S11 feet
0094 #    MF_plane_d0 = cms.untracked.vdouble(0, -900, 0),
0095 #    MF_plane_d1 = cms.untracked.vdouble(0, 0, 0),
0096 #    MF_plane_d2 = cms.untracked.vdouble(900, -900, 0),
0097
0098 #XY plane @ Z_SECTION
0099 #    MF_plane_d0 = cms.untracked.vdouble(-900, -900, Z_SECTION),
0100 #    MF_plane_d1 = cms.untracked.vdouble(-900, 900, Z_SECTION),
0101 #    MF_plane_d2 = cms.untracked.vdouble(900, -900, Z_SECTION),
0102
0103 #CASTOR detail
0104 #    MF_plane_d0 = cms.untracked.vdouble(-200, -200, -1519.49571429),
0105 #    MF_plane_d1 = cms.untracked.vdouble(-200, 200, -1519.49571429),
0106 #    MF_plane_d2 = cms.untracked.vdouble(200, -200, -1519.49571429),
0107
0108
0109     MF_plane_N  = cms.untracked.int32(1000), #Field map bins
0110 #    MF_plane_N2  = cms.untracked.int32(500),
0111
0112     MF_plane_draw_dir =  cms.untracked.int32(False)
0113 )
0114
0115 process.p = cms.Path(process.dump)
0116