Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:14:50

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Eras.Modifier_h2tb_cff import h2tb
0003 
0004 process = cms.Process("PROD", h2tb)
0005 
0006 process.load('FWCore.MessageService.MessageLogger_cfi')
0007 process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0008 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0009 process.load("Geometry.HcalTestBeamData.hcalDDDSimConstants_cff")
0010 process.load("Geometry.HcalTestBeamData.TB2007TestGeometryXML_cfi")
0011 process.load('GeneratorInterface.Core.generatorSmeared_cfi')
0012 process.load("Configuration.EventContent.EventContent_cff")
0013 process.load("SimG4Core.Application.g4SimHits_cfi")
0014 
0015 process.MessageLogger.cerr.FwkReport.reportEvery = 5
0016 if hasattr(process,'MessageLogger'):
0017     process.MessageLogger.HCalGeom=dict()
0018     process.MessageLogger.EcalGeom=dict()
0019     process.MessageLogger.VertexGenerator=dict()
0020 
0021 process.Timing = cms.Service("Timing")
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 beamPosition = -800.0
0029 process.common_beam_direction_parameters = cms.PSet(
0030     MaxEta       = cms.double(0.5655),
0031     MinEta       = cms.double(0.5655),
0032     MaxPhi       = cms.double(-0.1309),
0033     MinPhi       = cms.double(-0.1309),
0034     BeamPosition = cms.double(beamPosition)
0035 )
0036 
0037 process.source = cms.Source("EmptySource")
0038 
0039 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0040     PGunParameters = cms.PSet(
0041         process.common_beam_direction_parameters,
0042         MinE   = cms.double(9.99),
0043         MaxE   = cms.double(10.01),
0044         PartID = cms.vint32(211)
0045     ),
0046     Verbosity       = cms.untracked.int32(0),
0047     AddAntiParticle = cms.bool(False),
0048     firstRun        = cms.untracked.uint32(1)
0049 )
0050 
0051 from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import *
0052 process.VtxSmeared = cms.EDProducer("BeamProfileVtxGenerator",
0053     process.common_beam_direction_parameters,
0054     VtxSmearedCommon,
0055     BeamMeanX       = cms.double(0.0),
0056     BeamMeanY       = cms.double(0.0),
0057     BeamSigmaX      = cms.double(0.0001),
0058     BeamSigmaY      = cms.double(0.0001),
0059     Psi             = cms.double(999.9),
0060     GaussianProfile = cms.bool(False),
0061     BinX            = cms.int32(50),
0062     BinY            = cms.int32(50),
0063     File            = cms.string('beam.profile'),
0064     UseFile         = cms.bool(False),
0065     TimeOffset      = cms.double(0.)
0066 )
0067 
0068 process.maxEvents = cms.untracked.PSet(
0069     input = cms.untracked.int32(1)
0070 )
0071 
0072 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits)
0073 
0074 process.g4SimHits.NonBeamEvent = True
0075 process.g4SimHits.UseMagneticField = False
0076 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_FTFP_BERT_EML'
0077 process.g4SimHits.Physics.Region = 'HcalRegion'
0078 process.g4SimHits.Physics.DefaultCutValue = 1.
0079 
0080 process.g4SimHits.StackingAction.KillGamma = False
0081 process.g4SimHits.CaloSD.BeamPosition = beamPosition
0082 process.g4SimHits.ECalSD.UseBirkLaw = True
0083 process.g4SimHits.ECalSD.BirkL3Parametrization = True
0084 process.g4SimHits.ECalSD.BirkC1 = 0.033
0085 process.g4SimHits.ECalSD.BirkC2 = 0.0
0086 process.g4SimHits.ECalSD.SlopeLightYield = 0.05
0087 process.g4SimHits.HCalSD.UseBirkLaw = True
0088 process.g4SimHits.HCalSD.BirkC1 = 0.0052
0089 process.g4SimHits.HCalSD.BirkC2 = 0.142
0090 process.g4SimHits.HCalSD.BirkC3 = 1.75
0091 process.g4SimHits.HCalSD.UseLayerWt = False
0092 process.g4SimHits.HCalSD.WtFile     = ' '
0093 process.g4SimHits.HCalSD.UseShowerLibrary = False
0094 process.g4SimHits.HCalSD.TestNumberingScheme = True
0095 process.g4SimHits.HCalSD.UseHF   = False
0096 
0097 process.g4SimHits.HCalSD.ForTBHCAL = True
0098 process.g4SimHits.HCalSD.ForTBH2 = True
0099 process.g4SimHits.OnlySDs = ['CaloTrkProcessing',
0100                              'HcalTB06BeamDetector',
0101                              'EcalSensitiveDetector',
0102                              'HcalSensitiveDetector']
0103 process.g4SimHits.CaloTrkProcessing.TestBeam = True
0104 # Geant4 geometry check 
0105 process.g4SimHits.G4CheckOverlap.OutputBaseName = cms.string("hctb07")
0106 process.g4SimHits.G4CheckOverlap.OverlapFlag = cms.bool(True)
0107 process.g4SimHits.G4CheckOverlap.Tolerance  = cms.double(0.1)
0108 process.g4SimHits.G4CheckOverlap.Resolution = cms.int32(10000)
0109 process.g4SimHits.G4CheckOverlap.Depth      = cms.int32(-1)
0110 process.g4SimHits.G4CheckOverlap.RegionFlag = cms.bool(False)
0111 process.g4SimHits.G4CheckOverlap.NodeNames  = cms.vstring('TBHCal')
0112 process.g4SimHits.G4CheckOverlap.gdmlFlag   = cms.bool(False)
0113 process.g4SimHits.G4CheckOverlap.PVname     = ''
0114 process.g4SimHits.G4CheckOverlap.LVname     = ''
0115 process.g4SimHits.FileNameField   = ''
0116 process.g4SimHits.FileNameGDML    = ''
0117 process.g4SimHits.FileNameRegions = ''