Generator/python/BuToKstarJPsiToMuMu_forSTEAM_13TeV_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 """
0023 # This is the decay file for the decay B+ -> MU+ MU- K*+(-> Ks pi+)
0024 #
0025 # Descriptor: [B+ -> mu+ mu- {,gamma} {,gamma} (K*+ -> Ks pi+)]cc
0026 #
0027 # NickName:
0028 #
0029 # Physics: Includes radiative mode
0030 #
0031 # Tested: Yes
0032 # By:     K. Ulmer
0033 # Date:   2-26-08
0034 #
0035 Alias      MyB+        B+
0036 Alias      MyB-        B-
0037 ChargeConj MyB+        MyB-
0038 Alias      MyK*+       K*+
0039 Alias      MyK*-       K*-
0040 ChargeConj MyK*+       MyK*-
0041 Alias      MyK_S0      K_S0
0042 ChargeConj MyK_S0      MyK_S0
0043 #
0044 Decay MyB+
0045   1.000        MyK*+     mu+     mu-      BTOSLLBALL;
0046 Enddecay
0047 CDecay MyB-
0048 #
0049 Decay MyK*+
0050   1.000        MyK_S0    pi+              VSS;
0051 Enddecay
0052 CDecay MyK*-
0053 #
0054 Decay MyK_S0
0055   1.000        pi+       pi-              PHSP;
0056 Enddecay
0057 End
0058 """
0059                           ),
0060                 ),
0061                 parameterSets = cms.vstring('EvtGen130')
0062         ),
0063         PythiaParameters = cms.PSet(
0064         pythia8CommonSettingsBlock,
0065         pythia8CP5SettingsBlock,
0066         processParameters = cms.vstring(
0067             'HardQCD:gg2bbbar = on ',
0068             'HardQCD:qqbar2bbbar = on ',
0069             'HardQCD:hardbbbar = on',
0070             'PhaseSpace:pTHatMin = 20.',
0071             ),
0072         parameterSets = cms.vstring('pythia8CommonSettings',
0073                                     'pythia8CP5Settings',
0074                                     'processParameters',
0075                                     )
0076         )
0077                          )
0078
0079 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0080
0081 configurationMetadata = cms.untracked.PSet(
0082     version = cms.untracked.string('$Revision: 1.1 $'),
0083     name = cms.untracked.string('$Source: Configuration/Generator/python/BuToKstarMuMu_forSTEAM_13TeV_cfi.py $'),
0084     annotation = cms.untracked.string('Summer14: Pythia8+EvtGen130 generation of Bu --> K* Mu+Mu-, 13TeV, Tune CP5')
0085     )
0086
0087 ###########
0088 # Filters #
0089 ###########
0090 # Filter only pp events which produce a B+:
0091 bufilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(521))
0092
0093 # Filter on final state muons
0094 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
0095                              Status = cms.untracked.vint32(1, 1),
0096                              MinPt = cms.untracked.vdouble(2.8, 2.8),
0097                              MinP = cms.untracked.vdouble(2.8, 2.8),
0098                              MaxEta = cms.untracked.vdouble(2.3, 2.3),
0099                              MinEta = cms.untracked.vdouble(-2.3, -2.3),
0100                              ParticleID1 = cms.untracked.vint32(13,-13),
0101                              ParticleID2 = cms.untracked.vint32(13,-13)
0102                              )
0103
0104
0105 ProductionFilterSequence = cms.Sequence(generator*bufilter*mumugenfilter)