Configuration/test/gluglu_step1_GEN_SIM_2021.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 import random
0004 import math
0005
0006 from Configuration.StandardSequences.Eras import eras
0007 process = cms.Process('SIM',eras.Run3)
0008
0009 # import of standard configurations
0010 process.load('Configuration.StandardSequences.Services_cff')
0011 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0012 process.load('FWCore.MessageService.MessageLogger_cfi')
0013 process.load('Configuration.EventContent.EventContent_cff')
0014 process.load('Configuration.StandardSequences.MagneticField_cff')
0015 process.load('Configuration.StandardSequences.Generator_cff')
0016 process.load('IOMC.EventVertexGenerators.VtxSmearedHLLHC14TeV_cfi')
0017 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0018 process.load('Configuration.StandardSequences.SimIdeal_cff')
0019 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0020 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0021 process.load('Configuration.StandardSequences.GeometrySimDB_cff')
0022
0023 process.RandomNumberGeneratorService.generator.initialSeed = cms.untracked.uint32(random.randint(0,900000000))
0024
0025 process.load('SimG4Core.Application.g4SimHits_cfi')
0026 process.g4SimHits.LHCTransport = cms.bool(True)
0027
0028 nEvent_ = 100
0029 process.maxEvents = cms.untracked.PSet(
0030         input = cms.untracked.int32(nEvent_)
0031         )
0032
0033 process.source = cms.Source("EmptySource")
0034
0035 from Configuration.AlCa.GlobalTag import GlobalTag
0036 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase1_2022_realistic', '')
0037
0038 # generator
0039 process.generator = cms.EDFilter("ExhumeGeneratorFilter",
0040     ExhumeParameters = cms.PSet(
0041         AlphaEw = cms.double(0.0072974),
0042         B = cms.double(4.0),
0043         MuonMass = cms.double(0.1057),
0044         BottomMass = cms.double(4.65),
0045         CharmMass = cms.double(1.28),
0046         StrangeMass = cms.double(0.095),
0047         TauMass = cms.double(1.78),
0048         TopMass = cms.double(172.8),
0049         WMass = cms.double(80.38),
0050         ZMass = cms.double(91.187),
0051         HiggsMass = cms.double(125.1),
0052         HiggsVev = cms.double(246.0),
0053         LambdaQCD = cms.double(80.0),
0054         MinQt2 = cms.double(0.64),
0055         PDF = cms.double(11000),
0056         Rg = cms.double(1.2),
0057         Survive = cms.double(0.03)
0058     ),
0059     ExhumeProcess = cms.PSet(
0060         MassRangeHigh = cms.double(2000.0),
0061         MassRangeLow = cms.double(300.0),
0062         ProcessType = cms.string('GG'),
0063         ThetaMin = cms.double(0.3)
0064     ),
0065     PythiaParameters = cms.PSet(
0066         parameterSets = cms.vstring()
0067     ),
0068     comEnergy = cms.double(14000.0),
0069     maxEventsToPrint = cms.untracked.int32(2),
0070     pythiaHepMCVerbosity = cms.untracked.bool(False),
0071     pythiaPylistVerbosity = cms.untracked.int32(1)
0072 )
0073
0074 # Production Info
0075 process.configurationMetadata = cms.untracked.PSet(
0076     annotation = cms.untracked.string('GluGluTo2Jets_M_100_7TeV_exhume_cff.py nevts:1'),
0077     name = cms.untracked.string('Applications'),
0078     version = cms.untracked.string('$Revision: 1.19 $')
0079 )
0080
0081 process.ProductionFilterSequence = cms.Sequence(process.generator)
0082
0083 ############
0084 process.o1 = cms.OutputModule("PoolOutputModule",
0085         outputCommands = cms.untracked.vstring('keep *'),
0086         fileName = cms.untracked.string('GluGlu_step1_GEN_SIM_2021.root')
0087         )
0088
0089 process.generation_step = cms.Path(process.pgen)
0090 process.simulation_step = cms.Path(process.psim)
0091
0092 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0093 process.outpath = cms.EndPath(process.o1)
0094 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.outpath)
0095
0096 # filter all path with the production filter sequence
0097 for path in process.paths:
0098     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
0099