From 4b7fc51f0d56613926e719e4963fdc5d75cb1513 Mon Sep 17 00:00:00 2001
From: Christian Appelt <christian.appelt@cern.ch>
Date: Tue, 20 Sep 2022 09:02:10 +0200
Subject: [PATCH] Change HNL PIDs from 9900012, 9900014, 9900016 to 72, 74, 76
 according to PDG rule 11c.

---
 SM_HeavyN_NLO/parameters.py | 12 ++++++------
 SM_HeavyN_NLO/particles.py  |  6 +++---
 2 files changed, 9 insertions(+), 9 deletions(-)

diff --git a/SM_HeavyN_NLO/parameters.py b/SM_HeavyN_NLO/parameters.py
index 42e3261..855c12f 100644
--- a/SM_HeavyN_NLO/parameters.py
+++ b/SM_HeavyN_NLO/parameters.py
@@ -160,7 +160,7 @@ mN1 = Parameter(name = 'mN1',
                 value = 300.,
                 texname = '\\text{mN1}',
                 lhablock = 'MASS',
-                lhacode = [ 9900012 ])
+                lhacode = [ 72 ])
 
 mN2 = Parameter(name = 'mN2',
                 nature = 'external',
@@ -168,7 +168,7 @@ mN2 = Parameter(name = 'mN2',
                 value = 500.,
                 texname = '\\text{mN2}',
                 lhablock = 'MASS',
-                lhacode = [ 9900014 ])
+                lhacode = [ 74 ])
 
 mN3 = Parameter(name = 'mN3',
                 nature = 'external',
@@ -176,7 +176,7 @@ mN3 = Parameter(name = 'mN3',
                 value = 1000.,
                 texname = '\\text{mN3}',
                 lhablock = 'MASS',
-                lhacode = [ 9900016 ])
+                lhacode = [ 76 ])
 
 WZ = Parameter(name = 'WZ',
                nature = 'external',
@@ -216,7 +216,7 @@ WN1 = Parameter(name = 'WN1',
                 value = 0.303,
                 texname = '\\text{WN1}',
                 lhablock = 'DECAY',
-                lhacode = [ 9900012 ])
+                lhacode = [ 72 ])
 
 WN2 = Parameter(name = 'WN2',
                 nature = 'external',
@@ -224,7 +224,7 @@ WN2 = Parameter(name = 'WN2',
                 value = 1.5,
                 texname = '\\text{WN2}',
                 lhablock = 'DECAY',
-                lhacode = [ 9900014 ])
+                lhacode = [ 74 ])
 
 WN3 = Parameter(name = 'WN3',
                 nature = 'external',
@@ -232,7 +232,7 @@ WN3 = Parameter(name = 'WN3',
                 value = 12.3,
                 texname = '\\text{WN3}',
                 lhablock = 'DECAY',
-                lhacode = [ 9900016 ])
+                lhacode = [ 76 ])
 
 aEW = Parameter(name = 'aEW',
                 nature = 'internal',
diff --git a/SM_HeavyN_NLO/particles.py b/SM_HeavyN_NLO/particles.py
index 73ed05c..057421e 100644
--- a/SM_HeavyN_NLO/particles.py
+++ b/SM_HeavyN_NLO/particles.py
@@ -385,7 +385,7 @@ G__plus__ = Particle(pdg_code = 251,
 
 G__minus__ = G__plus__.anti()
 
-N1 = Particle(pdg_code = 9900012,
+N1 = Particle(pdg_code = 72,
               name = 'N1',
               antiname = 'N1',
               spin = 2,
@@ -399,7 +399,7 @@ N1 = Particle(pdg_code = 9900012,
               LeptonNumber = 0,
               Y = 0)
 
-N2 = Particle(pdg_code = 9900014,
+N2 = Particle(pdg_code = 74,
               name = 'N2',
               antiname = 'N2',
               spin = 2,
@@ -413,7 +413,7 @@ N2 = Particle(pdg_code = 9900014,
               LeptonNumber = 0,
               Y = 0)
 
-N3 = Particle(pdg_code = 9900016,
+N3 = Particle(pdg_code = 76,
               name = 'N3',
               antiname = 'N3',
               spin = 2,
-- 
GitLab