IAP GITLAB

Commit 5341caf3 authored by Xavier Garrido's avatar Xavier Garrido

add sibyll, qgsjet, gheisha ... directories

git-svn-id: https://devel-ik.fzk.de/svn/mc/epos/branches/root@998 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 177a66fa
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
Command for complie:
c f77 -o hijing aaa.f hijing.f hipyset1.35.f
program main
CHARACTER FRAME*8, PROJ*8, TARG*8
DIMENSION DNDPT(50),DNDY(50)
COMMON/HIPARNT/HIPR1(100), IHPR2(50), HINT1(100), IHNT2(50)
C....information of produced particles:
COMMON/HIMAIN1/NATT, EATT, JATT, NT, NP, N0, N01, N10, N11
COMMON/HIMAIN2/KATT(130000,4), PATT(130000,4)
C....information of produced partons:
COMMON/HIJJET1/NPJ(300), KFPJ(300,500), PJPX(300,500),
& PJPY(300,500), PJPZ(300,500), PJPE(300,500), PJPM(300,500),
& NTJ(300), KFTJ(300,500), PJTX(300,500), PJTY(300,500),
& PJTZ(300,500), PJTE(300,500), PJTM(300,500)
COMMON/HIJJET2/NSG, NJSG(900), IASG(900,3), K1SG(900,100),
& K2SG(900,100), PXSG(900,100), PYSG(900,100), PZSG(900,100),
& PESG(900,100), PMSG(900,100)
COMMON/HISTRNG/NFP(300,15), PP(300,15), NFT(300,15), PT(300,15)
C....initialize HIJING for PROJ+TARG collisions at c.m. energy of 20 GeV:
EFRM=20.0
FRAME='CMS'
PROJ='P'
TARG='A'
IAP=1
IZP=1
IAT=197
IZT=79
CALL HIJSET (EFRM, FRAME, PROJ, TARG, IAP, IZP, IAT, IZT)
C....generating 10 central events:
N_EVENT=2
BMIN=0.0
print*,"HIPR1(34),HIPR1(35)",HIPR1(34),HIPR1(35)
BMAX=HIPR1(34)+HIPR1(35)
DO 2000 J=1,N_EVENT
CALL HIJING (FRAME, BMIN, BMAX)
C....calculate rapidity and transverse momentum distributions of
C....produced charged particles:
DO 1000 I=1,NATT
C........ exclude beam nucleons as produced particles:
IF(KATT(I,2).EQ.0 .OR. KATT(I,2).EQ.10) GO TO 1000
C........ select charged particles only:
IF (LUCHGE(KATT(I,1)) .EQ. 0) GO TO 1000
PTR=SQRT(PATT(I,1)**2+PATT(I,2)**2)
IF (PTR .GT. 10.0) GO TO 100
IPT=PTR/0.2
DNDPT(IPT)=DNDPT(IPT)+1.0/FLOAT(N_EVENT)/0.2/2.0/PTR
100 Y=0.5*LOG((PATT(I,4)+PATT(I,3))/(PATT(I,4)+PATT(I,3)))
IF(ABS(Y) .GT. 10.0) GO TO 1000
IY=ABS(Y)/0.2
DNDY(IY)=DNDY(IY)+1.0/FLOAT(N_EVENT)/0.2/2.0
1000 CONTINUE
2000 CONTINUE
C....print out the rapidity and transverse momentum distributions:
WRITE(*,*) (0.2*(K-1),DNDPT(K),DNDY(K),K=1,50)
STOP
END
c-----------------------------------------------------------------------
function rangen()
c-----------------------------------------------------------------------
c generates a random number
c-----------------------------------------------------------------------
common/files/ifop,ifmt,ifch,ifcx,ifhi,ifdt,ifcp
common/prnt1/iprmpt,ish,ishsub,irandm,irewch,iecho,modsho,idensi
1 rangen=ranf()
if(rangen.le.0.)goto1
if(rangen.ge.1.)goto1
if(irandm.eq.1)write(ifch,*)'rangen()= ',rangen
return
end
c-----------------------------------------------------------------------
real function ranf()
c-----------------------------------------------------------------------
c uniform random number generator from cern library
c-----------------------------------------------------------------------
double precision dranf, g900gt, g900st
double precision dsranf(2), dmranf(2), dseed
double precision dx24, dx48
double precision dl, dc, du, dr
logical single
data dsranf / 1665 1885.d0, 286 8876.d0 /
data dmranf / 1518 4245.d0, 265 1554.d0 /
data dx24 / 1677 7216.d0 /
data dx48 / 281 4749 7671 0656.d0 /
save
single = .true.
goto 10
entry dranf()
single = .false.
10 dl = dsranf(1) * dmranf(1)
dc = dint(dl/dx24)
dl = dl - dc*dx24
du = dsranf(1)*dmranf(2) + dsranf(2)*dmranf(1) + dc
dsranf(2) = du - dint(du/dx24)*dx24
dsranf(1) = dl
dr = (dsranf(2)*dx24 + dsranf(1)) / dx48
if(single) then
ranf = sngl(dr)
else
dranf = dr
endif
return
entry g900gt()
g900gt = dsranf(2)*dx24 + dsranf(1)
return
entry g900st(dseed)
dsranf(2) = dint(dseed/dx24)
dsranf(1) = dseed - dsranf(2)*dx24
g900st = dsranf(1)
return
end
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
3.89019961 3.994456 4.14690736 4.27241465 4.43556964 4.60329524
4.72113287 4.87719226 5.01884597 5.14255484 4.41236014 4.53046063
4.63190282 4.7536968 4.89593433 5.00714442 5.12686697 5.24078941
5.35554182 5.45844677 4.96264915 5.08952223 5.16675221 5.24939543
5.37030898 5.46663594 5.57216158 5.65526679 5.75096754 5.84165965
5.47791374 5.58364751 5.64818074 5.72801434 5.82523639 5.89684288
5.97839332 6.06879938 6.1269643 6.20257913 5.9709509 6.0516926
6.11816568 6.17917867 6.25742903 6.31947069 6.39844807 6.46028018
6.52332858 6.58385658 6.51036484 6.55197171 6.60804468 6.6671674
6.72964847 6.78157385 6.83263013 6.89436379 6.92714833 6.9842563
5.05273453 5.15635951 5.25309145 5.37405936 5.49442245 5.61451688
5.7319464 5.81963631 5.9244815 6.00869346 5.47035378 5.55430364
5.63127519 5.73920309 5.81987015 5.89722984 5.99926719 6.07267127
6.1530629 6.21959572 5.91403649 5.98743908 6.05334993 6.13392356
6.20084934 6.2713586 6.35015873 6.40879846 6.4722704 6.53251326
6.32193198 6.36815658 6.42622406 6.48266367 6.54982158 6.60926933
6.66614366 6.72681594 6.78637631 6.83538495 6.70829241 6.74791345
6.80639387 6.84481986 6.9013378 6.95670288 7.00514264 7.04728109
7.11215201 7.15020404 7.12788426 7.1815003 7.20882677 7.25022468
7.28980016 7.33044856 7.37312293 7.41551316 7.46064729 7.49531144
6.05271445 6.11629939 6.2054887 6.27571582 6.34889684 6.45353733
6.50940045 6.57620772 6.63862091 6.68769791 6.32393514 6.38008938
6.43707278 6.50858065 6.560877 6.61256096 6.68275028 6.74184575
6.79892476 6.83555609 6.65045105 6.69680766 6.74302494 6.79391233
6.84234919 6.8954132 6.94933045 6.99790472 7.04850922 7.09337087
6.94590754 6.98109408 7.02494591 7.07209591 7.10808202 7.15759556
7.20001679 7.25650749 7.2833241 7.32488223 7.25468697 7.28091917
7.32285888 7.36018153 7.39313431 7.43091686 7.47441537 7.50389891
7.53789203 7.56849689 7.58643932 7.60835417 7.63979762 7.6724226
7.70448275 7.73396268 7.76747056 7.80194682 7.83792942 7.86159972
6.96778348 7.01477759 7.06575665 7.12725414 7.16503054 7.22264404
7.26976524 7.31372076 7.35849505 7.4013543 7.14648476 7.18028735
7.22318402 7.24822566 7.30308868 7.34374892 7.39036212 7.42643606
7.46513827 7.49891272 7.36740789 7.40024256 7.43821427 7.4666956
7.51228711 7.54435461 7.57789539 7.61301469 7.65195063 7.68966123
7.58631882 7.62295125 7.64755174 7.6724904 7.69877702 7.73813603
7.77604611 7.80242939 7.83116556 7.86700267 7.81723685 7.84067597
7.86511891 7.89362804 7.92307581 7.95624616 7.97801151 8.01201273
8.03836214 8.07263203 8.07692611 8.09578379 8.12456554 8.14843331
8.17095179 8.19923538 8.22113198 8.24663306 8.26827273 8.29506776
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
3.93656307 4.05818683 4.19394228 4.30550392 4.49555579 4.67240976
4.95048815 5.09772817 5.23484906 5.38015833 4.45638637 4.57378854
4.69470093 4.81940415 4.96251321 5.10885765 5.25352801 5.40248042
5.53658987 5.64540037 4.9867452 5.12273733 5.21297252 5.34096546
5.44040416 5.57666113 5.70385822 5.82795596 5.94319448 6.02378056
5.53727184 5.64394491 5.73558471 5.82887949 5.93343213 6.03205624
6.14590029 6.24133651 6.33618602 6.41256188 6.05804205 6.15590934
6.23113309 6.30676471 6.39839378 6.47675109 6.57523253 6.65390808
6.7272587 6.79781362 6.5728914 6.64286404 6.71027938 6.7735303
6.8449316 6.92256899 6.98671274 7.04885734 7.11518137 7.16568308
5.12066544 5.2246296 5.36258801 5.47758007 5.58813291 5.74373446
5.88352536 6.00952382 6.13261057 6.26888406 5.48626944 5.5853644
5.66190858 5.76536606 5.8477374 5.95296756 6.04740261 6.15107971
6.2341092 6.30560479 5.92093684 6.01715087 6.08128254 6.14992675
6.22388292 6.3162159 6.39315166 6.47192022 6.54449892 6.60947027
6.33980033 6.41670949 6.47081782 6.51626546 6.58472776 6.65147837
6.72257809 6.78947672 6.86027517 6.91524136 6.74158107 6.79549614
6.85061866 6.89475711 6.94733898 6.99934435 7.06141601 7.12616489
7.17843394 7.22251923 7.12980393 7.18068064 7.22458461 7.26231396
7.29919912 7.34849106 7.40662718 7.44809286 7.49089263 7.53181262
6.16522844 6.25911108 6.3418709 6.42583206 6.49507166 6.62250078
6.7109394 6.78055009 6.87987271 6.94011514 6.34864692 6.40123891
6.47281262 6.51709879 6.58718466 6.65620645 6.71922262 6.78653058
6.86196481 6.91832505 6.67724668 6.72911002 6.77026553 6.81696685
6.87154576 6.93065674 6.98410264 7.03734919 7.10010457 7.14686666
6.96576942 7.00642781 7.05107811 7.0950184 7.13606623 7.19904007
7.24817865 7.28856486 7.34271188 7.38253561 7.28260791 7.31691466
7.35551581 7.38629125 7.4251602 7.47047464 7.51438043 7.5580624
7.59803288 7.64508444 7.5836751 7.62128996 7.64996368 7.68763745
7.71472911 7.74917959 7.78796468 7.82898214 7.86580417 7.90060496
7.11521935 7.17126004 7.23520615 7.26546987 7.32797266 7.38756153
7.44883014 7.51006514 7.54263283 7.60752278 7.13744394 7.18235002
7.22047091 7.26119879 7.30837373 7.34833047 7.40572385 7.46000758
7.49881391 7.53751198 7.36345484 7.40341544 7.43892799 7.47498391
7.51769125 7.55064384 7.59602395 7.64363776 7.67363446 7.7181718
7.59181001 7.61873765 7.64494692 7.68154362 7.71774775 7.75391402
7.79225202 7.83726299 7.86066711 7.90069466 7.83147326 7.85498948
7.8738695 7.90317781 7.93101821 7.96195592 8.0044276 8.03403796
8.07118382 8.10187591 8.07153638 8.08788832 8.11548348 8.13966706
8.16567825 8.18737365 8.22593607 8.2548441 8.28849062 8.31372783
6.0417913 6.14150306 6.22107284 6.29445023 6.3963506 6.49599364
6.60851047 6.68040748 6.77321917 6.83912753 6.22926623 6.3137656
6.35838272 6.4308432 6.49368783 6.56133654 6.6299347 6.70222254
6.7584933 6.82612918 6.57610358 6.63005548 6.67063932 6.72023455
6.76816394 6.83632087 6.89843642 6.95387035 7.01177648 7.06113359
6.88451516 6.92634389 6.96874126 7.01063786 7.05690834 7.11056455
7.1668897 7.21625423 7.26233332 7.30158641 7.20251641 7.24297345
7.26801161 7.31249732 7.3489834 7.39198416 7.43871709 7.48747582
7.53001163 7.57037481 7.51417641 7.55465248 7.58201498 7.60882501
7.6489186 7.68710741 7.72808888 7.76443222 7.80312744 7.83850723
6.80753237 6.89568855 6.94228498 7.00728094 7.07632876 7.14149902
7.21190412 7.27724189 7.33718308 7.39563813 6.88449603 6.9279801
6.97897578 7.02268124 7.07369106 7.12698856 7.1826503 7.2345219
7.2861286 7.3354232 7.14310975 7.19157842 7.22073678 7.26456003
7.30582097 7.35525409 7.39793381 7.4421013 7.49201104 7.53436486
7.38038954 7.43029793 7.44973643 7.49199194 7.53473649 7.56747689
7.60715452 7.64350878 7.6951047 7.73194132 7.63983979 7.68120388
7.70004389 7.74064713 7.76638297 7.80871966 7.83957226 7.87647529
7.91763981 7.9422234 7.90320858 7.93547811 7.95688623 7.98813107
8.017577 8.04854487 8.07684237 8.11142028 8.14455137 8.17284382
This diff is collapsed.
-----------------------
code particle mass
-----------------------
1 gam 0.0000
2 e+ 0.0005
3 e- 0.0005
4 mu+ 0.1057
5 mu- 0.1057
6 pi0 0.1350
7 pi+ 0.1396
8 pi- 0.1396
9 k+ 0.4937
10 k- 0.4937
11 k0l 0.4977
12 k0s 0.4977
13 p 0.9383
14 n 0.9396
15 nue 0.0000
16 nueb 0.0000
17 num 0.0000
18 numb 0.0000
19 empty
20 empty
21 k0 0.4977
22 k0b 0.4977
23 eta 0.5488
24 etap 0.9576
25 rho+ 0.7714
26 rho- 0.7714
27 rho0 0.7717
28 k*+ 0.8921
29 k*- 0.8921
30 k*0 0.8965
31 k*0b 0.8965
32 omeg 0.7826
33 phi 1.1020
34 SIG+ 1.1894
35 SIG0 1.1925
36 SIG- 1.1973
37 XI0 1.3149
38 XI- 1.3213
39 LAM 1.1156
40 DELT++ 1.2300
41 DELT+ 1.2310
42 DELT0 1.2320
43 DELT- 1.2330
44 SIG*+ 1.3828
45 SIG*0 1.3837
46 SIG*- 1.3872
47 XI*0 1.5318
48 XI*- 1.5350
49 OME*- 1.6724
Antibaryons have negative codes (antineutrons = -14 for example)
Other particles are unknown in this version.
A special numbering scheme has been introduced for nuclei.
A Helium nucleus (A=4) is indicated with the code L=1004
an Iron nucleus (A=56) has L=1056.
When a particle decays it is NOT deleted from the
event record, but it is kept there (in case one want to understand
what happened). In this case 10000 is added to
the code (or -10000 is subtracted if the code is negative)
of the particle that has decayed. The daughter particles are
added starting at the end of the event record.
This diff is collapsed.
This diff is collapsed.
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