Geant4 Simulations for the Radon Electric Dipole Moment Search at

More documents

Recommendations

Info

1 % RacahW.m 2 function out = RacahW(a,b,c,d,e,f) 3 out = ((−1)ˆ(a+b+d+c))*Wigner6j(a,b,e,d,c,f); 1 % Rho.m 2 function out = Rho(k,ji,pops) 3 sum = 0; 4 J = abs(ji); 5 <strong>for</strong> mi=−J:1:J 6 mdex = mi+J; 7 M = mi; 8 pop = pops(mdex+1); 9 sum = sum + ((−1)ˆ(J−M))*ClebschGordan(J,M,J,−1*M,k,0)*pop; 10 end 11 out =((2*J+1)ˆ(1/2))*sum; 1 % Wigner3j.m 2 function out = Wigner3j(j1,j2,j3,m1,m2,m3) 3 % error checking ////////// 4 if ( 2*j1 ≠ floor(2*j1) || 2*j2 ≠ floor(2*j2) || 2*j3 ≠ ... floor(2*j3)... 5 || 2*m1 ≠ floor(2*m1) || 2*m2 ≠ floor(2*m2)... 6 || 2*m3 ≠ floor(2*m3) ) 7 error('All arguments must be integers or half−integers.'); 8 return; 9 end 10 if m1 + m2 + m3 ≠ 0 96
11 %warning('m1 + m2 + m3 must equal zero.'); 12 out = 0; 13 return; 14 end 15 if ( j1 + j2 + j3 ≠ floor(j1 + j2 + j3) ) 16 %warning('2*j1 and 2*m1 must have <strong>the</strong> same parity'); 17 out = 0; 18 return; 19 end 20 if j3 > j1 + j2 || j3 < abs(j1 − j2) 21 %warning('j3 is out of bounds.'); 22 out = 0; 23 return; 24 end 25 if abs(m1) > j1 26 %warning('m1 is out of bounds.'); 27 out = 0; 28 return; 29 end 30 if abs(m2) > j2 31 %warning('m2 is out of bounds.'); 32 out = 0; 33 return; 34 end 35 if abs(m3) > j3 36 %warning('m3 is out of bounds.'); 37 out = 0; 38 return; 39 end 40 % ////////// 97
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 and 104: 1 % LegendreP.m 2 function P = Lege
Page 105: 39 xx(i+1,j+1) = r(i+1,1)*sin((i*re
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?