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 #define BLANKET_SAVE
  10 #include "hiparnt.inc"
  11 #include "njet_hijing.inc"
  12 #include "bveg1_hijing.inc"
  13 #include "bveg2_hijing.inc"
  14 #include "bveg3_hijing.inc"
  15         EXTERNAL FJET,FJETRIG
  16         SAVE
  17 C
  18 c************************
  19 c       NCALL give the number of inner-iteration, ITMX
  20 C       gives the limit of out-iteration. Nprn is an option
  21 C       ( 1: print the integration process. 0: do not print)
  22 C
  23         NDIM=3
  24         IP_CRS=0
  25         CALL VEGAS(FJET,AVGI,SD,CHI2A)
  26         HINT1(14)=AVGI/2.5682
  27         IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
  28                 IP_CRS=1
  29                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  30                 HINT1(15)=AVGI/2.5682
  31         ENDIF
  32         IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
  33                 IP_CRS=2
  34                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  35                 HINT1(16)=AVGI/2.5682
  36         ENDIF
  37         IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
  38                 IP_CRS=3
  39                 CALL VEGAS(FJET,AVGI,SD,CHI2A)
  40                 HINT1(17)=AVGI/2.5682
  41         ENDIF
  42 C               ********Total inclusive jet cross section(Pt>P0)
  43 C
  44         IF(IHPR2(3).NE.0) THEN
  45            IP_CRS=0
  46            CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  47            HINT1(61)=AVGI/2.5682
  48            IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
  49               IP_CRS=1
  50               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  51               HINT1(62)=AVGI/2.5682
  52            ENDIF
  53            IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
  54               IP_CRS=2
  55               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  56               HINT1(63)=AVGI/2.5682
  57            ENDIF
  58            IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
  59               IP_CRS=3
  60               CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
  61               HINT1(64)=AVGI/2.5682
  62            ENDIF
  63         ENDIF
  64 C                       ********cross section of trigger jet
  65 C
  66         RETURN
  67         END