test/python/runHGC4_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Eras.Era_Phase2C11_cff import Phase2C11
0003
0004 process = cms.Process("PROD",Phase2C11)
0005 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007 process.load("Configuration.Geometry.GeometryExtended2026D77Reco_cff")
0008 process.load("Configuration.StandardSequences.MagneticField_cff")
0009 process.load("Configuration.EventContent.EventContent_cff")
0010 process.load('Configuration.StandardSequences.Generator_cff')
0011 process.load('Configuration.StandardSequences.SimIdeal_cff')
0012 process.load('FWCore.MessageService.MessageLogger_cfi')
0013 process.load('SimG4CMS.Calo.hgcalHitIdCheck_cff')
0014 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0015 from Configuration.AlCa.GlobalTag import GlobalTag
0016 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T21', '')
0017
0018 if hasattr(process,'MessageLogger'):
0019     process.MessageLogger.HGCalGeom=dict()
0020     process.MessageLogger.HGCSim=dict()
0021     process.MessageLogger.CaloSim=dict()
0022     process.MessageLogger.HitStudy=dict()
0023 #   process.MessageLogger.SimG4CoreGeometry=dict()
0024
0025 process.load("IOMC.RandomEngine.IOMC_cff")
0026 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0027 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0028 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0029
0030 process.Timing = cms.Service("Timing")
0031
0032 process.maxEvents = cms.untracked.PSet(
0033     input = cms.untracked.int32(5)
0034 )
0035
0036 process.source = cms.Source("EmptySource",
0037     firstRun        = cms.untracked.uint32(1),
0038     firstEvent      = cms.untracked.uint32(1)
0039 )
0040
0041 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0042     PGunParameters = cms.PSet(
0043         PartID = cms.vint32(211),
0044 #       MinEta = cms.double(2.95),
0045 #       MaxEta = cms.double(3.01),
0046 #       MinEta = cms.double(1.69),
0047 #       MaxEta = cms.double(2.32),
0048         MinEta = cms.double(1.50),
0049         MaxEta = cms.double(2.00),
0050         MinPhi = cms.double(-3.1415926),
0051         MaxPhi = cms.double(-1.5707963),
0052         MinE   = cms.double(100.00),
0053         MaxE   = cms.double(100.00)
0054     ),
0055     Verbosity       = cms.untracked.int32(0),
0056     AddAntiParticle = cms.bool(True)
0057 )
0058
0059 process.output = cms.OutputModule("PoolOutputModule",
0060     process.FEVTSIMEventContent,
0061     fileName = cms.untracked.string('hgcV14.root')
0062 )
0063
0064 process.generation_step = cms.Path(process.pgen)
0065 process.simulation_step = cms.Path(process.psim)
0066 process.analysis_step = cms.Path(process.hgcalHitIdCheckEE+process.hgcalHitIdCheckHEF+process.hgcalHitIdCheckHEB)
0067 process.out_step = cms.EndPath(process.output)
0068
0069 process.g4SimHits.Physics.type = 'SimG4Core/Physics/FTFP_BERT_EMM'
0070 process.g4SimHits.Physics.DefaultCutValue          = 0.1
0071 process.g4SimHits.HGCSD.CornerMinMask              = 3
0072 process.g4SimHits.CaloSD.UseFineCalo               = [2]
0073 process.g4SimHits.CaloTrkProcessing.UseFineCalo    = [2]
0074 process.g4SimHits.TrackingAction.UseFineCalo       = [2]
0075 process.g4SimHits.CaloSD.DoFineCalo                = True
0076 process.g4SimHits.CaloTrkProcessing.DoFineCalo     = True
0077 process.g4SimHits.TrackingAction.DoFineCalo        = True
0078 process.g4SimHits.CaloSD.EminFineTrack             = 1000.0
0079 process.g4SimHits.CaloTrkProcessing.EminFineTrack  = 1000.0
0080 process.g4SimHits.TrackingAction.EminFineTrack     = 1000.0
0081 process.hgcalHitIdCheckHEB.Verbosity               = 1
0082
0083 # Schedule definition
0084 process.schedule = cms.Schedule(process.generation_step,
0085                                 process.simulation_step,
0086                                 process.analysis_step,
0087                                 process.out_step
0088                                 )
0089
0090 # filter all path with the production filter sequence
0091 for path in process.paths:
0092         getattr(process,path)._seq = process.generator * getattr(process,path)._seq