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.

33 if L1 == L2 && ∆ == 0<br />

34 title( sprintf( 'Angular Profile Partially Aligned Nuclei For ...<br />

L%i Transition with Ji = %d and Jf = %d', L1,ji,jf ) );<br />

35 else<br />

36 title( sprintf( 'Angular Profile For Partially Aligned Nuclei ...<br />

For L%i, L%i (∆ = %d) Transition with Ji = %d and Jf = ...<br />

%d', L1,L2,∆,ji,jf) );<br />

37 end<br />

38 % Reference: Tables of Coefficients <strong>for</strong> Angular Distribution of<br />

39 % Gamma Rays From Aligned Nuclei, T. Yamazaki<br />

40 % Nuclear D<strong>at</strong>a A (1967), 3, 1−23<br />

1 % W Plot3DGammaAngularDistributions.m<br />

2 function [xy] = W Plot3DGammaAngularDistributions(ji,jf,L1,L2,∆,pops)<br />

3 minTheta=0;<br />

4 stepTheta=pi/180;<br />

5 maxTheta=pi;<br />

6 <strong>the</strong>ta =minTheta:stepTheta:maxTheta;<br />

7 d<strong>at</strong>a = W Plot2DGammaAngularDistributions(ji,jf,L1,L2,∆,minTheta,...<br />

8 stepTheta,maxTheta,pops);<br />

9 d<strong>at</strong>a3D=zeros(length(<strong>the</strong>ta),361);<br />

10 <strong>for</strong> i=1:1:length(<strong>the</strong>ta)<br />

11 <strong>for</strong> j=1:1:361<br />

12 d<strong>at</strong>a3D(i,j) = d<strong>at</strong>a(i,2);<br />

13 end<br />

14 end<br />

15 Plot3DSurface(d<strong>at</strong>a3D);<br />

16 % Reference: Tables of Coefficients <strong>for</strong> Angular Distribution of<br />

17 % Gamma Rays From Aligned Nuclei, T. Yamazaki<br />

101

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

Saved successfully!

Ooh no, something went wrong!