|
| 1 | +%% ---------------------- MIT_convertBlank2Struct.m --------------------- %% |
| 2 | +% DESCRIPTION |
| 3 | +% Converts matrix of blank values into a structure. |
| 4 | +% |
| 5 | +% CONSTANTS |
| 6 | +% None |
| 7 | +% |
| 8 | +% DEPENDENCIES |
| 9 | +% None |
| 10 | +% ----------------------------------------------------------------------- % |
| 11 | +% AUTHORS: Christine Y Chen [CYC], David McGee [DM] |
| 12 | +% Created: 06/12/15 |
| 13 | +% Last updated: 06/12/15 |
| 14 | +%% ---------------------- REVISION HISTORY ------------------------------ % |
| 15 | +% 08/13/15 [CYC] Changes reflected due to the new raw summary file |
| 16 | +% for U-tail having fewer columns. |
| 17 | +% 06/12/15 [CYC] Created by CYC. Based on Excel reduction sheet from |
| 18 | +% Dec 2014. |
| 19 | +% |
| 20 | +%% --------------------------------------------------------------------- %% |
| 21 | +% ------------------ Put blank values into structure ------------------ % |
| 22 | + |
| 23 | +blank.spk_wt = blank_spkWt; |
| 24 | + |
| 25 | +% ---------------- Read in Th data for the blank ------------------------ % |
| 26 | +blankTh_r.Th230_Th232 = blank_Thsamp_data(1); |
| 27 | +blankTh_r.Th230_Th232_err = blank_Thsamp_data(2); |
| 28 | +blankTh_r.Th232_Th229 = blank_Thsamp_data(3); |
| 29 | +blankTh_r.Th232_Th229_err = blank_Thsamp_data(4); |
| 30 | +blankTh_r.Th230_Th229 = blank_Thsamp_data(5); |
| 31 | +blankTh_r.Th230_Th229_err = blank_Thsamp_data(6); |
| 32 | +blankTh_r.Th232_V = blank_Thsamp_data(7); |
| 33 | +blankTh_r.Th232_V_err = blank_Thsamp_data(8); |
| 34 | +blankTh_r.Th230_cps = blank_Thsamp_data(9); |
| 35 | +blankTh_r.Th230_cps_err = blank_Thsamp_data(10); |
| 36 | +blankTh_r.Th229_V = blank_Thsamp_data(11); |
| 37 | +blankTh_r.Th229_V_err = blank_Thsamp_data(12); |
| 38 | +blankTh_r.U238_V = blank_Thsamp_data(13); |
| 39 | +blankTh_r.U238_V_err = blank_Thsamp_data(14); |
| 40 | +blankTh_r.Th230_bg = blank_Thsamp_data(15); |
| 41 | +blankTh_r.Th230_bg_err = blank_Thsamp_data(16); |
| 42 | + |
| 43 | +% ----------- Read in Th data for Standard #1 of blank ------------------ % |
| 44 | +stdBlnkTh1_r.Th230_Th232 = blank_Thstd1_data(1); |
| 45 | +stdBlnkTh1_r.Th230_Th232_err = blank_Thstd1_data(2); |
| 46 | +stdBlnkTh1_r.Th232_Th229 = blank_Thstd1_data(3); |
| 47 | +stdBlnkTh1_r.Th232_Th229_err = blank_Thstd1_data(4); |
| 48 | +stdBlnkTh1_r.Th230_Th229 = blank_Thstd1_data(5); |
| 49 | +stdBlnkTh1_r.Th230_Th229_err = blank_Thstd1_data(6); |
| 50 | +stdBlnkTh1_r.Th232_V = blank_Thstd1_data(7); |
| 51 | +stdBlnkTh1_r.Th232_V_err = blank_Thstd1_data(8); |
| 52 | +stdBlnkTh1_r.Th230_cps = blank_Thstd1_data(9); |
| 53 | +stdBlnkTh1_r.Th230_cps_err = blank_Thstd1_data(10); |
| 54 | +stdBlnkTh1_r.Th229_V = blank_Thstd1_data(11); |
| 55 | +stdBlnkTh1_r.Th229_V_err = blank_Thstd1_data(12); |
| 56 | +stdBlnkTh1_r.U238_V = blank_Thstd1_data(13); |
| 57 | +stdBlnkTh1_r.U238_V_err = blank_Thstd1_data(14); |
| 58 | +stdBlnkTh1_r.Th230_bg = blank_Thstd1_data(15); |
| 59 | +stdBlnkTh1_r.Th230_bg_err = blank_Thstd1_data(16); |
| 60 | + |
| 61 | +% ------------ Read in Th data for Standard #2 of blank------------------ % |
| 62 | +stdBlnkTh2_r.Th230_Th232 = blank_Thstd2_data(1); |
| 63 | +stdBlnkTh2_r.Th230_Th232_err = blank_Thstd2_data(2); |
| 64 | +stdBlnkTh2_r.Th232_Th229 = blank_Thstd2_data(3); |
| 65 | +stdBlnkTh2_r.Th232_Th229_err = blank_Thstd2_data(4); |
| 66 | +stdBlnkTh2_r.Th230_Th229 = blank_Thstd2_data(5); |
| 67 | +stdBlnkTh2_r.Th230_Th229_err = blank_Thstd2_data(6); |
| 68 | +stdBlnkTh2_r.Th232_V = blank_Thstd2_data(7); |
| 69 | +stdBlnkTh2_r.Th232_V_err = blank_Thstd2_data(8); |
| 70 | +stdBlnkTh2_r.Th230_cps = blank_Thstd2_data(9); |
| 71 | +stdBlnkTh2_r.Th230_cps_err = blank_Thstd2_data(10); |
| 72 | +stdBlnkTh2_r.Th229_V = blank_Thstd2_data(11); |
| 73 | +stdBlnkTh2_r.Th229_V_err = blank_Thstd2_data(12); |
| 74 | +stdBlnkTh2_r.U238_V = blank_Thstd2_data(13); |
| 75 | +stdBlnkTh2_r.U238_V_err = blank_Thstd2_data(14); |
| 76 | +stdBlnkTh2_r.Th230_bg = blank_Thstd2_data(15); |
| 77 | +stdBlnkTh2_r.Th230_bg_err = blank_Thstd2_data(16); |
| 78 | + |
| 79 | +% ---------------- Read in Standard 1 U of the blank -------------------- % |
| 80 | +stdBlnkU1_r.U.U238_U236 = blank_Ustd1_data(1); |
| 81 | +stdBlnkU1_r.U.U238_U236_err = blank_Ustd1_data(2); |
| 82 | +stdBlnkU1_r.U.U238_U235 = blank_Ustd1_data(3); |
| 83 | +stdBlnkU1_r.U.U238_U235_err = blank_Ustd1_data(4); |
| 84 | +stdBlnkU1_r.U.U236_U233 = blank_Ustd1_data(5); |
| 85 | +stdBlnkU1_r.U.U236_U233_err = blank_Ustd1_data(6); |
| 86 | +stdBlnkU1_r.U.U238_U234 = blank_Ustd1_data(7); |
| 87 | +stdBlnkU1_r.U.U238_U234_err = blank_Ustd1_data(8); |
| 88 | +stdBlnkU1_r.U.U238_V = blank_Ustd1_data(9); |
| 89 | +stdBlnkU1_r.U.U238_V_err = blank_Ustd1_data(10); |
| 90 | +stdBlnkU1_r.U.U236_V = blank_Ustd1_data(11); |
| 91 | +stdBlnkU1_r.U.U236_V_err = blank_Ustd1_data(12); |
| 92 | +stdBlnkU1_r.U.U234_cps = blank_Ustd1_data(13); |
| 93 | +stdBlnkU1_r.U.U234_cps_err = blank_Ustd1_data(14); |
| 94 | + |
| 95 | +% -------------- Read in Standard 1 U tail of the blank ----------------- % |
| 96 | +stdBlnkU1_r.Utail.U238_V = blank_Ustd1_tail_data(1); |
| 97 | +stdBlnkU1_r.Utail.U238_V_err = blank_Ustd1_tail_data(2); |
| 98 | +stdBlnkU1_r.Utail.U236_V = blank_Ustd1_tail_data(3); |
| 99 | +stdBlnkU1_r.Utail.U236_V_err = blank_Ustd1_tail_data(4); |
| 100 | +stdBlnkU1_r.Utail.U237_cps = blank_Ustd1_tail_data(5); |
| 101 | +stdBlnkU1_r.Utail.U237_err = blank_Ustd1_tail_data(6); |
| 102 | +stdBlnkU1_r.Utail.U2365_cps = blank_Ustd1_tail_data(7); |
| 103 | +stdBlnkU1_r.Utail.U2365_err = blank_Ustd1_tail_data(8); |
| 104 | +stdBlnkU1_r.Utail.U2345_cps = blank_Ustd1_tail_data(9); |
| 105 | +stdBlnkU1_r.Utail.U2345_err = blank_Ustd1_tail_data(10); |
| 106 | +stdBlnkU1_r.Utail.U2335_cps = blank_Ustd1_tail_data(11); |
| 107 | +stdBlnkU1_r.Utail.U2335_err = blank_Ustd1_tail_data(12); |
| 108 | + |
| 109 | +% ---------------- Read in Standard 2 U tail of the blank --------------- % |
| 110 | +stdBlnkU2_r.Utail.U238_V = blank_Ustd2_tail_data(1); |
| 111 | +stdBlnkU2_r.Utail.U238_V_err = blank_Ustd2_tail_data(2); |
| 112 | +stdBlnkU2_r.Utail.U236_V = blank_Ustd2_tail_data(3); |
| 113 | +stdBlnkU2_r.Utail.U236_V_err = blank_Ustd2_tail_data(4); |
| 114 | +stdBlnkU2_r.Utail.U237_cps = blank_Ustd2_tail_data(5); |
| 115 | +stdBlnkU2_r.Utail.U237_err = blank_Ustd2_tail_data(6); |
| 116 | +stdBlnkU2_r.Utail.U2365_cps = blank_Ustd2_tail_data(7); |
| 117 | +stdBlnkU2_r.Utail.U2365_err = blank_Ustd2_tail_data(8); |
| 118 | +stdBlnkU2_r.Utail.U2345_cps = blank_Ustd2_tail_data(9); |
| 119 | +stdBlnkU2_r.Utail.U2345_err = blank_Ustd2_tail_data(10); |
| 120 | +stdBlnkU2_r.Utail.U2335_cps = blank_Ustd2_tail_data(11); |
| 121 | +stdBlnkU2_r.Utail.U2335_err = blank_Ustd2_tail_data(12); |
| 122 | + |
| 123 | +% --------------------- Read in U data for the blank -------------------- % |
| 124 | +blankU_r.U.U238_U236 = blank_Usamp_data(1); |
| 125 | +blankU_r.U.U238_U236_err = blank_Usamp_data(2); |
| 126 | +blankU_r.U.U238_U235 = blank_Usamp_data(3); |
| 127 | +blankU_r.U.U238_U235_err = blank_Usamp_data(4); |
| 128 | +blankU_r.U.U236_U233 = blank_Usamp_data(5); |
| 129 | +blankU_r.U.U236_U233_err = blank_Usamp_data(6); |
| 130 | +blankU_r.U.U238_U234 = blank_Usamp_data(7); |
| 131 | +blankU_r.U.U238_U234_err = blank_Usamp_data(8); |
| 132 | +blankU_r.U.U238_V = blank_Usamp_data(9); |
| 133 | +blankU_r.U.U238_V_err = blank_Usamp_data(10); |
| 134 | +blankU_r.U.U236_V = blank_Usamp_data(11); |
| 135 | +blankU_r.U.U236_V_err = blank_Usamp_data(12); |
| 136 | +blankU_r.U.U234_cps = blank_Usamp_data(13); |
| 137 | +blankU_r.U.U234_cps_err = blank_Usamp_data(14); |
| 138 | + |
| 139 | +% -------------------- Read in U tail of the blank ---------------------- % |
| 140 | +blankU_r.Utail.U238_V = blank_Usamp_data(9); |
| 141 | +blankU_r.Utail.U238_V_err = blank_Usamp_data(10); |
| 142 | +blankU_r.Utail.U236_V = blank_Usamp_data(11); |
| 143 | +blankU_r.Utail.U236_V_err = blank_Usamp_data(12); |
| 144 | +%blankU_r.Utail.U237_cps = blank_Usamp_tail_data(5); |
| 145 | +blankU_r.Utail.U237_cps = ((blank_Ustd1_tail_data(5)/blank_Ustd1_tail_data(1)) + (blank_Ustd2_tail_data(5)/blank_Ustd2_tail_data(1)))/2 * blank_Usamp_data(9); |
| 146 | +%blankU_r.Utail.U237_err = blank_Usamp_tail_data(6); |
| 147 | +blankU_r.Utail.U237_err = ((blank_Ustd1_tail_data(6)/blank_Ustd1_tail_data(5)) + (blank_Ustd2_tail_data(6)/blank_Ustd2_tail_data(5)))/2 * blankU_r.Utail.U237_cps; |
| 148 | +%blankU_r.Utail.U2365_cps = blank_Usamp_tail_data(7); |
| 149 | +blankU_r.Utail.U2365_cps = ((blank_Ustd1_tail_data(7)/blank_Ustd1_tail_data(1)) + (blank_Ustd2_tail_data(7)/blank_Ustd2_tail_data(1)))/2 * blank_Usamp_data(9); |
| 150 | +%blankU_r.Utail.U2365_err = blank_Usamp_tail_data(8); |
| 151 | +blankU_r.Utail.U2365_err = ((blank_Ustd1_tail_data(8)/blank_Ustd1_tail_data(7)) + (blank_Ustd2_tail_data(8)/blank_Ustd2_tail_data(7)))/2 * blankU_r.Utail.U2365_cps; |
| 152 | +%blankU_r.Utail.U2345_cps = blank_Usamp_tail_data(9); |
| 153 | +blankU_r.Utail.U2345_cps = ((blank_Ustd1_tail_data(9)/blank_Ustd1_tail_data(1)) + (blank_Ustd2_tail_data(9)/blank_Ustd2_tail_data(1)))/2 * blank_Usamp_data(9); |
| 154 | +%blankU_r.Utail.U2345_err = blank_Usamp_tail_data(10); |
| 155 | +blankU_r.Utail.U2345_err = ((blank_Ustd1_tail_data(10)/blank_Ustd1_tail_data(9)) + (blank_Ustd2_tail_data(10)/blank_Ustd2_tail_data(9)))/2 * blankU_r.Utail.U2345_cps; |
| 156 | +%blankU_r.Utail.U2335_cps = blank_Usamp_tail_data(11); |
| 157 | +blankU_r.Utail.U2335_cps = ((blank_Ustd1_tail_data(11)/blank_Ustd1_tail_data(1)) + (blank_Ustd2_tail_data(11)/blank_Ustd2_tail_data(1)))/2 * blank_Usamp_data(9); |
| 158 | +%blankU_r.Utail.U2335_err = blank_Usamp_tail_data(12); |
| 159 | +blankU_r.Utail.U2335_err = ((blank_Ustd1_tail_data(i, 12)/blank_Ustd1_tail_data(i, 11)) + (blank_Ustd2_tail_data(i, 12)/blank_Ustd2_tail_data(11)))/2 * sampleU.Utail.U2335_cps; |
| 160 | + |
| 161 | + sampleU.Utail.U237_cps = ((Ustd1_tail(i, 5)/Ustd1_tail(i, 1)) + (Ustd2_tail(i, 5)/Ustd2_tail(i, 1)))/2 * Usamp(i, 9); |
| 162 | + sampleU.Utail.U237_err = ((Ustd1_tail(i, 6)/Ustd1_tail(i, 5)) + (Ustd2_tail(i, 6)/Ustd2_tail(i, 5)))/2 * sampleU.Utail.U237_cps; |
| 163 | + sampleU.Utail.U2365_cps = ((Ustd1_tail(i, 7)/Ustd1_tail(i, 1)) + (Ustd2_tail(i, 7)/Ustd2_tail(i, 1)))/2 * Usamp(i, 9); |
| 164 | + sampleU.Utail.U2365_err = ((Ustd1_tail(i, 8)/Ustd1_tail(i, 7)) + (Ustd2_tail(i, 8)/Ustd2_tail(i, 7)))/2 * sampleU.Utail.U2365_cps; |
| 165 | + sampleU.Utail.U2345_cps = ((Ustd1_tail(i, 9)/Ustd1_tail(i, 1)) + (Ustd2_tail(i, 9)/Ustd2_tail(i, 1)))/2 * Usamp(i, 9); |
| 166 | + sampleU.Utail.U2345_err = ((Ustd1_tail(i, 10)/Ustd1_tail(i, 9)) + (Ustd2_tail(i, 10)/Ustd2_tail(i, 9)))/2 * sampleU.Utail.U2345_cps; |
| 167 | + sampleU.Utail.U2335_cps = ((Ustd1_tail(i, 11)/Ustd1_tail(i, 1)) + (Ustd2_tail(i, 11)/Ustd2_tail(i, 1)))/2 * Usamp(i, 9); |
| 168 | + sampleU.Utail.U2335_err = ((Ustd1_tail(i, 12)/Ustd1_tail(i, 11)) + (Ustd2_tail(i, 12)/Ustd2_tail(i, 11)))/2 * sampleU.Utail.U2335_cps; |
| 169 | + |
| 170 | +% ----------------- Read in Standard 2 U of the blank-------------------- % |
| 171 | +stdBlnkU2_r.U.U238_U236 = blank_Ustd2_data(1); |
| 172 | +stdBlnkU2_r.U.U238_U236_err = blank_Ustd2_data(2); |
| 173 | +stdBlnkU2_r.U.U238_U235 = blank_Ustd2_data(3); |
| 174 | +stdBlnkU2_r.U.U238_U235_err = blank_Ustd2_data(4); |
| 175 | +stdBlnkU2_r.U.U236_U233 = blank_Ustd2_data(5); |
| 176 | +stdBlnkU2_r.U.U236_U233_err = blank_Ustd2_data(6); |
| 177 | +stdBlnkU2_r.U.U238_U234 = blank_Ustd2_data(7); |
| 178 | +stdBlnkU2_r.U.U238_U234_err = blank_Ustd2_data(8); |
| 179 | +stdBlnkU2_r.U.U238_V = blank_Ustd2_data(9); |
| 180 | +stdBlnkU2_r.U.U238_V_err = blank_Ustd2_data(10); |
| 181 | +stdBlnkU2_r.U.U236_V = blank_Ustd2_data(11); |
| 182 | +stdBlnkU2_r.U.U236_V_err = blank_Ustd2_data(12); |
| 183 | +stdBlnkU2_r.U.U234_cps = blank_Ustd2_data(13); |
| 184 | +stdBlnkU2_r.U.U234_cps_err = blank_Ustd2_data(14); |
| 185 | + |
| 186 | + |
| 187 | + |
| 188 | +% Clear blank data that is not in a structure |
| 189 | +clear blank_Thsamp_data blank_Thstd1_data blank_Thstd2_data blank_Usamp_data ... |
| 190 | + blank_Usamp_tail_data blank_Ustd1_data blank_Ustd1_tail_data blank_Ustd2_data blank_Ustd2_tail_data |
0 commit comments