test/python/runHGC4_cfg.py

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