Geant4 Simulations for the Radon Electric Dipole Moment Search at
Geant4 Simulations for the Radon Electric Dipole Moment Search at
Geant4 Simulations for the Radon Electric Dipole Moment Search at
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
44 return;<br />
45 end<br />
46 if abs(m) > j<br />
47 %warning('m is out of bounds.');<br />
48 cg = 0;<br />
49 return;<br />
50 end<br />
51 % //////////<br />
52 term1 = (((2*j+1)/factorial(j1+j2+j+1))*factorial(j2+j−j1)*...<br />
53 factorial(j+j1−j2)*factorial(j1+j2−j)*factorial(j1+m1)*...<br />
54 factorial(j1−m1)*factorial(j2+m2)*factorial(j2−m2)*...<br />
55 factorial(j+m)*factorial(j−m))ˆ(0.5);<br />
56 sum = 0;<br />
57 <strong>for</strong> k=0:1:99<br />
58 if ( (j1+j2−j−k < 0) || (j−j1−m2+k < 0) || (j−j2+m1+k < 0)...<br />
59 || (j1−m1−k < 0) || (j2+m2−k < 0) )<br />
60 else<br />
61 term = factorial(j1+j2−j−k)*factorial(j−j1−m2+k)*...<br />
62 factorial(j−j2+m1+k)*factorial(j1−m1−k)*...<br />
63 factorial(j2+m2−k)*factorial(k);<br />
64 if (mod(k,2) == 1)<br />
65 term = −1*term;<br />
66 end<br />
67 sum = sum + 1/term;<br />
68 end<br />
69 end<br />
70 cg =term1*sum;<br />
71 % Reference: An Effective Algorithm <strong>for</strong> Calcul<strong>at</strong>ion of <strong>the</strong> C.G.<br />
72 % Coefficients Liang Zuo, et. al.<br />
73 % J. Appl. Cryst. (1993). 26, 302−304<br />
91