| | 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 | #define BLANKET_SAVE |
| | 10 | #include "hiparnt.inc" |
| | 11 | #include "njet_hijing.inc" |
| | 12 | EXTERNAL FHIN,FTOT,FNJET,FTOTJET,FTOTRIG |
| | 13 | SAVE |
| | 14 | IF(HINT1(1).GE.10.0) CALL CRSJET |
| | 15 | C ********calculate jet cross section(in mb) |
| | 16 | C |
| | 17 | APHX1=HIPR1(6)*(IHNT2(1)**0.3333333-1.0) |
| | 18 | APHX2=HIPR1(6)*(IHNT2(3)**0.3333333-1.0) |
| | 19 | HINT1(11)=HINT1(14)-APHX1*HINT1(15) |
| | 20 | & -APHX2*HINT1(16)+APHX1*APHX2*HINT1(17) |
| | 21 | HINT1(10)=GAUSS1(FTOTJET,0.0,20.0,0.01) |
| | 22 | HINT1(12)=GAUSS1(FHIN,0.0,20.0,0.01) |
| | 23 | HINT1(13)=GAUSS1(FTOT,0.0,20.0,0.01) |
| | 24 | HINT1(60)=HINT1(61)-APHX1*HINT1(62) |
| | 25 | & -APHX2*HINT1(63)+APHX1*APHX2*HINT1(64) |
| | 26 | HINT1(59)=GAUSS1(FTOTRIG,0.0,20.0,0.01) |
| | 27 | IF(HINT1(59).EQ.0.0) HINT1(59)=HINT1(60) |
| | 28 | IF(HINT1(1).GE.10.0) Then |
| | 29 | DO 20 I=0,20 |
| | 30 | N=I |
| | 31 | HINT1(80+I)=GAUSS1(FNJET,0.0,20.0,0.01)/HINT1(12) |
| | 32 | 20 CONTINUE |
| | 33 | ENDIF |
| | 34 | HINT1(10)=HINT1(10)*HIPR1(31) |
| | 35 | HINT1(12)=HINT1(12)*HIPR1(31) |
| | 36 | HINT1(13)=HINT1(13)*HIPR1(31) |
| | 37 | HINT1(59)=HINT1(59)*HIPR1(31) |
| | 38 | C ********Total and Inel cross section are calculated |
| | 39 | C by Gaussian integration. |
| | 40 | IF(IHPR2(13).NE.0) THEN |
| | 41 | HIPR1(33)=1.36*(1.0+36.0/HINT1(1)**2) |
| | 42 | & *ALOG(0.6+0.1*HINT1(1)**2) |
| | 43 | HIPR1(33)=HIPR1(33)/HINT1(12) |
| | 44 | ENDIF |
| | 45 | C ********Parametrized cross section for single |
| | 46 | C diffractive reaction(Goulianos) |
| | 47 | RETURN |
| | 48 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental