HcalHits/test/run_HFV2_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("CaloTest")
0004 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0005
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007 process.load('Configuration.StandardSequences.Generator_cff')
0008
0009 #Magnetic Field
0010 process.load("Configuration.StandardSequences.MagneticField_cff")
0011
0012 #Geometry
0013 process.load("Validation.HcalHits.testGeometryPMTXML_cfi")
0014 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0015 process.load("Geometry.HcalCommonData.hcalDDConstants_cff")
0016
0017 # Calo geometry service model
0018 process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
0019
0020 process.load("Geometry.CaloEventSetup.CaloTopology_cfi")
0021
0022 process.load("SimG4Core.Application.g4SimHits_cfi")
0023
0024 process.load("Validation.HcalHits.HcalHitValidation_cfi")
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         CaloSim = cms.untracked.PSet(
0032             limit = cms.untracked.int32(0)
0033         ),
0034         G4cerr = cms.untracked.PSet(
0035             limit = cms.untracked.int32(0)
0036         ),
0037         G4cout = cms.untracked.PSet(
0038             limit = cms.untracked.int32(0)
0039         ),
0040         HFShower = cms.untracked.PSet(
0041             limit = cms.untracked.int32(0)
0042         ),
0043         HcalSim = cms.untracked.PSet(
0044             limit = cms.untracked.int32(0)
0045         ),
0046         ValidHcal = cms.untracked.PSet(
0047             limit = cms.untracked.int32(0)
0048         ),
0049         default = cms.untracked.PSet(
0050             limit = cms.untracked.int32(0)
0051         ),
0052         enable = cms.untracked.bool(True)
0053     )
0054 )
0055
0056 process.Timing = cms.Service("Timing")
0057
0058 process.load("IOMC.RandomEngine.IOMC_cff")
0059 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0060 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0061 process.rndmStore = cms.EDProducer("RandomEngineStateProducer")
0062
0063 process.maxEvents = cms.untracked.PSet(
0064     input = cms.untracked.int32(100)
0065 )
0066
0067 process.source = cms.Source("EmptySource")
0068
0069 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0070     PGunParameters = cms.PSet(
0071         PartID = cms.vint32(11),
0072         MinEta = cms.double(3.070),
0073         MaxEta = cms.double(3.071),
0074         MinPhi = cms.double(0.0872),
0075         MaxPhi = cms.double(0.0873),
0076         MinE   = cms.double(99.90),
0077         MaxE   = cms.double(100.1)
0078     ),
0079     AddAntiParticle = cms.bool(False),
0080     Verbosity       = cms.untracked.int32(0),
0081     firstRun        = cms.untracked.uint32(1)
0082 )
0083
0084 process.USER = cms.OutputModule("PoolOutputModule",
0085     fileName = cms.untracked.string('100_pi_bs.root')
0086 )
0087
0088 process.DQMStore = cms.Service("DQMStore",
0089     verbose = cms.untracked.int32(0)
0090 )
0091
0092 process.DQM = cms.Service("DQM",
0093     debug = cms.untracked.bool(False),
0094     publishFrequency = cms.untracked.double(5.0),
0095     collectorPort = cms.untracked.int32(9090),
0096     collectorHost = cms.untracked.string('')
0097 )
0098
0099 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.hcalHitValid)
0100 process.outpath = cms.EndPath(process.USER)
0101 process.VtxSmeared.SigmaX = 0.00001
0102 process.VtxSmeared.SigmaY = 0.00001
0103 process.VtxSmeared.SigmaZ = 0.00001
0104 process.g4SimHits.UseMagneticField = False
0105 process.g4SimHits.OnlySDs = ['EcalSensitiveDetector', 'CaloTrkProcessing', 'HcalSensitiveDetector']
0106 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_BERT_EML'
0107 process.g4SimHits.HCalSD.UseParametrize = True
0108 process.g4SimHits.HCalSD.UsePMTHits = True
0109 process.g4SimHits.HCalSD.BetaThreshold = 0.70
0110 process.g4SimHits.HFShower.TrackEM = True
0111 process.g4SimHits.G4Commands = ['/physics_engine/neutron/energyLimit 0 keV', '/physics_engine/neutron/timeLimit 0.01 ms']
0112 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0113     SimG4HcalValidation = cms.PSet(
0114         TimeLowLimit = cms.double(0.0),
0115         LabelNxNInfo = cms.untracked.string('HcalInfoNxN'),
0116         LabelLayerInfo = cms.untracked.string('HcalInfoLayer'),
0117         HcalHitThreshold = cms.double(1e-20),
0118         Phi0 = cms.double(0.3054),
0119         ConeSize = cms.double(0.5),
0120         InfoLevel = cms.int32(2),
0121         JetThreshold = cms.double(5.0),
0122         EcalHitThreshold = cms.double(1e-20),
0123         TimeUpLimit = cms.double(999.0),
0124         HcalClusterOnly = cms.bool(False),
0125         Eta0 = cms.double(0.3045),
0126         LabelJetsInfo = cms.untracked.string('HcalInfoJets'),
0127         Names = cms.vstring('HcalHits',
0128             'EcalHitsEB',
0129             'EcalHitsEE',
0130             'EcalHitsES'),
0131         HcalSampling = cms.bool(True)
0132     ),
0133     type = cms.string('SimG4HcalValidation')
0134 ))
0135 process.hcalHitValid.outputFile = '100_pi_bs_plots.root'
0136