[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / crsjet.F

* $Id$
C
C
C	THIS PROGRAM IS TO CALCULATE THE JET CROSS SECTION
C	THE INTEGRATION IS DONE BY USING VEGAS
C
	SUBROUTINE CRSJET
	IMPLICIT REAL*8(A-H,O-Z)
#define BLANKET_SAVE
#include "hiparnt.inc"
#include "njet_hijing.inc"
#include "bveg1_hijing.inc"
#include "bveg2_hijing.inc"
#include "bveg3_hijing.inc"
	EXTERNAL FJET,FJETRIG
	SAVE
C
c************************
c	NCALL give the number of inner-iteration, ITMX 
C       gives the limit of out-iteration. Nprn is an option
C       ( 1: print the integration process. 0: do not print)
C
	NDIM=3
	IP_CRS=0
	CALL VEGAS(FJET,AVGI,SD,CHI2A)
	HINT1(14)=AVGI/2.5682
	IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
		IP_CRS=1
		CALL VEGAS(FJET,AVGI,SD,CHI2A)
		HINT1(15)=AVGI/2.5682
	ENDIF
	IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
		IP_CRS=2
		CALL VEGAS(FJET,AVGI,SD,CHI2A)
		HINT1(16)=AVGI/2.5682
	ENDIF
	IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
		IP_CRS=3
		CALL VEGAS(FJET,AVGI,SD,CHI2A)
		HINT1(17)=AVGI/2.5682
	ENDIF
C		********Total inclusive jet cross section(Pt>P0) 
C
	IF(IHPR2(3).NE.0) THEN
	   IP_CRS=0
	   CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
	   HINT1(61)=AVGI/2.5682
	   IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
	      IP_CRS=1
	      CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
	      HINT1(62)=AVGI/2.5682
	   ENDIF
	   IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
	      IP_CRS=2
	      CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
	      HINT1(63)=AVGI/2.5682
	   ENDIF
	   IF(IHPR2(6).EQ.1.AND.IHNT2(1).GT.1.AND.IHNT2(3).GT.1) THEN
	      IP_CRS=3
	      CALL VEGAS(FJETRIG,AVGI,SD,CHI2A)
	      HINT1(64)=AVGI/2.5682
	   ENDIF
	ENDIF
C			********cross section of trigger jet
C
	RETURN
	END
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)
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
	17	C
	18	c************************
	19	c NCALL give the number of inner-iteration, ITMX
	20	C gives the limit of out-iteration. Nprn is an option
	21	C ( 1: print the integration process. 0: do not print)
	22	C
	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
	42	C ********Total inclusive jet cross section(Pt>P0)
	43	C
	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
	64	C ********cross section of trigger jet
	65	C
	66	RETURN
	67	END