diff --git a/src/conex_mod.F b/src/conex_mod.F
index 5a4ea78f7e89d5d4fb48d1dfe406f6541c303ace..bdbff367a16ee8bf9becfd1277111fe6575ca5b0 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