From 159ca648c1b6ae0a274d34f14bb5e87cbb2acf6a Mon Sep 17 00:00:00 2001
From: pierog <tanguy.pierog@kit.edu>
Date: Thu, 7 Oct 2021 15:40:23 +0200
Subject: [PATCH] add missing conex_extension par for pions and kaons
---
src/conex_mod.F | 49 ++++++++++++++++++++++++++++++++++++++++++++++++-
1 file changed, 48 insertions(+), 1 deletion(-)
diff --git a/src/conex_mod.F b/src/conex_mod.F
index 5a4ea78..bdbff36 100644
--- a/src/conex_mod.F
+++ b/src/conex_mod.F
@@ -49236,12 +49236,12 @@ C...Single nucleon (proton) case
#endif
CALL SIBYLL (13,IATARGET, ECM)
CALL DECSIB
- DO J=1,NP
c#ifdef CONEX_EXTENSIONS
c if (doResampling) then
c call resampleSIBYLL(NP, xselab, idprojxs, SQS, xsamproj)
c endif
c#endif
+ DO J=1,NP
LA = IABS(LLIST(J))
IF (LA .LT. 10000) THEN
NPA = NPA + 1
@@ -70106,6 +70106,30 @@ c taking slope parameterization from sigma_pp Donnie.-Landshoff
CALL SIG_RPP2014(2,1,SQS,SLOPE,SIGT,SIGEL,SIGINEL,RHO)
+#ifdef CONEX_EXTENSIONS
+
+* cross section as calculated in SIBYLL, but re-scaled (RU 02/2009)
+
+ ELSE IF(ICSPA.EQ.4) THEN
+
+ CALL SIB_SIGMA_HP(2,SQS,SIGT,SIGEL,SIGINEL,SIGDIF,SLOPE,RHO)
+ factMod=1.d0
+ npion=2
+ call modifiercx(factMod, e0*1.d3,npion)
+ sigt=factMod*sigt
+ sigdif(1)=factMod*sigdif(1)
+ sigdif(2)=factMod*sigdif(2)
+ sigdif(3)=factMod*sigdif(3)
+ siginel=factMod*siginel
+c rescaling variant A
+c sigel=factMod**2*sigel
+c rescaling variant B
+ sigel=factMod*sigel
+ slope=factMod*slope
+
+#endif
+c CONEX_EXTENSIONS
+
ENDIF
@@ -70261,6 +70285,29 @@ c taking slope parameterization from sigma_pp Donnie.-Landshoff
CALL SIG_RPP2014(3,1,SQS,SLOPE,SIGT,SIGEL,SIGINEL,RHO)
+ #ifdef CONEX_EXTENSIONS
+
+* cross section as calculated in SIBYLL, but re-scaled (RU 02/2009)
+
+ ELSE IF(ICSPA.EQ.4) THEN
+
+ CALL SIB_SIGMA_HP(3,SQS,SIGT,SIGEL,SIGINEL,SIGDIF,SLOPE,RHO)
+ factMod=1.d0
+ nkaon=4
+ call modifiercx(factMod, e0*1.d3, nkaon)
+ sigt=factMod*sigt
+ sigdif(1)=factMod*sigdif(1)
+ sigdif(2)=factMod*sigdif(2)
+ sigdif(3)=factMod*sigdif(3)
+ siginel=factMod*siginel
+c rescaling variant A
+c sigel=factMod**2*sigel
+c rescaling variant B
+ sigel=factMod*sigel
+ slope=factMod*slope
+
+#endif
+
ENDIF
RETURN
--
GitLab