Generator/python/B0ToKstarMuMu_14TeV_TuneCP5_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0004 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0005 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0006
0007 generator = cms.EDFilter("Pythia8GeneratorFilter",
0008     pythiaHepMCVerbosity = cms.untracked.bool(False),
0009     maxEventsToPrint = cms.untracked.int32(0),
0010     pythiaPylistVerbosity = cms.untracked.int32(0),
0011     filterEfficiency = cms.untracked.double(1.38e-3),
0012     crossSection = cms.untracked.double(540000000.),
0013     comEnergy = cms.double(14000.0),
0014     ExternalDecays = cms.PSet(
0015         EvtGen130 = cms.untracked.PSet(
0016             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
0017             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0018             user_decay_embedded= cms.vstring(
0019                 'Alias      MyB0        B0',
0020                 'Alias      Myanti-B0   anti-B0',
0021                 'ChargeConj MyB0        Myanti-B0',
0022                 'Alias      MyK*0       K*0',
0023                 'Alias      MyK*0bar    anti-K*0',
0024                 'ChargeConj MyK*0       MyK*0bar',
0025                 '#',
0026                 'Decay MyB0',
0027                 '1.000        MyK*0     mu+     mu-               BTOSLLBALL;',
0028                 'Enddecay',
0029                 'Decay Myanti-B0',
0030                 '1.000        MyK*0bar     mu+     mu-            BTOSLLBALL;',
0031                 'Enddecay',
0032                 '#',
0033                 'Decay MyK*0',
0034                 '1.000        K+        pi-                    VSS;',
0035                 'Enddecay',
0036                 'Decay MyK*0bar',
0037                 '1.000        K-        pi+                    VSS;',
0038                 'Enddecay ',
0039                 'End'
0040             ),
0041             list_forced_decays = cms.vstring('MyB0','Myanti-B0'),
0042             operates_on_particles = cms.vint32(),
0043         ),
0044         parameterSets = cms.vstring('EvtGen130')
0045     ),
0046     PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
0047         pythia8CP5SettingsBlock,
0048         processParameters = cms.vstring(
0049             "SoftQCD:nonDiffractive = on",
0050             'PTFilter:filter = on', # this turn on the filter

0051             'PTFilter:quarkToFilter = 5', # PDG id of q quark

0052             'PTFilter:scaleToFilter = 1.0'
0053         ),
0054         parameterSets = cms.vstring('pythia8CommonSettings',
0055             'pythia8CP5Settings',
0056             'processParameters',
0057         )
0058     )
0059 )
0060
0061 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0062
0063 configurationMetadata = cms.untracked.PSet(
0064     version = cms.untracked.string('$Revision: 1.0 $'),
0065     name = cms.untracked.string('$Source: /Configuration/Generator/python/BPH_BsMuMu_PhaseII_cfi.py $'),
0066     annotation = cms.untracked.string('PhaseII: Pythia8+EvtGen130 generation of B0 --> K* mu+mu-, 14TeV, Tune CP5')
0067 )
0068
0069 bfilter = cms.EDFilter(
0070     "PythiaFilter",
0071     MaxEta = cms.untracked.double(9999.),
0072     MinEta = cms.untracked.double(-9999.),
0073     ParticleID = cms.untracked.int32(511)  ## Bd

0074     )
0075
0076 decayfilter = cms.EDFilter(
0077     "PythiaDauVFilter",
0078     verbose         = cms.untracked.int32(1),
0079     NumberDaughters = cms.untracked.int32(3),
0080     ParticleID      = cms.untracked.int32(511),
0081     DaughterIDs     = cms.untracked.vint32(-13, 13, 313),  ## mu+, mu-, K*^0(892)

0082     MinPt           = cms.untracked.vdouble(2.5, 2.5, -1.),
0083     MinEta          = cms.untracked.vdouble(-2.9, -2.9, -9999.),
0084     MaxEta          = cms.untracked.vdouble( 2.9,  2.9,  9999.)
0085     )
0086
0087 kstarfilter = cms.EDFilter(
0088     "PythiaDauVFilter",
0089     verbose         = cms.untracked.int32(1),
0090     NumberDaughters = cms.untracked.int32(2),
0091     MotherID        = cms.untracked.int32(511),  ## Bd

0092     ParticleID      = cms.untracked.int32(313),  ## K*^0(892)

0093     DaughterIDs     = cms.untracked.vint32(321, -211), ## K+, pi-

0094     MinPt           = cms.untracked.vdouble(0.4, 0.4),
0095     MinEta          = cms.untracked.vdouble(-4.1, -4.1),
0096     MaxEta          = cms.untracked.vdouble( 4.1,  4.1)
0097     )
0098
0099 ProductionFilterSequence = cms.Sequence(generator*bfilter*decayfilter*kstarfilter)