Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:13:58

0001 # The configuration for comparison of Pythia8 tau decays with Taula tau decays (Py8_Z2tautau_tauola_cfg.py)
0002 # hadronization, MPI settings, PDF settings to be eventually updated in both
0003 import FWCore.ParameterSet.Config as cms
0004 
0005 process = cms.Process("TEST")
0006 
0007 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0008 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0009 
0010 process.source = cms.Source("EmptySource",
0011      firstLuminosityBlock = cms.untracked.uint32(1),
0012      numberEventsInLuminosityBlock = cms.untracked.uint32(200),
0013      firstEvent = cms.untracked.uint32(1),
0014      firstRun = cms.untracked.uint32(1),
0015      numberEventsInRun = cms.untracked.uint32(1000)
0016 
0017 )
0018 
0019 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0020     maxEventsToPrint = cms.untracked.int32(1),
0021     pythiaPylistVerbosity = cms.untracked.int32(1),
0022     filterEfficiency = cms.untracked.double(1.0),
0023     pythiaHepMCVerbosity = cms.untracked.bool(True),
0024     comEnergy = cms.double(10000.),
0025         
0026     PythiaParameters = cms.PSet(
0027         py8UESettings = cms.vstring( 'ParticleDecays:limitTau0 = on',  # - decay unstable particles
0028                      'ParticleDecays:tau0Max   = 10.', # - for which ctau < 10 mm
0029                      'HadronLevel:Hadronize = on',
0030                      'PartonLevel:MPI = on'
0031     ),
0032         py8ProcessSettings = cms.vstring( 'WeakSingleBoson:ffbar2gmZ = on'
0033     ),
0034     py8ZDecaySettings = cms.vstring(  '23:onMode = off', # turn OFF all Z decays
0035                                       #'23:onIfAny = 13', # turn ON Z->mumu
0036                       '23:onIfAny = 15'  # turn ON Z->tautau
0037                       #'23:onIfAny = 5'
0038     ),
0039     py8TauDecaySettings = cms.vstring('15:onMode = off', # turn OFF all tau decays
0040                                       '15:onIfAny = 11', # turn ON tau -> e
0041                       '15:onIfAny = 13'  # turn ON tau -> mu
0042     ),
0043         parameterSets = cms.vstring(   'py8UESettings',
0044                                 'py8ProcessSettings',
0045                     'py8ZDecaySettings')
0046     )
0047 )
0048 
0049 # The following three lines reduce the clutter of repeated printouts
0050 # of the same exception message.
0051 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0052 process.MessageLogger = cms.Service("MessageLogger",
0053     cerr = cms.untracked.PSet(
0054         enable = cms.untracked.bool(False)
0055     ),
0056     cout = cms.untracked.PSet(
0057         default = cms.untracked.PSet(
0058             limit = cms.untracked.int32(2)
0059         ),
0060         enable = cms.untracked.bool(True)
0061     )
0062 )
0063 
0064 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0065     generator = cms.PSet(
0066         initialSeed = cms.untracked.uint32(123456789),
0067         engineName = cms.untracked.string('HepJamesRandom')
0068     )
0069 )
0070 
0071 process.maxEvents = cms.untracked.PSet(
0072     input = cms.untracked.int32(1000)
0073 )
0074 
0075 process.GEN = cms.OutputModule("PoolOutputModule",
0076     fileName = cms.untracked.string('Py8_Z2tautau.root')
0077 )
0078 
0079 #
0080 # needed for the BasicGenTester (analyzer)
0081 #
0082 #process.TFileService = cms.Service("TFileService",
0083 #        fileName = cms.string("BasicGenTester_Py8_Ztautau.root")
0084 #)
0085 # the analyzer itself  
0086 #
0087 #process.BasicGenTest = cms.EDAnalyzer( "BasicGenTester",
0088 #        NPartForHisto = cms.untracked.int32(1000),
0089 #        PtMaxForHisto = cms.untracked.double(5.0)
0090 #)
0091 #process.p1 = cms.Path( process.BasicGenTest )
0092 
0093 process.p = cms.Path(process.generator)
0094 process.outpath = cms.EndPath(process.GEN)
0095 
0096 process.schedule = cms.Schedule(process.p, process.outpath)
0097 #process.schedule = cms.Schedule(process.p, process.p1)