d3/d14/Main__HIJING_8f_source.html

      program test

      CHARACTER frame*8,proj*8,targ*8

      common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)

      SAVE  /hiparnt/

      common/himain1/natt,eatt,jatt,nt,np,n0,n01,n10,n11

      SAVE  /himain1/

      common/himain2/katt(130000,4),patt(130000,4)

      SAVE  /himain2/

      common/ranseed/nseed

      SAVE  /ranseed/

C

C

C....initialize HIJING for Au+Au collisions at c.m. energy of 200 GeV:

c        FRAME="CMS"

c        PROJ="A"

c TARG="A"

        iap=1

  izp=1

  iat=197

  izt=79

        efrm=200.0   !GeV

        n_event=100

        ihpr2(10)=0

c       ---------------------------------------------------------

c        initialize HIJING for geant4

  ihnt2(1)=iap

  ihnt2(2)=izp

  ihnt2(3)=iat

  ihnt2(4)=izt

  ihnt2(5)=211   !pi+

  ihnt2(6)=0

cc     --------------Ap>1----------------------


        hint1(8)=max(ulmass(2112),ulmass(2212))

        hint1(9)=hint1(8)

        write(*,*)'mass at the beginning',hint1(8), hint1(9)

cc---------------------------------------------

        if(ihnt2(5).ne.0) then

        hint1(8)=ulmass(ihnt2(5))

        endif


c        write(*,*)'modified mass',HINT1(8)

c

c

c         ranseed=1097569630

c         print *,"seed",ranseed

c         call sseed(ranseed)


c         do i=1, 10

c         write(*,*)'ran',NSEED, RAN(NSEED), RLU(NSEED)

c         end do

C

C

C*** initialize HIJING

c        CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)

         CALL hijset(efrm,iap,izp,iat,izt)

c         CALL HIJSET(EFRM)


          write(*,*)'Sjet=', hint1(11),'mb','Stot=',hint1(13),'mb'

C

C....set BMIN=0 and BMAX=0.0 for central interaction

      bmin=0.0

      bmax=0.0    !HIPR1(34)+HIPR1(35)

C....generating N_EVENT events of central AA interaction:

      DO 200 ie=1,n_event

c         CALL HIJING(FRAME,BMIN,BMAX)

         CALL hijing(bmin,bmax)

         WRITE(*,*) ie,natt,eatt

 200  continue

      stop

      END


c FUNCTION RAN(NSEED)

c RAN=RLU(NSEED)

c RETURN

c END