ExternalDecays/test/Py8_Z2tautau_tauolacxx_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004
0005 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0006 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0007
0008 process.source = cms.Source("EmptySource")
0009
0010 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0011     maxEventsToPrint = cms.untracked.int32(0),
0012     pythiaPylistVerbosity = cms.untracked.int32(1),
0013     filterEfficiency = cms.untracked.double(1.0),
0014     pythiaHepMCVerbosity = cms.untracked.bool(True),
0015     comEnergy = cms.double(7000.),
0016
0017     ExternalDecays = cms.PSet(
0018         Tauola = cms.untracked.PSet(
0019          UseTauolaPolarization = cms.bool(True),
0020          InputCards = cms.PSet
0021          (
0022             pjak1 = cms.int32(1), # electron mode
0023         pjak2 = cms.int32(2), # muon mode
0024         mdtau = cms.int32(0)  # does NOT matter at this point...
0025          )
0026     ),
0027         parameterSets = cms.vstring('Tauola')
0028     ),
0029
0030     PythiaParameters = cms.PSet(
0031         py8GeneralSettings = cms.vstring(  'Main:timesAllowErrors    = 10000',
0032                                            'ParticleDecays:limitTau0 = on',
0033                                        'ParticleDecays:tauMax = 10'
0034         ),
0035         py8UESettings = cms.vstring(       'Tune:pp 2',
0036                                        'Tune:ee 3'
0037     ),
0038         py8SpecialSettings   = cms.vstring('HadronLevel:Hadronize = off',     # mstp(111)=1
0039                            #' ParticleDecays:limitTau0 = on',  # mstj(22)=2 - decay unstable particles
0040                            # 'ParticleDecays:tau0Max   = 10.', # parj(71)=10.- for which ctau < 10 mm
0041                        'PartonLevel:ISR = off',
0042                        'PartonLevel:FSR = off',
0043                        'PartonLevel:MPI = off'           # mstp(81)=0
0044                        # for pp intractions, SpaceShower things are not important,
0045                        # but TimeShower settings are relevat
0046                        'TimeShower:QEDshowerByL = off',
0047                        'TimeShower:QEDshowerByQ = off',
0048                        'SpaceShower:QEDshowerByL = off',
0049                        'SpaceShower:QEDshowerByQ = off',
0050                        'PhaseSpace:pTHatMin = 5.',
0051                        'PhaseSpace:pTHatMax = 150.',
0052                        'StringZ:usePetersonB = on', # these 2 together ==
0053                        'StringZ:usePetersonC = on'  # mstj(11)=3
0054     ),
0055         py8ProcessSettings = cms.vstring(  # 'HiggsSM:gg2H = on',
0056                        # 'HiggsSM:ff2Hff(t:ZZ) = on',
0057                        # 'HiggsSM:ff2Hff(t:WW) = on'
0058                        # 'WeakDoubleBoson:ffbar2ZW = on'
0059                        'WeakSingleBoson:ffbar2gmZ = on' # msel=11
0060     ),
0061     py8ZDecaySettings = cms.vstring(   # '24:onMode = off', '24:onIfAny = 15',
0062                        # '25:onMode = off', '25:onIfAny = 23',
0063                                        '23:onMode = off', # turn OFF all Z decays
0064                                        # '23:onIfAny = 13', # turn ON Z->mumu
0065                        '23:onIfAny = 15'  # turn ON Z->tautau
0066     ),
0067         parameterSets = cms.vstring('py8GeneralSettings',
0068                                 'py8UESettings',
0069                                 'py8ProcessSettings',
0070                     'py8ZDecaySettings'
0071                     )
0072     )
0073 )
0074
0075 # The following three lines reduce the clutter of repeated printouts
0076 # of the same exception message.
0077 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0078
0079 process.MessageLogger.cerr.enableStatistics = False
0080
0081
0082 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0083     generator = cms.PSet(
0084         initialSeed = cms.untracked.uint32(123456789),
0085         engineName = cms.untracked.string('HepJamesRandom')
0086     )
0087 )
0088
0089 process.maxEvents = cms.untracked.PSet(
0090     input = cms.untracked.int32(20)
0091 )
0092
0093 process.GEN = cms.OutputModule("PoolOutputModule",
0094     fileName = cms.untracked.string('Py8_Z2tautau.root')
0095 )
0096
0097 process.p = cms.Path(process.generator)
0098 process.outpath = cms.EndPath(process.GEN)
0099
0100 process.schedule = cms.Schedule(process.p, process.outpath)
0101 #process.schedule = cms.Schedule(process.p, process.p1)
0102