Adding macros to create Calibration objects
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hijcrs.F
CommitLineData
e74335a4 1* $Id$
2C
3C
4C
5C
6 SUBROUTINE HIJCRS
7C THIS IS TO CALCULATE THE CROSS SECTIONS OF JET PRODUCTION AND
8C THE TOTAL INELASTIC CROSS SECTIONS.
bc676b8e 9#define BLANKET_SAVE
e74335a4 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
15C ********calculate jet cross section(in mb)
16C
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)
38C ********Total and Inel cross section are calculated
39C 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
45C ********Parametrized cross section for single
46C diffractive reaction(Goulianos)
47 RETURN
48 END