Generator/python/B0ToJpsiK0s_JMM_Filter_DGamma0_TuneCP5_13p6TeV-pythia8-evtgen_cfi.py

0001 # Found 139 output events for 5000 input events.
0002 # Filter efficiency = 0.0278
0003 # Timing = 0.684644 sec/event
0004 # Event size = 581.0 kB/event
0005
0006
0007 import FWCore.ParameterSet.Config as cms
0008 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0009 from Configuration.Generator.MCTunesRun3ECM13p6TeV.PythiaCP5Settings_cfi import *
0010 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0011
0012 _generator = cms.EDFilter("Pythia8GeneratorFilter",
0013                          comEnergy = cms.double(13600.0),
0014                          crossSection = cms.untracked.double(54000000000),
0015                          filterEfficiency = cms.untracked.double(3.0e-4),
0016                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0017                          maxEventsToPrint = cms.untracked.int32(0),
0018                          pythiaPylistVerbosity = cms.untracked.int32(0),
0019 #                          maxEventsToPrint = cms.untracked.int32(1),
0020 #                          pythiaPylistVerbosity = cms.untracked.int32(12),
0021                          ExternalDecays = cms.PSet(
0022         EvtGen130 = cms.untracked.PSet(
0023 ##            decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010_NOLONGLIFE.DEC'),
0024             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0025 #            particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0026             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0027             #user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Bs_Jpsiphi.dec'),
0028             user_decay_embedded= cms.vstring(
0029 '#',
0030 '# Particles updated from PDG2018 https://journals.aps.org/prd/abstract/10.1103/PhysRevD.98.030001',
0031 'Particle   pi+           1.3957061e-01   0.0000000e+00',
0032 'Particle   pi-           1.3957061e-01   0.0000000e+00',
0033 'Particle   mu+           1.0565837e-01   0.0000000e+00', ## id 13
0034 'Particle   mu-           1.0565837e-01   0.0000000e+00',
0035 'Particle   K+            4.9367700e-01   0.0000000e+00', ## id 321
0036 'Particle   K-            4.9367700e-01   0.0000000e+00',
0037 'Particle   p+            9.3827203e-01   0.0000000e+00', ## id 2212
0038 'Particle   anti-p-       9.3827203e-01   0.0000000e+00',
0039 'Particle   K_S0          4.9761100e-01   0.0000000e+00', ## id 310
0040 'Particle   K*+           8.9176000e-01   5.0300000e-02',
0041 'Particle   K*-           8.9176000e-01   5.0300000e-02',
0042 'Particle   K*0           8.9555000e-01   4.7300000e-02',
0043 'Particle   anti-K*0      8.9555000e-01   4.7300000e-02',
0044 'Particle   rho0          7.7526000e-01   1.4910000e-01',
0045 'Particle   phi           1.0194610e+00   4.2490000e-03',
0046 'Particle   Lambda0       1.1156830e+00   0.0000000e+00', ## id 3122
0047 'Particle   anti-Lambda0  1.1156830e+00   0.0000000e+00',
0048 'Particle   Sigma0        1.1926420e+00   8.8947595e-06', ## id 3212
0049 'Particle   anti-Sigma0   1.1926420e+00   8.8947595e-06',
0050 'Particle   B-            5.2793200e+00   0.0000000e+00',
0051 'Particle   B+            5.2793200e+00   0.0000000e+00',
0052 'Particle   B0            5.2796300e+00   0.0000000e+00', ## id 511
0053 'Particle   anti-B0       5.2796300e+00   0.0000000e+00',
0054 'Particle   B_s0          5.3668900e+00   0.0000000e+00',
0055 'Particle   anti-B_s0     5.3668900e+00   0.0000000e+00',
0056 'Particle   J/psi         3.0969000e+00   9.2900006e-05', ## id 443
0057 'Particle   psi(2S)       3.6860970e+00   2.9400000e-04', ## id 100443
0058 'Particle   Lambda_b0     5.6196000e+00   0.0000000e-04', ## id 5122,
0059 'Particle anti-Lambda_b0  5.6196000e+00   0.0000000e-04',
0060 '#',
0061 '#',
0062 '#',
0063 'Alias       MyB        B0',
0064 'Alias       Myanti-B   anti-B0',
0065 'ChargeConj  Myanti-B   MyB',
0066 '#',
0067 'Alias       Mypsi      J/psi',
0068 'ChargeConj  Mypsi      Mypsi',
0069 '#',
0070 'Alias       MyK0s      K_S0',
0071 'ChargeConj  MyK0s      MyK0s',
0072 '#',
0073 'Decay Mypsi',
0074 '1.000       mu+    mu-        PHOTOS VLL;',
0075 'Enddecay',
0076 '#',
0077 'Decay MyB',
0078 '1.000       Mypsi  MyK0s      PHSP;',
0079 'Enddecay',
0080 'CDecay Myanti-B',
0081 'End'
0082 ),
0083             list_forced_decays = cms.vstring('MyB','Myanti-B'),
0084             operates_on_particles = cms.vint32(),
0085             convertPythiaCodes = cms.untracked.bool(False)
0086             ),
0087         parameterSets = cms.vstring('EvtGen130')
0088         ),
0089         PythiaParameters = cms.PSet(
0090         pythia8CommonSettingsBlock,
0091         pythia8CP5SettingsBlock,
0092         processParameters = cms.vstring(
0093             "SoftQCD:nonDiffractive = on",
0094             "511:m0=5.27963",       ## changing also lambda_b0 mass in pythia
0095             'PTFilter:filter = on', # this turn on the filter
0096             'PTFilter:quarkToFilter = 5', # PDG id of q quark
0097             'PTFilter:scaleToFilter = 1.0'),
0098         parameterSets = cms.vstring(
0099             'pythia8CommonSettings',
0100             'pythia8CP5Settings',
0101             'processParameters',
0102         )
0103     )
0104 )
0105 _generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0106
0107 from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
0108 generator = ExternalGeneratorFilter(_generator)
0109
0110 ###########
0111 # Filters #
0112 ###########
0113
0114 #
0115 bfilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(511))
0116
0117 # psifilter = cms.EDFilter("PythiaFilter",
0118 #         MotherID        = cms.untracked.int32(5122),
0119 #         ParticleID      = cms.untracked.int32(443)
0120 # #         MinPt           = cms.untracked.double(4.95),
0121 # #         MinEta          = cms.untracked.double(-3.0),
0122 # #         MaxEta          = cms.untracked.double( 3.0)
0123 # )
0124 psifilter = cms.EDFilter("PythiaDauVFilter",
0125         verbose         = cms.untracked.int32(0),
0126         NumberDaughters = cms.untracked.int32(2),
0127         MotherID        = cms.untracked.int32(511),
0128         ParticleID      = cms.untracked.int32(443),
0129         DaughterIDs     = cms.untracked.vint32(13, -13),
0130         MinPt           = cms.untracked.vdouble(3., 3.),
0131         MinEta          = cms.untracked.vdouble(-2.5, -2.5),
0132         MaxEta          = cms.untracked.vdouble(2.5, 2.5)
0133 )
0134 decayfilter = cms.EDFilter("PythiaDauVFilter",
0135         verbose         = cms.untracked.int32(0),
0136         NumberDaughters = cms.untracked.int32(2),
0137         MotherID        = cms.untracked.int32(0),
0138         ParticleID      = cms.untracked.int32(511),
0139         DaughterIDs     = cms.untracked.vint32(443, 310),
0140         MinPt           = cms.untracked.vdouble(5, 0.5),
0141         MinEta          = cms.untracked.vdouble(-99999, -3),
0142         MaxEta          = cms.untracked.vdouble( 99999,  3)
0143 )
0144
0145
0146 # ProductionFilterSequence = cms.Sequence(generator*lbfilter*psifilter)
0147 ProductionFilterSequence = cms.Sequence(generator*bfilter*decayfilter*psifilter)