test/python/runWithGun_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PROD")
0004 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0005
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007 process.load('Configuration/StandardSequences/Generator_cff')
0008
0009 process.load("Geometry.CMSCommonData.cmsExtendedGeometryXML_cfi")
0010 process.load("Geometry.TrackerNumberingBuilder.trackerNumberingGeometry_cfi")
0011 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0012 process.load("Geometry.HcalCommonData.hcalDDConstants_cff")
0013 process.load("Geometry.MuonNumbering.muonGeometryConstants_cff")
0014 process.load("Geometry.MuonNumbering.muonOffsetESProducer_cff")
0015
0016 process.load("Configuration.StandardSequences.MagneticField_cff")
0017
0018 process.load("Configuration.EventContent.EventContent_cff")
0019
0020 process.load("SimG4Core.Application.g4SimHits_cfi")
0021
0022 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0023 from Configuration.AlCa.autoCond import autoCond
0024 process.GlobalTag.globaltag = autoCond['run1_mc']
0025
0026 process.MessageLogger = cms.Service("MessageLogger",
0027     cerr = cms.untracked.PSet(
0028         enable = cms.untracked.bool(False)
0029     ),
0030     cout = cms.untracked.PSet(
0031         DEBUG = cms.untracked.PSet(
0032             limit = cms.untracked.int32(0)
0033         ),
0034         G4cerr = cms.untracked.PSet(
0035             limit = cms.untracked.int32(-1)
0036         ),
0037         G4cout = cms.untracked.PSet(
0038             limit = cms.untracked.int32(-1)
0039         ),
0040         INFO = cms.untracked.PSet(
0041             limit = cms.untracked.int32(-1)
0042         ),
0043         SimTrackManager = cms.untracked.PSet(
0044             limit = cms.untracked.int32(-1)
0045         ),
0046         enable = cms.untracked.bool(True)
0047     )
0048 )
0049
0050 process.load("IOMC.RandomEngine.IOMC_cff")
0051 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0052 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0053 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0054
0055 process.maxEvents = cms.untracked.PSet(
0056     input = cms.untracked.int32(2)
0057 )
0058
0059 process.source = cms.Source("EmptySource",
0060     firstRun        = cms.untracked.uint32(1),
0061     firstEvent      = cms.untracked.uint32(1)
0062 )
0063
0064 process.generator = cms.EDProducer("FlatRandomPtGunProducer",
0065     PGunParameters = cms.PSet(
0066         PartID = cms.vint32(211),
0067         MinEta = cms.double(-6.5),
0068         MaxEta = cms.double(-5.4),
0069         MinPhi = cms.double(-3.14159265359),
0070         MaxPhi = cms.double(3.14159265359),
0071         MinPt  = cms.double(20.),
0072         MaxPt  = cms.double(20.)
0073     ),
0074     Verbosity       = cms.untracked.int32(0),
0075     AddAntiParticle = cms.bool(False)
0076 )
0077
0078 process.o1 = cms.OutputModule("PoolOutputModule",
0079     process.FEVTSIMEventContent,
0080     fileName = cms.untracked.string('simevent_QGSP_BERT_EML.root')
0081 )
0082
0083 process.Timing = cms.Service("Timing")
0084
0085 process.SimpleMemoryCheck = cms.Service("SimpleMemoryCheck",
0086     oncePerEventMode = cms.untracked.bool(True),
0087     showMallocInfo = cms.untracked.bool(True),
0088     dump = cms.untracked.bool(True),
0089     ignoreTotal = cms.untracked.int32(1)
0090 )
0091
0092 process.common_maximum_timex = cms.PSet(
0093     MaxTrackTime  = cms.double(1000.0),
0094     MaxTrackTimeForward = cms.double(2000.0), # ns
0095     MaxTimeNames  = cms.vstring(),
0096     MaxTrackTimes = cms.vdouble(),
0097     MaxZCentralCMS = cms.double(50.0), # m
0098     DeadRegions   = cms.vstring('QuadRegion','InterimRegion'),
0099     CriticalEnergyForVacuum = cms.double(2.0),   # MeV
0100     CriticalDensity         = cms.double(1e-15)  # g/cm3
0101 )
0102 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits)
0103 process.outpath = cms.EndPath(process.o1)
0104 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_BERT_EML'
0105 process.g4SimHits.ECalSD.StoreSecondary = True
0106 process.g4SimHits.CaloTrkProcessing.PutHistory = True
0107 process.g4SimHits.StackingAction = cms.PSet(
0108    process.common_heavy_suppression,
0109    process.common_maximum_timex,        # need to be localy redefined
0110    KillDeltaRay  = cms.bool(False),
0111    TrackNeutrino = cms.bool(False),
0112    KillHeavy     = cms.bool(False),
0113    KillGamma     = cms.bool(True),
0114    GammaThreshold = cms.double(0.0001), ## (MeV)
0115    SaveFirstLevelSecondary = cms.untracked.bool(True),
0116    SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0117    SavePrimaryDecayProductsAndConversionsInCalo    = cms.untracked.bool(True),
0118    SavePrimaryDecayProductsAndConversionsInMuon    = cms.untracked.bool(True),
0119    RusRoGammaEnergyLimit  = cms.double(5.0), ## (MeV)
0120    RusRoEcalGamma         = cms.double(0.3),
0121    RusRoHcalGamma         = cms.double(0.3),
0122    RusRoMuonIronGamma     = cms.double(0.3),
0123    RusRoPreShowerGamma    = cms.double(0.3),
0124    RusRoCastorGamma       = cms.double(0.3),
0125    RusRoWorldGamma        = cms.double(0.3),
0126    RusRoNeutronEnergyLimit  = cms.double(10.0), ## (MeV)
0127    RusRoEcalNeutron         = cms.double(0.1),
0128    RusRoHcalNeutron         = cms.double(0.1),
0129    RusRoMuonIronNeutron     = cms.double(0.1),
0130    RusRoPreShowerNeutron    = cms.double(0.1),
0131    RusRoCastorNeutron       = cms.double(0.1),
0132    RusRoWorldNeutron        = cms.double(0.1),
0133    RusRoProtonEnergyLimit  = cms.double(0.0),
0134    RusRoEcalProton         = cms.double(1.0),
0135    RusRoHcalProton         = cms.double(1.0),
0136    RusRoMuonIronProton     = cms.double(1.0),
0137    RusRoPreShowerProton    = cms.double(1.0),
0138    RusRoCastorProton       = cms.double(1.0),
0139    RusRoWorldProton        = cms.double(1.0)
0140 )
0141 process.g4SimHits.SteppingAction = cms.PSet(
0142    process.common_maximum_timex, # need to be localy redefined
0143    MaxNumberOfSteps        = cms.int32(50000),
0144    EkinNames               = cms.vstring(),
0145    EkinThresholds          = cms.vdouble(),
0146    EkinParticles           = cms.vstring(),
0147    Verbosity               = cms.untracked.int32(1)
0148 )
0149 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0150     CheckForHighEtPhotons = cms.untracked.bool(False),
0151     EventMin     = cms.untracked.int32(0),
0152     EventMax     = cms.untracked.int32(0),
0153     EventStep    = cms.untracked.int32(1),
0154     TrackMin     = cms.untracked.int32(0),
0155     TrackMax     = cms.untracked.int32(999999999),
0156     TrackStep    = cms.untracked.int32(1),
0157     PDGids       = cms.untracked.vint32(),
0158     VerboseLevel = cms.untracked.int32(2),
0159     G4Verbose    = cms.untracked.bool(True),
0160     DEBUG        = cms.untracked.bool(False),
0161     type         = cms.string('TrackingVerboseAction')
0162 ))
0163