hijing/examples/example_1.f

0001 C****************************************************************************
0002 C
0003 C
0004 C
0005 C       The following is an example program for calling HIJING. one should
0006 C       include all the common blocks and the data values which are listed
0007 C       below in his own program.
0008 C***************************************************************************
0009
0010         CHARACTER FRAME*8,PROJ*8,TARG*8
0011         COMMON/HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11
0012         SAVE  /HIMAIN1/
0013
0014         COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4)
0015         SAVE  /HIMAIN2/
0016
0017 C         ********information of produced particles
0018 C
0019         COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),PJPY(300,500)
0020      &               ,PJPZ(300,500),PJPE(300,500),PJPM(300,500)
0021      &               ,NTJ(300),KFTJ(300,500),PJTX(300,500),PJTY(300,500)
0022      &               ,PJTZ(300,500),PJTE(300,500),PJTM(300,500)
0023         SAVE  /HIJJET1/
0024
0025         COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100)
0026      &         ,K2SG(900,100),PXSG(900,100),PYSG(900,100),PZSG(900,100)
0027      &         ,PESG(900,100),PMSG(900,100)
0028         SAVE  /HIJJET2/
0029
0030 C         ********information of produced partons
0031
0032       COMMON/RANSEED/NSEED                                     ! Uzhi
0033       SAVE  /RANSEED/                                          ! Uzhi
0034       NSEED=0                                                  ! Uzhi
0035
0036       EFRM =200.0
0037       FRAME='CMS'
0038       PROJ ='P'    ! 'A'                                       ! Uzhi
0039       TARG ='P'    ! 'A'                                       ! Uzhi
0040       IAP  =1      ! 197                                       ! Uzhi
0041       IZP  =1      !  79                                       ! Uzhi
0042       IAT  =1      ! 197                                       ! Uzhi
0043       IZT  =1      !  79                                       ! Uzhi
0044       N_events=100                                             ! Uzhi
0045
0046 C   Simulation of PP-interactions at \sqrt{s}=200 GeV          ! Uzhi
0047
0048
0049       CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
0050 C                       ********Initialize HIJING
0051
0052       WRITE(6,*)' Simulation of interactions with'             ! Uzhi
0053       WRITE(6,*)                                               ! Uzhi
0054       WRITE(6,*)' Proj = ',PROJ,' and Targ = ',TARG            ! Uzhi
0055       WRITE(6,*)' IAP  =',IAP  ,'            IAT  =',IAT       ! Uzhi
0056       WRITE(6,*)' IZP  =',IZP  ,'            IZT  =',IZT       ! Uzhi
0057       WRITE(6,*)                                               ! Uzhi
0058       WRITE(6,*)' Reference frame -   ',FRAME                  ! Uzhi
0059       WRITE(6,*)' ENERGY            ',EFRM,' GeV'              ! Uzhi
0060       WRITE(6,*)' Number of generated events -',N_events       ! Uzhi
0061       WRITE(6,*)                                               ! Uzhi
0062
0063       BMIN=0.0
0064       BMAX=0.0
0065       DO 2000 I_event=1,N_events
0066
0067         WRITE(6,*)' Event # ',I_event,' ------------------'    ! Uzhi
0068
0069         CALL HIJING(FRAME,BMIN,BMAX)
0070 C
0071         WRITE(6,*)' Multiplicity of produced particles - ',NATT! Uzhi
0072         write(6,*)                                             ! Uzhi
0073         WRITE(6,*)'   ID Charge  Mass (GeV)  Px           Py', ! Uzhi
0074      &            '           Pz    (GeV/c)'                   ! Uzhi
0075         WRITE(6,*)' ----------------------------------------'  ! Uzhi
0076
0077
0078         DO 1000 I=1,NATT
0079                                                                ! Uzhi
0080 C         IF(LUCHGE(KATT(I,1)).NE.0) THEN                      ! Uzhi
0081 C           this select charged particles only
0082 C                                  !information of produced particles
0083 C                                  !is stored in common blocks HIMAIN1 and
0084 C                                       HIMAIN2
0085 C         ENDIF                                                ! Uzhi
0086
0087           ICH=LUCHGE(KATT(I,1))/3                              ! Uzhi
0088
0089           Amass=ULMASS(KATT(I,1))                              ! Uzhi
0090
0091           WRITE(6,900) KATT(I,1),ICH,Amass,                    ! Uzhi
0092      &                 PATT(I,1),PATT(I,2),PATT(I,3)           ! Uzhi
0093  900      FORMAT(1X,I6,I4,2x,F7.3,3(2X,E11.4))                 ! Uzhi
0094
0095 1000    CONTINUE
0096
0097         PAUSE                                                  ! Uzhi
0098
0099 2000    CONTINUE
0100         STOP
0101         END
0102
0103       FUNCTION RAN(NSEED)                                      ! Uzhi
0104       RAN=RLU(NSEED)                                           ! Uzhi
0105       RETURN                                                   ! Uzhi
0106       END                                                      ! Uzhi