test/python/runHGCalTBTest_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('SIM')
0004
0005 # import of standard configurations
0006 process.load("FWCore.MessageService.MessageLogger_cfi")
0007 process.load('Configuration.StandardSequences.Services_cff')
0008 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0009 process.load('Configuration.EventContent.EventContent_cff')
0010 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0011 process.load('Geometry.HGCalCommonData.testTB181V1XML_cfi')
0012 process.load('Geometry.HGCalCommonData.hgcalNumberingInitialization_cfi')
0013 process.load('Geometry.HGCalCommonData.hgcalParametersInitialization_cfi')
0014 process.load('Configuration.StandardSequences.MagneticField_0T_cff')
0015 process.load('Configuration.StandardSequences.Generator_cff')
0016 process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_cfi')
0017 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0018 process.load('Configuration.StandardSequences.SimIdeal_cff')
0019 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0020 process.load('SimG4CMS.HGCalTestBeam.HGCalTBAnalyzer_cfi')
0021 process.load('SimG4CMS.HGCalTestBeam.HGCalTBCheckGunPosition_cfi')
0022 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0023
0024 process.maxEvents = cms.untracked.PSet(
0025     input = cms.untracked.int32(1)
0026 )
0027
0028 if 'MessageLogger' in process.__dict__:
0029     process.MessageLogger.HGCalGeom=dict()
0030     process.MessageLogger.SimG4CoreGeometry=dict()
0031     process.MessageLogger.HGCSim=dict()
0032
0033 # Input source
0034 process.source = cms.Source("EmptySource")
0035
0036 process.options = cms.untracked.PSet(
0037 )
0038
0039 # Production Info
0040 process.configurationMetadata = cms.untracked.PSet(
0041     annotation = cms.untracked.string('SingleMuonE200_cfi nevts:10'),
0042     name = cms.untracked.string('Applications'),
0043     version = cms.untracked.string('$Revision: 1.19 $')
0044 )
0045
0046 # Additional output definition
0047 process.TFileService = cms.Service("TFileService",
0048                                    fileName = cms.string('TBGenSim.root')
0049 )
0050
0051 # Other statements
0052 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0053 from Configuration.AlCa.GlobalTag import GlobalTag
0054 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0055
0056 process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
0057     AddAntiParticle = cms.bool(False),
0058     PGunParameters = cms.PSet(
0059         MinE = cms.double(99.99),
0060         MaxE = cms.double(100.01),
0061         MinTheta = cms.double(0.0),
0062         MaxTheta = cms.double(0.0),
0063         MinPhi = cms.double(-3.14159265359),
0064         MaxPhi = cms.double(3.14159265359),
0065         PartID = cms.vint32(211)
0066     ),
0067     Verbosity = cms.untracked.int32(0),
0068     firstRun = cms.untracked.uint32(1),
0069     psethack = cms.string('single pion E 100')
0070 )
0071
0072 process.VtxSmeared.MinZ                 = -800.0
0073 process.VtxSmeared.MaxZ                 = -800.0
0074 process.VtxSmeared.MinX                 = -7.5
0075 process.VtxSmeared.MaxX                 =  7.5
0076 process.VtxSmeared.MinY                 = -7.5
0077 process.VtxSmeared.MaxY                 =  7.5
0078 process.g4SimHits.HGCSD.RejectMouseBite = True
0079 process.g4SimHits.HGCSD.RotatedWafer    = True
0080 process.HGCalTBAnalyzer.doDigis         = False
0081 process.HGCalTBAnalyzer.doRecHits       = False
0082 process.HGCalTBAnalyzer.useFH           = True
0083 process.HGCalTBAnalyzer.useEE           = False
0084 process.HGCalTBAnalyzer.useBeam         = False
0085 process.HGCalTBAnalyzer.zFrontEE        = 1110.0
0086 process.HGCalTBAnalyzer.zFrontFH        = 1176.5
0087 process.HGCalTBAnalyzer.zFrontFH        = 1307.5
0088 process.HGCalTBAnalyzer.maxDepth        = 39
0089 process.HGCalTBAnalyzer.deltaZ          = 26.2
0090 process.HGCalTBAnalyzer.zFirst          = 22.8
0091
0092 # Path and EndPath definitions
0093 process.generation_step = cms.Path(process.pgen)
0094 process.simulation_step = cms.Path(process.psim)
0095 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0096 process.gunfilter_step  = cms.Path(process.HGCalTBCheckGunPostion)
0097 process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
0098 process.endjob_step = cms.EndPath(process.endOfProcess)
0099
0100 # Schedule definition
0101 process.schedule = cms.Schedule(process.generation_step,
0102                                 process.genfiltersummary_step,
0103                                 process.simulation_step,
0104                                 process.gunfilter_step,
0105                                 process.analysis_step,
0106                                 process.endjob_step
0107 )
0108
0109 # filter all path with the production filter sequence
0110 for path in process.paths:
0111         getattr(process,path)._seq = process.generator * getattr(process,path)._seq