Generator/python/Chib1PToUpsilon1SGamma_MuFilter_TuneCP5_14TeV-pythia8_evtgen_cfi.py

0001 # 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
0002 #
0003 import FWCore.ParameterSet.Config as cms
0004 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0005 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0006
0007 generator = cms.EDFilter("Pythia8GeneratorFilter",
0008     pythiaPylistVerbosity = cms.untracked.int32(0),
0009     pythiaHepMCVerbosity = cms.untracked.bool(False),
0010     maxEventsToPrint = cms.untracked.int32(0),
0011     comEnergy = cms.double(14000.0),
0012     ExternalDecays = cms.PSet(
0013         EvtGen130 = cms.untracked.PSet(
0014             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
0015             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
0016             list_forced_decays = cms.vstring('mychi_b1','mychi_b2'),        # will force one at the time
0017             operates_on_particles = cms.vint32(20553,555),                  # we care just about our signal particles
0018             convertPythiaCodes = cms.untracked.bool(False),
0019             user_decay_embedded= cms.vstring(
0020 """
0021 Particle Upsilon 9.4603000 0.00005402
0022 Particle chi_b1 9.8927800 0.00000
0023 Particle chi_b2 9.9122100 0.00000
0024 Alias myUpsilon Upsilon
0025 Alias mychi_b1 chi_b1
0026 Alias mychi_b2 chi_b2
0027 Decay myUpsilon
0028 1.0   mu+  mu-          PHOTOS  VLL;
0029 Enddecay
0030 Decay mychi_b1
0031 1.0   gamma  myUpsilon  HELAMP 1. 0. 1. 0. -1. 0. -1. 0.;
0032 Enddecay
0033 Decay mychi_b2
0034 1.0   gamma  myUpsilon  PHSP;
0035 Enddecay
0036 End
0037 """
0038             )
0039             ),
0040         parameterSets = cms.vstring('EvtGen130')
0041     ),
0042     PythiaParameters = cms.PSet(
0043         pythia8CommonSettingsBlock,
0044         pythia8CP5SettingsBlock,
0045         processParameters = cms.vstring(
0046             'Bottomonium:states(3PJ) = 20553,555',
0047             'Bottomonium:O(3PJ)[3P0(1)] = 0.085,0.085',
0048             'Bottomonium:O(3PJ)[3S1(8)] = 0.04,0.04',
0049             'Bottomonium:gg2bbbar(3PJ)[3PJ(1)]g = on,on',
0050             'Bottomonium:qg2bbbar(3PJ)[3PJ(1)]q = on,on',
0051             'Bottomonium:qqbar2bbbar(3PJ)[3PJ(1)]g = on,on',
0052             'Bottomonium:gg2bbbar(3PJ)[3S1(8)]g = on,on',
0053             'Bottomonium:qg2bbbar(3PJ)[3S1(8)]q = on,on',
0054             'Bottomonium:qqbar2bbbar(3PJ)[3S1(8)]g = on,on',
0055             'PhaseSpace:pTHatMin = 2.',
0056             '20553:m0 = 9.89278',
0057             '555:m0 = 9.91221',
0058             '20553:onMode = off',
0059             '555:onMode = off'
0060             ),
0061         parameterSets = cms.vstring('pythia8CommonSettings',
0062                                     'pythia8CP5Settings',
0063                                     'processParameters',
0064                                     )
0065     )
0066 )
0067
0068 # We will filter for chi_b1, and chi_b2, first on the ID, then in the mass, this will constraint the photon daughter
0069
0070 pwaveIDfilter = cms.EDFilter("MCSingleParticleFilter",
0071     ParticleID = cms.untracked.vint32(20553,555),
0072     MinPt = cms.untracked.vdouble(0.0, 0.0),
0073     MinEta = cms.untracked.vdouble(-9., -9.),
0074     MaxEta = cms.untracked.vdouble(9., 9.),
0075     Status = cms.untracked.vint32(2, 2)
0076 )
0077
0078 pwaveMassfilter = cms.EDFilter("MCParticlePairFilter",
0079     Status = cms.untracked.vint32(2, 1),
0080     MinPt = cms.untracked.vdouble(7.9, 0.2),
0081     MaxEta = cms.untracked.vdouble(1.6, 1.6),
0082     MinEta = cms.untracked.vdouble(-1.6, -1.6),
0083     ParticleCharge = cms.untracked.int32(0),
0084     MinP = cms.untracked.vdouble(0.,0.),
0085     ParticleID1 = cms.untracked.vint32(553),
0086     ParticleID2 = cms.untracked.vint32(22),
0087     MinInvMass = cms.untracked.double(9.88),
0088     MaxInvMass = cms.untracked.double(9.93),
0089 )
0090
0091 # Next two muon filter are derived from muon reconstruction
0092
0093 muminusfilter = cms.EDFilter("PythiaDauVFilter",
0094     MotherID = cms.untracked.int32(0),
0095     MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
0096     ParticleID = cms.untracked.int32(553),
0097     ChargeConjugation = cms.untracked.bool(False),
0098     MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
0099     MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
0100     NumberDaughters = cms.untracked.int32(1),
0101     DaughterIDs = cms.untracked.vint32(-13, -13, -13)
0102 )
0103
0104 muplusfilter = cms.EDFilter("PythiaDauVFilter",
0105     MotherID = cms.untracked.int32(0),
0106     MinPt = cms.untracked.vdouble(2.5, 2.5, 3.5),
0107     ParticleID = cms.untracked.int32(553),
0108     ChargeConjugation = cms.untracked.bool(False),
0109     MinEta = cms.untracked.vdouble(1.2, -1.6, -1.2),
0110     MaxEta = cms.untracked.vdouble(1.6, -1.2, 1.2),
0111     NumberDaughters = cms.untracked.int32(1),
0112     DaughterIDs = cms.untracked.vint32(13, 13, 13)
0113 )
0114
0115 ProductionFilterSequence = cms.Sequence(generator*pwaveIDfilter*pwaveMassfilter*muminusfilter*muplusfilter)