ExternalDecays/test/TauSpinnerWTest.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 process.load("GeneratorInterface.ExternalDecays.TauSpinner_cfi")
0015
0016 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0017                                                    generator = cms.PSet(initialSeed = cms.untracked.uint32(123456789),
0018                                                                         engineName = cms.untracked.string('HepJamesRandom')
0019                                                                         ),
0020                                                    TauSpinnerGen  = cms.PSet(initialSeed = cms.untracked.uint32(123456789),
0021                                                                              engineName = cms.untracked.string('HepJamesRandom')
0022
0023                                                                              )
0024                                                    )
0025 process.randomEngineStateProducer = cms.EDProducer("RandomEngineStateProducer")
0026
0027 process.Timing=cms.Service("Timing",
0028                            summaryOnly=cms.untracked.bool(True))
0029
0030
0031 # The following three lines reduce the clutter of repeated printouts
0032 # of the same exception message.
0033 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0034
0035 process.MessageLogger.cerr.enableStatistics = False
0036
0037
0038 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(100))
0039
0040 process.source = cms.Source("EmptySource")
0041
0042 from Configuration.Generator.PythiaUESettings_cfi import *
0043 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0044                                  pythiaHepMCVerbosity = cms.untracked.bool(True),
0045                                  maxEventsToPrint = cms.untracked.int32(0),
0046                                  pythiaPylistVerbosity = cms.untracked.int32(1),
0047                                      # this shows how to turn ON some of the general Py6 printouts, like banner...
0048                                  ## --> displayPythiaBanner = cms.untracked.bool(True),
0049                                  ## --> displayPythiaCards = cms.untracked.bool(True),
0050                                  comEnergy = cms.double(7000.0),
0051
0052                                  ExternalDecays = cms.PSet(
0053     Tauola = cms.untracked.PSet(
0054     # these settings below exemplfy how to use "native" Tauola approach:
0055     # one MUST set mdtau=1 !!! then pjak1 & pjak2 will translate into
0056     # native variables jak1 & jak2 (jak1/jak2=4 means that both tau's
0057     # decay into the rho-mode
0058     #
0059     UseTauolaPolarization = cms.bool(True),
0060     InputCards = cms.PSet
0061     (
0062     pjak1 = cms.int32(0),
0063     pjak2 = cms.int32(0),
0064     mdtau = cms.int32(0) # generic tau decays
0065     ## mdtau = cms.int32(240) # (any) tau -> nu pi+-
0066     )
0067     #           TauolaDefaultInputCards,
0068     #          TauolaPolar
0069     ),
0070     parameterSets = cms.vstring('Tauola')
0071     ),
0072                                  PythiaParameters = cms.PSet(
0073     pythiaUESettingsBlock,
0074     processParameters = cms.vstring('MSEL        = 0    !User defined processes',
0075                                     'MSUB(2)     = 1    !W production',
0076                                     'MDME(190,1) = 0    !W decay into dbar u',
0077                                     'MDME(191,1) = 0    !W decay into dbar c',
0078                                     'MDME(192,1) = 0    !W decay into dbar t',
0079                                     'MDME(194,1) = 0    !W decay into sbar u',
0080                                     'MDME(195,1) = 0    !W decay into sbar c',
0081                                     'MDME(196,1) = 0    !W decay into sbar t',
0082                                     'MDME(198,1) = 0    !W decay into bbar u',
0083                                     'MDME(199,1) = 0    !W decay into bbar c',
0084                                     'MDME(200,1) = 0    !W decay into bbar t',
0085                                     'MDME(205,1) = 0    !W decay into bbar tp',
0086                                     'MDME(206,1) = 0    !W decay into e+ nu_e',
0087                                     'MDME(207,1) = 0    !W decay into mu+ nu_mu',
0088                                     'MDME(208,1) = 1    !W decay into tau+ nu_tau'),
0089     # This is a vector of ParameterSet names to be read, in this order
0090     parameterSets = cms.vstring('pythiaUESettings',
0091                                 'processParameters')
0092     )
0093                                  )
0094
0095 # Produce PDF weights (maximum is 3)
0096 process.pdfWeights = cms.EDProducer("PdfWeightProducer",
0097                                     # Fix POWHEG if buggy (this PDF set will also appear on output,
0098                                     # so only two more PDF sets can be added in PdfSetNames if not "")
0099                                     #FixPOWHEG = cms.untracked.string("cteq66.LHgrid"),
0100                                     #GenTag = cms.untracked.InputTag("genParticles"),
0101                                     PdfInfoTag = cms.untracked.InputTag("VtxSmeared"),
0102                                     PdfSetNames = cms.untracked.vstring(
0103     #    "cteq66.LHgrid"
0104     #    , "MRST2006nnlo.LHgrid" ,
0105         "MSTW2008nnlo90cl.LHgrid"
0106         )
0107                                     )
0108
0109
0110 process.p1 = cms.Path( process.TauSpinnerGen )
0111
0112 process.GEN = cms.OutputModule("PoolOutputModule",
0113                                fileName = cms.untracked.string('Test_Py6_W2TauNu_Tauola.root')
0114                                )
0115
0116 process.p = cms.Path(process.generator)
0117
0118 process.outpath = cms.EndPath(process.GEN)
0119 process.p1 = cms.Path(process.randomEngineStateProducer*process.TauSpinnerGen)
0120 process.schedule = cms.Schedule(process.p, process.p1, process.outpath)
0121