Generator/python/UpsMM_8TeV_TuneCUETP8M1_cfi.py

0001 # fragment from https://github.com/alberto-sanchez/my-genproductions/blob/master/py8_UpsilonMM_FSR_13TeV_cfi.py
0002 # Updated to CUEP8M1 tune by Ian M. Nugent
0003 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0004 from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
0005 import FWCore.ParameterSet.Config as cms
0006 generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
0007                          pythiaPylistVerbosity = cms.untracked.int32(0),
0008                          filterEfficiency = cms.untracked.double(0.0757),
0009                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0010                          comEnergy = cms.double(8000.0),
0011                          crossSection = cms.untracked.double(13775390),
0012                          maxEventsToPrint = cms.untracked.int32(0),
0013                          PythiaParameters = cms.PSet(
0014         pythia8CommonSettingsBlock,
0015         pythia8CUEP8M1SettingsBlock,
0016         processParameters = cms.vstring(
0017             'Bottomonium:all = on', # Quarkonia, MSEL=62
0018             'ParticleDecays:allowPhotonRadiation = on', # Turn on QED FSR
0019             'StringFlav:mesonBvector = 3.000', # relative production vector/pseudoscalar for charm mesons - needs work
0020             'StringFlav:mesonBL1S1J0 = 0.072', # relative scalar production (L=1,S=1,J=0)/pseudoscalar for charm mesons
0021             'StringFlav:mesonBL1S0J1 = 3.000', # relative pseudovector production (L=1,S=0,J=1)/pseudoscalar for charm mesons
0022             'StringFlav:mesonBL1S1J1 = 0.216', # relative pseudovector production (L=1,S=1,J=1)/pseudoscalar for charm mesons
0023             'StringFlav:mesonBL1S1J2 = 0.000', # relative tensor production (L=1,S=1,J=2)/pseudoscalar for charm mesons
0024             '553:onMode = off', # Turn off Upsilon decays
0025             '553:onIfMatch = 13 -13' # just let Upsilon -> mu+ mu-
0026             ),
0027         parameterSets = cms.vstring('pythia8CommonSettings',
0028                                     'pythia8CUEP8M1Settings',
0029                                     'processParameters',
0030                                     )
0031         )
0032                          )
0033
0034 oniafilter = cms.EDFilter("PythiaFilter",
0035                           Status = cms.untracked.int32(2),
0036                           MaxEta = cms.untracked.double(1000.0),
0037                           MinEta = cms.untracked.double(-1000.0),
0038                           MinPt = cms.untracked.double(0.0),
0039                           ParticleID = cms.untracked.int32(553)
0040 )
0041
0042 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
0043                              Status = cms.untracked.vint32(1, 1),
0044                              MinPt = cms.untracked.vdouble(0.5, 0.5),
0045                              MinP = cms.untracked.vdouble(2.7, 2.7),
0046                              MaxEta = cms.untracked.vdouble(2.5, 2.5),
0047                              MinEta = cms.untracked.vdouble(-2.5, -2.5),
0048                              MinInvMass = cms.untracked.double(5.0),
0049                              MaxInvMass = cms.untracked.double(20.0),
0050                              ParticleCharge = cms.untracked.int32(-1),
0051                              ParticleID1 = cms.untracked.vint32(13),
0052                              ParticleID2 = cms.untracked.vint32(13)
0053                              )
0054
0055 ProductionFilterSequence = cms.Sequence(generator*oniafilter*mumugenfilter)
0056
0057