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

Susy gen fragment for trigger validation #35203

Merged
merged 1 commit into from Sep 9, 2021
Merged

Susy gen fragment for trigger validation #35203

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
276 changes: 276 additions & 0 deletions ...n/Generator/python/SMS-Higgsino_mN2-170_mC1-160_mN1-150_HT60_TuneCP5_13TeV_pythia8_cfi.py
View file
Open in desktop
@@ -0,0 +1,276 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *

SLHA_TABLE = """
Block MODSEL # Model selection
1 0 # Generic MSSM
Block MINPAR # Input parameters
3 1.00000000E+00 # tanb at m_Z
#
Block SMINPUTS # SM parameters
1 1.27931277E+02 # alpha_em^-1(MZ)^MSbar
2 1.16639000E-05 # G_mu [GeV^-2]
3 1.17200000E-01 # alpha_s(MZ)^MSbar
4 9.11876000E+01 # m_Z(pole)
5 4.20000000E+00 # m_b(m_b), MSbar
6 1.74300000E+02 # m_t(pole)
7 1.77700000E+00 # m_tau(pole)
Block MASS # Mass spectrum
# PDG code mass particle
24 8.04009772E+01 # W+
25 10.0E+4 # h0
35 10.0E+4 # H0
36 10.0E+4 # A0
37 10.0E+4 # H+
1000001 10.0E+4 # ~d_L
2000001 10.0E+4 # ~d_R
1000002 10.0E+4 # ~u_L
2000002 10.0E+4 # ~u_R
1000003 10.0E+4 # ~s_L
2000003 10.0E+4 # ~s_R
1000004 10.0E+4 # ~c_L
2000004 10.0E+4 # ~c_R
1000005 10.0E+4 # ~b_1
2000005 10.0E+4 # ~b_2
1000006 10.0E+4 # ~t_1
2000006 10.0E+4 # ~t_2
1000011 10.0E+4 # ~e_L-
2000011 10.0E+4 # ~e_R-
1000012 10.0E+4 # ~nu_eL
1000013 10.0E+4 # ~mu_L-
2000013 10.0E+4 # ~mu_R-
1000014 10.0E+4 # ~nu_muL
1000015 10.0E+4 # ~tau_1-
2000015 10.0E+4 # ~tau_2-
1000016 10.0E+4 # ~nu_tauL
1000021 10.0E+4 # ~g
1000022 150 # ~chi_10
1000023 170 # ~chi_20
1000025 10.0E+4 # ~chi_30
1000035 10.0E+4 # ~chi_40
1000024 160 # ~chi_1+
1000037 10.0E+4 # ~chi_2+
#
BLOCK NMIX # Neutralino Mixing Matrix
1 1 0.00000000E+00 # N_11
1 2 0.00000000E+00 # N_12
1 3 0.70710678E+00 # N_13
1 4 0.70710678E+00 # N_14
2 1 0.00000000E+00 # N_21
2 2 0.00000000E+00 # N_22
2 3 0.70710678E+00 # N_23
2 4 -0.70710678E+00 # N_24
3 1 0.00000000E+00 # N_31
3 2 1.00000000E+00 # N_32
3 3 0.00000000E+00 # N_33
3 4 0.00000000E+00 # N_34
4 1 1.00000000E+00 # N_41
4 2 0.00000000E+00 # N_42
4 3 0.00000000E+00 # N_43
4 4 0.00000000E+00 # N_44
#
BLOCK UMIX # Chargino Mixing Matrix U
1 1 0.00000000E+00 # U_11
1 2 1.00000000E+00 # U_12
2 1 1.00000000E+00 # U_21
2 2 0.00000000E+00 # U_22
#
BLOCK VMIX # Chargino Mixing Matrix V
1 1 0.00000000E+00 # V_11
1 2 1.00000000E+00 # V_12
2 1 1.00000000E+00 # V_21
2 2 0.00000000E+00 # V_22
#
BLOCK STOPMIX # Stop Mixing Matrix
1 1 1.0 # cos(theta_t)
1 2 0.0 # sin(theta_t)
2 1 0.0 # -sin(theta_t)
2 2 1.0 # cos(theta_t)
#
BLOCK SBOTMIX # Sbottom Mixing Matrix
1 1 9.66726392E-01 # cos(theta_b)
1 2 2.55812594E-01 # sin(theta_b)
2 1 -2.55812594E-01 # -sin(theta_b)
2 2 9.66726392E-01 # cos(theta_b)
#
BLOCK STAUMIX # Stau Mixing Matrix
1 1 4.51419848E-01 # cos(theta_tau)
1 2 8.92311672E-01 # sin(theta_tau)
2 1 -8.92311672E-01 # -sin(theta_tau)
2 2 4.51419848E-01 # cos(theta_tau)
#
BLOCK ALPHA # Higgs mixing
-1.13676047E-01 # Mixing angle in the neutral Higgs boson sector
#
BLOCK HMIX Q= 2.90528802E+02 # DRbar Higgs Parameters
1 3.05599351E+02 # mu(Q)MSSM
#
BLOCK AU Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_u(Q) DRbar
2 2 0.00000000E+00 # A_c(Q) DRbar
3 3 -4.46245994E+02 # A_t(Q) DRbar
#
BLOCK AD Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_d(Q) DRbar
2 2 0.00000000E+00 # A_s(Q) DRbar
3 3 -8.28806503E+02 # A_b(Q) DRbar
#
BLOCK AE Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_e(Q) DRbar
2 2 0.00000000E+00 # A_mu(Q) DRbar
3 3 -4.92306701E+02 # A_tau(Q) DRbar
#
BLOCK MSOFT Q= 2.90528802E+02 # The soft SUSY breaking masses at the scale Q
1 6.39136864E+01 # M_1(Q)
2 1.22006983E+02 # M_2(Q)
3 3.90619532E+02 # M_3(Q)
21 4.42860395E+04 # mH1^2(Q)
22 -9.76585434E+04 # mH2^2(Q)
31 2.26648170E+02 # meL(Q)
32 2.26648170E+02 # mmuL(Q)
33 2.24355944E+02 # mtauL(Q)
34 2.08394096E+02 # meR(Q)
35 2.08394096E+02 # mmuR(Q)
36 2.03337218E+02 # mtauR(Q)
41 4.08594291E+02 # mqL1(Q)
42 4.08594291E+02 # mqL2(Q)
43 3.46134575E+02 # mqL3(Q)
44 3.98943379E+02 # muR(Q)
45 3.98943379E+02 # mcR(Q)
46 2.58021672E+02 # mtR(Q)
47 3.95211849E+02 # mdR(Q)
48 3.95211849E+02 # msR(Q)
49 3.90320031E+02 # mbR(Q)
#
#
#
# =================
# |The decay table|
# =================
#
# - The QCD corrections to the decays gluino -> squark + quark
# squark -> gaugino + quark_prime
# squark -> squark_prime + Higgs
# squark -> gluino + quark
# are included.
#
# - The multi-body decays for the inos, stops and sbottoms are included.
#
# - The loop induced decays for the gluino, neutralinos and stops
# are included.
#
# - The SUSY decays of the top quark are included.
#
#
#
# PDG Width
DECAY 1000022 0.00000000E+00 # neutralino1 decays
DECAY 1000021 1.00000000E+00 # gluino decays
# BR NDA ID1 ID2
0.000000E+00 3 1000022 1 -1
0.000000E+00 3 1000022 2 -2
0.000000E+00 3 1000022 3 -3
0.000000E+00 3 1000022 4 -4
0.000000E+00 3 1000022 5 -5
1.000000E+00 3 1000022 6 -6
#
# PDG Width
DECAY 1000006 0.00000000E+00 # stop1 decays
DECAY 2000006 0.00000000E+00 # stop2 decays
DECAY 1000005 0.00000000E+00 # sbottom1 decays
DECAY 2000005 0.00000000E+00 # sbottom2 decays
#
# PDG Width
DECAY 1000002 1.00000000E+00 # sup_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 2 # BR(~u_L -> ~chi_10 u)
#
# PDG Width
DECAY 2000002 1.00000000E+00 # sup_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 2 # BR(~u_R -> ~chi_10 u)
#
# PDG Width
DECAY 1000001 1.00000000E+00 # sdown_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 1 # BR(~d_L -> ~chi_10 d)
#
# PDG Width
DECAY 2000001 1.00000000E+00 # sdown_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 1 # BR(~d_R -> ~chi_10 d)
#
# PDG Width
DECAY 1000004 1.00000000E+00 # scharm_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 4 # BR(~c_L -> ~chi_10 c)
#
# PDG Width
DECAY 2000004 1.00000000E+00 # scharm_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 4 # BR(~c_R -> ~chi_10 c)
#
# PDG Width
DECAY 1000003 1.00000000E+00 # sstrange_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 3 # BR(~s_L -> ~chi_10 s)
#
# PDG Width
DECAY 2000003 1.00000000E+00 # sstrange_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 3 # BR(~s_R -> ~chi_10 s)
#
# PDG Width
DECAY 1000011 0.00000000E+00 # selectron_L decays
DECAY 2000011 0.00000000E+00 # selectron_R decays
DECAY 1000013 0.00000000E+00 # smuon_L decays
DECAY 2000013 0.00000000E+00 # smuon_R decays
DECAY 1000015 0.00000000E+00 # stau_1 decays
DECAY 2000015 0.00000000E+00 # stau_2 decays
#
# PDG Width
DECAY 1000012 0.00000000E+00 # snu_elL decays
DECAY 1000014 0.00000000E+00 # snu_muL decays
DECAY 1000016 0.00000000E+00 # snu_tauL decays
# PDG Width
DECAY 1000024 2.02592183E-05 # chargino1+ decays
# BR NDA ID1 ID2 ID3
0.33333333E+00 3 1000022 -11 12 # BR(~chi_1+ -> ~chi_10 e+ nu_e)
0.33333333E+00 3 1000022 -13 14 # BR(~chi_1+ -> ~chi_10 mu+ nu_mu)
0.33333333E+00 3 1000022 -15 16 # BR(~chi_1+ -> ~chi_10 tau+ nu_tau)
# PDG Width
DECAY 1000037 0.00000000E+00 # chargino2+ decays
DECAY 1000023 1.00000000E-05 # neutralino2 decays
0.33333333E+00 3 1000022 -11 11 # BR(~chi_1+ -> ~chi_10 e+ nu_e)
0.33333333E+00 3 1000022 -13 13 # BR(~chi_1+ -> ~chi_10 mu+ nu_mu)
0.33333333E+00 3 1000022 -15 15 # BR(~chi_1+ -> ~chi_10 tau+ nu_tau)
DECAY 1000025 0.00000000E+00 # neutralino3 decays
DECAY 1000035 0.00000000E+00 # neutralino4 decays
"""

generator = cms.EDFilter("Pythia8GeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(13000),
SLHATableForPythia8 = cms.string(SLHA_TABLE),
maxEventsToPrint = cms.untracked.int32(0),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'SUSY:all = off',
'SUSY:qqbar2chi0chi0 = on',
'SUSY:qqbar2chi+-chi0 = on',
'SUSY:qqbar2chi+chi- = on',
'PhaseSpace:pTHatMin = 60.',
),
parameterSets = cms.vstring(
'pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters'
)
)
)