Geant4 Simulations for the Radon Electric Dipole Moment Search at

Appendix B
MATLABCode: γ-rayangulardistributions
for any transition and degree of orientation.
1 % Alpha.m
2 function out = Alpha(k,ji,pops)
3 rho = Rho(k,ji,pops);
4 b = B(k,ji);
5 % zero checking //////////
6 if (rho == 0 || b == 0)
7 out = 0;
8 return;
9 end
10 % //////////
11 out = (rho)/(b);
1 % Amax.m
2 function out = Amax(k,ji,L1,L2,jf,∆)
3 if ∆ == 0
4 out = (1/(1+∆ˆ2))*(f(k,jf,L1,L1,ji));
5 else
88