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