HIJING/hijing1_36/hijcrs.F

   1 * $Id$
   2 C
   3 C
   4 C
   5 C
   6         SUBROUTINE HIJCRS
   7 C       THIS IS TO CALCULATE THE CROSS SECTIONS OF JET PRODUCTION AND
   8 C       THE TOTAL INELASTIC CROSS SECTIONS.
   9 #define BLANKET_SAVE
  10 #include "hiparnt.inc"
  11 #include "njet_hijing.inc"
  12         EXTERNAL FHIN,FTOT,FNJET,FTOTJET,FTOTRIG
  13         SAVE
  14         IF(HINT1(1).GE.10.0) CALL CRSJET
  15 C                       ********calculate jet cross section(in mb)
  16 C
  17         APHX1=HIPR1(6)*(IHNT2(1)**0.3333333-1.0)
  18         APHX2=HIPR1(6)*(IHNT2(3)**0.3333333-1.0)
  19         HINT1(11)=HINT1(14)-APHX1*HINT1(15)
  20      &                  -APHX2*HINT1(16)+APHX1*APHX2*HINT1(17)
  21         HINT1(10)=GAUSS1(FTOTJET,0.0,20.0,0.01)
  22         HINT1(12)=GAUSS1(FHIN,0.0,20.0,0.01)
  23         HINT1(13)=GAUSS1(FTOT,0.0,20.0,0.01)
  24         HINT1(60)=HINT1(61)-APHX1*HINT1(62)
  25      &                  -APHX2*HINT1(63)+APHX1*APHX2*HINT1(64)
  26         HINT1(59)=GAUSS1(FTOTRIG,0.0,20.0,0.01)
  27         IF(HINT1(59).EQ.0.0) HINT1(59)=HINT1(60)
  28         IF(HINT1(1).GE.10.0) Then
  29            DO 20 I=0,20
  30               N=I
  31               HINT1(80+I)=GAUSS1(FNJET,0.0,20.0,0.01)/HINT1(12)
  32  20        CONTINUE
  33         ENDIF
  34         HINT1(10)=HINT1(10)*HIPR1(31)
  35         HINT1(12)=HINT1(12)*HIPR1(31)
  36         HINT1(13)=HINT1(13)*HIPR1(31)
  37         HINT1(59)=HINT1(59)*HIPR1(31)
  38 C               ********Total and Inel cross section are calculated
  39 C                       by Gaussian integration.
  40         IF(IHPR2(13).NE.0) THEN
  41         HIPR1(33)=1.36*(1.0+36.0/HINT1(1)**2)
  42      &             *ALOG(0.6+0.1*HINT1(1)**2)
  43         HIPR1(33)=HIPR1(33)/HINT1(12)
  44         ENDIF
  45 C               ********Parametrized cross section for single
  46 C                       diffractive reaction(Goulianos)
  47         RETURN
  48         END