HIJING/hijing1_36/crsjet.F

   1 * $Id$
   2 C
   3 C
   4 C       THIS PROGRAM IS TO CALCULATE THE JET CROSS SECTION
   5 C       THE INTEGRATION IS DONE BY USING VEGAS
   6 C
   7         SUBROUTINE CRSJET
   8         IMPLICIT REAL*8(A-H,O-Z)
   9 #include "hiparnt.inc"
  10 #include "njet_hijing.inc"
  11 #include "bveg1_hijing.inc"
  12 #include "bveg2_hijing.inc"
  13 #include "bveg3_hijing.inc"
  14         EXTERNAL FJET,FJETRIG
  15         SAVE
  16 C
  17 c************************
  18 c       NCALL give the number of inner-iteration, ITMX
  19 C       gives the limit of out-iteration. Nprn is an option
  20 C       ( 1: print the integration process. 0: do not print)
  21 C
  22         NDIM=3
  23         IP_CRS=0
  24         CALL VEGAS(FJET,AVGI,SD,CHI2A)
  25         HINT1(14)=AVGI/2.5682
  26         IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
  27                 IP_CRS=1
  28                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  29                 HINT1(15)=AVGI/2.5682
  30         ENDIF
  31         IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
  32                 IP_CRS=2
  33                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  34                 HINT1(16)=AVGI/2.5682
  35         ENDIF
  36         IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
  37                 IP_CRS=3
  38                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  39                 HINT1(17)=AVGI/2.5682
  40         ENDIF
  41 C               ********Total inclusive jet cross section(Pt>P0)
  42 C
  43         IF(IHPR2(3).NE.0) THEN
  44            IP_CRS=0
  45            CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  46            HINT1(61)=AVGI/2.5682
  47            IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
  48               IP_CRS=1
  49               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  50               HINT1(62)=AVGI/2.5682
  51            ENDIF
  52            IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
  53               IP_CRS=2
  54               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  55               HINT1(63)=AVGI/2.5682
  56            ENDIF
  57            IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
  58               IP_CRS=3
  59               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  60               HINT1(64)=AVGI/2.5682
  61            ENDIF
  62         ENDIF
  63 C                       ********cross section of trigger jet
  64 C
  65         RETURN
  66         END