Geant4 Simulations for the Radon Electric Dipole Moment Search at

3 minTheta,stepTheta,maxTheta,pops)
4 Ji = abs(ji);
5 Jf = abs(jf);
6 Alpha2 = Alpha(2,Ji,pops);
7 Alpha4 = Alpha(4,Ji,pops);
8 A2max = Amax(2,Ji,L1,L2,Jf,∆);
9 A4max = Amax(4,Ji,L1,L2,Jf,∆);
10 A2 = Alpha2*A2max;
11 A4 = Alpha4*A4max;
12 sum = 0;
13 theta =minTheta:stepTheta:maxTheta;
14 xy=zeros(length(theta),2);
15 for i=1:1:length(theta)
16 xy(i,1) = theta(i);
17 xy(i,2) = 1 + A2*LegendreP(2,cos(theta(i))) +...
18 A4*LegendreP(4,cos(theta(i)));
19 sum = sum + sin(theta(i)).*xy(i,2)*stepTheta;
20 end
21 disp( sprintf( 'Amax%d = %d', 2,A2max ) );
22 disp( sprintf( 'Amax%d = %d', 4,A4max ) );
23 disp( sprintf( 'MyAlpha%d = %d', 2,Alpha2 ) );
24 disp( sprintf( 'MyAlpha%d = %d', 4,Alpha4 ) );
25 disp( sprintf( 'Integral = %d', sum ) );
26 colors={'−+k';'−−or';':*b';'−.xk';'−sr';'−−db';':ˆk';'−.vr';'−>b'};
27 plot(xy(:,1),xy(:,2),char(colors(1)))
28 hold on
29 set(gca,'XTick',0:pi/4:pi)
30 set(gca,'XTickLabel',{'0','pi/4','pi/2','3*pi/2','pi'})
31 xlabel('\theta')
32 ylabel('W')
100