Generator/python/DStarToD0Pi_D0ToKsPiPi_inclusive_SoftQCD_TuneCP5_13p6TeV-pythia8-evtgen.py

0001 # Found 20 output events for 1500 input events.
0002 # Filter efficiency = 0.013333
0003 # Timing = 0.179750 sec/event
0004 # Event size = 636.9 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     pythiaPylistVerbosity = cms.untracked.int32(0),
0014     pythiaHepMCVerbosity = cms.untracked.bool(False),
0015     maxEventsToPrint = cms.untracked.int32(0),
0016     comEnergy = cms.double(13600.0),
0017     ExternalDecays = cms.PSet(
0018         EvtGen130 = cms.untracked.PSet(
0019             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0020             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0021             list_forced_decays = cms.vstring('MyDStar+','MyDStar-'),
0022             operates_on_particles = cms.vint32(),
0023             convertPythiaCodes = cms.untracked.bool(False),
0024             user_decay_embedded= cms.vstring(
0025 """
0026 #
0027 # Particles updated from PDG2018 https://journals.aps.org/prd/abstract/10.1103/PhysRevD.98.030001
0028 Particle   pi+         1.3957061e-01   0.0000000e+00
0029 Particle   pi-         1.3957061e-01   0.0000000e+00
0030 Particle   K_S0        4.9761100e-01   0.0000000e+00
0031 Particle   B-          5.2793200e+00   0.0000000e+00
0032 Particle   B+          5.2793200e+00   0.0000000e+00
0033 Particle   B0          5.2796300e+00   0.0000000e+00
0034 Particle   anti-B0     5.2796300e+00   0.0000000e+00
0035 Particle   B_s0        5.3668900e+00   0.0000000e+00
0036 Particle   anti-B_s0   5.3668900e+00   0.0000000e+00
0037 Particle   phi         1.0194610e+00   4.2490000e-03
0038 Particle   D+          1.8696500e+00   0.0000000e+00
0039 Particle   D-          1.8696500e+00   0.0000000e+00
0040 Particle   D0          1.8648300e+00   0.0000000e+00
0041 Particle   D*+         2.0102600e+00   0.0000834e-05
0042 Particle   D*-         2.0102600e+00   0.0000834e-05
0043 Particle   K+          4.9367700e-01   0.0000000e+00
0044 Particle   K-          4.9367700e-01   0.0000000e+00
0045
0046 #
0047 Alias      MyDStar+    D*+
0048 Alias      MyDStar-    D*-
0049 ChargeConj MyDStar-    MyDStar+
0050 #
0051 Alias      MyD0        D0
0052 Alias      MyAntiD0    anti-D0
0053 ChargeConj MyAntiD0    MyD0
0054 #
0055 Decay MyDStar+
0056 1.000  MyD0  pi+    VSS;
0057 Enddecay
0058 CDecay MyDStar-
0059 #
0060 Decay MyD0
0061 1.000  K_S0  pi+  pi-    PHSP;
0062 Enddecay
0063 CDecay MyAntiD0
0064 #
0065 End
0066 """
0067             )
0068         ),
0069         parameterSets = cms.vstring('EvtGen130')
0070     ),
0071     PythiaParameters = cms.PSet(
0072         pythia8CommonSettingsBlock,
0073         pythia8CP5SettingsBlock,
0074         processParameters = cms.vstring('SoftQCD:nonDiffractive = on',
0075                                         'PTFilter:filter = on', # this turn on the filter
0076                                         'PTFilter:quarkToFilter = 4', # PDG id of q quark
0077                                         'PTFilter:scaleToFilter = 2.0'
0078             ),
0079         parameterSets = cms.vstring('pythia8CommonSettings',
0080                                     'pythia8CP5Settings',
0081                                     'processParameters',
0082                                     )
0083     )
0084 )
0085
0086 _generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0087
0088 from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
0089 generator = ExternalGeneratorFilter(_generator)
0090
0091 ###### Filters ##########
0092 decayfilter = cms.EDFilter(
0093     "PythiaDauVFilter",
0094     verbose         = cms.untracked.int32(1),
0095     NumberDaughters = cms.untracked.int32(2),
0096     ParticleID      = cms.untracked.int32(413),  # DStar+
0097     DaughterIDs     = cms.untracked.vint32(421, 211), # D0 and pi+
0098     MinPt           = cms.untracked.vdouble(3.5, 0.1),
0099     MinEta          = cms.untracked.vdouble(-3., -3.),
0100     MaxEta          = cms.untracked.vdouble( 3.,  3.)
0101 )
0102
0103 D0filter = cms.EDFilter(
0104     "PythiaDauVFilter",
0105     verbose         = cms.untracked.int32(1),
0106     NumberDaughters = cms.untracked.int32(2),
0107     MotherID        = cms.untracked.int32(413), # DStar+
0108     ParticleID      = cms.untracked.int32(421),  # D0
0109     DaughterIDs     = cms.untracked.vint32(310, 211, -211), # K0s pi+ pi-
0110     MinPt           = cms.untracked.vdouble(0.5, 0.1, 0.1),
0111     MinEta          = cms.untracked.vdouble(-3., -3., -3.),
0112     MaxEta          = cms.untracked.vdouble( 3.,  3.,  3.)
0113 )
0114
0115
0116 ProductionFilterSequence = cms.Sequence(generator*decayfilter*D0filter)