Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:52:37

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0003 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0004 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0005 
0006 generator = cms.EDFilter("Pythia8GeneratorFilter",
0007                          comEnergy = cms.double(13000.0),
0008                          crossSection = cms.untracked.double(54000000000),
0009                          filterEfficiency = cms.untracked.double(3.0e-4),
0010                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0011                          maxEventsToPrint = cms.untracked.int32(0),
0012                          pythiaPylistVerbosity = cms.untracked.int32(0),
0013                          ExternalDecays = cms.PSet(
0014                          #using alternative name for decayer
0015                          EvtGen1 = cms.untracked.PSet(
0016                         #uses latest evt and decay tables from evtgen 
0017                             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_NOLONGLIFE.DEC'),
0018                             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0019                             convertPythiaCodes = cms.untracked.bool(False),
0020                             #here we will use the user.dec store in the release
0021                             user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Bs_mumu.dec'),
0022                             list_forced_decays = cms.vstring('MyB_s0','Myanti-B_s0'),
0023                             operates_on_particles = cms.vint32()
0024                          ),
0025                          parameterSets = cms.vstring('EvtGen1')
0026                         ),
0027                         PythiaParameters = cms.PSet(
0028                            pythia8CommonSettingsBlock,
0029                            pythia8CP5SettingsBlock,
0030                            processParameters = cms.vstring(
0031                               #filter of a b-quark before hadronizing, and use a better data-like process
0032                               'PTFilter:filter = on',
0033                               'PTFilter:quarkToFilter = 5',
0034                               'PTFilter:scaleToFilter = 1.0',
0035                               'SoftQCD:nonDiffractive = on',
0036                            ),
0037                            parameterSets = cms.vstring('pythia8CommonSettings',
0038                                     'pythia8CP5Settings',
0039                                     'processParameters',
0040                            )
0041                         )
0042 )
0043 
0044 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0045 
0046 MuMuFilter = cms.EDFilter("MCParticlePairFilter",
0047                           Status = cms.untracked.vint32(1, 1),
0048                           MinPt = cms.untracked.vdouble(3., 3.),
0049                           MaxEta = cms.untracked.vdouble(2.5, 2.5),
0050                           MinEta = cms.untracked.vdouble(-2.5, -2.5),
0051                           ParticleCharge = cms.untracked.int32(-1),
0052                           ParticleID1 = cms.untracked.vint32(13,-13),
0053                           )
0054 
0055 # -- Require Muon from Bs
0056 MuFilter = cms.EDFilter("PythiaFilter",
0057                         Status = cms.untracked.int32(1),
0058                         MotherID = cms.untracked.int32(531),
0059                         MinPt = cms.untracked.double(3.),
0060                         ParticleID = cms.untracked.int32(13),
0061                         MaxEta = cms.untracked.double(2.5),
0062                         MinEta = cms.untracked.double(-2.5)
0063                         )
0064 
0065 ProductionFilterSequence = cms.Sequence(generator*MuMuFilter*MuFilter)