HGCalTestBeam/test/HGCalTBCERN181Oct1_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Eras.Modifier_hgcaltb_cff import hgcaltb
0003
0004 process = cms.Process('SIM', hgcaltb)
0005
0006 # import of standard configurations
0007 process.load("FWCore.MessageService.MessageLogger_cfi")
0008 process.load('Configuration.StandardSequences.Services_cff')
0009 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0010 process.load('Configuration.EventContent.EventContent_cff')
0011 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0012 process.load('SimG4CMS.HGCalTestBeam.HGCalTB181Oct1XML_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.HcalCommonData.caloSimulationParameters_cff')
0017 process.load('Geometry.CaloEventSetup.HGCalTBTopology_cfi')
0018 process.load('Geometry.HGCalGeometry.HGCalTBGeometryESProducer_cfi')
0019 process.load('Configuration.StandardSequences.MagneticField_0T_cff')
0020 process.load('Configuration.StandardSequences.Generator_cff')
0021 process.load('IOMC.EventVertexGenerators.VtxSmearedFlat_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.HGCalTBAnalyzer_cfi')
0027 process.load('SimG4CMS.HGCalTestBeam.HGCalTBCheckGunPosition_cfi')
0028
0029 process.maxEvents = cms.untracked.PSet(
0030     input = cms.untracked.int32(10)
0031 )
0032
0033 if 'MessageLogger' in process.__dict__:
0034     process.MessageLogger.HGCalGeom=dict()
0035     process.MessageLogger.HGCSim=dict()
0036     process.MessageLogger.HcalSim=dict()
0037     process.MessageLogger.HcalTB06BeamSD=dict()
0038 #    process.MessageLogger.SimG4CoreGeometry=dict()
0039
0040 # Input source
0041 process.source = cms.Source("EmptySource")
0042
0043 process.options = cms.untracked.PSet(
0044 )
0045
0046 # Production Info
0047 process.configurationMetadata = cms.untracked.PSet(
0048     annotation = cms.untracked.string('SingleMuonE200_cfi nevts:10'),
0049     name = cms.untracked.string('Applications'),
0050     version = cms.untracked.string('$Revision: 1.19 $')
0051 )
0052
0053 # Output definition
0054
0055 process.FEVTDEBUGoutput = cms.OutputModule("PoolOutputModule",
0056     SelectEvents = cms.untracked.PSet(
0057         SelectEvents = cms.vstring('generation_step')
0058     ),
0059     dataset = cms.untracked.PSet(
0060         dataTier = cms.untracked.string('GEN-SIM'),
0061         filterName = cms.untracked.string('')
0062     ),
0063     fileName = cms.untracked.string('file:TBGenSim181Oct1.root'),
0064     outputCommands = process.FEVTDEBUGEventContent.outputCommands,
0065     splitLevel = cms.untracked.int32(0)
0066 )
0067
0068 # Additional output definition
0069 process.TFileService = cms.Service("TFileService",
0070                                    fileName = cms.string('TBGenSim.root')
0071                                    )
0072
0073 # Other statements
0074 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0075 from Configuration.AlCa.GlobalTag import GlobalTag
0076 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0077
0078 process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
0079     AddAntiParticle = cms.bool(False),
0080     PGunParameters = cms.PSet(
0081         MinE = cms.double(99.99),
0082         MaxE = cms.double(100.01),
0083         MinTheta = cms.double(0.0),
0084         MaxTheta = cms.double(0.0),
0085         MinPhi = cms.double(-3.14159265359),
0086         MaxPhi = cms.double(3.14159265359),
0087         PartID = cms.vint32(11)
0088     ),
0089     Verbosity = cms.untracked.int32(0),
0090     firstRun = cms.untracked.uint32(1),
0091     psethack = cms.string('single muon E 100')
0092 )
0093 process.VtxSmeared.MinZ                 = -800.0
0094 process.VtxSmeared.MaxZ                 = -800.0
0095 process.VtxSmeared.MinX                 = -7.5
0096 process.VtxSmeared.MaxX                 =  7.5
0097 process.VtxSmeared.MinY                 = -7.5
0098 process.VtxSmeared.MaxY                 =  7.5
0099 process.g4SimHits.HGCSD.RejectMouseBite = True
0100 process.g4SimHits.HGCSD.RotatedWafer    = True
0101 process.g4SimHits.G4CheckOverlap.OutputBaseName = "2018"
0102 process.g4SimHits.G4CheckOverlap.gdmlFlag       = True
0103 process.g4SimHits.FileNameGDML                  = "TBHGCal181Oct.gdml"
0104 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0105         HGCalTBPassive = cms.PSet(
0106                     LVNames = cms.vstring('HGCalEE','HGCalHE','HGCalAH', 'HGCalBeam', 'CMSE'),
0107                     MotherName = cms.string('OCMS'),
0108                     IfDD4hep = cms.bool(False),
0109                 ),
0110         type = cms.string('HGCalTBPassive'),
0111         )
0112                        )
0113 process.HGCalTBAnalyzer.doDigis         = False
0114 process.HGCalTBAnalyzer.doRecHits       = False
0115 process.HGCalTBAnalyzer.useFH           = True
0116 process.HGCalTBAnalyzer.useBH           = True
0117 process.HGCalTBAnalyzer.useBeam         = True
0118 process.HGCalTBAnalyzer.zFrontEE        = 1110.0
0119 process.HGCalTBAnalyzer.zFrontFH        = 1176.5
0120 process.HGCalTBAnalyzer.zFrontFH        = 1307.5
0121 process.HGCalTBAnalyzer.maxDepth        = 39
0122 process.HGCalTBAnalyzer.deltaZ          = 26.2
0123 process.HGCalTBAnalyzer.zFirst          = 22.8
0124 process.HGCalTBAnalyzer.doPassive       = True
0125 process.HGCalTBAnalyzer.doPassiveEE     = True
0126 process.HGCalTBAnalyzer.doPassiveHE     = True
0127 process.HGCalTBAnalyzer.doPassiveBH     = True
0128
0129 # Path and EndPath definitions
0130 process.generation_step = cms.Path(process.pgen)
0131 process.gunfilter_step  = cms.Path(process.HGCalTBCheckGunPostion)
0132 process.simulation_step = cms.Path(process.psim)
0133 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0134 process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
0135 process.endjob_step = cms.EndPath(process.endOfProcess)
0136 process.FEVTDEBUGoutput_step = cms.EndPath(process.FEVTDEBUGoutput)
0137
0138 # Schedule definition
0139 process.schedule = cms.Schedule(process.generation_step,
0140                 process.genfiltersummary_step,
0141                 process.simulation_step,
0142                 process.gunfilter_step,
0143                     process.analysis_step,
0144                 process.endjob_step,
0145                                 process.FEVTDEBUGoutput_step
0146                 )
0147 # filter all path with the production filter sequence
0148 for path in process.paths:
0149     getattr(process,path)._seq = process.generator * getattr(process,path)._seq
0150
0151