test/python/runHF1_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PROD")
0004 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0005 process.load('FWCore.MessageService.MessageLogger_cfi')
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007 process.load("Geometry.CMSCommonData.cmsExtendedGeometryHFParametrizeXML_cfi")
0008 process.load("Geometry.TrackerNumberingBuilder.trackerNumberingGeometry_cff")
0009 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0010 process.load("Geometry.HcalCommonData.hcalDDConstants_cff")
0011 process.load("Geometry.MuonNumbering.muonGeometryConstants_cff")
0012 process.load("Geometry.MuonNumbering.muonOffsetESProducer_cff")
0013 process.load("Configuration.StandardSequences.MagneticField_cff")
0014 process.load("Configuration.EventContent.EventContent_cff")
0015 process.load('Configuration.StandardSequences.Generator_cff')
0016 process.load('Configuration.StandardSequences.SimIdeal_cff')
0017 process.load("SimG4CMS.Calo.HFPMTHitAnalyzer_cfi")
0018
0019 if 'MessageLogger' in process.__dict__:
0020     process.MessageLogger.G4cerr=dict()
0021     process.MessageLogger.G4cout=dict()
0022     process.MessageLogger.HFShower=dict()
0023
0024 process.load("IOMC.RandomEngine.IOMC_cff")
0025 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0026 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0027 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0028
0029 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0030 from Configuration.AlCa.autoCond import autoCond
0031 process.GlobalTag.globaltag = autoCond['run1_mc']
0032
0033 process.Timing = cms.Service("Timing")
0034
0035 process.maxEvents = cms.untracked.PSet(
0036     input = cms.untracked.int32(200)
0037 )
0038
0039 process.source = cms.Source("EmptySource",
0040     firstRun        = cms.untracked.uint32(1),
0041     firstEvent      = cms.untracked.uint32(1)
0042 )
0043
0044 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0045     PGunParameters = cms.PSet(
0046         PartID = cms.vint32(211),
0047         MinEta = cms.double(3.25),
0048         MaxEta = cms.double(3.50),
0049         MinPhi = cms.double(-3.1415926),
0050         MaxPhi = cms.double(3.1415926),
0051         MinE   = cms.double(1000.00),
0052         MaxE   = cms.double(1000.00)
0053     ),
0054     Verbosity       = cms.untracked.int32(0),
0055     AddAntiParticle = cms.bool(False)
0056 )
0057
0058 process.output = cms.OutputModule("PoolOutputModule",
0059     process.FEVTSIMEventContent,
0060     fileName = cms.untracked.string('simevent.root')
0061 )
0062
0063 process.TFileService = cms.Service("TFileService",
0064     fileName = cms.string('HFParametrizeGFlash.root')
0065 )
0066
0067 process.common_maximum_timex = cms.PSet(
0068     MaxTrackTime  = cms.double(500.0),
0069     MaxTimeNames  = cms.vstring(),
0070     MaxTrackTimes = cms.vdouble()
0071 )
0072
0073 process.generation_step = cms.Path(process.pgen)
0074 process.simulation_step = cms.Path(process.psim)
0075 process.analysis_step   = cms.Path(process.HFPMTHitAnalyzer)
0076 process.out_step = cms.EndPath(process.output)
0077
0078 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_FTFP_BERT_EML'
0079 process.g4SimHits.Physics.DefaultCutValue   = 0.1
0080 process.g4SimHits.HCalSD.UseShowerLibrary   = False
0081 process.g4SimHits.HCalSD.UseParametrize     = True
0082 process.g4SimHits.HCalSD.UsePMTHits         = True
0083 process.g4SimHits.HCalSD.UseFibreBundleHits = True
0084 process.g4SimHits.HFShower.UseShowerLibrary = False
0085 process.g4SimHits.HFShower.UseHFGflash      = True
0086 process.g4SimHits.HFShower.TrackEM          = False
0087 process.g4SimHits.HFShower.EminLibrary      = 0.0
0088 process.g4SimHits.HFShower.HFShowerBlock.OnlyLong = cms.bool(True)
0089 process.g4SimHits.HFShower.Aperture         = 0.33
0090 process.g4SimHits.HFShower.ParametrizeLast  = True
0091 # Schedule definition
0092 process.schedule = cms.Schedule(process.generation_step,
0093                                 process.simulation_step,
0094                                 process.analysis_step,
0095                                 process.out_step
0096                                 )
0097
0098 # filter all path with the production filter sequence
0099 for path in process.paths:
0100         getattr(process,path)._seq = process.generator * getattr(process,path)._seq