test/python/testHcalTBParameterDDD_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PROD")
0004
0005 process.load('FWCore.MessageService.MessageLogger_cfi')
0006 process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0007 process.load("SimG4CMS.HcalTestBeam.TB2004GeometryXML_cfi")
0008 #process.load("SimG4CMS.HcalTestBeam.TB2002GeometryXML_cfi")
0009 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0010 process.load("Geometry.HcalTestBeamData.hcalDDDSimConstants_cff")
0011 process.load('GeneratorInterface.Core.generatorSmeared_cfi')
0012 process.load("Configuration.EventContent.EventContent_cff")
0013
0014 process.maxEvents = cms.untracked.PSet(
0015     input = cms.untracked.int32(1)
0016 )
0017
0018 process.MessageLogger.cerr.FwkReport.reportEvery = 5
0019 if hasattr(process,'MessageLogger'):
0020     process.MessageLogger.HCalGeom=dict()
0021     process.MessageLogger.EcalGeom=dict()
0022     process.MessageLogger.HcalSim=dict()
0023     process.MessageLogger.HcalTBSim=dict()
0024     process.MessageLogger.EcalSim=dict()
0025     process.MessageLogger.CaloSim=dict()
0026
0027 process.load("IOMC.RandomEngine.IOMC_cff")
0028 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0029 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0030 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0031
0032 process.common_beam_direction_parameters = cms.PSet(
0033     MinEta       = cms.double(0.5655),
0034     MaxEta       = cms.double(0.5655),
0035     MinPhi       = cms.double(-0.1309),
0036     MaxPhi       = cms.double(-0.1309),
0037     BeamPosition = cms.double(-521.5)
0038 )
0039
0040 from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import *
0041 process.VtxSmeared = cms.EDProducer("BeamProfileVtxGenerator",
0042     process.common_beam_direction_parameters,
0043     VtxSmearedCommon,
0044     BeamMeanX       = cms.double(0.0),
0045     BeamMeanY       = cms.double(0.0),
0046     BeamSigmaX      = cms.double(0.0001),
0047     BeamSigmaY      = cms.double(0.0001),
0048     Psi             = cms.double(999.9),
0049     GaussianProfile = cms.bool(False),
0050     BinX            = cms.int32(50),
0051     BinY            = cms.int32(50),
0052     File            = cms.string('beam.profile'),
0053     UseFile         = cms.bool(False),
0054     TimeOffset      = cms.double(0.)
0055 )
0056
0057 process.source = cms.Source("EmptySource")
0058
0059 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0060     PGunParameters = cms.PSet(
0061         process.common_beam_direction_parameters,
0062         MinE   = cms.double(9.99),
0063         MaxE   = cms.double(10.01),
0064         PartID = cms.vint32(211)
0065     ),
0066     Verbosity       = cms.untracked.int32(0),
0067     AddAntiParticle = cms.bool(False),
0068     firstRun        = cms.untracked.uint32(1)
0069 )
0070
0071 process.Timing = cms.Service("Timing")
0072
0073 process.testH4 = cms.EDAnalyzer("HcalTBParameterTester",
0074                                 Name = cms.untracked.string("EcalHitsEB"),
0075                                 Mode = cms.untracked.int32(1)
0076 )
0077
0078 process.testH2EE = cms.EDAnalyzer("HcalTBParameterTester",
0079                                   Name = cms.untracked.string("EcalHitsEB"),
0080                                   Mode = cms.untracked.int32(0)
0081 )
0082
0083 process.testH2HC = cms.EDAnalyzer("HcalTBParameterTester",
0084                                   Name = cms.untracked.string("HcalHits"),
0085                                   Mode = cms.untracked.int32(0)
0086 )
0087
0088 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.testH4)
0089 #process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.testH2EE*process.testH2HC)