Generator/python/BuToKstarJPsiToMuMu_forSTEAM_13TeV_ExtGen_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0003 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0004 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0005
0006 _generator = cms.EDFilter("Pythia8GeneratorFilter",
0007                          pythiaPylistVerbosity = cms.untracked.int32(0),
0008                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0009                          comEnergy = cms.double(13000.0),
0010                          maxEventsToPrint = cms.untracked.int32(0),
0011                          ExternalDecays = cms.PSet(
0012                          EvtGen130 = cms.untracked.PSet(
0013                             #uses latest evt and decay tables from evtgen
0014                             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_NOLONGLIFE.DEC'),
0015                             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0016                             convertPythiaCodes = cms.untracked.bool(False),
0017                             #user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Bu_Kstarmumu_Kspi.dec'),
0018                             #content was dump in the embed string below. This should test this feature.
0019                             list_forced_decays = cms.vstring('MyB+','MyB-'),
0020                             operates_on_particles = cms.vint32(),
0021                             user_decay_embedded= cms.vstring(
0022                                 '# This is the decay file for the decay B+ -> MU+ MU- K*+(-> Ks pi+)',
0023                                 '#',
0024                                 '# Descriptor: [B+ -> mu+ mu- {,gamma} {,gamma} (K*+ -> Ks pi+)]cc',
0025                                 '#',
0026                                 '# NickName:',
0027                                 '#',
0028                                 '# Physics: Includes radiative mode',
0029                                 '#',
0030                                 '# Tested: Yes',
0031                                 '# By:     K. Ulmer',
0032                                 '# Date:   2-26-08',
0033                                 '#',
0034                                 'Alias      MyB+        B+',
0035                                 'Alias      MyB-        B-',
0036                                 'ChargeConj MyB+        MyB-',
0037                                 'Alias      MyK*+       K*+',
0038                                 'Alias      MyK*-       K*-',
0039                                 'ChargeConj MyK*+       MyK*-',
0040                                 'Alias      MyK_S0      K_S0',
0041                                 'ChargeConj MyK_S0      MyK_S0',
0042                                 '#',
0043                                 'Decay MyB+',
0044                                 '  1.000        MyK*+     mu+     mu-      BTOSLLBALL;',
0045                                 'Enddecay',
0046                                 'CDecay MyB-',
0047                                 '#',
0048                                 'Decay MyK*+',
0049                                 '  1.000        MyK_S0    pi+              VSS;',
0050                                 'Enddecay',
0051                                 'CDecay MyK*-',
0052                                 '#',
0053                                 'Decay MyK_S0',
0054                                 '  1.000        pi+       pi-              PHSP;',
0055                                 'Enddecay',
0056                                 'End'
0057                           ),
0058                 ),
0059                 parameterSets = cms.vstring('EvtGen130')
0060         ),
0061         PythiaParameters = cms.PSet(
0062         pythia8CommonSettingsBlock,
0063         pythia8CP5SettingsBlock,
0064         processParameters = cms.vstring(
0065             'HardQCD:gg2bbbar = on ',
0066             'HardQCD:qqbar2bbbar = on ',
0067             'HardQCD:hardbbbar = on',
0068             'PhaseSpace:pTHatMin = 20.',
0069             ),
0070         parameterSets = cms.vstring('pythia8CommonSettings',
0071                                     'pythia8CP5Settings',
0072                                     'processParameters',
0073                                     )
0074         )
0075                          )
0076
0077 from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
0078 generator = ExternalGeneratorFilter(_generator)
0079
0080 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0081
0082 configurationMetadata = cms.untracked.PSet(
0083     version = cms.untracked.string('$Revision: 1.1 $'),
0084     name = cms.untracked.string('$Source: Configuration/Generator/python/BuToKstarMuMu_forSTEAM_13TeV_cfi.py $'),
0085     annotation = cms.untracked.string('Summer14: Pythia8+EvtGen130 generation of Bu --> K* Mu+Mu-, 13TeV, Tune CP5')
0086     )
0087
0088 ###########
0089 # Filters #
0090 ###########
0091 # Filter only pp events which produce a B+:
0092 bufilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(521))
0093
0094 # Filter on final state muons
0095 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
0096                              Status = cms.untracked.vint32(1, 1),
0097                              MinPt = cms.untracked.vdouble(2.8, 2.8),
0098                              MinP = cms.untracked.vdouble(2.8, 2.8),
0099                              MaxEta = cms.untracked.vdouble(2.3, 2.3),
0100                              MinEta = cms.untracked.vdouble(-2.3, -2.3),
0101                              ParticleID1 = cms.untracked.vint32(13,-13),
0102                              ParticleID2 = cms.untracked.vint32(13,-13)
0103                              )
0104
0105
0106 ProductionFilterSequence = cms.Sequence(generator*bufilter*mumugenfilter)