Merge pull request #4903 from heppye/CMSSW_7_2_X

BPH cfi for STEAM

Configuration/Generator/python/ BsToMuMu_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,55 @@
+# Only for trigger study
+import FWCore.ParameterSet.Config as cms
+from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
+generator = cms.EDFilter(
+"Pythia6GeneratorFilter",
+comEnergy = cms.double(13000.0),
+crossSection = cms.untracked.double(54000000000),
+filterEfficiency = cms.untracked.double(3.0e-4),
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+maxEventsToPrint = cms.untracked.int32(0),
+pythiaPylistVerbosity = cms.untracked.int32(0),
+ExternalDecays = cms.PSet(
+EvtGen = cms.untracked.PSet(
+operates_on_particles = cms.vint32( 0 ), # 0 (zero) means default list (hardcoded)
+# you can put here the list of particles (PDG IDs)
+# that you want decayed by EvtGen
+use_default_decay = cms.untracked.bool(False),
+decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
+particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
+user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bs_mumu.dec'),
+list_forced_decays = cms.vstring('MyB_s0',
+'Myanti-B_s0'),
+),
+parameterSets = cms.vstring('EvtGen')
+),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+bbbarSettings = cms.vstring('MSEL = 1'),
+# This is a vector of ParameterSet names to be read, in this order
+parameterSets = cms.vstring(
+'pythiaUESettings',
+'bbbarSettings')
+)
+)
+MuMuFilter = cms.EDFilter("MCParticlePairFilter",
+Status = cms.untracked.vint32(1, 1),
+MinPt = cms.untracked.vdouble(3., 3.),
+MaxEta = cms.untracked.vdouble(2.5, 2.5),
+MinEta = cms.untracked.vdouble(-2.5, -2.5),
+ParticleCharge = cms.untracked.int32(-1),
+MaxInvMass = cms.untracked.double(5.5),
+MinInvMass = cms.untracked.double(5.0),
+ParticleID1 = cms.untracked.vint32(13),
+ParticleID2 = cms.untracked.vint32(13)
+)
+# -- Require Muon from Bs
+MuFilter = cms.EDFilter("PythiaFilter",
+Status = cms.untracked.int32(1),
+MotherID = cms.untracked.int32(531),
+MinPt = cms.untracked.double(3.),
+ParticleID = cms.untracked.int32(13),
+MaxEta = cms.untracked.double(2.5),
+MinEta = cms.untracked.double(-2.5)
+)
+ProductionFilterSequence = cms.Sequence(generator*MuMuFilter*MuFilter)

Configuration/Generator/python/ EtaBToJpsiJpsi_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,55 @@
+# Only for Trigger Study
+import FWCore.ParameterSet.Config as cms
+from Configuration.Generator.PythiaUESettings_cfi import *
+generator = cms.EDFilter("Pythia6GeneratorFilter",
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+maxEventsToPrint = cms.untracked.int32(3),
+pythiaPylistVerbosity = cms.untracked.int32(1),
+displayPythiaCards = cms.untracked.bool(False),
+comEnergy = cms.double(13000.0),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+pythiaEtab = cms.vstring(
+'MSEL=0',
+'MSUB(471)=1',
+'MDME(1521,1)=0',
+'KFDP(1520,1)=443',
+'KFDP(1520,2)=443',
+'443:ALLOFF',
+'443:ONIFMATCH 13 13',
+'333:ALLOFF',
+'333:ONIFMATCH 321 321',
+'PMAS(C10551,1)=7.0'),
+# This is a vector of ParameterSet names to be read, in this order
+parameterSets = cms.vstring(
+'pythiaUESettings',
+'pythiaEtab')
+)
+)
+etafilter = cms.EDFilter(
+"PythiaFilter",
+MaxEta = cms.untracked.double(9999.),
+MinEta = cms.untracked.double(-9999.),
+ParticleID = cms.untracked.int32(10551)
+)
+jpsifilter = cms.EDFilter(
+"PythiaDauVFilter",
+verbose = cms.untracked.int32(0),
+NumberDaughters = cms.untracked.int32(2),
+MotherID = cms.untracked.int32(10551),
+ParticleID = cms.untracked.int32(443),
+DaughterIDs = cms.untracked.vint32(13, -13),
+MinP = cms.untracked.vdouble(2.5,2.5),
+MinEta = cms.untracked.vdouble(-2.6, -2.6),
+MaxEta = cms.untracked.vdouble( 2.6, 2.6)
+)
+#mumugenfilter = cms.EDFilter(
+# "MCParticlePairFilter",
+# Status = cms.untracked.vint32(1, 1),
+# MaxEta = cms.untracked.vdouble(2.6, 2.6),
+# MinEta = cms.untracked.vdouble(-2.6, -2.6),
+# MinP = cms.untracked.vdouble(2.5, 2.5),
+# ParticleID1 = cms.untracked.vint32(13),
+# ParticleID2 = cms.untracked.vint32(-13)
+# )
+ProductionFilterSequence = cms.Sequence(generator*etafilter*jpsifilter)

