Pythia6Interface/test/Py6GenFilter_Photos_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, for testing
0006 #          newly implemented PhotosInterface - which is NOT yet
0007 #          released via ExternalDecayDeriver
0008 #
0009
0010 process = cms.Process("TEST")
0011 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0012 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0013 #process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0014
0015 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0016     generator = cms.PSet(
0017         initialSeed = cms.untracked.uint32(123456789),
0018         engineName = cms.untracked.string('HepJamesRandom')
0019     )
0020 )
0021
0022
0023 # The following three lines reduce the clutter of repeated printouts
0024 # of the same exception message.
0025 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0026
0027
0028
0029
0030 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(20))
0031
0032 process.source = cms.Source("EmptySource")
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     comEnergy = cms.double(10000.0),
0042
0043     ExternalDecays = cms.PSet(
0044         Photos = cms.untracked.PSet(),
0045         parameterSets = cms.vstring( "Photos" )
0046     ),
0047
0048     PythiaParameters = cms.PSet(
0049
0050         pythiaSimpleSettings = cms.vstring(
0051        'PMAS(5,1)=4.8 ! b quark mass',
0052            'PMAS(6,1)=172.3 ! t quark mass',
0053            'MSTP(61)=0 ! turn off initial state radiation',
0054            'mstj(41)=1' # per Steve M., instead of mstp(71)... btw, shoult it be 0 or 1 ?
0055            #'MSTP(71)=0 ! final state radiation (or final-state showers ?)'
0056
0057     ),
0058     pythiaSpecialSettings = cms.vstring(
0059         'PMAS(25,1)=190.0        !mass of Higgs',
0060             '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)',
0061             'MSTJ(11)=3              !Choice of the fragmentation function',
0062             'MSTJ(41)=1              !Switch off Pythia QED bremsshtrahlung',
0063             'MSTP(51)=7              !structure function chosen',
0064             'MSTP(61)=0              ! no initial-state showers',
0065             'MSTP(71)=0              ! no final-state showers',
0066             'MSTP(81)=0              ! no multiple interactions',
0067             'MSTP(111)=0             ! no hadronization',
0068             'MSTU(21)=1              !Check on possible errors during program execution',
0069         # these 4 below are irrelevant if the hadronization is off (mstp(111)=0)
0070             'PARP(82)=1.9            !pt cutoff for multiparton interactions',
0071             'PARP(83)=0.5            !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch',
0072             'PARP(84)=0.4            !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch',
0073             'PARP(90)=0.16           !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch',
0074         #........................................
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     ),
0080     pythiaProcessSettings = cms.vstring(
0081             'MSUB(102)=1             !ggH',
0082             'MSUB(123)=1             !ZZ fusion to H',
0083             'MSUB(124)=1             !WW fusion to H',
0084     ),
0085     pythiaHZZDecays = cms.vstring(
0086             'MDME(174,1)=0           !Z decay into d dbar',
0087             'MDME(175,1)=0           !Z decay into u ubar',
0088             'MDME(176,1)=0           !Z decay into s sbar',
0089             'MDME(177,1)=0           !Z decay into c cbar',
0090             # set it to 4 for the 1st Z to go to b bbar...
0091         'MDME(178,1)=0           !Z decay into b bbar',
0092         # ............................................
0093             'MDME(179,1)=0           !Z decay into t tbar',
0094             'MDME(182,1)=0           !Z decay into e- e+',
0095             'MDME(183,1)=0           !Z decay into nu_e nu_ebar',
0096             'MDME(184,1)=4           !Z decay into mu- mu+',
0097             'MDME(185,1)=0           !Z decay into nu_mu nu_mubar',
0098         # ... and this one to 5, for the 2nd Z to goe to tau- tau+
0099             'MDME(186,1)=5           !Z decay into tau- tau+',
0100         # ............................................
0101             'MDME(187,1)=0           !Z decay into nu_tau nu_taubar',
0102             'MDME(210,1)=0           !Higgs decay into dd',
0103             'MDME(211,1)=0           !Higgs decay into uu',
0104             'MDME(212,1)=0           !Higgs decay into ss',
0105             'MDME(213,1)=0           !Higgs decay into cc',
0106             'MDME(214,1)=0           !Higgs decay into bb',
0107             'MDME(215,1)=0           !Higgs decay into tt',
0108             'MDME(216,1)=0           !Higgs decay into',
0109             'MDME(217,1)=0           !Higgs decay into Higgs decay',
0110             'MDME(218,1)=0           !Higgs decay into e nu e',
0111             'MDME(219,1)=0           !Higgs decay into mu nu mu',
0112             'MDME(220,1)=0           !Higgs decay into tau nu tau',
0113             'MDME(221,1)=0           !Higgs decay into Higgs decay',
0114             'MDME(222,1)=0           !Higgs decay into g g',
0115             'MDME(223,1)=0           !Higgs decay into gam gam',
0116             'MDME(224,1)=0           !Higgs decay into gam Z',
0117             'MDME(225,1)=1           !Higgs decay into Z Z',
0118             'MDME(226,1)=0           !Higgs decay into W W'
0119         ),
0120         pythiaTauL = cms.vstring(
0121            "mdme(89,1)=1", # tau -> e
0122            "mdme(90,1)=1", # tau -> mu
0123            # all other tau decays OFF
0124            "mdme(91,1)=0",
0125            "mdme(92,1)=0",
0126            "mdme(93,1)=0",
0127            "mdme(94,1)=0",
0128            "mdme(95,1)=0",
0129            "mdme(96,1)=0",
0130            "mdme(97,1)=0",
0131            "mdme(98,1)=0",
0132            "mdme(99,1)=0",
0133            "mdme(100,1)=0",
0134            "mdme(101,1)=0",
0135            "mdme(102,1)=0",
0136            "mdme(103,1)=0",
0137            "mdme(104,1)=0",
0138            "mdme(105,1)=0",
0139            "mdme(106,1)=0",
0140            "mdme(107,1)=0",
0141            "mdme(108,1)=0",
0142            "mdme(109,1)=0",
0143            "mdme(110,1)=0",
0144            "mdme(111,1)=0",
0145            "mdme(112,1)=0",
0146            "mdme(113,1)=0",
0147            "mdme(114,1)=0",
0148            "mdme(115,1)=0",
0149            "mdme(116,1)=0",
0150            "mdme(117,1)=0",
0151            "mdme(118,1)=0",
0152            "mdme(119,1)=0",
0153            "mdme(120,1)=0",
0154            "mdme(121,1)=0",
0155            "mdme(122,1)=0",
0156            "mdme(123,1)=0",
0157            "mdme(124,1)=0",
0158            "mdme(125,1)=0",
0159            "mdme(126,1)=0",
0160            "mdme(127,1)=0",
0161            "mdme(128,1)=0",
0162            "mdme(129,1)=0",
0163            "mdme(130,1)=0",
0164            "mdme(131,1)=0",
0165            "mdme(132,1)=0",
0166            "mdme(133,1)=0",
0167            "mdme(134,1)=0",
0168            "mdme(135,1)=0",
0169            "mdme(136,1)=0",
0170            "mdme(137,1)=0",
0171            "mdme(138,1)=0",
0172            "mdme(139,1)=0",
0173            "mdme(140,1)=0",
0174            "mdme(141,1)=0",
0175            "mdme(142,1)=0"
0176         ),
0177         # This is a vector of ParameterSet names to be read, in this order
0178         parameterSets = cms.vstring(
0179         'pythiaSpecialSettings',
0180         ##'pythiaSimpleSettings',
0181         'pythiaProcessSettings',
0182             'pythiaHZZDecays', 'pythiaTauL' )
0183     )
0184 )
0185
0186 process.GEN = cms.OutputModule("PoolOutputModule",
0187     fileName = cms.untracked.string('TestHZZleptons.root')
0188 )
0189
0190 process.p = cms.Path(process.generator)
0191 process.outpath = cms.EndPath(process.GEN)
0192
0193 process.schedule = cms.Schedule(process.p, process.outpath)