PhotosInterface/test/DYtauola_Tune4C_13TeV_pythia8.py

0001 # A script to produce Z (Z -> tau tau) events using pythia8 and with no FSR
0002 # To decay tau tauola++ is called. tauola++ in turn calls photos++ internally for QED corrections
0003
0004 import FWCore.ParameterSet.Config as cms
0005
0006 process = cms.Process('VALIDATION')
0007
0008 # import of standard configurations
0009 process.load('Configuration.StandardSequences.Services_cff')
0010 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0011 process.load('FWCore.MessageService.MessageLogger_cfi')
0012 process.load('Configuration.EventContent.EventContent_cff')
0013 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0014 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0015 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0016 process.load('Configuration.StandardSequences.Generator_cff')
0017 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
0018 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0019 process.load('Configuration.StandardSequences.Validation_cff')
0020 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0021 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0022
0023 process.maxEvents = cms.untracked.PSet(
0024     input = cms.untracked.int32(3000)
0025 )
0026
0027 # Input source
0028 process.source = cms.Source("EmptySource")
0029
0030 process.options = cms.untracked.PSet(
0031
0032 )
0033
0034 # Production Info
0035 process.configurationMetadata = cms.untracked.PSet(
0036     version = cms.untracked.string('$Revision: 1.19 $'),
0037     annotation = cms.untracked.string('Configuration/Generator/python/ThirteenTeV/WToMuNu_M_500_Tune4C_13TeV_pythia8_cfi.py nevts:100000'),
0038     name = cms.untracked.string('Applications')
0039 )
0040
0041 # Output definition
0042
0043 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0044     splitLevel = cms.untracked.int32(0),
0045     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0046     outputCommands = process.RAWSIMEventContent.outputCommands,
0047     fileName = cms.untracked.string('DYnorad_Tune4C_13TeV_pythia8_cfi_py_GEN_VALIDATION.root'),
0048     dataset = cms.untracked.PSet(
0049         filterName = cms.untracked.string(''),
0050         dataTier = cms.untracked.string('GEN')
0051     ),
0052     SelectEvents = cms.untracked.PSet(
0053         SelectEvents = cms.vstring('generation_step')
0054     )
0055 )
0056
0057 # Additional output definition
0058
0059 # Other statements
0060 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0061 process.mix.playback = True
0062 process.mix.digitizers = cms.PSet()
0063 for a in process.aliases: delattr(process, a)
0064 from Configuration.AlCa.GlobalTag import GlobalTag
0065 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:startup', '')
0066
0067 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0068     ExternalDecays = cms.PSet(
0069         Tauola = cms.untracked.PSet(
0070             UseTauolaPolarization = cms.bool(True),
0071             InputCards = cms.PSet(
0072                 mdtau = cms.int32(0),
0073                 pjak2 = cms.int32(3),
0074                 pjak1 = cms.int32(3)
0075             )
0076         ),
0077         parameterSets = cms.vstring('Tauola')
0078     ),
0079     UseExternalGenerators = cms.untracked.bool(True),
0080     pythiaPylistVerbosity = cms.untracked.int32(1),
0081     filterEfficiency = cms.untracked.double(1),
0082     pythiaHepMCVerbosity = cms.untracked.bool(False),
0083     comEnergy = cms.double(13000.0),
0084     crossSection = cms.untracked.double(0.2188),
0085     maxEventsToPrint = cms.untracked.int32(1),
0086     PythiaParameters = cms.PSet(
0087         processParameters = cms.vstring('Main:timesAllowErrors = 10000',
0088                     'ParticleDecays:limitTau0 = on',
0089                     'ParticleDecays:tauMax = 10',
0090                     'Tune:ee 3',
0091                     'Tune:pp 5',
0092                     'PartonLevel:FSR = off',
0093                     'WeakSingleBoson:ffbar2gmZ = on',
0094                     '23:onMode = off',
0095                     '23:onIfAny = 15',
0096                     '23:mMin = 50.',
0097                     '23:mMax = 120.'
0098                     ),
0099         parameterSets = cms.vstring('processParameters')
0100     )
0101 )
0102
0103 process.ProductionFilterSequence = cms.Sequence(process.generator)
0104
0105 # Path and EndPath definitions
0106 process.generation_step = cms.Path(process.pgen)
0107 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0108 process.validation_step = cms.EndPath(process.genstepfilter+process.genvalid_all)
0109 process.endjob_step = cms.EndPath(process.endOfProcess)
0110 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0111
0112 # Schedule definition
0113 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.validation_step,process.endjob_step,process.RAWSIMoutput_step)
0114 # filter all path with the production filter sequence
0115 for path in process.paths:
0116     getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq
0117
0118 # customisation of the process.
0119
0120 # Automatic addition of the customisation function from SimGeneral.MixingModule.fullMixCustomize_cff
0121 from SimGeneral.MixingModule.fullMixCustomize_cff import setCrossingFrameOn
0122
0123 #call to customisation function setCrossingFrameOn imported from SimGeneral.MixingModule.fullMixCustomize_cff
0124 process = setCrossingFrameOn(process)
0125
0126 # End of customisation functions