]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HIJING/hijing1_36/crsjet.F
Coding conventions
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / crsjet.F
CommitLineData
e74335a4 1* $Id$
2C
3C
4C THIS PROGRAM IS TO CALCULATE THE JET CROSS SECTION
5C THE INTEGRATION IS DONE BY USING VEGAS
6C
7 SUBROUTINE CRSJET
8 IMPLICIT REAL*8(A-H,O-Z)
bc676b8e 9#define BLANKET_SAVE
e74335a4 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
17C
18c************************
19c NCALL give the number of inner-iteration, ITMX
20C gives the limit of out-iteration. Nprn is an option
21C ( 1: print the integration process. 0: do not print)
22C
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
42C ********Total inclusive jet cross section(Pt>P0)
43C
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
64C ********cross section of trigger jet
65C
66 RETURN
67 END