11.01.2015 Views

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

41 out = ((−1)ˆ(j1−j2−m3))/((2*j3+1)ˆ(1/2))*...<br />

42 ClebschGordan(j1,m1,j2,m2,j3,−m3);<br />

1 % Wigner6j.m<br />

2 function sum = Wigner6j(J1,J2,J3,J4,J5,J6)<br />

3 % error checking //////////<br />

4 if J3 > J1 + J2 || J3 < abs(J1 − J2)<br />

5 warning('first J3 triange condition not s<strong>at</strong>isfied. J3 > J1 + ...<br />

J2<br />

|| J3 < abs(J1 − J2)');<br />

6 sum = 0;<br />

7 return;<br />

8 end<br />

9 if J3 > J4 + J5 || J3 < abs(J4 − J5)<br />

10 warning('second J3 triange condition not s<strong>at</strong>isfied. J3 > J4 + ...<br />

J5<br />

|| J3 < abs(J4 − J5)');<br />

11 sum = 0;<br />

12 return;<br />

13 end<br />

14 if J6 > J2 + J4 || J6 < abs(J2 − J4)<br />

15 warning('first J6 triange condition not s<strong>at</strong>isfied. J6 > J2 + ...<br />

J4<br />

|| J6 < abs(J2 − J4)');<br />

16 sum = 0;<br />

17 return;<br />

18 end<br />

19 if J6 > J1 + J5 || J6 < abs(J1 − J5)<br />

20 warning('second J6 triange condition not s<strong>at</strong>isfied. J6 > J1 + ...<br />

J5<br />

|| J6 < abs(J1 − J5)');<br />

21 sum = 0;<br />

22 return;<br />

98

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!