IAP GITLAB

Commit 96ff6faf authored by Felix Riehn's avatar Felix Riehn

updated sibyll2.3c to version from corsika release

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@6430 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent f13ecdcb
...@@ -7,7 +7,7 @@ C SSSSSS IIIIIII BBBBB YY LLLLLLL LLLLLLL ...@@ -7,7 +7,7 @@ C SSSSSS IIIIIII BBBBB YY LLLLLLL LLLLLLL
C======================================================================= C=======================================================================
C Code for SIBYLL: hadronic interaction Monte Carlo event generator C Code for SIBYLL: hadronic interaction Monte Carlo event generator
C======================================================================= C=======================================================================
C Version 2.3c01 (Jun-01-2017, modified Aug-06-2017) C Version 2.3c01 (Jun-01-2017, modified Sept-05-2017)
C C
C with CHARM production C with CHARM production
C C
...@@ -33,6 +33,7 @@ C stanev@bartol.udel.edu ...@@ -33,6 +33,7 @@ C stanev@bartol.udel.edu
C C
C differences to Sibyll 2.3 include: C differences to Sibyll 2.3 include:
C C
C *extend eta' and phi decay to include prompt muons
C *allow initialization from file C *allow initialization from file
C *explicit hyperon production in 2-particle fireballs C *explicit hyperon production in 2-particle fireballs
C *forbid baryon pair production immediately next to leading baryon C *forbid baryon pair production immediately next to leading baryon
...@@ -465,10 +466,10 @@ C----------------------------------------------------------------------- ...@@ -465,10 +466,10 @@ C-----------------------------------------------------------------------
* /,' ','| |', * /,' ','| |',
* /,' ','| Publication to be cited when using this program: |', * /,' ','| Publication to be cited when using this program: |',
* /,' ','| Eun-Joo AHN et al., Phys.Rev. D80 (2009) 094003 |', * /,' ','| Eun-Joo AHN et al., Phys.Rev. D80 (2009) 094003 |',
* /,' ','| F. RIEHN et al., Proc. 34th Int. Cosmic Ray Conf.|', * /,' ','| F. RIEHN et al., Proc. 35th Int. Cosmic Ray Conf.|',
* /,' ','| The Hague, The Netherlands, cont. 1313 (2015) |', * /,' ','| Bexco, Busan, Korea, cont. 301 (2017) |',
* /,' ','| |', * /,' ','| |',
* /,' ','| last modifications: R. Engel (08/06/2017) |', * /,' ','| last modifications: F. Riehn (09/18/2017) |',
* /,' ','====================================================', * /,' ','====================================================',
* /) * /)
...@@ -1446,9 +1447,11 @@ C---------------------------------------------------------------- ...@@ -1446,9 +1447,11 @@ C----------------------------------------------------------------
DO J=1,99 DO J=1,99
IA = ISIB_PID2PDG( j ) IA = ISIB_PID2PDG( j )
IF(IA.ne.0)THEN IF(IA.ne.0)THEN
ISIBPDG2PIDIA=ISIB_PDG2PID( IA ) ISIBPDG2PIDIA=ISIB_PDG2PID( IA )
WRITE (LUN,120) NAMP(J), J, IA, ISIBPDG2PIDIA, AM(J) ELSE
WRITE(LUN,*) 'PDG conversion not found!'
ENDIF ENDIF
WRITE (LUN,120) NAMP(J), J, IA, ISIBPDG2PIDIA, AM(J)
ENDDO ENDDO
120 FORMAT(4X,A6,4X,I4,7X,I7,8X,I4,5X,F9.3) 120 FORMAT(4X,A6,4X,I4,7X,I7,8X,I4,5X,F9.3)
...@@ -6747,14 +6750,25 @@ C----------------------------------------------------------------------- ...@@ -6747,14 +6750,25 @@ C-----------------------------------------------------------------------
CHARACTER*6 NAMP CHARACTER*6 NAMP
COMMON /S_CNAM/ NAMP (0:99) COMMON /S_CNAM/ NAMP (0:99)
SAVE SAVE
c CBR contains the normed sum of the branching ratios of the decay channels
c indexed by IDB, i.e. a particle with 4 decay channels will have the entries
c [B1/Btot, (B1+B2)/Btot, (B1+B2+B3)/Btot, 1.]
DATA CBR /3*1.D0,0.D0,1.D0,1.D0,0.6354D0,0.8422D0,0.8981D0, DATA CBR /3*1.D0,0.D0,1.D0,1.D0,0.6354D0,0.8422D0,0.8981D0,
+ 0.9157D0,0.9492D0,1.D0,0.6354D0,0.8422D0,0.8981D0,0.9157D0, + 0.9157D0,0.9492D0,1.D0,0.6354D0,0.8422D0,0.8981D0,0.9157D0,
+ 0.9492D0,1.D0,0.1965D0,0.3224D0,0.4579D0,0.5934D0,0.7967D0,1.D0, + 0.9492D0,1.D0,0.1965D0,0.3224D0,0.4579D0,0.5934D0,0.7967D0,1.D0,
+ 0.6925D0,1.D0,3*0.D0,0.5D0,1.D0,0.5D0,1.D0,0.3941D0,0.7197D0, + 0.6925D0,1.D0,3*0.D0,0.5D0,1.D0,0.5D0,1.D0,
+ 0.9470D0,0.9930D0,1.D0,0.D0,0.4460D0,0.6530D0,0.9470D0,0.9770D0, + 0.3941D0,0.7197D0,0.9470D0,0.9930D0,1.D0, ! eta
+ 0.9980D0,4*1.D0,0.6670D0,1.D0,9*0.D0,0.6670D0,1.D0,0.6670D0,1.D0, + 0.4285D0,0.7193D0,0.9487D0,0.9750D0,0.9973D0,0.9999D0,1.D0, ! eta'
+ 0.6670D0,1.D0,0.8940D0,0.9830D0,1.D0,0.4930D0,0.8340D0,0.9870D0, + 3*1.D0, ! rho-mesons
+ 1.D0,0.5160D0,5*1.D0,0.6410D0,2*1.D0,0.67D0,1.D0,0.33D0,2*1.D0, + 0.6670D0,1.D0, ! K*+
+ 0.4894D0,0.8317D0,0.9850D0,0.9981D0,0.9994D0,0.9997D0,1.D0, ! phi(1020)
+ 2*0.D0, ! (empty)
+ 0.6670D0,1.D0, ! K*-
+ 0.6670D0,1.D0, ! K*0
+ 0.6670D0,1.D0, ! K*0 bar
+ 0.8940D0,0.9830D0,1.D0, ! omega
+ 4*0.D0, ! (empty)
+ 0.5160D0,5*1.D0,0.6410D0,2*1.D0,0.67D0,1.D0,0.33D0,2*1.D0,
+ 0.88D0,0.94D0,1.D0,0.88D0,0.94D0,1.D0,0.88D0,0.94D0,1.D0,0.33D0, + 0.88D0,0.94D0,1.D0,0.88D0,0.94D0,1.D0,0.88D0,0.94D0,1.D0,0.33D0,
+ 1.D0,0.67D0,1.D0,0.678D0,0.914D0,1.D0,0.217D0,0.398D0,0.506D0, + 1.D0,0.67D0,1.D0,0.678D0,0.914D0,1.D0,0.217D0,0.398D0,0.506D0,
+ 0.595D0,0.684D0,0.768D0,0.852D0,0.923D0,0.976D0,1.D0,0.217D0, + 0.595D0,0.684D0,0.768D0,0.852D0,0.923D0,0.976D0,1.D0,0.217D0,
...@@ -6815,13 +6829,22 @@ C----------------------------------------------------------------------- ...@@ -6815,13 +6829,22 @@ C-----------------------------------------------------------------------
& 0.001295D0,0.00155D0,8.281D-05,9.801D-05,0.D0,0.D0,0.09D0, & 0.001295D0,0.00155D0,8.281D-05,9.801D-05,0.D0,0.D0,0.09D0,
& 0.01D0,0.09D0,0.01D0,6*0.D0,0.01D0,0.01D0,0.01D0,4*0.0729D0, & 0.01D0,0.09D0,0.01D0,6*0.D0,0.01D0,0.01D0,0.01D0,4*0.0729D0,
& 32*0.D0/ & 32*0.D0/
c IDB is the index to the branching ratios (CBR) and decay channels (KDEC).
c always indicates the first decay channel
DATA IDB / DATA IDB /
+ 0,0,0,1,2,3,5,6,7,13,19,25,8*0,30,32,34,40,46,47,48,49,60,62, + 0,0,0,1,2, ! leptons
+ 64,66,69,73,75,76,77,78,79,81,82,84,86,87,90,93,96,98,100, + 3,5,6,7,13,19,25, ! pions and kaons
+ 0,224,228,232,239,4*0, ! < 59 + 8*0,30,32,34,39,46,47,48,49,60,62,64,66,51, !69, ! meson resonances
+ 103,113,246,248,250, 252,254,256,258,3*0, + 73,75,76,77,78,79,81,82,84,86,87,90,93,96,98,100, ! baryons : Sibyll 2.1
+ 123,134,145,204,214,200,202,151,154,157,159,0, + 0,224,228,232,239,4*0, ! Nucleon resonaces
+ 161,164,165,166,167,175,179,4*0,189,190,191,192,194,196 / + 103,113,246,248,250, 252,254,256,258,3*0,
+ 123,134,145,204,214,200,202,151,154,157,159,0,
+ 161,164,165,166,167,175,179,4*0,189,190,191,192,194,196 /
c KDEC contains decay channels, format is [ND, MAT, LL(1:4)]
c where ND is the number of particles in the final state (max 4)
C MAT is 0, 1 for semi-leptonic (weak decay) or not
c (adds primitive matrix element)
c LL(1:4) are the particle ids of the final state particles
DATA KDEC / DATA KDEC /
+ 3,1,15,2,18,0,3,1,16,3,17,0,2,0,1,1,8*0,2,0,4,17,0,0,2,0,5,18,0, + 3,1,15,2,18,0,3,1,16,3,17,0,2,0,1,1,8*0,2,0,4,17,0,0,2,0,5,18,0,
+ 0,2,0,4,17,0,0,2,0,7,6,0,0,3,0,7,7,8,0,3,0,7,6,6,0,3,1,17,4,6,0, + 0,2,0,4,17,0,0,2,0,7,6,0,0,3,0,7,7,8,0,3,0,7,6,6,0,3,1,17,4,6,0,
...@@ -6829,13 +6852,22 @@ C----------------------------------------------------------------------- ...@@ -6829,13 +6852,22 @@ C-----------------------------------------------------------------------
+ 1,18,5,6,0,3,1,16,3,6,0,3,0,6,6,6,0,3,0,7,8,6,0,3,1,18,5,7,0,3, + 1,18,5,6,0,3,1,16,3,6,0,3,0,6,6,6,0,3,0,7,8,6,0,3,1,18,5,7,0,3,
+ 1,17,4,8,0,3,1,16,3,7,0,3,1,15,2,8,0,2,0,7,8,0,0,2,0,6,6,20*0,1, + 1,17,4,8,0,3,1,16,3,7,0,3,1,15,2,8,0,2,0,7,8,0,0,2,0,6,6,20*0,1,
+ 0,11,3*0,1,0,12,0,0,0,1,0,11,0,0,0,1,0,12,0,0,0,2,0,1,1,0,0,3,0, + 0,11,3*0,1,0,12,0,0,0,1,0,11,0,0,0,1,0,12,0,0,0,2,0,1,1,0,0,3,0,
+ 6,6,6,0,3,0,7,8,6,0,3,0,1,7,8,0,3,0,1,3,2,7*0,3,0,7,8,23,0,3,0,6 + 6,6,6,0,3,0,7,8,6,0,3,0,1,7,8,0,3,0,1,3,2,0,
+ ,6,23,0,2,0,1,27,0,0,2,0,1,32,0,0,2,0,1,1,0,0,3,0,6,6,6,0,2,0,7, + 3,0,7,8,23,0, 3,0,6,6,23,0, 2,0,1,27,0,0, 2,0,1,32,0,0, ! eta'
+ 6,0,0,2,0,8,6,0,0,2,0,7,8,0,0,2,0,21,7,0,0,2,0,9,6,0,0,54*0,2,0, + 2,0,1,1,0,0, 3,0,6,6,6,0, 3,0,1,4,5,0, ! eta'
+ 22,8,0,0,2,0,10,6,0,0,2,0,9,8,0,0,2,0,21,6,0,0,2,0,10,7,0,0, + 2,0,7,6,0,0, ! rho+
+ 2,0,22,6,0,0,3,0,7,8,6,0,2,0,1,6,0,0,2,0,7,8,0,0,2,0,9,10,0, + 2,0,8,6,0,0, ! rho-
+ 0,2,0,11,12,0,0,3,0,7, + 2,0,7,8,0,0, ! rho0
+ 8,6,0,2,0,1,23,0,0,2,0,13,6,0,0,2,0,14,7,0,0,2,0,39,1,0,0,2, + 2,0,21,7,0,0, 2,0,9,6,0,0, ! K*+
+ 2,0,9,10,0,0, 2,0,11,12,0,0, 3,0,7,8,6,0, 2,0,1,23,0,0, ! phi(1020)
+ 2,0,1,6,0,0, 2,0,2,3,0,0, 2,0,4,5,0,0, ! phi(1020)
+ 12*0,
+ 2,0,22,8,0,0, 2,0,10,6,0,0, ! K*-
+ 2,0,9,8,0,0, 2,0,21,6,0,0, ! K*0
+ 2,0,10,7,0,0, 2,0,22,6,0,0, ! K*0 bar
+ 3,0,7,8,6,0, 2,0,1,6,0,0, 2,0,7,8,0,0, ! omega
+ 24*0,
+ 2,0,13,6,0,0,2,0,14,7,0,0,2,0,39,1,0,0,2, ! baryons
+ 0,14,8,0,0,2,0,39,6,0,0,2,0,39,8,0,0,2,0,13,8,0,0,2,0, + 0,14,8,0,0,2,0,39,6,0,0,2,0,39,8,0,0,2,0,13,8,0,0,2,0,
+ 14,6,0,0,2,0,13,7,0,0,2,0,13,6, + 14,6,0,0,2,0,13,7,0,0,2,0,13,6,
+ 0,0,2,0,14,7,0,0,2,0,13,8,0,0,2,0,14,6,0,0,2,0,14,8,0,0,2,0, + 0,0,2,0,14,7,0,0,2,0,13,8,0,0,2,0,14,6,0,0,2,0,14,8,0,0,2,0,
......
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