ExternalDecays/test/Py6GenFilter_W2TauNu_Tauola_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 process.randomEngineStateProducer = cms.EDProducer("RandomEngineStateProducer")
0022
0023 process.Timing=cms.Service("Timing",
0024     summaryOnly=cms.untracked.bool(True))
0025
0026
0027 # The following three lines reduce the clutter of repeated printouts
0028 # of the same exception message.
0029 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0030
0031 process.MessageLogger.cerr.enableStatistics = False
0032
0033
0034 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(5))
0035
0036 process.source = cms.Source("EmptySource")
0037
0038 from Configuration.Generator.PythiaUESettings_cfi import *
0039 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0040     pythiaHepMCVerbosity = cms.untracked.bool(True),
0041     maxEventsToPrint = cms.untracked.int32(0),
0042     pythiaPylistVerbosity = cms.untracked.int32(1),
0043     # this shows how to turn ON some of the general Py6 printouts, like banner...
0044     ## --> displayPythiaBanner = cms.untracked.bool(True),
0045     ## --> displayPythiaCards = cms.untracked.bool(True),
0046     comEnergy = cms.double(7000.0),
0047
0048     ExternalDecays = cms.PSet(
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(0) # generic tau decays
0061                 ## mdtau = cms.int32(240) # (any) tau -> nu pi+-
0062              )
0063 #           TauolaDefaultInputCards,
0064 #          TauolaPolar
0065         ),
0066         parameterSets = cms.vstring('Tauola')
0067     ),
0068     PythiaParameters = cms.PSet(
0069         pythiaUESettingsBlock,
0070         processParameters = cms.vstring('MSEL        = 0    !User defined processes',
0071                                         'MSUB(2)     = 1    !W production',
0072                                         'MDME(190,1) = 0    !W decay into dbar u',
0073                                         'MDME(191,1) = 0    !W decay into dbar c',
0074                                         'MDME(192,1) = 0    !W decay into dbar t',
0075                                         'MDME(194,1) = 0    !W decay into sbar u',
0076                                         'MDME(195,1) = 0    !W decay into sbar c',
0077                                         'MDME(196,1) = 0    !W decay into sbar t',
0078                                         'MDME(198,1) = 0    !W decay into bbar u',
0079                                         'MDME(199,1) = 0    !W decay into bbar c',
0080                                         'MDME(200,1) = 0    !W decay into bbar t',
0081                                         'MDME(205,1) = 0    !W decay into bbar tp',
0082                                         'MDME(206,1) = 0    !W decay into e+ nu_e',
0083                                         'MDME(207,1) = 0    !W decay into mu+ nu_mu',
0084                                         'MDME(208,1) = 1    !W decay into tau+ nu_tau'),
0085         # This is a vector of ParameterSet names to be read, in this order
0086         parameterSets = cms.vstring('pythiaUESettings',
0087                                     'processParameters')
0088     )
0089 )
0090
0091 process.GEN = cms.OutputModule("PoolOutputModule",
0092     fileName = cms.untracked.string('Test_Py6_W2TauNu_Tauola.root')
0093 )
0094
0095 process.p = cms.Path(process.generator)
0096
0097 process.outpath = cms.EndPath(process.GEN)
0098 process.p1 = cms.Path(process.randomEngineStateProducer)
0099 # process.schedule = cms.Schedule(process.p, process.p1, process.outpath)
0100 process.schedule = cms.Schedule(process.p)