Pythia6Interface/test/Py6GenFilter_ElectronProtonInitialState_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0005 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0006 #process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0007
0008 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0009     generator = cms.PSet(
0010         initialSeed = cms.untracked.uint32(123456789),
0011         engineName = cms.untracked.string('HepJamesRandom')
0012     )
0013 )
0014
0015
0016 # The following three lines reduce the clutter of repeated printouts
0017 # of the same exception message.
0018 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0019
0020 process.MessageLogger.cerr.enableStatistics = False
0021
0022
0023 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(20))
0024
0025 ###process.source = cms.Source("EmptySource")
0026 process.source = cms.Source("EmptySource",
0027      firstLuminosityBlock = cms.untracked.uint32(1),
0028      numberEventsInLuminosityBlock = cms.untracked.uint32(5),
0029      firstEvent = cms.untracked.uint32(1),
0030      firstRun = cms.untracked.uint32(1),
0031      numberEventsInRun = cms.untracked.uint32(10)
0032 )
0033
0034 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0035     pythiaHepMCVerbosity = cms.untracked.bool(True),
0036     maxEventsToPrint = cms.untracked.int32(3),
0037     pythiaPylistVerbosity = cms.untracked.int32(1),
0038     # this shows how to turn ON some of the general Py6 printouts, like banner...
0039     # --> displayPythiaBanner = cms.untracked.bool(True),
0040     displayPythiaCards = cms.untracked.bool(True),
0041     ElectronProtonInitialState = cms.PSet(
0042        electronMomentum = cms.double(400.),
0043        protonMomentum   = cms.double(-400.)
0044     ),
0045     comEnergy = cms.double(0.0), # set to 0.0 as it's NOT used with this initial state
0046                                  # the parameter but must be present (by Fwk)
0047     PythiaParameters = cms.PSet(
0048
0049         pythiaHZZ4tau = cms.vstring('PMAS(25,1)=190.0        !mass of Higgs',
0050             'MSEL=0                  !(D=1) to select between full user control (0, then use MSUB) and some preprogrammed alternative: QCD hight pT processes (1, then ISUB=11, 12, 13, 28, 53, 68), QCD low pT processes (2, then ISUB=11, 12, 13, 28, 53, 68, 91, 92, 94, 95)',
0051             'MSTJ(11)=3              !Choice of the fragmentation function',
0052             'MSTJ(41)=3              !Switch off Pythia QED bremsshtrahlung off leptons (only)',
0053             'MSTP(51)=7              !structure function chosen',
0054             'MSTP(61)=0              ! no initial-state showers',
0055             'MSTP(71)=0              ! no final-state showers',
0056             'MSTP(81)=0              ! no multiple interactions',
0057             'MSTP(111)=0             ! no hadronization',
0058             'MSTU(21)=1              !Check on possible errors during program execution',
0059             'MSUB(102)=1             !ggH',
0060             'MSUB(123)=1             !ZZ fusion to H',
0061             'MSUB(124)=1             !WW fusion to H',
0062             'PARP(82)=1.9            !pt cutoff for multiparton interactions',
0063             'PARP(83)=0.5            !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch',
0064             'PARP(84)=0.4            !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch',
0065             'PARP(90)=0.16           !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch',
0066             'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch',
0067             'CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0068             'CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0069             'CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0070             'MDME(174,1)=0           !Z decay into d dbar',
0071             'MDME(175,1)=0           !Z decay into u ubar',
0072             'MDME(176,1)=0           !Z decay into s sbar',
0073             'MDME(177,1)=0           !Z decay into c cbar',
0074             'MDME(178,1)=0           !Z decay into b bbar',
0075             'MDME(179,1)=0           !Z decay into t tbar',
0076             'MDME(182,1)=0           !Z decay into e- e+',
0077             'MDME(183,1)=0           !Z decay into nu_e nu_ebar',
0078             'MDME(184,1)=0           !Z decay into mu- mu+',
0079             'MDME(185,1)=0           !Z decay into nu_mu nu_mubar',
0080             'MDME(186,1)=1           !Z decay into tau- tau+',
0081             'MDME(187,1)=0           !Z decay into nu_tau nu_taubar',
0082             'MDME(210,1)=0           !Higgs decay into dd',
0083             'MDME(211,1)=0           !Higgs decay into uu',
0084             'MDME(212,1)=0           !Higgs decay into ss',
0085             'MDME(213,1)=0           !Higgs decay into cc',
0086             'MDME(214,1)=0           !Higgs decay into bb',
0087             'MDME(215,1)=0           !Higgs decay into tt',
0088             'MDME(216,1)=0           !Higgs decay into',
0089             'MDME(217,1)=0           !Higgs decay into Higgs decay',
0090             'MDME(218,1)=0           !Higgs decay into e nu e',
0091             'MDME(219,1)=0           !Higgs decay into mu nu mu',
0092             'MDME(220,1)=0           !Higgs decay into tau nu tau',
0093             'MDME(221,1)=0           !Higgs decay into Higgs decay',
0094             'MDME(222,1)=0           !Higgs decay into g g',
0095             'MDME(223,1)=0           !Higgs decay into gam gam',
0096             'MDME(224,1)=0           !Higgs decay into gam Z',
0097             'MDME(225,1)=1           !Higgs decay into Z Z',
0098             'MDME(226,1)=0           !Higgs decay into W W',
0099             'MSTP(128)=0             !dec.prods out of doc section, point at parents in the main section'),
0100
0101         # This is a vector of ParameterSet names to be read, in this order
0102         parameterSets = cms.vstring(
0103             'pythiaHZZ4tau')
0104     )
0105 )
0106
0107 process.GEN = cms.OutputModule("PoolOutputModule",
0108     fileName = cms.untracked.string('Py6_HZZ4tau.root')
0109 )
0110
0111 process.p = cms.Path(process.generator)
0112 process.outpath = cms.EndPath(process.GEN)
0113
0114 process.schedule = cms.Schedule(process.p, process.outpath)