Generator/python/HTo2LongLivedTo4bMH-125_MFF-12_CTau-900mm_TuneCP5_14TeV_pythia8_cfi.py

0001 #CROSS_SECTION = 1 # pb
0002 #MASS_HIGGS = 125 # in GeV
0003 #MASS_HIGGS = 125 # in GeV
0004 #WIDTH_HIGGS   = 0.027*MASS_HIGGS # Same as default for id=35
0005 #MASS_X = 12 # in GeV
0006 #MASS_X = 12 # in GeV
0007 #CTAU_X = 900 # in mm
0008 #WIDTH_X = 2.1925219866666667e-16 # in mm
0009
0010 import FWCore.ParameterSet.Config as cms
0011
0012 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0013 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0014 from Configuration.Generator.PSweightsPythia.PythiaPSweightsSettings_cfi import *
0015
0016 generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
0017     pythiaPylistVerbosity = cms.untracked.int32(1),
0018     filterEfficiency = cms.untracked.double(1),
0019     pythiaHepMCVerbosity = cms.untracked.bool(False),
0020     comEnergy = cms.double(13600.),
0021     crossSection = cms.untracked.double(CROSS_SECTION),
0022     maxEventsToPrint = cms.untracked.int32(10),
0023     PythiaParameters = cms.PSet(
0024         pythia8CommonSettingsBlock,
0025 #        pythia8CUEP8M1SettingsBlock, # Old PYTHIA tune
0026         pythia8CP5SettingsBlock,
0027         pythia8PSweightsSettingsBlock,
0028         processParameters = cms.vstring(
0029             "Higgs:useBSM = on",
0030             "HiggsBSM:all = off",
0031 # Gluon-fusion production only
0032 #            "HiggsBSM:ffbar2H2 = on",
0033             "HiggsBSM:gg2H2 = on",
0034 # Long-lived scalar decaying to nunu
0035             "6000111:new = LL_nunu LLbar_nunu 1 0 0",
0036             "6000111:m0 = %s" % MASS_X,
0037             "6000111:mWidth = %s" % WIDTH_X,
0038             "6000111:tau0 = %s" % CTAU_X,
0039             "6000111:isResonance = on",
0040             "6000111:mayDecay = on",
0041             "6000111:oneChannel = 1  0.333 100 12 -12",
0042             "6000111:addChannel = 1  0.333 100 14 -14",
0043             "6000111:addChannel = 1  0.333 100 16 -16",
0044 # Long-lived scalar decaying only to bb
0045             "6000113:new = LL_b LLbar_b 1 0 0",
0046             "6000113:m0 = %s" % MASS_X,
0047             "6000113:mWidth = %s" % WIDTH_X,
0048             "6000113:tau0 = %s" % CTAU_X,
0049             "6000113:isResonance = on",
0050             "6000113:mayDecay = on",
0051             "6000113:oneChannel = 1 1.0 100 5 -5",
0052 # Shut down H0 decays to ordinary particles and A0
0053             "35:m0 = %s" % MASS_HIGGS,
0054             "35:mWidth = %s" % WIDTH_HIGGS,
0055             "35:2:bRatio = 0.0",
0056             "35:3:bRatio = 0.0",
0057             "35:4:bRatio = 0.0",
0058             "35:5:bRatio = 0.0",
0059             "35:7:bRatio = 0.0",
0060             "35:8:bRatio = 0.0",
0061 # Keep small coupling of H0 to gluons in order to be produced
0062 #            "35:9:bRatio = 0.0",
0063             "35:10:bRatio= 0.0",
0064             "35:12:bRatio= 0.0",
0065             "35:13:bRatio= 0.0",
0066             "35:15:bRatio= 0.0", # h0(25)h0(25); would be open at high m(H0)
0067             "35:18:bRatio= 0.0", # Z0(23)A0(36); would be open at high m(H0)
0068             "35:36:bRatio= 0.0", # A0(36)A0(36); would be open at high m(H0)
0069             "35:2:meMode = 100",
0070             "35:3:meMode = 100",
0071             "35:4:meMode = 100",
0072             "35:5:meMode = 100",
0073             "35:7:meMode = 100",
0074             "35:8:meMode = 100",
0075 #            "35:9:meMode = 100",
0076             "35:10:meMode= 100",
0077             "35:12:meMode= 100",
0078             "35:13:meMode= 100",
0079             "35:15:meMode= 100",
0080             "35:18:meMode= 100",
0081             "35:20:meMode= 100",
0082
0083 # Enable H0-->X(mumu)X(mumu) decay
0084             "35:addChannel = 1 1. 100 6000113 6000113",# Enable H0-->X(mumu)X(jetjet) decay
0085 #            "35:addChannel = 1 1. 100 6000111 6000113",
0086             "35:onMode = off",
0087 #`           "35:onIfAny = 6000111 6000113"
0088             "35:onIfAny = 6000113 6000113"
0089         ),
0090
0091         parameterSets = cms.vstring(
0092             'pythia8CommonSettings',
0093 #            'pythia8CUEP8M1Settings', #Old Pythia Tune
0094             'pythia8CP5Settings',
0095             'pythia8PSweightsSettings',
0096             'processParameters'
0097         )
0098     )
0099 )
0100
0101 ProductionFilterSequence = cms.Sequence(generator)