example_1.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file example_1.f

C****************************************************************************
C
C
C
C The following is an example program for calling HIJING. one should
C include all the common blocks and the data values which are listed
C below in his own program.
C***************************************************************************

        CHARACTER frame*8,proj*8,targ*8
  common/himain1/ natt,eatt,jatt,nt,np,n0,n01,n10,n11
        SAVE  /himain1/

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

C         ********information of produced particles
C
  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)
        SAVE  /hijjet1/

        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)
        SAVE  /hijjet2/

C         ********information of produced partons

      common/ranseed/nseed                                     ! Uzhi
      SAVE  /ranseed/                                          ! Uzhi
      nseed=0                                                  ! Uzhi

      efrm =200.0
      frame='CMS'
      proj ='P'    ! 'A'                                       ! Uzhi
      targ ='P'    ! 'A'                                       ! Uzhi
      iap  =1      ! 197                                       ! Uzhi
      izp  =1      !  79                                       ! Uzhi
      iat  =1      ! 197                                       ! Uzhi
      izt  =1      !  79                                       ! Uzhi
      n_events=100                                             ! Uzhi

C   Simulation of PP-interactions at \sqrt{s}=200 GeV          ! Uzhi


      CALL hijset(efrm,frame,proj,targ,iap,izp,iat,izt)
C     ********Initialize HIJING

      WRITE(6,*)' Simulation of interactions with'             ! Uzhi
      WRITE(6,*)                                               ! Uzhi
      WRITE(6,*)' Proj = ',proj,' and Targ = ',targ            ! Uzhi
      WRITE(6,*)' IAP  =',iap  ,'            IAT  =',iat       ! Uzhi
      WRITE(6,*)' IZP  =',izp  ,'            IZT  =',izt       ! Uzhi
      WRITE(6,*)                                               ! Uzhi
      WRITE(6,*)' Reference frame -   ',frame                  ! Uzhi
      WRITE(6,*)' ENERGY            ',efrm,' GeV'              ! Uzhi
      WRITE(6,*)' Number of generated events -',n_events       ! Uzhi
      WRITE(6,*)                                               ! Uzhi

      bmin=0.0
      bmax=0.0
      DO 2000 i_event=1,n_events

        WRITE(6,*)' Event # ',i_event,' ------------------'    ! Uzhi

        CALL hijing(frame,bmin,bmax)
C
        WRITE(6,*)' Multiplicity of produced particles - ',natt! Uzhi
        write(6,*)                                             ! Uzhi
        WRITE(6,*)'   ID Charge  Mass (GeV)  Px           Py', ! Uzhi
     &            '           Pz    (GeV/c)'                   ! Uzhi
        WRITE(6,*)' ----------------------------------------'  ! Uzhi


        DO 1000 i=1,natt
                                                               ! Uzhi
C         IF(LUCHGE(KATT(I,1)).NE.0) THEN                      ! Uzhi
C           this select charged particles only
C                            !information of produced particles
C                !is stored in common blocks HIMAIN1 and
C         HIMAIN2
C         ENDIF                                                ! Uzhi

          ich=luchge(katt(i,1))/3                              ! Uzhi

          amass=ulmass(katt(i,1))                              ! Uzhi

          WRITE(6,900) katt(i,1),ich,amass,                    ! Uzhi
     &                 patt(i,1),patt(i,2),patt(i,3)           ! Uzhi
 900      FORMAT(1x,i6,i4,2x,f7.3,3(2x,e11.4))                 ! Uzhi

1000    CONTINUE

        pause                                                  ! Uzhi

2000  CONTINUE
  stop
  END

      FUNCTION ran(NSEED)                                      ! Uzhi
      ran=rlu(nseed)                                           ! Uzhi
      RETURN                                                   ! Uzhi
      END                                                      ! Uzhi