test/python/runHF_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.cmsHFPMTFibreXML_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_38T_cff")
0014 process.load("Configuration.EventContent.EventContent_cff")
0015 process.load("SimG4Core.Application.g4SimHits_cfi")
0016 process.load("SimG4CMS.Calo.HFPMTHitAnalyzer_cfi")
0017
0018 if 'MessageLogger' in process.__dict__:
0019     process.MessageLogger.G4cerr=dict()
0020     process.MessageLogger.G4cout=dict()
0021     process.MessageLogger.HFShower=dict()
0022
0023 process.load("IOMC.RandomEngine.IOMC_cff")
0024 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0025 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0026 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0027
0028 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0029 from Configuration.AlCa.GlobalTag import GlobalTag
0030 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run1_mc', '')
0031
0032 process.Timing = cms.Service("Timing")
0033
0034 process.maxEvents = cms.untracked.PSet(
0035     input = cms.untracked.int32(200)
0036 )
0037
0038 process.source = cms.Source("EmptySource",
0039     firstRun        = cms.untracked.uint32(1),
0040     firstEvent      = cms.untracked.uint32(1)
0041 )
0042
0043 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0044     PGunParameters = cms.PSet(
0045         PartID = cms.vint32(13),
0046         MinEta = cms.double(3.25),
0047         MaxEta = cms.double(3.50),
0048         MinPhi = cms.double(-3.1415926),
0049         MaxPhi = cms.double(3.1415926),
0050         MinE   = cms.double(1000.00),
0051         MaxE   = cms.double(1000.00)
0052     ),
0053     Verbosity       = cms.untracked.int32(0),
0054     AddAntiParticle = cms.bool(False)
0055 )
0056
0057 process.o1 = cms.OutputModule("PoolOutputModule",
0058     process.FEVTSIMEventContent,
0059     fileName = cms.untracked.string('simevent.root')
0060 )
0061
0062 process.TFileService = cms.Service("TFileService",
0063     fileName = cms.string('HFPMTFibre1.root')
0064 )
0065
0066 process.common_maximum_timex = cms.PSet(
0067     MaxTrackTime            = cms.double(500.0),
0068     MaxTrackTimeForward     = cms.double(2000.0), # ns
0069     MaxTimeNames            = cms.vstring(),
0070     MaxTrackTimes           = cms.vdouble(),
0071     MaxZCentralCMS          = cms.double(50.0), # m
0072     DeadRegions             = cms.vstring(),
0073     CriticalEnergyForVacuum = cms.double(2.0),
0074     CriticalDensity         = cms.double(1e-15)
0075 )
0076
0077 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.g4SimHits*process.HFPMTHitAnalyzer)
0078 process.outpath = cms.EndPath(process.o1)
0079 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_BERT_EML'
0080 process.g4SimHits.Physics.DefaultCutValue   = 0.1
0081 process.g4SimHits.HCalSD.UseShowerLibrary   = False
0082 process.g4SimHits.HCalSD.UseParametrize     = True
0083 process.g4SimHits.HCalSD.UsePMTHits         = True
0084 process.g4SimHits.HCalSD.UseFibreBundleHits = True
0085 process.g4SimHits.HFShower.UseShowerLibrary = False
0086 process.g4SimHits.HFShower.UseHFGflash      = True
0087 process.g4SimHits.HFShower.TrackEM          = False
0088 process.g4SimHits.HFShower.HFShowerBlock.OnlyLong = cms.bool(True)
0089 process.g4SimHits.HFShower.EminLibrary      = 0.0
0090 process.g4SimHits.HFShower.ApplyFiducialCut = True
0091 process.g4SimHits.StackingAction = cms.PSet(
0092     process.common_heavy_suppression,
0093     process.common_maximum_timex,
0094     KillDeltaRay  = cms.bool(False),
0095     TrackNeutrino = cms.bool(False),
0096     KillHeavy     = cms.bool(False),
0097     KillGamma     = cms.bool(True),
0098     GammaThreshold = cms.double(0.0001), ## (MeV)
0099     SaveFirstLevelSecondary = cms.untracked.bool(True),
0100     SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0101     SavePrimaryDecayProductsAndConversionsInCalo    = cms.untracked.bool(True),
0102     SavePrimaryDecayProductsAndConversionsInMuon    = cms.untracked.bool(True),
0103     RusRoGammaEnergyLimit  = cms.double(5.0), ## (MeV)
0104     RusRoEcalGamma         = cms.double(0.3),
0105     RusRoHcalGamma         = cms.double(0.3),
0106     RusRoMuonIronGamma     = cms.double(0.3),
0107     RusRoPreShowerGamma    = cms.double(0.3),
0108     RusRoCastorGamma       = cms.double(0.3),
0109     RusRoWorldGamma        = cms.double(0.3),
0110     RusRoNeutronEnergyLimit  = cms.double(10.0), ## (MeV)
0111     RusRoEcalNeutron         = cms.double(0.1),
0112     RusRoHcalNeutron         = cms.double(0.1),
0113     RusRoMuonIronNeutron     = cms.double(0.1),
0114     RusRoPreShowerNeutron    = cms.double(0.1),
0115     RusRoCastorNeutron       = cms.double(0.1),
0116     RusRoWorldNeutron        = cms.double(0.1),
0117     RusRoProtonEnergyLimit  = cms.double(0.0),
0118     RusRoEcalProton         = cms.double(1.0),
0119     RusRoHcalProton         = cms.double(1.0),
0120     RusRoMuonIronProton     = cms.double(1.0),
0121     RusRoPreShowerProton    = cms.double(1.0),
0122     RusRoCastorProton       = cms.double(1.0),
0123     RusRoWorldProton        = cms.double(1.0)
0124 )
0125 process.g4SimHits.SteppingAction = cms.PSet(
0126     process.common_maximum_timex,
0127     MaxNumberOfSteps        = cms.int32(50000),
0128     EkinNames               = cms.vstring(),
0129     EkinThresholds          = cms.vdouble(),
0130     EkinParticles           = cms.vstring(),
0131     Verbosity               = cms.untracked.int32(2)
0132 )
0133 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0134     CheckForHighEtPhotons = cms.untracked.bool(False),
0135     TrackMin     = cms.untracked.int32(0),
0136     TrackMax     = cms.untracked.int32(99999999),
0137     TrackStep    = cms.untracked.int32(1),
0138     EventMin     = cms.untracked.int32(0),
0139     EventMax     = cms.untracked.int32(0),
0140     EventStep    = cms.untracked.int32(1),
0141     PDGids       = cms.untracked.vint32(),
0142     VerboseLevel = cms.untracked.int32(0),
0143     G4Verbose    = cms.untracked.bool(True),
0144     DEBUG        = cms.untracked.bool(False),
0145     type         = cms.string('TrackingVerboseAction')
0146 ))
0147