" non-diffractive xsec = 17.0318 mb with nd factor = 0.3142."),
Run=cms.vstring(
" (run){", " CORE_SCALE TTBar;", " METS_BBAR_MODE 5;",
" NJET:=2; LJET:=2; QCUT:=20.;",
" ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;",
" OL_PREFIX={0};".format(
os.environ['CMS_OPENLOOPS_PREFIX'] if 'CMS_OPENLOOPS_PREFIX' in
os.environ else ''), " LOOPGEN:=OpenLoops;",
" MI_HANDLER=Amisic;", " NLO_SMEAR_THRESHOLD 1;",
" NLO_SMEAR_POWER 2;", " HARD_DECAYS On;",
" HARD_SPIN_CORRELATIONS=1;", " SOFT_SPIN_CORRELATIONS=1;",
" PDF_LIBRARY LHAPDFSherpa;",
" PDF_SET NNPDF31_nnlo_hessian_pdfas;", " USE_PDF_ALPHAS=1;",
" BEAM_1=2212; BEAM_ENERGY_1=6500;",
" BEAM_2=2212; BEAM_ENERGY_2=6500;",
" STABLE[6] 0; WIDTH[6] 0; STABLE[24] 0;",
" EXCLUSIVE_CLUSTER_MODE 1;", " HEPMC_TREE_LIKE=1;",
" PRETTY_PRINT=Off;", "}(run)", " (processes){",
" Process : 93 93 -> 6 -6 93{NJET};",
" Order (*,0); CKKW sqr(QCUT/E_CMS);",
" NLO_QCD_Mode MC@NLO {LJET};", " ME_Generator Amegic {LJET};",
" RS_ME_Generator Comix {LJET};",
" Loop_Generator LOOPGEN {LJET};",
" Integration_Error 0.05 {3,4};", " End process;", "}(processes)"),
))
from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)
ProductionFilterSequence = cms.Sequence(generator)
'MultipartonInteractions:ecmPow=0.25208',
'MultipartonInteractions:expPow=1.6'),
pythia8CommonSettings=cms.vstring(
'Tune:preferLHAPDF = 2', 'Main:timesAllowErrors = 10000',
'Check:epTolErr = 0.01', 'Beams:setProductionScalesFromLHEF = off',
'SLHA:minMassSM = 1000.', 'ParticleDecays:limitTau0 = on',
'ParticleDecays:tau0Max = 10',
'ParticleDecays:allowPhotonRadiation = on')),
comEnergy=cms.double(13000.0),
maxEventsToPrint=cms.untracked.int32(0),
pythiaHepMCVerbosity=cms.untracked.bool(False),
pythiaPylistVerbosity=cms.untracked.int32(0))
from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
process.generator = ExternalGeneratorFilter(_generator)
# The line below removes messages like "particle not recognized by pythia"
from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
process.generator.PythiaParameters.processParameters.extend(
EvtGenExtraParticles)
process.bufilter = cms.EDFilter("PythiaDauVFilter",
DaughterIDs=cms.untracked.vint32(20443, 321),
MaxEta=cms.untracked.vdouble(9999.0, 2.5),
MinEta=cms.untracked.vdouble(-9999.0, -2.5),
MinPt=cms.untracked.vdouble(0.0, 0.4),
MotherID=cms.untracked.int32(0),
NumberDaughters=cms.untracked.int32(2),
ParticleID=cms.untracked.int32(521))
" NJET:=2; LJET:=2; QCUT:=20.;",
" ME_SIGNAL_GENERATOR Comix Amegic LOOPGEN;",
" OL_PREFIX={0} ".format(os.environ['CMS_OPENLOOPS_PREFIX']),
" LOOPGEN:=OpenLoops;", " MI_HANDLER=Amisic;",
" NLO_SMEAR_THRESHOLD 1;", " NLO_SMEAR_POWER 2;",
" HARD_DECAYS On;", " HARD_SPIN_CORRELATIONS=1;",
" SOFT_SPIN_CORRELATIONS=1;", " PDF_LIBRARY LHAPDFSherpa;",
" PDF_SET NNPDF31_nnlo_hessian_pdfas;", " USE_PDF_ALPHAS=1;",
" BEAM_1=2212; BEAM_ENERGY_1=6500;",
" BEAM_2=2212; BEAM_ENERGY_2=6500;",
" STABLE[6] 0; WIDTH[6] 0; STABLE[24] 0;",
" EXCLUSIVE_CLUSTER_MODE 1;", " HEPMC_TREE_LIKE=1;",
" PRETTY_PRINT=Off;", "}(run)", " (processes){",
" Process : 93 93 -> 6 -6 93{NJET};",
" Order (*,0); CKKW sqr(QCUT/E_CMS);",
" NLO_QCD_Mode MC@NLO {LJET};", " ME_Generator Amegic {LJET};",
" RS_ME_Generator Comix {LJET};",
" Loop_Generator LOOPGEN {LJET};",
" Integration_Error 0.05 {3,4};", " End process;", "}(processes)"),
))
from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(
_generator,
_external_process_waitTime_=cms.untracked.uint32(
1200), # increase the max waitTime for Sherpa
)
ProductionFilterSequence = cms.Sequence(generator)
"Pythia8GeneratorFilter",
pythiaHepMCVerbosity=cms.untracked.bool(False),
maxEventsToPrint=cms.untracked.int32(0),
pythiaPylistVerbosity=cms.untracked.int32(0),
filterEfficiency=cms.untracked.double(1.0),
comEnergy=cms.double(8000.0),
ExternalDecays=cms.PSet(Tauola=cms.untracked.PSet(
TauolaPolar,
InputCards=cms.PSet(pjak1=cms.int32(0),
pjak2=cms.int32(0),
mdtau=cms.int32(214)),
),
parameterSets=cms.vstring('Tauola')),
PythiaParameters=cms.PSet(pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters=cms.vstring(
'WeakSingleBoson:ffbar2gmZ = on',
'23:onMode = off',
'23:onIfAny = 15',
'PhaseSpace:mHatMin = 75.',
),
parameterSets=cms.vstring(
'pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters',
)))
from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(
_generator, _external_process_components_=cms.vstring("HepPDTESSource"))
process.RandomNumberGeneratorService.gen2 = process.RandomNumberGeneratorService.generator.clone(
)
process.gen1 = cms.EDFilter(
"Pythia8GeneratorFilter",
comEnergy=cms.double(7000.),
PythiaParameters=cms.PSet(pythia8_example02=cms.vstring(
'HardQCD:all = on', 'PhaseSpace:pTHatMin = 20.'),
parameterSets=cms.vstring('pythia8_example02')))
_pythia8 = cms.EDFilter(
"Pythia8GeneratorFilter",
comEnergy=cms.double(7000.),
PythiaParameters=cms.PSet(pythia8_example02=cms.vstring(
'HardQCD:all = on', 'PhaseSpace:pTHatMin = 20.'),
parameterSets=cms.vstring('pythia8_example02')))
from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
process.gen2 = ExternalGeneratorFilter(_pythia8)
process.sleeper = cms.EDProducer("timestudy::SleepingProducer",
ivalue=cms.int32(1),
consumes=cms.VInputTag(),
eventTimes=cms.vdouble(1.))
process.compare = cms.EDAnalyzer("CompareGeneratorResultsAnalyzer",
module1=cms.untracked.string("gen1"),
module2=cms.untracked.string("gen2"),
allowXSecDifferences=cms.untracked.bool(True))
process.p = cms.Path(process.sleeper + process.gen1 + process.gen2 +
process.compare)