HGCalTestBeam/test/HGCalTBGenSimMB_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('SIM')
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('SimGeneral.MixingModule.mixNoPU_cfi')
0011 process.load('SimG4CMS.HGCalTestBeam.HGCalTB160Module16XML_cfi')
0012 process.load('Geometry.HGCalTBCommonData.hgcalTBNumberingInitialization_cfi')
0013 process.load('Geometry.HGCalTBCommonData.hgcalTBParametersInitialization_cfi')
0014 process.load('Configuration.StandardSequences.MagneticField_0T_cff')
0015 process.load('Configuration.StandardSequences.Generator_cff')
0016 process.load('IOMC.EventVertexGenerators.VtxSmearedGauss_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('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0021 process.load('SimG4CMS.HGCalTestBeam.HGCalTBAnalyzer_cfi')
0022 process.load('SimG4CMS.HGCalTestBeam.hgcalTBMBFNAL_cfi')
0023
0024 process.maxEvents = cms.untracked.PSet(
0025     input = cms.untracked.int32(100)
0026 )
0027
0028 if 'MessageLogger' in process.__dict__:
0029     process.MessageLogger.HGCalGeom=dict()
0030     process.MessageLogger.HGCSim=dict()
0031
0032 # Input source
0033 process.source = cms.Source("EmptySource")
0034
0035 process.options = cms.untracked.PSet(
0036 )
0037
0038 # Production Info
0039 process.configurationMetadata = cms.untracked.PSet(
0040     annotation = cms.untracked.string('SingleNeutrinoE100_cfi nevts:100'),
0041     name = cms.untracked.string('Applications'),
0042     version = cms.untracked.string('$Revision: 1.19 $')
0043 )
0044
0045 # Output definition
0046
0047 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0048     SelectEvents = cms.untracked.PSet(
0049         SelectEvents = cms.vstring('generation_step')
0050     ),
0051     dataset = cms.untracked.PSet(
0052         dataTier = cms.untracked.string('GEN-SIM'),
0053         filterName = cms.untracked.string('')
0054     ),
0055     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0056     fileName = cms.untracked.string('file:gensim.root'),
0057     outputCommands = process.RAWSIMEventContent.outputCommands,
0058     splitLevel = cms.untracked.int32(0)
0059 )
0060
0061 # Additional output definition
0062 process.TFileService = cms.Service("TFileService",
0063                                    fileName = cms.string('TBGenSim.root')
0064                                    )
0065
0066 # Other statements
0067 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0068 from Configuration.AlCa.GlobalTag import GlobalTag
0069 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0070
0071 process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
0072     AddAntiParticle = cms.bool(False),
0073     PGunParameters = cms.PSet(
0074         MinE = cms.double(99.99),
0075         MaxE = cms.double(100.01),
0076         MinTheta = cms.double(0.0),
0077         MaxTheta = cms.double(0.0),
0078         MinPhi = cms.double(-3.14159265359),
0079         MaxPhi = cms.double(3.14159265359),
0080         PartID = cms.vint32(14)
0081     ),
0082     Verbosity = cms.untracked.int32(0),
0083     firstRun = cms.untracked.uint32(1),
0084     psethack = cms.string('single neutrino E 100')
0085 )
0086 process.VtxSmeared.MeanZ = -3600.0
0087 process.VtxSmeared.SigmaZ = 0
0088 process.HGCalTBAnalyzer.doDigis = False
0089 process.HGCalTBAnalyzer.doRecHits = False
0090 process.g4SimHits.StackingAction.TrackNeutrino = True
0091 process.g4SimHits.OnlySDs = ['HGCSensitiveDetector',
0092                              'HGCalTB1601SensitiveDetector']
0093 process.g4SimHits.HGCSD.Detectors = 1
0094 process.g4SimHits.HGCSD.UseDetector = 1
0095 # Path and EndPath definitions
0096 process.generation_step = cms.Path(process.pgen)
0097 process.simulation_step = cms.Path(process.psim)
0098 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0099 process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
0100 process.endjob_step = cms.EndPath(process.endOfProcess)
0101 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0102
0103 # Schedule definition
0104 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.analysis_step,process.endjob_step,process.RAWSIMoutput_step)
0105 # filter all path with the production filter sequence
0106 for path in process.paths:
0107     getattr(process,path)._seq = process.generator * getattr(process,path)._seq
0108
0109