Generator/python/LbToJpsiLambda_JMM_Filter_DGamma0_TuneCP5_13p6TeV-pythia8-evtgen_cfi.py

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