4 C THIS PROGRAM IS TO CALCULATE THE JET CROSS SECTION
5 C THE INTEGRATION IS DONE BY USING VEGAS
8 IMPLICIT REAL*8(A-H,O-Z)
10 #include "hiparnt.inc"
11 #include "njet_hijing.inc"
12 #include "bveg1_hijing.inc"
13 #include "bveg2_hijing.inc"
14 #include "bveg3_hijing.inc"
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)
25 CALL VEGAS(FJET,AVGI,SD,CHI2A)
27 IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
29 CALL VEGAS(FJET,AVGI,SD,CHI2A)
32 IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
34 CALL VEGAS(FJET,AVGI,SD,CHI2A)
37 IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
39 CALL VEGAS(FJET,AVGI,SD,CHI2A)
42 C ********Total inclusive jet cross section(Pt>P0)
44 IF(IHPR2(3).NE.0) THEN
46 CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
48 IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
50 CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
53 IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
55 CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
58 IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
60 CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
64 C ********cross section of trigger jet