Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:30:12

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 process = cms.Process("HFShowerLib")
0004 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0005 #process.load("Geometry.HcalCommonData.hcalforwardshowerLong_cfi")
0006 process.load("SimG4CMS.ShowerLibraryProducer.hcalforwardshower_cfi")
0007 process.load("Geometry.HcalCommonData.hcalDDConstants_cff")
0008 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0009 process.load('Configuration.StandardSequences.Generator_cff')
0010 process.load('Configuration.StandardSequences.SimIdeal_cff')
0011 process.load("Configuration.StandardSequences.MagneticField_cff")
0012 process.load("Configuration.StandardSequences.Services_cff")
0013 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0014 
0015 if 'MessageLogger' in process.__dict__:
0016     process.MessageLogger.FiberSim=dict()
0017     process.MessageLogger.FlatThetaGun=dict()
0018     process.MessageLogger.HFShower=dict()
0019     process.MessageLogger.HcalForwardLib=dict()
0020     process.MessageLogger.HCalGeom=dict()
0021     process.MessageLogger.SensitiveDetector=dict()
0022 
0023 process.RandomNumberGeneratorService.generator.initialSeed = 12345
0024 
0025 from IOMC.RandomEngine.RandomServiceHelper import RandomNumberServiceHelper
0026 randSvc = RandomNumberServiceHelper(process.RandomNumberGeneratorService)
0027 randSvc.populate()
0028 
0029 process.Timing = cms.Service("Timing")
0030 
0031 process.maxEvents = cms.untracked.PSet(
0032     input = cms.untracked.int32(1)
0033 )
0034 
0035 process.source = cms.Source("EmptySource")
0036 
0037 process.generator = cms.EDProducer("FlatRandomEThetaGunProducer",
0038     PGunParameters = cms.PSet(
0039         PartID   = cms.vint32(11),
0040         #MinTheta = cms.double(-1.145762838),
0041         #MaxTheta = cms.double(1.145762838),
0042         MinTheta = cms.double(0.019997),
0043         MaxTheta = cms.double(0.019997),
0044         MinPhi   = cms.double(3.14500926),
0045         MaxPhi   = cms.double(3.14500926),
0046         MinE     = cms.double(100.0),
0047         MaxE     = cms.double(100.0)
0048     ),
0049     Verbosity = cms.untracked.int32(2),
0050     AddAntiParticle = cms.bool(False),
0051     firstRun = cms.untracked.uint32(1)
0052 )
0053 
0054 process.o1 = cms.OutputModule("PoolOutputModule",
0055     fileName = cms.untracked.string('simevent_50GeVElec.root')
0056 )
0057 
0058 process.TFileService = cms.Service("TFileService",
0059     fileName = cms.string('hfShowerLibSimu_extended2_50GeVElec_deneme.root')
0060 )
0061 
0062 # for GEN produced since 760pre6, for older GEN - just "":
0063 process.VtxSmeared.src = cms.InputTag("generator", "unsmeared")
0064 process.generatorSmeared = cms.EDProducer("GeneratorSmearedProducer")
0065 process.g4SimHits.Generator.HepMCProductLabel = cms.InputTag('VtxSmeared')
0066 
0067 process.p1 = cms.Path(
0068  process.generator *
0069  process.VtxSmeared *
0070  process.generatorSmeared *
0071  process.g4SimHits
0072 )
0073 process.outpath = cms.EndPath(process.o1)
0074 
0075 process.g4SimHits.HCalSD.UseShowerLibrary = True
0076 process.g4SimHits.HCalSD.UseParametrize = False
0077 process.g4SimHits.HCalSD.UsePMTHits = False
0078 process.g4SimHits.HCalSD.UseFibreBundleHits = False
0079 process.g4SimHits.OnlySDs = ['CaloTrkProcessing',
0080                              'HFWedgeSensitiveDetector',
0081                              'FiberSensitiveDetector',
0082                              'HFChamberSensitiveDetector',
0083                              'HcalSensitiveDetector']
0084 
0085 process.g4SimHits.HFShower.UseShowerLibrary= True
0086 process.g4SimHits.HFShower.UseR7600UPMT    = True
0087 process.g4SimHits.HFShower.UseHFGflash = False
0088 process.g4SimHits.HFShower.ApplyFiducialCut = False
0089 process.g4SimHits.UseMagneticField = False
0090 
0091 process.g4SimHits.NonBeamEvent = True
0092 process.g4SimHits.Generator.ApplyPCuts   = False
0093 process.g4SimHits.Generator.ApplyEtaCuts = False
0094 process.g4SimHits.Physics.type = 'SimG4Core/Physics/FTFP_BERT_EMM'
0095 process.g4SimHits.Physics.DefaultCutValue = 0.1
0096 process.g4SimHits.G4Commands = ['/tracking/verbose 1']
0097 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0098     HFShowerLibraryProducer = cms.PSet(
0099         Names = cms.vstring('FibreHits', 
0100             'ChamberHits', 
0101             'WedgeHits')
0102     ),
0103     type = cms.string('HcalForwardAnalysis')
0104 ))
0105 
0106