Geant4 Simulations for the Radon Electric Dipole Moment Search at

More documents

Recommendations

Info

9 for j=0:1:360 10 sum = sum + data(i+1,j+1); 11 end 12 end 13 for i=0:1:180 14 for j=0:1:360 15 out(i+1) = out(i+1) + data(i+1,j+1)/sum; 16 end 17 end 18 indexTheta = 1; 19 for i=0:1:180 20 if(i == indexTheta*resTheta && i ≠180) 21 outRes(indexTheta) = outRes(indexTheta) + out(i+1); 22 indexTheta = indexTheta + 1; 23 end 24 outRes(indexTheta) = outRes(indexTheta) + out(i+1); 25 end 26 for i=0:1:(180/resTheta)−1 27 for j=0:1:(360/resPhi)−1 28 surfarea = 1; 29 resSurface(i+1,j+1) = outRes(i+1)/surfarea; 30 end 31 end 32 r = resSurface; 33 xx = zeros(size(resSurface)); 34 yy = zeros(size(resSurface)); 35 zz = zeros(size(resSurface)); 36 for i=0:1:((180/resTheta)−1) 37 for j=0:1:(360/resPhi) 38 if(j == (360/resPhi)) % make it the same as theta = 0; 94
39 xx(i+1,j+1) = r(i+1,1)*sin((i*resTheta*pi/180.0)+... 40 ((resTheta/2)*pi/180.0))*cos((0*pi/180.0)); 41 yy(i+1,j+1) = r(i+1,1)*sin((i*resTheta*pi/180.0)+... 42 ((resTheta/2)*pi/180.0))*sin((0*pi/180.0)); 43 zz(i+1,j+1) = r(i+1,1)*cos((i*resTheta*pi/180.0)+... 44 ((resTheta/2)*pi/180.0)); 45 else 46 xx(i+1,j+1) = r(i+1,j+1)*sin((i*resTheta*pi/180.0)+... 47 ((resTheta/2)*pi/180.0))*cos((j*resPhi*pi/180.0)); 48 yy(i+1,j+1) = r(i+1,j+1)*sin((i*resTheta*pi/180.0)+... 49 ((resTheta/2)*pi/180.0))*sin((j*resPhi*pi/180.0)); 50 zz(i+1,j+1) = r(i+1,j+1)*cos((i*resTheta*pi/180.0)+... 51 ((resTheta/2)*pi/180.0)); 52 end 53 end 54 end 55 clf 56 h1 = surf(xx,yy,zz); 57 set(h1,'facealpha',0.85); 58 colormap cool 59 light 60 lighting phong 61 axis tight equal off 62 % xlabel('Normalized Units') 63 % ylabel('Normalized Units') 64 % zlabel('Normalized Units') 65 view(40,30) 66 camzoom(1.5) 67 fh = figure(1); % returns the handle to the figure object 68 set(fh, 'color', 'white'); % sets the color to white 95
Page 1 and 2:
GEANT4 SIMULATIONS FOR THE RADON EL
Page 3 and 4:
Dictated but not read. i
Page 5 and 6:
Contents Acknowledgements ii 1 Intr
Page 7 and 8:
4 Results 60 4.1 γ-Ray Spectroscop
Page 9 and 10:
List of Figures 1.1 An illustration
Page 11 and 12:
Chapter 1 Introduction Thesearchfor
Page 13 and 14:
ˆP ր ˆT ց Figure 1.1: An illust
Page 15 and 16:
1.1.1 CP Violation Following the ob
Page 17 and 18:
Figure 1.2: Experimental upper limi
Page 19 and 20:
Figure 1.3: Illustration of the dou
Page 21 and 22:
Chapter 2 RnEDM Experiment at TRIUM
Page 23 and 24:
Figure 2.1: A schematic of the ISAC
Page 25 and 26:
a measurement cell. The RnEDM exper
Page 27 and 28:
NMR techniques, and supplementing s
Page 29 and 30:
Occupied Orbitals Fine Structure Hy
Page 31 and 32:
energy levels are separated accordi
Page 33 and 34:
(σ + ) along the axis of the B fie
Page 35 and 36:
atom and absorbed by another. Radia
Page 37 and 38:
anaverageof120kHzphotopeakcountrate
Page 39 and 40:
Figure 2.7: One GRIFFIN/TIGRESS HPG
Page 41 and 42:
(a) One GRIFFIN detector in the HPG
Page 43 and 44:
3.1 Introduction 3.1.1 Previous Wor
Page 45 and 46:
energies, branching ratios, emissio
Page 47 and 48:
used around the measurement cell in
Page 49 and 50:
3.2.5 Simulated Data Thecodewaswrit
Page 51 and 52:
lower its total energy. Internal co
Page 53 and 54: 10000 1000 2 χ red = 1.21 t 1/2 =
Page 55 and 56: as [47] I(E) = G 2π 3 7 c 6 | M i
Page 57 and 58: was to use an average X-ray energy
Page 59 and 60: 1.5 1.4 1.3 1.2 P = 100% P = 75% P
Page 61 and 62: where (a b c d | e f) are Clebsch-G
Page 63 and 64: 2 2 1.5 J i = 5/2, J f = 5/2 J i =
Page 65 and 66: of radius 10 cm, allowing the GRIFF
Page 67 and 68: (a) Screenshot showing the detector
Page 69 and 70: 3.6.2 Output Data and Sort Codes Th
Page 71 and 72: (user-defined option), secondly toa
Page 73 and 74: 30 Absolute Efficiency (%) 25 20 15
Page 75 and 76: 11 10 9 8 7 6 5 4 3 2 1 0 1e+06 Cou
Page 77 and 78: 4.2.1 Multipolarity Effects The sha
Page 79 and 80: Figure 4.7: The time projections an
Page 81 and 82: distribution, this average intensit
Page 83 and 84: Table 4.1: The fit results for the
Page 85 and 86: 12 10 linear regression: y = 0.2243
Page 87 and 88: 15 Linear Counts (e+04) 10 5 0 1e+0
Page 89 and 90: Chapter 5 Conclusions and Future Di
Page 91 and 92: simulated given sufficient parent a
Page 93 and 94: Appendix A Figure A.1: Simulated de
Page 95 and 96: Figure A.3: Simulated decay scheme
Page 97 and 98: Figure A.5: Simulated decay scheme
Page 99 and 100: 6 out = (1/(1+∆ˆ2))*(f(k,jf,L1,L
Page 101 and 102: 44 return; 45 end 46 if abs(m) > j
Page 103: 1 % LegendreP.m 2 function P = Lege
Page 107 and 108: 11 %warning('m1 + m2 + m3 must equa
Page 109 and 110: 23 end 24 % ////////// 25 j1 = J1;
Page 111 and 112: 33 if L1 == L2 && ∆ == 0 34 title
Page 113 and 114: Bibliography [1] J. M. Pendlebury a
Page 115 and 116: [25] P. G. Bricault, M. Dombsky, Je
Page 117: [47] RichardA.Dunlap. An introducti
show all

Geant4 Simulations for the Radon Electric Dipole Moment Search at

Create successful ePaper yourself

Delete template?

Save as template?