Generator/python/BuToJPsiPrimeKToJPsiPiPiK_14TeV_TuneCP5_pythia8_cfi.py

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(14000.0),
0008                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0009                          maxEventsToPrint = cms.untracked.int32(0),
0010                          pythiaPylistVerbosity = cms.untracked.int32(0),
0011                          ExternalDecays = cms.PSet(
0012         EvtGen130 = cms.untracked.PSet(
0013             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0014             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0015             user_decay_embedded= cms.vstring(
0016 """## my decay
0017 Alias      MyB+    B+
0018 Alias      MyB-   B-
0019 ChargeConj MyB-   MyB+
0020 #
0021 Alias      MyJpsi J/psi
0022 ChargeConj MyJpsi MyJpsi
0023 #
0024 # Alias      Myrho0       rho0
0025 # ChargeConj Myrho0       Myrho0
0026 #
0027 #
0028 Alias       Mypsi(2S)  psi(2S)
0029 ChargeConj  Mypsi(2S)  Mypsi(2S)
0030 #
0031 Decay MyJpsi
0032 1.000  mu+       mu-                     PHOTOS VLL;
0033 Enddecay
0034 #
0035 Decay MyB+
0036 1.000     Mypsi(2S) K+         PHSP;
0037 Enddecay
0038 #
0039 Decay MyB-
0040 1.000     Mypsi(2S) K-         PHSP;
0041 Enddecay
0042 #
0043 Decay Mypsi(2S)
0044 1.000   MyJpsi  pi+   pi-         VVPIPI;
0045 Enddecay
0046 #
0047 # Decay Myrho0
0048 # 1.000   pi+   pi-                  PHSP;
0049 # Enddecay
0050 #
0051 End"""
0052      ),
0053             list_forced_decays = cms.vstring('MyB+','MyB-'),
0054             operates_on_particles = cms.vint32(),
0055             convertPythiaCodes = cms.untracked.bool(False)
0056             ),
0057         parameterSets = cms.vstring('EvtGen130')
0058         ),
0059         PythiaParameters = cms.PSet(
0060         pythia8CommonSettingsBlock,
0061         pythia8CP5SettingsBlock,
0062         processParameters = cms.vstring(
0063             "SoftQCD:nonDiffractive = on",
0064             'PTFilter:filter = on', # this turn on the filter
0065             'PTFilter:quarkToFilter = 5', # PDG id of q quark
0066             'PTFilter:scaleToFilter = 1.0'),
0067         parameterSets = cms.vstring(
0068             'pythia8CommonSettings',
0069             'pythia8CP5Settings',
0070             'processParameters',
0071         )
0072     )
0073 )
0074
0075
0076
0077 ###########
0078 # Filters #
0079 ###########
0080 bfilter = cms.EDFilter("PythiaFilter",
0081         MaxEta = cms.untracked.double(9999.),
0082         MinEta = cms.untracked.double(-9999.),
0083         ParticleID = cms.untracked.int32(521)
0084 )
0085
0086 jpsifilter = cms.EDFilter("PythiaDauVFilter",
0087         verbose         = cms.untracked.int32(0),
0088         NumberDaughters = cms.untracked.int32(2),
0089         MotherID        = cms.untracked.int32(100443),
0090         ParticleID      = cms.untracked.int32(443),
0091         DaughterIDs     = cms.untracked.vint32(13, -13),
0092         MinPt           = cms.untracked.vdouble(1.5, 1.5),
0093         MinEta          = cms.untracked.vdouble(-2.5, -2.5),
0094         MaxEta          = cms.untracked.vdouble( 2.5,  2.5)
0095 )
0096
0097 xxxfilter = cms.EDFilter("PythiaDauVFilter",
0098         verbose         = cms.untracked.int32(0),
0099         NumberDaughters = cms.untracked.int32(3),
0100         MotherID        = cms.untracked.int32(521),
0101         ParticleID      = cms.untracked.int32(100443),
0102         DaughterIDs     = cms.untracked.vint32(443, 211, -211),
0103         MinPt           = cms.untracked.vdouble(5, 0.0, 0.0),
0104         MinEta          = cms.untracked.vdouble(-9999, -9999, -9999),
0105         MaxEta          = cms.untracked.vdouble( 9999,  9999,  9999)
0106 )
0107
0108 decayfilter = cms.EDFilter("PythiaDauVFilter",
0109         verbose         = cms.untracked.int32(0),
0110         NumberDaughters = cms.untracked.int32(2),
0111         MotherID        = cms.untracked.int32(0),
0112         ParticleID      = cms.untracked.int32(521),
0113         DaughterIDs     = cms.untracked.vint32(100443, 321),
0114         MinPt           = cms.untracked.vdouble(5., 0.0),
0115         MinEta          = cms.untracked.vdouble(-9999., -9999.),
0116         MaxEta          = cms.untracked.vdouble( 9999.,  9999.)
0117 )
0118
0119 ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*xxxfilter*decayfilter)