Configuration/Generator/python/BdToMuMu_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,55 @@
+# Only for Trigger Study
+import FWCore.ParameterSet.Config as cms
+from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
+generator = cms.EDFilter(
+"Pythia6GeneratorFilter",
+comEnergy = cms.double(13000.0),
+crossSection = cms.untracked.double(54000000000),
+filterEfficiency = cms.untracked.double(0.0013),
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+maxEventsToPrint = cms.untracked.int32(0),
+pythiaPylistVerbosity = cms.untracked.int32(0),
+ExternalDecays = cms.PSet(
+EvtGen = cms.untracked.PSet(
+operates_on_particles = cms.vint32( 0 ), # 0 (zero) means default list (hardcoded)
+# you can put here the list of particles (PDG IDs)
+# that you want decayed by EvtGen
+use_default_decay = cms.untracked.bool(False),
+decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
+particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
+user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bd_mumu.dec'),
+list_forced_decays = cms.vstring('MyB0',
+'Myanti-B0'),
+),
+parameterSets = cms.vstring('EvtGen')
+),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+bbbarSettings = cms.vstring('MSEL = 1'),
+# This is a vector of ParameterSet names to be read, in this order
+parameterSets = cms.vstring(
+'pythiaUESettings',
+'bbbarSettings')
+)
+)
+MuMuFilter = cms.EDFilter("MCParticlePairFilter",
+Status = cms.untracked.vint32(1, 1),
+MinPt = cms.untracked.vdouble(3., 3.),
+MaxEta = cms.untracked.vdouble(2.5, 2.5),
+MinEta = cms.untracked.vdouble(-2.5, -2.5),
+ParticleCharge = cms.untracked.int32(-1),
+MaxInvMass = cms.untracked.double(5.5),
+MinInvMass = cms.untracked.double(5.0),
+ParticleID1 = cms.untracked.vint32(13),
+ParticleID2 = cms.untracked.vint32(13)
+)
+# -- Require Muon from Bd
+MuFilter = cms.EDFilter("PythiaFilter",
+Status = cms.untracked.int32(1),
+MotherID = cms.untracked.int32(511),
+MinPt = cms.untracked.double(3.),
+ParticleID = cms.untracked.int32(13),
+MaxEta = cms.untracked.double(2.5),
+MinEta = cms.untracked.double(-2.5)
+)
+ProductionFilterSequence = cms.Sequence(generator*MuMuFilter*MuFilter)

Configuration/Generator/python/BuToKstarMuMu_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,42 @@
+# Only for Trigger Study
+import FWCore.ParameterSet.Config as cms
+configurationMetadata = cms.untracked.PSet(
+version = cms.untracked.string('$Revision: 1.1 $'),
+name = cms.untracked.string('$Source: /local/reps/CMSSW/CMSSW/Configuration/GenProduction/python/EightTeV/PYTHIA6_BuToKstarMuMu_EtaPtFilter_TuneZ2star_8TeV_cff.py,v $'),
+annotation = cms.untracked.string('Summer11: Pythia6+EvtGen generation of Bu --> K* Mu+Mu-, 8TeV, D6T tune'))
+from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
+generator = cms.EDFilter("Pythia6GeneratorFilter",
+pythiaPylistVerbosity = cms.untracked.int32(0),
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+comEnergy = cms.double(13000.0),
+crossSection = cms.untracked.double(54000000000), # Given by PYTHIA after running
+filterEfficiency = cms.untracked.double(0.004), # Given by PYTHIA after running
+maxEventsToPrint = cms.untracked.int32(0),
+ExternalDecays = cms.PSet(
+EvtGen = cms.untracked.PSet(
+operates_on_particles = cms.vint32(0),
+use_default_decay = cms.untracked.bool(False),
+decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY.DEC'),
+particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
+user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bu_Kstarmumu_Kspi.dec'),
+list_forced_decays = cms.vstring('MyB+','MyB-')),
+parameterSets = cms.vstring('EvtGen')),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+bbbarSettings = cms.vstring('MSEL = 1'),
+parameterSets = cms.vstring('pythiaUESettings','bbbarSettings')))
+###########
+# Filters #
+###########
+# Filter only pp events which produce a B+:
+bufilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(521))
+# Filter on final state muons
+mumugenfilter = cms.EDFilter("MCParticlePairFilter",
+Status = cms.untracked.vint32(1, 1),
+MinPt = cms.untracked.vdouble(2.8, 2.8),
+MaxEta = cms.untracked.vdouble(2.3, 2.3),
+MinEta = cms.untracked.vdouble(-2.3, -2.3),
+ParticleCharge = cms.untracked.int32(-1),
+ParticleID1 = cms.untracked.vint32(13),
+ParticleID2 = cms.untracked.vint32(13))
+ProductionFilterSequence = cms.Sequence(generator*bufilter*mumugenfilter)

