hijcrs.f

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

C
C
C
C
        SUBROUTINE hijcrs
C       THIS IS TO CALCULATE THE CROSS SECTIONS OF JET PRODUCTION AND
C       THE TOTAL INELASTIC CROSS SECTIONS.
        common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
        SAVE  /hiparnt/
        common/njet/n,ip_crs
        SAVE  /njet/
        EXTERNAL fhin,ftot,fnjet,ftotjet,ftotrig
        IF(hint1(1).GE.10.0) CALL crsjet
C                       ********calculate jet cross section(in mb)
C
        aphx1=hipr1(6)*(ihnt2(1)**0.3333333-1.0)
        aphx2=hipr1(6)*(ihnt2(3)**0.3333333-1.0)
        hint1(11)=hint1(14)-aphx1*hint1(15)
     &                  -aphx2*hint1(16)+aphx1*aphx2*hint1(17)
        hint1(10)=gauss1(ftotjet,0.0,20.0,0.01)
        hint1(12)=gauss1(fhin,0.0,20.0,0.01)
        hint1(13)=gauss1(ftot,0.0,20.0,0.01)
        hint1(60)=hint1(61)-aphx1*hint1(62)
     &                  -aphx2*hint1(63)+aphx1*aphx2*hint1(64)
        hint1(59)=gauss1(ftotrig,0.0,20.0,0.01)
        IF(hint1(59).EQ.0.0) hint1(59)=hint1(60)
        IF(hint1(1).GE.10.0) Then
           DO 20 i=0,20
              n=i
              hint1(80+i)=gauss1(fnjet,0.0,20.0,0.01)/hint1(12)
 20        CONTINUE
        ENDIF
        hint1(10)=hint1(10)*hipr1(31)
        hint1(12)=hint1(12)*hipr1(31)
        hint1(13)=hint1(13)*hipr1(31)
        hint1(59)=hint1(59)*hipr1(31)
C               ********Total and Inel cross section are calculated
C                       by Gaussian integration.
        IF(ihpr2(13).NE.0) THEN
        hipr1(33)=1.36*(1.0+36.0/hint1(1)**2)
     &             *alog(0.6+0.1*hint1(1)**2)
        hipr1(33)=hipr1(33)/hint1(12)
        ENDIF
C               ********Parametrized cross section for single
C                       diffractive reaction(Goulianos)
        RETURN
        END