]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | C | |
6 | SUBROUTINE HIJCRS | |
7 | C THIS IS TO CALCULATE THE CROSS SECTIONS OF JET PRODUCTION AND | |
8 | C THE TOTAL INELASTIC CROSS SECTIONS. | |
9 | #include "hiparnt.inc" | |
10 | #include "njet_hijing.inc" | |
11 | EXTERNAL FHIN,FTOT,FNJET,FTOTJET,FTOTRIG | |
12 | SAVE | |
13 | IF(HINT1(1).GE.10.0) CALL CRSJET | |
14 | C ********calculate jet cross section(in mb) | |
15 | C | |
16 | APHX1=HIPR1(6)*(IHNT2(1)**0.3333333-1.0) | |
17 | APHX2=HIPR1(6)*(IHNT2(3)**0.3333333-1.0) | |
18 | HINT1(11)=HINT1(14)-APHX1*HINT1(15) | |
19 | & -APHX2*HINT1(16)+APHX1*APHX2*HINT1(17) | |
20 | HINT1(10)=GAUSS1(FTOTJET,0.0,20.0,0.01) | |
21 | HINT1(12)=GAUSS1(FHIN,0.0,20.0,0.01) | |
22 | HINT1(13)=GAUSS1(FTOT,0.0,20.0,0.01) | |
23 | HINT1(60)=HINT1(61)-APHX1*HINT1(62) | |
24 | & -APHX2*HINT1(63)+APHX1*APHX2*HINT1(64) | |
25 | HINT1(59)=GAUSS1(FTOTRIG,0.0,20.0,0.01) | |
26 | IF(HINT1(59).EQ.0.0) HINT1(59)=HINT1(60) | |
27 | IF(HINT1(1).GE.10.0) Then | |
28 | DO 20 I=0,20 | |
29 | N=I | |
30 | HINT1(80+I)=GAUSS1(FNJET,0.0,20.0,0.01)/HINT1(12) | |
31 | 20 CONTINUE | |
32 | ENDIF | |
33 | HINT1(10)=HINT1(10)*HIPR1(31) | |
34 | HINT1(12)=HINT1(12)*HIPR1(31) | |
35 | HINT1(13)=HINT1(13)*HIPR1(31) | |
36 | HINT1(59)=HINT1(59)*HIPR1(31) | |
37 | C ********Total and Inel cross section are calculated | |
38 | C by Gaussian integration. | |
39 | IF(IHPR2(13).NE.0) THEN | |
40 | HIPR1(33)=1.36*(1.0+36.0/HINT1(1)**2) | |
41 | & *ALOG(0.6+0.1*HINT1(1)**2) | |
42 | HIPR1(33)=HIPR1(33)/HINT1(12) | |
43 | ENDIF | |
44 | C ********Parametrized cross section for single | |
45 | C diffractive reaction(Goulianos) | |
46 | RETURN | |
47 | END |