Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 23:30:12

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 
0021 
0022 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(5))
0023 
0024 process.source = cms.Source("EmptySource")
0025 
0026 from Configuration.Generator.PythiaUESettings_cfi import *
0027 from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
0028 
0029 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0030     pythiaHepMCVerbosity = cms.untracked.bool(False),
0031     maxEventsToPrint = cms.untracked.int32(3),
0032     pythiaPylistVerbosity = cms.untracked.int32(1),
0033     # this shows how to turn ON some of the general Py6 printouts, like banner...
0034     displayPythiaBanner = cms.untracked.bool(True),
0035     displayPythiaCards = cms.untracked.bool(True),
0036     comEnergy = cms.double(10000.0),
0037     ExternalDecays = cms.PSet(
0038         Tauola = cms.untracked.PSet(
0039 # these settings below exemplfy how to use "native" Tauola approach:
0040 # one MUST set mdtau=1 !!! then pjak1 & pjak2 will translate into
0041 # native variables jak1 & jak2 (jak1/jak2=4 means that both tau's
0042 # decay into the rho-mode
0043 #
0044          UseTauolaPolarization = cms.bool(True),
0045          InputCards = cms.PSet
0046          ( 
0047             pjak1 = cms.int32(0),
0048         pjak2 = cms.int32(0), 
0049         mdtau = cms.int32(214) 
0050          )
0051 #           TauolaDefaultInputCards,
0052 #      TauolaPolar
0053     ),
0054         parameterSets = cms.vstring('Tauola')
0055     ),
0056     PythiaParameters = cms.PSet(
0057 
0058         pythiaHZZ4tau = cms.vstring('PMAS(25,1)=190.0        !mass of Higgs', 
0059             '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)', 
0060             'MSTJ(11)=3              !Choice of the fragmentation function', 
0061             'MSTJ(41)=1              !Switch off Pythia QED bremsshtrahlung', 
0062             'MSTP(51)=7              !structure function chosen', 
0063             'MSTP(61)=0              ! no initial-state showers', 
0064             'MSTP(71)=0              ! no final-state showers', 
0065             'MSTP(81)=0              ! no multiple interactions', 
0066             'MSTP(111)=0             ! no hadronization', 
0067             'MSTU(21)=1              !Check on possible errors during program execution', 
0068             'MSUB(102)=1             !ggH', 
0069             'MSUB(123)=1             !ZZ fusion to H', 
0070             'MSUB(124)=1             !WW fusion to H', 
0071             'PARP(82)=1.9            !pt cutoff for multiparton interactions', 
0072             'PARP(83)=0.5            !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch', 
0073             'PARP(84)=0.4            !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch', 
0074             'PARP(90)=0.16           !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch', 
0075             'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch', 
0076             'CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
0077             'CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
0078             'CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
0079             'MDME(174,1)=0           !Z decay into d dbar', 
0080             'MDME(175,1)=0           !Z decay into u ubar', 
0081             'MDME(176,1)=0           !Z decay into s sbar', 
0082             'MDME(177,1)=0           !Z decay into c cbar', 
0083             'MDME(178,1)=0           !Z decay into b bbar', 
0084             'MDME(179,1)=0           !Z decay into t tbar', 
0085             'MDME(182,1)=0           !Z decay into e- e+', 
0086             'MDME(183,1)=0           !Z decay into nu_e nu_ebar', 
0087             'MDME(184,1)=0           !Z decay into mu- mu+', 
0088             'MDME(185,1)=0           !Z decay into nu_mu nu_mubar', 
0089             'MDME(186,1)=1           !Z decay into tau- tau+', 
0090             'MDME(187,1)=0           !Z decay into nu_tau nu_taubar', 
0091             'MDME(210,1)=0           !Higgs decay into dd', 
0092             'MDME(211,1)=0           !Higgs decay into uu', 
0093             'MDME(212,1)=0           !Higgs decay into ss', 
0094             'MDME(213,1)=0           !Higgs decay into cc', 
0095             'MDME(214,1)=0           !Higgs decay into bb', 
0096             'MDME(215,1)=0           !Higgs decay into tt', 
0097             'MDME(216,1)=0           !Higgs decay into', 
0098             'MDME(217,1)=0           !Higgs decay into Higgs decay', 
0099             'MDME(218,1)=0           !Higgs decay into e nu e', 
0100             'MDME(219,1)=0           !Higgs decay into mu nu mu', 
0101             'MDME(220,1)=0           !Higgs decay into tau nu tau', 
0102             'MDME(221,1)=0           !Higgs decay into Higgs decay', 
0103             'MDME(222,1)=0           !Higgs decay into g g', 
0104             'MDME(223,1)=0           !Higgs decay into gam gam', 
0105             'MDME(224,1)=0           !Higgs decay into gam Z', 
0106             'MDME(225,1)=1           !Higgs decay into Z Z', 
0107             'MDME(226,1)=0           !Higgs decay into W W', 
0108             'MSTP(128)=0             !dec.prods out of doc section, point at parents in the main section'),
0109 
0110         # This is a vector of ParameterSet names to be read, in this order
0111         parameterSets = cms.vstring( 
0112         ##'pythiaUESettings', 
0113             'pythiaHZZ4tau')
0114     )
0115 )
0116 
0117 process.writer = cms.EDAnalyzer("HepMCEventWriter",
0118                                 hepMCProduct = cms.InputTag("VtxSmeared"))
0119 
0120 process.p = cms.Path(process.generator)
0121 process.outpath = cms.EndPath(process.writer)
0122 
0123 process.schedule = cms.Schedule(process.p, process.outpath)