Adding macros to create Calibration objects
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hijwds.F
CommitLineData
e74335a4 1* $Id$
2C*************************
3C
4C
5C
6C
7C ********************************************************
8C ************************ WOOD-SAX
9 SUBROUTINE HIJWDS(IA,IDH,XHIGH)
10C SETS UP HISTOGRAM IDH WITH RADII FOR
11C NUCLEUS IA DISTRIBUTED ACCORDING TO THREE PARAM WOOD SAXON
bc676b8e 12#define BLANKET_SAVE
e74335a4 13#include "hiparnt.inc"
14#include "wood_hijing.inc"
15 DIMENSION IAA(20),RR(20),DD(20),WW(20),RMS(20)
16 EXTERNAL RWDSAX,WDSAX
3cc37230 17
e74335a4 18C
19C PARAMETERS OF SPECIAL NUCLEI FROM ATOMIC DATA AND NUC DATA TABLES
20C VOL 14, 5-6 1974
3cc37230 21 DATA IAA/2,4,12,16,27,32,40,56,63,93,184,197,208,7*0./
22 DATA RR/0.01,.964,2.355,2.608,2.84,3.458,3.766,3.971,4.214,
23 1 4.87,6.51,6.38,6.624,7*0./
24 DATA DD/0.5882,.322,.522,.513,.569,.61,.586,.5935,.586,.573,
25 1 .535,.535,.549,7*0./
26 DATA WW/0.0,.517,-0.149,-0.051,0.,-0.208,-0.161,13*0./
27 DATA RMS/2.11,1.71,2.46,2.73,3.05,3.247,3.482,3.737,3.925,4.31,
28 1 5.42,5.33,5.521,7*0./
e74335a4 29 SAVE
30C
31 A=IA
32C
33C ********SET WOOD-SAX PARAMS FIRST AS IN DATE ET AL
34 D=0.54
35C ********D IS WOOD SAX DIFFUSE PARAM IN FM
36 R=1.19*A**(1./3.) - 1.61*A**(-1./3.)
37C ********R IS RADIUS PARAM
38 W=0.
39C ********W IS The third of three WOOD-SAX PARAM
40C
41C ********CHECK TABLE FOR SPECIAL CASES
3cc37230 42 DO 10 I=1,13
e74335a4 43 IF (IA.EQ.IAA(I)) THEN
44 R=RR(I)
45 D=DD(I)
46 W=WW(I)
47 RS=RMS(I)
48 END IF
4910 CONTINUE
50C ********FNORM is the normalize factor
51 FNORM=1.0
52 XLOW=0.
53 XHIGH=R+ 12.*D
54 IF (W.LT.-0.01) THEN
55 IF (XHIGH.GT.R/SQRT(ABS(W))) XHIGH=R/SQRT(ABS(W))
56 END IF
57 FGAUS=GAUSS1(RWDSAX,XLOW,XHIGH,0.001)
58 FNORM=1./FGAUS
59C
60 IF (IDH.EQ.1) THEN
61 HINT1(72)=R
62 HINT1(73)=D
63 HINT1(74)=W
64 HINT1(75)=FNORM/4.0/HIPR1(40)
65 ELSE IF (IDH.EQ.2) THEN
66 HINT1(76)=R
67 HINT1(77)=D
68 HINT1(78)=W
69 HINT1(79)=FNORM/4.0/HIPR1(40)
70 ENDIF
71C
72C NOW SET UP HBOOK FUNCTIONS IDH FOR R**2*RHO(R)
73C THESE HISTOGRAMS ARE USED TO GENERATE RANDOM RADII
74 CALL HIFUN(IDH,XLOW,XHIGH,RWDSAX)
75 RETURN
76 END