HGCalTestBeam/test/HGCalTBGenSimDigiReco_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('SIMDIGIRECO')
0004
0005 # import of standard configurations
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0008 process.load('FWCore.MessageService.MessageLogger_cfi')
0009 process.load('Configuration.EventContent.EventContent_cff')
0010 process.load('FWCore.MessageService.MessageLogger_cfi')
0011 process.load('Configuration.EventContent.EventContent_cff')
0012 process.load('SimG4CMS.HGCalTestBeam.HGCalTB160Module4XML_cfi')
0013 process.load('Geometry.HGCalTBCommonData.hgcalTBNumberingInitialization_cfi')
0014 process.load('Geometry.HGCalTBCommonData.hgcalTBParametersInitialization_cfi')
0015 process.load('Geometry.HcalTestBeamData.hcalTB06Parameters_cff')
0016 process.load('Geometry.CaloEventSetup.HGCalTBTopology_cfi')
0017 process.load('Geometry.HGCalGeometry.HGCalTBGeometryESProducer_cfi')
0018 process.load('Configuration.StandardSequences.MagneticField_0T_cff')
0019 process.load('Configuration.StandardSequences.Generator_cff')
0020 process.load('GeneratorInterface.Core.generatorSmeared_cfi')
0021 process.load('IOMC.EventVertexGenerators.VtxSmearedGauss_cfi')
0022 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0023 process.load('Configuration.StandardSequences.SimIdeal_cff')
0024 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0025 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0026 process.load('SimG4CMS.HGCalTestBeam.DigiHGCalTB160_cff')
0027 process.load('RecoLocalCalo.HGCalRecProducers.HGCalLocalRecoTestBeamSequence_cff')
0028 process.load('SimG4CMS.HGCalTestBeam.HGCalTBAnalyzer_cfi')
0029
0030 process.maxEvents = cms.untracked.PSet(
0031     input = cms.untracked.int32(1000)
0032 )
0033
0034 # Input source
0035 process.source = cms.Source("EmptySource")
0036
0037 process.options = cms.untracked.PSet(
0038
0039 )
0040
0041 # Production Info
0042 process.configurationMetadata = cms.untracked.PSet(
0043     version = cms.untracked.string('$Revision: 1.20 $'),
0044     annotation = cms.untracked.string('SingleElectronPt10_cfi nevts:10'),
0045     name = cms.untracked.string('Applications')
0046 )
0047
0048 # Output definition
0049
0050 process.FEVTDEBUGHLToutput = cms.OutputModule("PoolOutputModule",
0051     splitLevel = cms.untracked.int32(0),
0052     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0053 #    outputCommands = process.FEVTDEBUGHLTEventContent.outputCommands,
0054     outputCommands = cms.untracked.vstring('keep *'),
0055     fileName = cms.untracked.string('file:gensimdigireco.root'),
0056     dataset = cms.untracked.PSet(
0057         filterName = cms.untracked.string(''),
0058         dataTier = cms.untracked.string('GEN-SIM-DIGI-RECO')
0059     ),
0060     SelectEvents = cms.untracked.PSet(
0061         SelectEvents = cms.vstring('generation_step')
0062     )
0063 )
0064
0065 # Additional output definition
0066 process.TFileService = cms.Service("TFileService",
0067                                    fileName = cms.string('TBProt4MGenSimDigiReco.root')
0068                                    )
0069
0070 # Other statements
0071 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0072 from Configuration.AlCa.GlobalTag import GlobalTag
0073 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2017_design', '')
0074
0075 process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
0076     AddAntiParticle = cms.bool(True),
0077     PGunParameters = cms.PSet(
0078         MinE = cms.double(120.0),
0079         MaxE = cms.double(120.0),
0080         MinTheta = cms.double(0.0),
0081         MaxTheta = cms.double(0.0),
0082         MinPhi = cms.double(-3.14159265359),
0083         MaxPhi = cms.double(3.14159265359),
0084         PartID = cms.vint32(2212)
0085     ),
0086     Verbosity = cms.untracked.int32(0),
0087     firstRun = cms.untracked.uint32(1),
0088     psethack = cms.string('single muon E 100')
0089 )
0090 process.VtxSmeared.MeanX  = 0.0
0091 process.VtxSmeared.SigmaX = 0.55
0092 process.VtxSmeared.MeanY  = 0.0
0093 process.VtxSmeared.SigmaY = 0.65
0094 process.VtxSmeared.MeanZ  = -3500.0
0095 process.VtxSmeared.SigmaZ = 0
0096 process.HGCalUncalibRecHit.HGCHEFConfig.isSiFE = False
0097 process.g4SimHits.OnlySDs = ['AHcalSensitiveDetector',
0098                              'HGCSensitiveDetector',
0099                              'HGCalTB1601SensitiveDetector',
0100                              'HcalTB06BeamDetector']
0101
0102 # Path and EndPath definitions
0103 process.generation_step = cms.Path(process.pgen)
0104 process.simulation_step = cms.Path(process.psim)
0105 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0106 process.digitisation_step = cms.Path(process.mix)
0107 process.reconstruction_step = cms.Path(process.HGCalLocalRecoTestBeamSequence)
0108 process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
0109 process.endjob_step = cms.EndPath(process.endOfProcess)
0110 process.FEVTDEBUGHLToutput_step = cms.EndPath(process.FEVTDEBUGHLToutput)
0111
0112 # Schedule definition
0113 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.digitisation_step,process.reconstruction_step,process.analysis_step,process.endjob_step,process.FEVTDEBUGHLToutput_step)
0114 # filter all path with the production filter sequence
0115 for path in process.paths:
0116         getattr(process,path)._seq = process.generator * getattr(process,path)._seq
0117
0118 for label, prod in process.producers_().items():
0119         if prod.type_() == "OscarMTProducer":
0120             # ugly hack
0121             prod.__dict__['_TypedParameterizable__type'] = "OscarProducer"