Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

add genfragment for Chi_b -> Y(1S)Gamma and Chi_c -> J/psiGamma #35330

Merged
merged 5 commits into from Nov 2, 2021
Merged

add genfragment for Chi_b -> Y(1S)Gamma and Chi_c -> J/psiGamma #35330

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
fd3db7f
add genfragment for Chi_b -> Y(1S)Gamma and Chi_c -> J/psiGamma
mtosi Sep 18, 2021
4b1f1d3
PhaseSpace:pTHatMin
mtosi Sep 23, 2021
1acd8fa
Rename Chib1PToUpsilon1SGamma_TuneCP5_14TeV-pythia8-evtgen_cfi.py to …
mtosi Oct 25, 2021
b009538
Rename Chib1PToUpsilon1SGamma_MuFilter_TuneCP5_14TeV-pythia8-evtgen_c…
mtosi Oct 26, 2021
79c1503
Rename ChicToJpsiGamma_MuFilter_TuneCP5_14TeV_pythia8-evtgen_cfi.py t…
mtosi Oct 26, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
115 changes: 115 additions & 0 deletions ...tion/Generator/python/Chib1PToUpsilon1SGamma_MuFilter_TuneCP5_14TeV-pythia8_evtgen_cfi.py
View file
Open in desktop
@@ -0,0 +1,115 @@
# cfg file for Chi_{b1,b2}(1P) -> Upsilon(1S) gamma, we are neglecting Chi_b0. Masses and widths are matched between pythia, evtgen and PDG 2016
#
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *

generator = cms.EDFilter("Pythia8GeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
comEnergy = cms.double(14000.0),
ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
list_forced_decays = cms.vstring('mychi_b1','mychi_b2'), # will force one at the time
operates_on_particles = cms.vint32(20553,555), # we care just about our signal particles
convertPythiaCodes = cms.untracked.bool(False),
user_decay_embedded= cms.vstring(
"""
Particle Upsilon 9.4603000 0.00005402
Particle chi_b1 9.8927800 0.00000
Particle chi_b2 9.9122100 0.00000
Alias myUpsilon Upsilon
Alias mychi_b1 chi_b1
Alias mychi_b2 chi_b2
Decay myUpsilon
1.0 mu+ mu- PHOTOS VLL;
Enddecay
Decay mychi_b1
1.0 gamma myUpsilon HELAMP 1. 0. 1. 0. -1. 0. -1. 0.;
Enddecay
Decay mychi_b2
1.0 gamma myUpsilon PHSP;
Enddecay
End
"""
)
),
parameterSets = cms.vstring('EvtGen130')
),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'Bottomonium:states(3PJ) = 20553,555',
'Bottomonium:O(3PJ)[3P0(1)] = 0.085,0.085',
'Bottomonium:O(3PJ)[3S1(8)] = 0.04,0.04',
'Bottomonium:gg2bbbar(3PJ)[3PJ(1)]g = on,on',
'Bottomonium:qg2bbbar(3PJ)[3PJ(1)]q = on,on',
'Bottomonium:qqbar2bbbar(3PJ)[3PJ(1)]g = on,on',
'Bottomonium:gg2bbbar(3PJ)[3S1(8)]g = on,on',
'Bottomonium:qg2bbbar(3PJ)[3S1(8)]q = on,on',
'Bottomonium:qqbar2bbbar(3PJ)[3S1(8)]g = on,on',
'PhaseSpace:pTHatMin = 2.',
'20553:m0 = 9.89278',
'555:m0 = 9.91221',
'20553:onMode = off',
'555:onMode = off'
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)
)
)

# We will filter for chi_b1, and chi_b2, first on the ID, then in the mass, this will constraint the photon daughter

