1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("TQAF")
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
process.source = cms.Source("EmptySource")
process.generator = cms.EDFilter("Pythia8GeneratorFilter",
maxEventsToPrint = cms.untracked.int32(1),
pythiaPylistVerbosity = cms.untracked.int32(1),
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(True),
comEnergy = cms.double(13000.),
PythiaParameters = cms.PSet(
processParameters = cms.vstring('Top:all = on'),
parameterSets = cms.vstring('processParameters')
)
)
process.genParticles = cms.EDProducer("GenParticleProducer",
abortOnUnknownPDGCode = cms.untracked.bool(False),
saveBarCodes = cms.untracked.bool(True),
src = cms.InputTag("generator:unsmeared")
)
## define maximal number of events to loop over
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(100)
)
## configure process options
process.options = cms.untracked.PSet(
allowUnscheduled = cms.untracked.bool(True),
wantSummary = cms.untracked.bool(True)
)
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load("GeneratorInterface.RivetInterface.genParticles2HepMC_cfi")
process.load("GeneratorInterface.RivetInterface.particleLevel_cfi")
process.particleLevel.src = cms.InputTag("genParticles2HepMC:unsmeared")
process.path = cms.Path(process.generator*process.genParticles*process.genParticles2HepMC*process.particleLevel)
process.out = cms.OutputModule("PoolOutputModule",
fileName = cms.untracked.string("particleLevel.root"),
outputCommands = cms.untracked.vstring(
"drop *",
"keep *_genParticles_*_*",
"keep *_particleLevel_*_*",
),
)
process.outPath = cms.EndPath(process.out)
|