IAP GITLAB

Commit f13ecdcb authored by Tanguy Pierog's avatar Tanguy Pierog

correction in the interface (no effect ?)


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@6249 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent d58187e4
......@@ -7,7 +7,7 @@ C SSSSSS IIIIIII BBBBB YY LLLLLLL LLLLLLL
C=======================================================================
C Code for SIBYLL: hadronic interaction Monte Carlo event generator
C=======================================================================
C Version 2.3c (Jun-01-2017)
C Version 2.3c01 (Jun-01-2017, modified Aug-06-2017)
C
C with CHARM production
C
......@@ -468,7 +468,7 @@ C-----------------------------------------------------------------------
* /,' ','| F. RIEHN et al., Proc. 34th Int. Cosmic Ray Conf.|',
* /,' ','| The Hague, The Netherlands, cont. 1313 (2015) |',
* /,' ','| |',
* /,' ','| last modifications: F. Riehn (06/01/2017) |',
* /,' ','| last modifications: R. Engel (08/06/2017) |',
* /,' ','====================================================',
* /)
......@@ -8406,8 +8406,10 @@ C--------------------------------------------------------------------
c internal types
INTEGER LL,LCON,LRES,LRES1,NTRYS,NRJECT,LA,N1,IREJ,I,J,IFLA,
& IFL1,IFL2,IFBAD,NPI,IRES,LA1,JQQ,JQTOT,K,JQR
DOUBLE PRECISION PD,BE,EMIN,EMIN2,PCHEX,PRES,DELTAE,SQS_0,
& IFL1,IFL2,IFBAD,NPI,IRES,LA1,JQQ,JQTOT,K,JQR,
& KB_0,IAT_0
DOUBLE PRECISION PD,BE,EMIN,EMIN2,PCHEX,PRES,DELTAE,
& SQS_0,S_0,PTmin_0,XMIN_0,ZMIN_0,
& PAR1_def,PAR24_def,PAR53_def,GA,BEP,S_RNDM,AV,GASDEV,PCXG,
& XI1,XI2,XSMR !,FERMI
DIMENSION LL(10), PD(10,5), BE(3), LCON(6:99),LRES1(6:99)
......@@ -8475,11 +8477,23 @@ C... pomeron-hadron scattering (pi0 is used instead of pomeron)
IPFLAG= IPFLAG*100
c create subevent
SQS_0 = SQS
S_0 = S
PTmin_0 = PTmin
XMIN_0 = XMIN
ZMIN_0 = ZMIN
KB_0 = KB
IAT_0 = IAT
CALL INI_EVENT(P0(5),L0,6,0)
c create L0 - pi0 interaction, pi0(pid=6) target
CALL SIB_NDIFF(L0, 1, P0(5), 0, IREJ) ! ori
c restore main event
SQS = SQS_0
S = S_0
PTmin = PTmin_0
XMIN = XMIN_0
ZMIN = ZMIN_0
KB = KB_0
IAT = IAT_0
IF(IREJ.NE.0) THEN
NP = N1-1
GOTO 50
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment