Generator/python/LbToJpsiXiK0sPi_JMM_Filter_DGamma0_TuneCP5_13p6TeV-pythia8-evtgen_cfi.py

0001 # Found 48 output events for 8000 input events.
0002 # Filter efficiency = 0.006
0003 # Timing = 0.260728 sec/event
0004 # Event size = 546.5 kB/event
0005
0006 import FWCore.ParameterSet.Config as cms
0007 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0008 from Configuration.Generator.MCTunesRun3ECM13p6TeV.PythiaCP5Settings_cfi import *
0009 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0010
0011 _generator = cms.EDFilter("Pythia8GeneratorFilter",
0012                          comEnergy = cms.double(13600.0),
0013                          crossSection = cms.untracked.double(54000000000),
0014                          filterEfficiency = cms.untracked.double(3.0e-4),
0015                          #pythiaHepMCVerbosity = cms.untracked.bool(False),
0016                          maxEventsToPrint = cms.untracked.int32(0),
0017                          pythiaPylistVerbosity = cms.untracked.int32(0),
0018                          ExternalDecays = cms.PSet(
0019         EvtGen130 = cms.untracked.PSet(
0020             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0021             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0022             user_decay_embedded= cms.vstring(
0023 '#',
0024 '# Particles updated from PDG2020 published in Prog. Theor. Exp. Phys. 2020, 083C01 (2020)',
0025 'Particle   pi+           1.3957039e-01   0.0000000e+00',
0026 'Particle   pi-           1.3957039e-01   0.0000000e+00',
0027 'Particle   Xi-           1.3217100e+00   0.0000000e+00', ## id 3312
0028 'Particle   anti-Xi+      1.3217100e+00   0.0000000e+00', ## id -3312
0029 'Particle   Lambda0       1.1156830e+00   0.0000000e+00', ## id 3122
0030 'Particle   anti-Lambda0  1.1156830e+00   0.0000000e+00', ## id -3122
0031 'Particle   Lambda_b0     5.6196000e+00   0.0000000e+00', ## id 5122
0032 'Particle   anti-Lambda_b0 5.6196000e+00   0.0000000e+00', ## id -5122
0033 'Particle   K+            4.9367700e-01   0.0000000e+00', ## id 321
0034 'Particle   K-            4.9367700e-01   0.0000000e+00',
0035 'Particle   K_S0          4.9761100e-01   0.0000000e+00', ## id 310
0036 'Particle   J/psi         3.0969000e+00   9.2900000e-05',
0037 '#',
0038 'Alias      MyLambda_b0    Lambda_b0',
0039 'Alias      Myanti-Lambda_b0   anti-Lambda_b0',
0040 'ChargeConj Myanti-Lambda_b0   MyLambda_b0',
0041 'Alias      MyJ/psi  J/psi',
0042 'ChargeConj MyJ/psi  MyJ/psi',
0043 '#',
0044 'Alias       MyK0s      K_S0',
0045 'ChargeConj  MyK0s      MyK0s',
0046 '#',
0047 'Decay MyLambda_b0',
0048 '1.000      MyJ/psi     Xi-         MyK0s    pi+  PHSP;',
0049 'Enddecay',
0050 'Decay Myanti-Lambda_b0',
0051 '1.000      MyJ/psi     anti-Xi+    MyK0s    pi-  PHSP;',
0052 'Enddecay',
0053 '#',
0054 'Decay MyJ/psi',
0055 '1.000         mu+         mu-          PHOTOS VLL;',
0056 'Enddecay',
0057 '#',
0058 'End'),
0059             list_forced_decays = cms.vstring('MyLambda_b0','Myanti-Lambda_b0'),
0060             operates_on_particles = cms.vint32(),
0061             convertPythiaCodes = cms.untracked.bool(False)
0062             ),
0063         parameterSets = cms.vstring('EvtGen130')
0064         ),
0065         PythiaParameters = cms.PSet(
0066         pythia8CommonSettingsBlock,
0067         pythia8CP5SettingsBlock,
0068         processParameters = cms.vstring(
0069             "SoftQCD:nonDiffractive = on",
0070             "5122:m0=5.61960",     ## changing also Lambda_b0 mass in pythia
0071             'PTFilter:filter = on', # this turn on the filter
0072             'PTFilter:quarkToFilter = 5', # PDG id of q quark
0073             'PTFilter:scaleToFilter = 1.0'),
0074         parameterSets = cms.vstring(
0075             'pythia8CommonSettings',
0076             'pythia8CP5Settings',
0077             'processParameters',
0078         )
0079     )
0080 )
0081 _generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0082
0083 from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
0084 generator = ExternalGeneratorFilter(_generator)
0085
0086 ###########
0087 # Filters #
0088 ###########
0089 # Filter only pp events which produce a Lambdab:
0090 lbfilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(5122))
0091
0092 decayfilter = cms.EDFilter(
0093         "PythiaDauVFilter",
0094   verbose         = cms.untracked.int32(0),
0095   NumberDaughters = cms.untracked.int32(3),
0096   MotherID        = cms.untracked.int32(0),
0097   ParticleID      = cms.untracked.int32(5122),
0098   DaughterIDs     = cms.untracked.vint32(443, 3312, 310, 211),
0099   MinPt           = cms.untracked.vdouble(2.5, 0.7, 0.5, 0.1),
0100   MinEta          = cms.untracked.vdouble(-9999., -3.0, -3.0, -3.0),
0101   MaxEta          = cms.untracked.vdouble( 9999.,  3.0,  3.0,  3.0)
0102         )
0103
0104 jpsifilter = cms.EDFilter("PythiaDauVFilter",
0105   verbose         = cms.untracked.int32(0),
0106   NumberDaughters = cms.untracked.int32(2),
0107   MotherID        = cms.untracked.int32(5122),
0108   ParticleID      = cms.untracked.int32(443),
0109   DaughterIDs     = cms.untracked.vint32(13, -13),
0110   MinPt           = cms.untracked.vdouble(1., 1.),
0111   MinEta          = cms.untracked.vdouble(-3., -3.),
0112   MaxEta          = cms.untracked.vdouble( 3.,  3.)
0113           )
0114
0115 ProductionFilterSequence = cms.Sequence(generator*lbfilter*decayfilter*jpsifilter)