Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
example_1.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file example_1.f
1 C****************************************************************************
2 C
3 C
4 C
5 C The following is an example program for calling HIJING. one should
6 C include all the common blocks and the data values which are listed
7 C below in his own program.
8 C***************************************************************************
9 
10  CHARACTER frame*8,proj*8,targ*8
11  common/himain1/ natt,eatt,jatt,nt,np,n0,n01,n10,n11
12  SAVE /himain1/
13 
14  common/himain2/katt(130000,4),patt(130000,4)
15  SAVE /himain2/
16 
17 C ********information of produced particles
18 C
19  common/hijjet1/npj(300),kfpj(300,500),pjpx(300,500),pjpy(300,500)
20  & ,pjpz(300,500),pjpe(300,500),pjpm(300,500)
21  & ,ntj(300),kftj(300,500),pjtx(300,500),pjty(300,500)
22  & ,pjtz(300,500),pjte(300,500),pjtm(300,500)
23  SAVE /hijjet1/
24 
25  common/hijjet2/nsg,njsg(900),iasg(900,3),k1sg(900,100)
26  & ,k2sg(900,100),pxsg(900,100),pysg(900,100),pzsg(900,100)
27  & ,pesg(900,100),pmsg(900,100)
28  SAVE /hijjet2/
29 
30 C ********information of produced partons
31 
32  common/ranseed/nseed ! Uzhi
33  SAVE /ranseed/ ! Uzhi
34  nseed=0 ! Uzhi
35 
36  efrm =200.0
37  frame='CMS'
38  proj ='P' ! 'A' ! Uzhi
39  targ ='P' ! 'A' ! Uzhi
40  iap =1 ! 197 ! Uzhi
41  izp =1 ! 79 ! Uzhi
42  iat =1 ! 197 ! Uzhi
43  izt =1 ! 79 ! Uzhi
44  n_events=100 ! Uzhi
45 
46 C Simulation of PP-interactions at \sqrt{s}=200 GeV ! Uzhi
47 
48 
49  CALL hijset(efrm,frame,proj,targ,iap,izp,iat,izt)
50 C ********Initialize HIJING
51 
52  WRITE(6,*)' Simulation of interactions with' ! Uzhi
53  WRITE(6,*) ! Uzhi
54  WRITE(6,*)' Proj = ',proj,' and Targ = ',targ ! Uzhi
55  WRITE(6,*)' IAP =',iap ,' IAT =',iat ! Uzhi
56  WRITE(6,*)' IZP =',izp ,' IZT =',izt ! Uzhi
57  WRITE(6,*) ! Uzhi
58  WRITE(6,*)' Reference frame - ',frame ! Uzhi
59  WRITE(6,*)' ENERGY ',efrm,' GeV' ! Uzhi
60  WRITE(6,*)' Number of generated events -',n_events ! Uzhi
61  WRITE(6,*) ! Uzhi
62 
63  bmin=0.0
64  bmax=0.0
65  DO 2000 i_event=1,n_events
66 
67  WRITE(6,*)' Event # ',i_event,' ------------------' ! Uzhi
68 
69  CALL hijing(frame,bmin,bmax)
70 C
71  WRITE(6,*)' Multiplicity of produced particles - ',natt! Uzhi
72  write(6,*) ! Uzhi
73  WRITE(6,*)' ID Charge Mass (GeV) Px Py', ! Uzhi
74  & ' Pz (GeV/c)' ! Uzhi
75  WRITE(6,*)' ----------------------------------------' ! Uzhi
76 
77 
78  DO 1000 i=1,natt
79  ! Uzhi
80 C IF(LUCHGE(KATT(I,1)).NE.0) THEN ! Uzhi
81 C this select charged particles only
82 C !information of produced particles
83 C !is stored in common blocks HIMAIN1 and
84 C HIMAIN2
85 C ENDIF ! Uzhi
86 
87  ich=luchge(katt(i,1))/3 ! Uzhi
88 
89  amass=ulmass(katt(i,1)) ! Uzhi
90 
91  WRITE(6,900) katt(i,1),ich,amass, ! Uzhi
92  & patt(i,1),patt(i,2),patt(i,3) ! Uzhi
93  900 FORMAT(1x,i6,i4,2x,f7.3,3(2x,e11.4)) ! Uzhi
94 
95 1000 CONTINUE
96 
97  pause ! Uzhi
98 
99 2000 CONTINUE
100  stop
101  END
102 
103  FUNCTION ran(NSEED) ! Uzhi
104  ran=rlu(nseed) ! Uzhi
105  RETURN ! Uzhi
106  END ! Uzhi