Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
atlas-physics
pmg
MC Job Options
Commits
acce6453
Commit
acce6453
authored
Nov 11, 2021
by
lderamo
Browse files
Adding new RPV JO names [skip modfiles]
parent
18cc3a94
Pipeline
#3229277
failed with stages
in 53 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
506xxx/506179/MadGraphControl_SimplifiedModel_RPV_LLP.py
View file @
acce6453
...
...
@@ -21,37 +21,37 @@ print gentype, decaytype, mass_stopgluino, masses, decayflavor, lifetimeString,
#Decay flavour
decays
[
'1000006'
]
=
"""DECAY 1000006 5.69449678E+00 # stop1 decays
# BR NDA ID1 ID2
1.000 2 1000022 6
1.000 2 1000022 6
"""
if
decayflavor
==
"rpvHF"
:
if
decayflavor
==
"rpvHF"
:
decays
[
'1000021'
]
=
"""DECAY 1000021 7.40992706E-02 # gluino decays
# BR NDA ID1 ID2 ID3
1.000 3 1000022 6 -6
1.000 3 1000022 6 -6
"""
elif
decayflavor
==
"rpvLF"
:
decays
[
'1000021'
]
=
"""DECAY 1000021 7.40992706E-02 # gluino decays
# BR NDA ID1 ID2 ID3
0.250 3 1000022 1 -1
0.250 3 1000022 2 -2
0.250 3 1000022 3 -3
0.250 3 1000022 4 -4
0.250 3 1000022 1 -1
0.250 3 1000022 2 -2
0.250 3 1000022 3 -3
0.250 3 1000022 4 -4
"""
elif
decayflavor
==
"rpvdem"
:
decays
[
'1000021'
]
=
"""DECAY 1000021 7.40992706E-02 # gluino decays
# BR NDA ID1 ID2 ID3
0.166 3 1000022 1 -1
0.166 3 1000022 2 -2
0.166 3 1000022 3 -3
0.166 3 1000022 4 -4
0.166 3 1000022 5 -5
0.166 3 1000022 6 -6
0.166 3 1000022 1 -1
0.166 3 1000022 2 -2
0.166 3 1000022 3 -3
0.166 3 1000022 4 -4
0.166 3 1000022 5 -5
0.166 3 1000022 6 -6
"""
decayWidthStr
=
'%e'
%
decayWidth
decayStringHeader
=
'DECAY 1000022 '
header
=
decayStringHeader
+
decayWidthStr
header
=
decayStringHeader
+
decayWidthStr
evgenLog
.
info
(
'lifetime of 1000022 is set to %s ns'
%
neutralinoLifetime
)
if
decayflavor
==
"rpvHF"
:
branchingRatios
=
'''
...
...
@@ -69,7 +69,7 @@ elif decayflavor == "rpvLF":
#'''
elif
decayflavor
==
"rpvdem"
:
branchingRatios
=
'''
# BR NDA ID1 ID2 ID3
# BR NDA ID1 ID2 ID3
0.05555556 3 2 1 3
0.05555556 3 2 1 5
0.05555556 3 2 3 5
...
...
@@ -117,7 +117,7 @@ elif gentype == "TT":
'''
njets
=
2
evgenConfig
.
contact
=
[
"rebecca.carney@cern.ch"
,
"jmontejo@cern.ch"
]
evgenConfig
.
contact
=
[
"rebecca.carney@cern.ch"
,
"jmontejo@cern.ch"
,
"lderamo@cern.ch"
]
if
gentype
==
"GG"
:
evgenLog
.
info
(
'Registered generation of gluino grid'
)
...
...
@@ -137,7 +137,7 @@ run_settings['event_norm']='average' # see AGENE-1725
if
njets
>
0
:
evgenLog
.
info
(
'Using Merging:Process = guess'
)
genSeq
.
Pythia8
.
Commands
+=
[
"Merging:Process = guess"
]
genSeq
.
Pythia8
.
Commands
+=
[
"Merging:Process = guess"
]
genSeq
.
Pythia8
.
UserHooks
+=
[
'JetMergingaMCatNLO'
]
testSeq
.
TestHepMC
.
MaxVtxDisp
=
1e8
#in mm
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment