Pythia6Interface/test/ZtautauLeptonsNoFSRPtGun_cfg.py

0001 # The following comments couldn't be translated into the new config version:
0002
0003 # Single tau(+decays) ptgun
0004
0005 import FWCore.ParameterSet.Config as cms
0006
0007 process = cms.Process("Gen")
0008 # this example configuration offers some minimum
0009 # annotation, to help users get through; please
0010 # don't hesitate to read through the comments
0011 # use MessageLogger to redirect/suppress multiple
0012 # service messages coming from the system
0013 #
0014 # in this config below, we use the replace option to make
0015 # the logger let out messages of severity ERROR (INFO level
0016 # will be suppressed), and we want to limit the number to 10
0017 #
0018 process.load("FWCore.MessageService.MessageLogger_cfi")
0019
0020 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0021
0022 # Event output
0023 process.load("Configuration.EventContent.EventContent_cff")
0024
0025 process.maxEvents = cms.untracked.PSet(
0026     input = cms.untracked.int32(5)
0027 )
0028
0029 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0030     generator = cms.PSet(
0031         initialSeed = cms.untracked.uint32(123456789),
0032         engineName = cms.untracked.string('HepJamesRandom')
0033     )
0034 )
0035
0036
0037 process.source = cms.Source("EmptySource")
0038
0039 process.generator = cms.EDProducer("Pythia6PtGun",
0040     maxEventsToPrint = cms.untracked.int32(5),
0041     pythiaPylistVerbosity = cms.untracked.int32(1),
0042     pythiaHepMCVerbosity = cms.untracked.bool(True),
0043     PGunParameters = cms.PSet(
0044         ParticleID = cms.vint32(23),
0045         AddAntiParticle = cms.bool(False),
0046         MinPhi = cms.double(-3.14159265359),
0047         MaxPhi = cms.double(3.14159265359),
0048         MinPt = cms.double(50.0),
0049         MaxPt = cms.double(50.0001),
0050         MinEta = cms.double(-2.4),
0051         MaxEta = cms.double(2.4)
0052     ),
0053     PythiaParameters = cms.PSet(
0054         pythiaDefault = cms.vstring(
0055        'PMAS(5,1)=4.8 ! b quark mass',
0056            'PMAS(6,1)=172.3 ! t quark mass',
0057            'MSTP(61)=0 ! initial state radiation',
0058        'mstj(41)=1' # per Steve M., instead of mstp(71)... btw, shoult it be 0 or 1 ?
0059            #'MSTP(71)=0 !final state radiation'
0060     ),
0061     pythiaZtoMuons = cms.vstring(
0062        "MDME(174,1)=0",          # !Z decay into d dbar,
0063            "MDME(175,1)=0",          # !Z decay into u ubar,
0064            "MDME(176,1)=0",          # !Z decay into s sbar,
0065            "MDME(177,1)=0",          # !Z decay into c cbar,
0066            "MDME(178,1)=0",          # !Z decay into b bbar,
0067            "MDME(179,1)=0",          # !Z decay into t tbar,
0068            "MDME(182,1)=0",          # !Z decay into e- e+,
0069            "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0070            "MDME(184,1)=0",          # !Z decay into mu- mu+,
0071            "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0072            "MDME(186,1)=1",          # !Z decay into tau- tau+,
0073            "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0074         ),
0075     pythiaTauL = cms.vstring(
0076        "mdme(89,1)=1", # tau -> e
0077        "mdme(90,1)=1", # tau -> mu
0078        # all other tau decays OFF
0079        "mdme(91,1)=0",
0080        "mdme(92,1)=0",
0081        "mdme(93,1)=0",
0082        "mdme(94,1)=0",
0083        "mdme(95,1)=0",
0084        "mdme(96,1)=0",
0085        "mdme(97,1)=0",
0086        "mdme(98,1)=0",
0087        "mdme(99,1)=0",
0088        "mdme(100,1)=0",
0089        "mdme(101,1)=0",
0090        "mdme(102,1)=0",
0091        "mdme(103,1)=0",
0092        "mdme(104,1)=0",
0093        "mdme(105,1)=0",
0094        "mdme(106,1)=0",
0095        "mdme(107,1)=0",
0096        "mdme(108,1)=0",
0097        "mdme(109,1)=0",
0098        "mdme(110,1)=0",
0099        "mdme(111,1)=0",
0100        "mdme(112,1)=0",
0101        "mdme(113,1)=0",
0102        "mdme(114,1)=0",
0103        "mdme(115,1)=0",
0104        "mdme(116,1)=0",
0105        "mdme(117,1)=0",
0106        "mdme(118,1)=0",
0107        "mdme(119,1)=0",
0108        "mdme(120,1)=0",
0109        "mdme(121,1)=0",
0110        "mdme(122,1)=0",
0111        "mdme(123,1)=0",
0112        "mdme(124,1)=0",
0113        "mdme(125,1)=0",
0114        "mdme(126,1)=0",
0115        "mdme(127,1)=0",
0116        "mdme(128,1)=0",
0117        "mdme(129,1)=0",
0118        "mdme(130,1)=0",
0119        "mdme(131,1)=0",
0120        "mdme(132,1)=0",
0121        "mdme(133,1)=0",
0122        "mdme(134,1)=0",
0123        "mdme(135,1)=0",
0124        "mdme(136,1)=0",
0125        "mdme(137,1)=0",
0126        "mdme(138,1)=0",
0127        "mdme(139,1)=0",
0128        "mdme(140,1)=0",
0129        "mdme(141,1)=0",
0130        "mdme(142,1)=0"
0131     ),
0132         parameterSets = cms.vstring(
0133             'pythiaDefault',
0134             'pythiaZtoMuons', 'pythiaTauL'
0135         )
0136     )
0137 )
0138
0139 process.FEVT = cms.OutputModule("PoolOutputModule",
0140     process.FEVTSIMEventContent,
0141     fileName = cms.untracked.string('gen_singleTau.root')
0142 )
0143
0144 process.p = cms.Path(process.generator)
0145 process.outpath = cms.EndPath(process.FEVT)
0146 process.schedule = cms.Schedule(process.p,process.outpath)
0147
0148