pwaveIDfilter = cms.EDFilter("MCSingleParticleFilter",
ParticleID = cms.untracked.vint32(20553,555),
MinPt = cms.untracked.vdouble(0.0, 0.0),
MinEta = cms.untracked.vdouble(-9., -9.),
MaxEta = cms.untracked.vdouble(9., 9.),
Status = cms.untracked.vint32(2, 2)
)

pwaveMassfilter = cms.EDFilter("MCParticlePairFilter",
Status = cms.untracked.vint32(2, 1),
MinPt = cms.untracked.vdouble(7.9, 0.2),
MaxEta = cms.untracked.vdouble(1.6, 1.6),
MinEta = cms.untracked.vdouble(-1.6, -1.6),
ParticleCharge = cms.untracked.int32(0),
MinP = cms.untracked.vdouble(0.,0.),
ParticleID1 = cms.untracked.vint32(553),
ParticleID2 = cms.untracked.vint32(22),
MinInvMass = cms.untracked.double(9.88),
MaxInvMass = cms.untracked.double(9.93),
)

# Next two muon filter are derived from muon reconstruction

muminusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
ParticleID = cms.untracked.int32(553),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(-13, -13, -13)
)

muplusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
ParticleID = cms.untracked.int32(553),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(13, 13, 13)
)

ProductionFilterSequence = cms.Sequence(generator*pwaveIDfilter*pwaveMassfilter*muminusfilter*muplusfilter)
79 changes: 79 additions & 0 deletions Configuration/Generator/python/ChicToJpsiGamma_MuFilter_TuneCP5_14TeV_pythia8_evtgen_cfi.py
View file
Open in desktop
@@ -0,0 +1,79 @@
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from Configuration.Generator.Pythia8CommonSettings_cfi import *

generator = cms.EDFilter("Pythia8GeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
comEnergy = cms.double(14000.0),
ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
operates_on_particles = cms.vint32(20443,445), # we care just about our signal particles
convertPythiaCodes = cms.untracked.bool(False),
user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Onia_chic_jpsigamma.dec'),
list_forced_decays = cms.vstring('Mychi_c1','Mychi_c2'),
),
parameterSets = cms.vstring('EvtGen130')
),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
# generate just the needed and nothing else
'Charmonium:states(3PJ) = 20443,445',
'Charmonium:O(3PJ)[3P0(1)] = 0.05,0.05',
'Charmonium:O(3PJ)[3S1(8)] = 0.0031,0.0031',
'Charmonium:gg2ccbar(3PJ)[3PJ(1)]g = on,on',
'Charmonium:qg2ccbar(3PJ)[3PJ(1)]q = on,on',
'Charmonium:qqbar2ccbar(3PJ)[3PJ(1)]g = on,on',
'Charmonium:gg2ccbar(3PJ)[3S1(8)]g = on,on',
'Charmonium:qg2ccbar(3PJ)[3S1(8)]q = on,on',
'Charmonium:qqbar2ccbar(3PJ)[3S1(8)]g = on,on',
#
'PhaseSpace:pTHatMin = 10.' # (filter efficiency 4.671e-02) be aware of this ckin(3) equivalent
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)
)
)

# Next two muon filter are derived from muon reconstruction
oniafilter = cms.EDFilter("MCParticlePairFilter",
Status = cms.untracked.vint32(2, 1),
MinPt = cms.untracked.vdouble(16., 0.2),
MaxEta = cms.untracked.vdouble(1.6, 1.6),
MinEta = cms.untracked.vdouble(-1.6, -1.6),
ParticleCharge = cms.untracked.int32(0),
MinP = cms.untracked.vdouble(0.,0.),
ParticleID1 = cms.untracked.vint32(443),
ParticleID2 = cms.untracked.vint32(22)
)

muminusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
ParticleID = cms.untracked.int32(443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(-13, -13, -13)
)

muplusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
ParticleID = cms.untracked.int32(443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(13, 13, 13)
)

ProductionFilterSequence = cms.Sequence(generator*oniafilter*muminusfilter*muplusfilter)