Changes between Version 13 and Version 14 of pSPSS

Timestamp:

Oct 18, 2022, 10:30:28 AM (2 years ago)

Author:

Jan Hajer

Comment:

--

pSPSS

@@ -3 +3 @@
 The motivation and implementation of the `pSPSS` is discussed in 
 
-S. Antusch, J. Hajer, and J. Rosskopp. ‘Simulating heavy neutrino-antineutrino oscillations at colliders’ (2022).
+S. Antusch, J. Hajer, and J. Rosskopp. ‘Simulating lepton number violation caused by heavy neutrino-antineutrino oscillations at colliders’ (2022).
 
 == Model description ==
@@ -39 +39 @@
 
 {{{
+mass_splitting = param_card.get_value('PSPSS', 2)
+damping = param_card.get_value('PSPSS', 6)
 for event in lhe:
     leptonnumber = 0
@@ -51 +53 @@
         id = particle.pid
         width = param_card['decay'].get((abs(id),)).value
-        mass splitting = param_card.get_value('PSPSS', 2)
         if width:
             if id in [8000011, 8000012]:
                 tau0 = random.expovariate(width / cst)
-                if 0.5 * (1 + math.cos(mass splitting * tau0 / cst)) >= random.random():
+                if 0.5 * (1 + math.exp(=damping)*math.cos(mass_splitting * tau0 / cst)) >= random.random():
                     write_event = (leptonnumber == 0)
                 else: