RivetInterface/test/runRivetHTXS_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 process = cms.Process("runRivetAnalysis")
0003
0004 process.options   = cms.untracked.PSet(
0005     allowUnscheduled = cms.untracked.bool(False)
0006 )
0007
0008 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
0009 process.source = cms.Source("PoolSource",  fileNames = cms.untracked.vstring(
0010 # compare AOD and MINIAOD
0011 #'/store/mc/RunIISpring16MiniAODv2/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUgenV6_pythia8/MINIAODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/10000/2C7F3153-393B-E611-9323-0CC47AA98A3A.root'
0012 #'/store/mc/RunIISpring16reHLT80/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUgenV6_pythia8/AODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/10000/52E144D7-793A-E611-B70F-0025904A8ECC.root',
0013 #'/store/mc/RunIISpring16reHLT80/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUgenV6_pythia8/AODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/10000/D650466F-A13A-E611-AA11-0CC47A13CD56.root'
0014
0015 # just run some MINIAOD
0016 '/store/mc/RunIIAutumn18MiniAOD/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUGenV7011_pythia8/MINIAODSIM/102X_upgrade2018_realistic_v15-v2/80000/6F4F411E-8111-684D-827D-B5962A0CB94F.root',
0017 '/store/mc/RunIIAutumn18MiniAOD/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUGenV7011_pythia8/MINIAODSIM/102X_upgrade2018_realistic_v15-v2/80000/49CB36B2-E124-2249-A0F8-CE867CF4F8A6.root',
0018 '/store/mc/RunIIAutumn18MiniAOD/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUGenV7011_pythia8/MINIAODSIM/102X_upgrade2018_realistic_v15-v2/80000/79F49EC1-42B4-3349-A268-59510E899BCC.root',
0019 '/store/mc/RunIIAutumn18MiniAOD/GluGluHToZZTo4L_M125_13TeV_powheg2_JHUGenV7011_pythia8/MINIAODSIM/102X_upgrade2018_realistic_v15-v2/80000/D65A4D51-2E80-AD41-B50D-E4083BA2A668.root',
0020
0021 ),
0022 )
0023
0024 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0025
0026 process.rivetProducerHTXS = cms.EDProducer('HTXSRivetProducer',
0027   HepMCCollection = cms.InputTag('myGenerator','unsmeared'),
0028   LHERunInfo = cms.InputTag('externalLHEProducer'),
0029   #ProductionMode = cms.string('GGF'),
0030   ProductionMode = cms.string('AUTO'),
0031 )
0032
0033 #MINIAOD
0034 process.mergedGenParticles = cms.EDProducer("MergedGenParticleProducer",
0035     inputPruned = cms.InputTag("prunedGenParticles"),
0036     inputPacked = cms.InputTag("packedGenParticles"),
0037 )
0038 process.myGenerator = cms.EDProducer("GenParticles2HepMCConverter",
0039     genParticles = cms.InputTag("mergedGenParticles"),
0040     genEventInfo = cms.InputTag("generator"),
0041 )
0042 process.p = cms.Path(process.mergedGenParticles*process.myGenerator*process.rivetProducerHTXS)
0043
0044 # # AOD
0045 #process.myGenerator = cms.EDProducer("GenParticles2HepMCConverterHTXS",
0046 #    genParticles = cms.InputTag("genParticles"),
0047 #    genEventInfo = cms.InputTag("generator"),
0048 #)
0049 #process.p = cms.Path(process.myGenerator*process.rivetProducerHTXS)
0050
0051 process.out = cms.OutputModule("PoolOutputModule",
0052     outputCommands = cms.untracked.vstring('drop *','keep *_rivetProducerHTXS_*_runRivetAnalysis','keep *_generator_*_*','keep *_externalLHEProducer_*_*'),
0053     fileName = cms.untracked.string('testHTXSRivet_ggH4l_MINIAOD_100k.root')
0054 )
0055 process.o = cms.EndPath( process.out )