Pythia6Interface/test/Py6GenFilter_EvtGenTest_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #
0004 # WARNING: This is NOT an example for users -
0005 #          it's my private (JY) "development" cfg,
0006 #          for testing several ideas and/or suspects.
0007 #          Do NOT use it, especially if you're just
0008 #          learning things !!!
0009 #
0010
0011 process = cms.Process("TEST")
0012 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0013 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0014 #process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0015
0016 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0017     generator = cms.PSet(
0018         initialSeed = cms.untracked.uint32(123456789),
0019         engineName = cms.untracked.string('HepJamesRandom')
0020     )
0021 )
0022
0023
0024 # The following three lines reduce the clutter of repeated printouts
0025 # of the same exception message.
0026 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0027
0028
0029
0030
0031 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(5))
0032
0033 process.source = cms.Source("EmptySource")
0034
0035 from Configuration.Generator.PythiaUESettings_cfi import *
0036 from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
0037
0038 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0039     pythiaHepMCVerbosity = cms.untracked.bool(False),
0040     maxEventsToPrint = cms.untracked.int32(3),
0041     pythiaPylistVerbosity = cms.untracked.int32(1),
0042     # this shows how to turn ON some of the general Py6 printouts, like banner...
0043     ## --> displayPythiaBanner = cms.untracked.bool(True),
0044     ## --> displayPythiaCards = cms.untracked.bool(True),
0045     comEnergy = cms.double(10000.0),
0046
0047     ExternalDecays = cms.PSet(
0048
0049         Tauola = cms.untracked.PSet(
0050 # these settings below exemplfy how to use "native" Tauola approach:
0051 # one MUST set mdtau=1 !!! then pjak1 & pjak2 will translate into
0052 # native variables jak1 & jak2 (jak1/jak2=4 means that both tau's
0053 # decay into the rho-mode
0054 #
0055          UseTauolaPolarization = cms.bool(True),
0056          InputCards = cms.PSet
0057          (
0058             pjak1 = cms.int32(0),
0059         pjak2 = cms.int32(0),
0060         mdtau = cms.int32(202)
0061          )
0062 #           TauolaDefaultInputCards,
0063 #      TauolaPolar
0064     ),
0065
0066         EvtGen = cms.untracked.PSet(
0067              operates_on_particles = cms.vint32( 0 ), # 0 (zero) means default list (hardcoded)
0068                                                       # you can put here the list of particles (PDG IDs)
0069                                                       # that you want decayed by EvtGen
0070          use_internal_pythia = cms.untracked.bool(False),
0071          use_default_decay = cms.untracked.bool(True),
0072              decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
0073              # decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY.DEC'),
0074              particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
0075              user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Onia_mumu.dec'), # irrelevat in this case
0076              list_forced_decays = cms.vstring(),
0077              ),
0078
0079
0080
0081         parameterSets = cms.vstring( 'Tauola', # do NOT use these 2 together !
0082                                  'EvtGen'
0083                    )
0084     ),
0085
0086     PythiaParameters = cms.PSet(
0087
0088         pythiaGeneralSettings = cms.vstring(
0089         'PMAS(25,1)=190.0        !mass of Higgs',
0090             'MSTP(128)=0             !dec.prods out of doc section, point at parents in the main section',
0091             '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)',
0092             'MSTJ(11)=3              !Choice of the fragmentation function',
0093             #'MSTJ(41)=1              !Switch off Pythia QED bremsshtrahlung',
0094             'MSTP(51)=7              !structure function chosen',
0095             #'MSTP(61)=0              ! no initial-state showers',
0096             #'MSTP(71)=0              ! no final-state showers',
0097             #'MSTP(81)=0              ! no multiple interactions',
0098             #'MSTP(111)=0             ! no hadronization',
0099             'MSTU(21)=1              !Check on possible errors during program execution',
0100             'MSUB(102)=1             !ggH',
0101             'MSUB(123)=1             !ZZ fusion to H',
0102             'MSUB(124)=1             !WW fusion to H',
0103         # these 4 below are irrelevant if the hadronization is off (mstp(111)=0)
0104             'PARP(82)=1.9            !pt cutoff for multiparton interactions',
0105             'PARP(83)=0.5            !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch',
0106             'PARP(84)=0.4            !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch',
0107             'PARP(90)=0.16           !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch',
0108         #........................................
0109             'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch',
0110             'CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0111             'CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0112             'CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch'
0113     ),
0114     pythiaHZZDecays = cms.vstring(
0115             'MDME(174,1)=0           !Z decay into d dbar',
0116             'MDME(175,1)=0           !Z decay into u ubar',
0117             'MDME(176,1)=0           !Z decay into s sbar',
0118             'MDME(177,1)=0           !Z decay into c cbar',
0119             # set it to 4 for the 1st Z to go to b bbar...
0120         'MDME(178,1)=4           !Z decay into b bbar',
0121         # ............................................
0122             'MDME(179,1)=0           !Z decay into t tbar',
0123             'MDME(182,1)=0           !Z decay into e- e+',
0124             'MDME(183,1)=0           !Z decay into nu_e nu_ebar',
0125             'MDME(184,1)=0           !Z decay into mu- mu+',
0126             'MDME(185,1)=0           !Z decay into nu_mu nu_mubar',
0127         # ... and this one to 5, for the 2nd Z to goe to tau- tau+
0128             'MDME(186,1)=5           !Z decay into tau- tau+',
0129         # ............................................
0130             'MDME(187,1)=0           !Z decay into nu_tau nu_taubar',
0131             'MDME(210,1)=0           !Higgs decay into dd',
0132             'MDME(211,1)=0           !Higgs decay into uu',
0133             'MDME(212,1)=0           !Higgs decay into ss',
0134             'MDME(213,1)=0           !Higgs decay into cc',
0135             'MDME(214,1)=0           !Higgs decay into bb',
0136             'MDME(215,1)=0           !Higgs decay into tt',
0137             'MDME(216,1)=0           !Higgs decay into',
0138             'MDME(217,1)=0           !Higgs decay into Higgs decay',
0139             'MDME(218,1)=0           !Higgs decay into e nu e',
0140             'MDME(219,1)=0           !Higgs decay into mu nu mu',
0141             'MDME(220,1)=0           !Higgs decay into tau nu tau',
0142             'MDME(221,1)=0           !Higgs decay into Higgs decay',
0143             'MDME(222,1)=0           !Higgs decay into g g',
0144             'MDME(223,1)=0           !Higgs decay into gam gam',
0145             'MDME(224,1)=0           !Higgs decay into gam Z',
0146             'MDME(225,1)=1           !Higgs decay into Z Z',
0147             'MDME(226,1)=0           !Higgs decay into W W'
0148          ),
0149         # This is a vector of ParameterSet names to be read, in this order
0150         parameterSets = cms.vstring(
0151         ##'pythiaUESettings',
0152         'pythiaGeneralSettings',
0153             'pythiaHZZDecays')
0154     )
0155 )
0156
0157 process.GEN = cms.OutputModule("PoolOutputModule",
0158     fileName = cms.untracked.string('TestHZZ4tau.root')
0159 )
0160
0161 process.p = cms.Path(process.generator)
0162 process.outpath = cms.EndPath(process.GEN)
0163
0164 process.schedule = cms.Schedule(process.p, process.outpath)