ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
hijwds.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file hijwds.f
1 C*************************
2 C
3 C
4 C
5 C
6 C ********************************************************
7 C ************************ WOOD-SAX
8  SUBROUTINE hijwds(IA,IDH,XHIGH)
9 C SETS UP HISTOGRAM IDH WITH RADII FOR
10 C NUCLEUS IA DISTRIBUTED ACCORDING TO THREE PARAM WOOD SAXON
11  common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
12  SAVE /hiparnt/
13  common/wood/r,d,fnorm,w
14  SAVE /wood/
15  dimension iaa(20),rr(20),dd(20),ww(20),rms(20)
16  EXTERNAL rwdsax,wdsax
17 C
18 C PARAMETERS OF SPECIAL NUCLEI FROM ATOMIC DATA AND NUC DATA TABLES
19 C VOL 14, 5-6 1974
20  DATA iaa/2,4,12,16,27,32,40,56,63,93,184,197,208,7*0./
21  DATA rr/0.01,.964,2.355,2.608,2.84,3.458,3.766,3.971,4.214,
22  1 4.87,6.51,6.38,6.624,7*0./
23  DATA dd/0.5882,.322,.522,.513,.569,.61,.586,.5935,.586,.573,
24  1 .535,.535,.549,7*0./
25  DATA ww/0.0,.517,-0.149,-0.051,0.,-0.208,-0.161,13*0./
26  DATA rms/2.11,1.71,2.46,2.73,3.05,3.247,3.482,3.737,3.925,4.31,
27  1 5.42,5.33,5.521,7*0./
28 C
29  a=ia
30 C
31 C ********SET WOOD-SAX PARAMS FIRST AS IN DATE ET AL
32  d=0.54
33 C ********D IS WOOD SAX DIFFUSE PARAM IN FM
34  r=1.19*a**(1./3.) - 1.61*a**(-1./3.)
35 C ********R IS RADIUS PARAM
36  w=0.
37 C ********W IS The third of three WOOD-SAX PARAM
38 C
39 C ********CHECK TABLE FOR SPECIAL CASES
40  DO 10 i=1,13
41  IF (ia.EQ.iaa(i)) THEN
42  r=rr(i)
43  d=dd(i)
44  w=ww(i)
45  rs=rms(i)
46  END IF
47 10 CONTINUE
48 C ********FNORM is the normalize factor
49  fnorm=1.0
50  xlow=0.
51  xhigh=r+ 12.*d
52  IF (w.LT.-0.01) THEN
53  IF (xhigh.GT.r/sqrt(abs(w))) xhigh=r/sqrt(abs(w))
54  END IF
55  fgaus=gauss1(rwdsax,xlow,xhigh,0.001)
56  fnorm=1./fgaus
57 C
58  IF (idh.EQ.1) THEN
59  hint1(72)=r
60  hint1(73)=d
61  hint1(74)=w
62  hint1(75)=fnorm/4.0/hipr1(40)
63  ELSE IF (idh.EQ.2) THEN
64  hint1(76)=r
65  hint1(77)=d
66  hint1(78)=w
67  hint1(79)=fnorm/4.0/hipr1(40)
68  ENDIF
69 C
70 C NOW SET UP HBOOK FUNCTIONS IDH FOR R**2*RHO(R)
71 C THESE HISTOGRAMS ARE USED TO GENERATE RANDOM RADII
72  CALL hifun(idh,xlow,xhigh,rwdsax)
73  RETURN
74  END