Configuration/Generator/python/EtaBToUpsilon1SUpsilon1S_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,46 @@
+# Only for Trigger Study
+import FWCore.ParameterSet.Config as cms
+from Configuration.Generator.PythiaUESettings_cfi import *
+generator = cms.EDFilter("Pythia6GeneratorFilter",
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+maxEventsToPrint = cms.untracked.int32(3),
+pythiaPylistVerbosity = cms.untracked.int32(1),
+displayPythiaCards = cms.untracked.bool(False),
+comEnergy = cms.double(13000.0),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+pythiaEtab = cms.vstring(
+'MSEL=0',
+'MSUB(471)=1',
+'MDME(1521,1)=0',
+'KFDP(1520,1)=553',
+'KFDP(1520,2)=553',
+'553:ALLOFF',
+'553:ONIFMATCH 13 13',
+'333:ALLOFF',
+'333:ONIFMATCH 321 321',
+'PMAS(C10551,1)=20.0'),
+# This is a vector of ParameterSet names to be read, in this order
+parameterSets = cms.vstring(
+'pythiaUESettings',
+'pythiaEtab')
+)
+)
+etafilter = cms.EDFilter(
+"PythiaFilter",
+MaxEta = cms.untracked.double(9999.),
+MinEta = cms.untracked.double(-9999.),
+ParticleID = cms.untracked.int32(10551)
+)
+upsilonfilter = cms.EDFilter(
+"PythiaDauVFilter",
+verbose = cms.untracked.int32(0),
+NumberDaughters = cms.untracked.int32(2),
+MotherID = cms.untracked.int32(10551),
+ParticleID = cms.untracked.int32(553),
+DaughterIDs = cms.untracked.vint32(13, -13),
+MinP = cms.untracked.vdouble(2.5,2.5),
+MinEta = cms.untracked.vdouble(-2.6, -2.6),
+MaxEta = cms.untracked.vdouble( 2.6, 2.6)
+)
+ProductionFilterSequence = cms.Sequence(generator*etafilter*upsilonfilter)

Configuration/Generator/python/JpsiMuMu_forSTEAM_13TeV_cfi.py

