]>
Commit | Line | Data |
---|---|---|
e74335a4 | 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 |