@@ -0,0 +1,82 @@
+# Only for Trigger Study
+import FWCore.ParameterSet.Config as cms
+from Configuration.Generator.PythiaUEZ2starSettings_cfi import *
+generator = cms.EDFilter("Pythia6GeneratorFilter",
+pythiaPylistVerbosity = cms.untracked.int32(0),
+filterEfficiency = cms.untracked.double(0.15),
+pythiaHepMCVerbosity = cms.untracked.bool(False),
+crossSection = cms.untracked.double(1712000.0),
+comEnergy = cms.double(13000.0),
+maxEventsToPrint = cms.untracked.int32(0),
+PythiaParameters = cms.PSet(
+pythiaUESettingsBlock,
+processParameters = cms.vstring(
+'MSEL=61 ! Quarkonia',
+'CKIN(3)=10. ! Min pthard',
+'CKIN(4)=-1. ! Max pthard',
+'MDME(858,1) = 0 ! 0.060200 e- e+',
+'MDME(859,1) = 1 ! 0.060100 mu- mu+',
+'MDME(860,1) = 0 ! 0.879700 rndmflav rndmflavbar',
+'MSTP(142)=2 ! turns on the PYEVWT Pt re-weighting routine',
+'PARJ(13)=0.750 ! probability that a c or b meson has S=1',
+'PARJ(14)=0.162 ! probability that a meson with S=0 is produced with L=1, J=1',
+'PARJ(15)=0.018 ! probability that a meson with S=1 is produced with L=1, J=0',
+'PARJ(16)=0.054 ! probability that a meson with S=1 is produced with L=1, J=1',
+'MSTP(145)=0 ! choice of polarization',
+'MSTP(146)=0 ! choice of polarization frame ONLY when mstp(145)=1',
+'MSTP(147)=0 ! particular helicity or density matrix component when mstp(145)=1',
+'MSTP(148)=1 ! possibility to allow for final-state shower evolution, extreme case !',
+'MSTP(149)=1 ! if mstp(148)=1, it determines the kinematics of the QQ~3S1(8)->QQ~3S1(8)+g branching',
+'PARP(141)=1.16 ! New values for COM matrix elements',
+'PARP(142)=0.0119 ! New values for COM matrix elements',
+'PARP(143)=0.01 ! New values for COM matrix elements',
+'PARP(144)=0.01 ! New values for COM matrix elements',
+'PARP(145)=0.05 ! New values for COM matrix elements',
+'PARP(146)=9.28 ! New values for COM matrix elements',
+'PARP(147)=0.15 ! New values for COM matrix elements',
+'PARP(148)=0.02 ! New values for COM matrix elements',
+'PARP(149)=0.02 ! New values for COM matrix elements',
+'PARP(150)=0.085 ! New values for COM matrix elements',
+'BRAT(859)=1.000 ! J/psi->mu+mu-',
+'BRAT(861)=0.000 ! chi_2c->J/psi gamma',
+'BRAT(862)=0.798 ! chi_2c->rndmflav rndmflavbar',
+'BRAT(1501)=0.000 ! chi_0c->J/psi gamma',
+'BRAT(1502)=0.987 ! chi_0c->rndmflav rndmflavbar',
+'BRAT(1555)=0.000 ! chi_1c->J/psi gamma',
+'BRAT(1556)=0.644 ! chi_1c->rndmflav rndmflavbar',
+'BRAT(1569)=0.186600 ! psi(2S) -> rndmflav rndmflavbar',
+'BRAT(1570)=0.000 ! psi(2S) ->J/psi pi+ pi-',
+'BRAT(1571)=0.000 ! psi(2S) ->J/psi pi0 pi0',
+'BRAT(1572)=0.000 ! psi(2S) ->J/psi eta',
+'BRAT(1573)=0.000 ! psi(2S) ->J/psi pi0'
+),
+# This is a vector of ParameterSet names to be read, in this order
+parameterSets = cms.vstring('pythiaUESettings',
+'processParameters',
+'CSAParameters'),
+CSAParameters = cms.vstring('CSAMODE = 6 ! cross-section reweighted quarkonia')
+)
+)
+oniafilter = cms.EDFilter("PythiaFilter",
+Status = cms.untracked.int32(2),
+MaxEta = cms.untracked.double(1000.0),
+MinEta = cms.untracked.double(-1000.0),
+MinPt = cms.untracked.double(0.0),
+ParticleID = cms.untracked.int32(443)
+)
+mumugenfilter = cms.EDFilter("MCParticlePairFilter",
+Status = cms.untracked.vint32(1, 1),
+MinPt = cms.untracked.vdouble(0.5, 0.5),
+MinP = cms.untracked.vdouble(2.7, 2.7),
+MaxEta = cms.untracked.vdouble(2.5, 2.5),
+MinEta = cms.untracked.vdouble(-2.5, -2.5),
+ParticleCharge = cms.untracked.int32(-1),
+ParticleID1 = cms.untracked.vint32(13),
+ParticleID2 = cms.untracked.vint32(13)
+)
+mugenfilter = cms.EDFilter("MCSingleParticleFilter",
+Status = cms.untracked.vint32(1,1),
+MinPt = cms.untracked.vdouble(10.0,10.0),
+ParticleID = cms.untracked.vint32(13,-13),
+)
+ProductionFilterSequence = cms.Sequence(generator*oniafilter*mumugenfilter*~mugenfilter)