polarizationThermalNoisePlots

0001 d=read_arcSouth('15-Jan-2013:12:00:00','15-Jan-2013:13:30:00') % get the polarization BU/K
0002 
0003 %first we get the Q and U
0004 Q=mean(d.antenna0.roach1.Q(:,7:61),2);
0005 U=mean(d.antenna0.roach1.U(:,7:61),2);
0006 P=sqrt(Q.^2 + U.^2);
0007 Qon=mean(d.antenna0.roach1.Q(1.305e5:1.31e5,7:61),2)
0008 Uon=mean(d.antenna0.roach1.U(1.305e5:1.31e5,7:61),2)
0009 
0010 Qoff=mean(d.antenna0.roach1.Q(1.205e5:1.21e5,7:61),2)
0011 Uoff=mean(d.antenna0.roach1.U(1.205e5:1.21e5,7:61),2)
0012 
0013 ND1=0.9 % this is the total intensity noise diode in ch1
0014 ND2=0.6 % this is the total intensity noise diode in ch2
0015 NDTot= sqrt(ND1.^2 + ND2.^2)
0016 Pscale=0.0045 % scaling factor thatis adjustable to get things correct such that PND = NDtot
0017 
0018 QK=Q.*Pscale;
0019 UK=U.*Pscale;
0020 
0021 
0022 indextt=[1.205e5:1.21e5,7:61]; % index of the 100Hz Samples
0023 indextt2=[1.205e5:1.21e5,7:61]; % index of the downsampled 1Hz Samples
0024 
0025 QND=mean(Qon.*Pscale)-mean(Qoff.*Pscale)
0026 UND=mean(Uon.*Pscale)-mean(Uoff.*Pscale)
0027 PND = sqrt((QND.^2+UND.^2)) % This should be the same as NDtot
0028 QonK=mean(Pscale.*d.antenna0.roach1.Q(1.305e5:1.31e5,7:61),2);
0029 UonK=mean(Pscale.*d.antenna0.roach1.U(1.305e5:1.31e5,7:61),2);
0030 
0031 QoffK=mean(Pscale.*d.antenna0.roach1.Q(1.205e5:1.21e5,7:61),2);
0032 UoffK=mean(Pscale.*d.antenna0.roach1.U(1.205e5:1.21e5,7:61),2);
0033 
0034 
0035 PBUK = PND./1.081
0036 
0037 stdQ=std(QoffK)
0038 stdU=(std(UoffK))
0039 
0040 plot(QoffK)
0041 hold all
0042 plot(UoffK)
0043 
0044 
0045 
0046 stdP = sqrt(stdQ.^2 + stdU.^2)
0047 stdPK=stdP./PBUK
0048 
0049 subplot(311)
0050 plot(Q)
0051 title('Q')
0052 subplot(312)
0053 plot(U)
0054 title('Q')
0055 subplot(313)
0056 plot(P)
0057 title('P=sqrt(Q^2+U^2)')
0058 
0059 d=read_arcSouth('15-Jan-2013:12:30:00','15-Jan-2013:13:00:00') 
0060 
0061 %balanced load
0062 LL1=d.antenna0.roach1.LL-d.antenna0.roach1.load2;
0063 RR1=d.antenna0.roach1.RR-d.antenna0.roach1.load1;
0064 LL2=d.antenna0.roach2.LL-d.antenna0.roach2.load2;
0065 RR2=d.antenna0.roach2.RR-d.antenna0.roach2.load1;
0066 Q=d.antenna0.roach1.Q;
0067 U=d.antenna0.roach1.U;
0068 % to combine both use this
0069 %LLcomb = [LL1 LL2];
0070 %RRcomb = [RR1 RR2];
0071 
0072 LLcomb = [LL1];
0073 RRcomb = [RR1];
0074 Qcomb = [Q];
0075 Ucomb=[U];
0076 
0077 LLmean=mean(LLcomb,2);
0078 RRmean=mean(RRcomb,2);
0079 Imean=mean([LLcomb RRcomb],2);
0080 Qmean=mean(Qcomb,2);
0081 Umean=mean(Ucomb,2)
0082 
0083 indextt=[1:12000]; % index of the 100Hz Samples
0084 indextt2=[1000:1200]; % index of the downsampled 1Hz Samples
0085 
0086 
0087 
0088 BU2Kelvin=2500;
0089 BU2KelvinPol=220 % scaling factor thatis adjustable to get things correct such that PND = NDtot
0090 
0091 LLmeanK=LLmean./BU2Kelvin;
0092 RRmeanK=RRmean./BU2Kelvin;
0093 ImeanK=Imean./BU2Kelvin;
0094 QmeanK=Qmean./BU2KelvinPol;
0095 UmeanK=Umean./BU2KelvinPol;
0096 
0097 stdRa=std(RRmeanK(indextt))
0098 stdLa=std(LLmeanK(indextt))
0099 
0100 stdI=std(ImeanK(indextt))
0101 stqQ=std(QmeanK(indextt))
0102 stqU=std(UmeanK(indextt))
0103 
0104 
0105 
0106 
0107 
0108 %Ok now what about at different load temperatures
0109 d_12=read_arcSouth('15-Jan-2013:12:30:00','15-Jan-2013:12:35:00') % @ 12K both sky and load
0110 d_16=read_arcSouth('15-Jan-2013:14:20:00','15-Jan-2013:14:25:00') % @ 16K both sky and load
0111 d_18=read_arcSouth('15-Jan-2013:15:00:00','15-Jan-2013:15:05:00') % @ 18K both sky and load
0112 d_26=read_arcSouth('19-Jan-2013:12:00:00','19-Jan-2013:12:05:00') % @ 26K both sky and load
0113 
0114 Q_12=Pscale.*d_12.antenna0.roach1.Q(:,7:61);
0115 U_12=Pscale.*d_12.antenna0.roach1.U(:,7:61);
0116 Q_16=Pscale.*d_16.antenna0.roach1.Q(:,7:61);
0117 U_16=Pscale.*d_16.antenna0.roach1.U(:,7:61);
0118 Q_18=Pscale.*d_18.antenna0.roach1.Q(:,7:61);
0119 U_18=Pscale.*d_18.antenna0.roach1.U(:,7:61);
0120 Q_26=Pscale.*d_26.antenna0.roach1.Q(:,7:61);
0121 U_26=Pscale.*d_26.antenna0.roach1.U(:,7:61);
0122 
0123 Q_12b=Pscale.*d_12.antenna0.roach2.Q(:,7:61);
0124 U_12b=Pscale.*d_12.antenna0.roach2.U(:,7:61);
0125 Q_16b=Pscale.*d_16.antenna0.roach2.Q(:,7:61);
0126 U_16b=Pscale.*d_16.antenna0.roach2.U(:,7:61);
0127 Q_18b=Pscale.*d_18.antenna0.roach2.Q(:,7:61);
0128 U_18b=Pscale.*d_18.antenna0.roach2.U(:,7:61);
0129 Q_26b=Pscale.*d_26.antenna0.roach2.Q(:,7:61);
0130 U_26b=Pscale.*d_26.antenna0.roach2.U(:,7:61);
0131 
0132 Q_12comb=[Q_12 Q_12b];
0133 Q_16comb=[Q_16 Q_16b];
0134 Q_18comb=[Q_18 Q_18b];
0135 Q_26comb=[Q_26 Q_26b];
0136 U_12comb=[U_12 U_12b];
0137 U_16comb=[U_16 U_16b];
0138 U_18comb=[U_18 U_18b];
0139 U_26comb=[U_26 U_26b];
0140 
0141 Q_12K=mean(Q_12,2);
0142 U_12K=mean(U_12,2);
0143 Q_16K=mean(Q_16,2);
0144 U_16K=mean(U_16,2);
0145 Q_18K=mean(Q_18,2);
0146 U_18K=mean(U_18,2);
0147 Q_26K=mean(Q_26,2);
0148 U_26K=mean(U_26,2);
0149 
0150 Q_12K2=mean(Q_12b,2);
0151 U_12K2=mean(U_12b,2);
0152 Q_16K2=mean(Q_16b,2);
0153 U_16K2=mean(U_16b,2);
0154 Q_18K2=mean(Q_18b,2);
0155 U_18K2=mean(U_18b,2);
0156 Q_26K2=mean(Q_26b,2);
0157 U_26K2=mean(U_26b,2);
0158 
0159 Q_12Kcomb=mean(Q_12comb,2);
0160 U_12Kcomb=mean(U_12comb,2);
0161 Q_16Kcomb=mean(Q_16comb,2);
0162 U_16Kcomb=mean(U_16comb,2);
0163 Q_18Kcomb=mean(Q_18comb,2);
0164 U_18Kcomb=mean(U_18comb,2);
0165 Q_26Kcomb=mean(Q_26comb,2);
0166 U_26Kcomb=mean(U_26comb,2);
0167 
0168 
0169 plot(Q_12K)
0170 hold all
0171 plot(Q_12K2)
0172 plot(Q_12Kcomb)
0173 
0174 
0175 
0176 stdQ_12K = std(Q_12K)
0177 stdU_12K = std(U_12K)
0178 stdQ_16K = std(Q_16K)
0179 stdU_16K = std(U_16K)
0180 stdQ_18K = std(Q_18K)
0181 stdU_18K = std(U_18K)
0182 stdQ_26K = std(Q_26K)
0183 stdU_26K = std(U_26K)
0184 stdQ_12K2 = std(Q_12K2)
0185 stdU_12K2 = std(U_12K2)
0186 stdQ_16K2 = std(Q_16K2)
0187 stdU_16K2 = std(U_16K2)
0188 stdQ_18K2 = std(Q_18K2)
0189 stdU_18K2 = std(U_18K2)
0190 stdQ_26K2 = std(Q_26K2)
0191 stdU_26K2 = std(U_26K2)
0192 
0193 stdQComb_12K=std(Q_12Kcomb)
0194 stdUComb_12K=std(U_12Kcomb)
0195 stdQComb_16K=std(Q_16Kcomb)
0196 stdUComb_16K=std(U_16Kcomb)
0197 stdQComb_18K=std(Q_18Kcomb)
0198 stdUComb_18K=std(U_18Kcomb)
0199 stdQComb_26K=std(Q_26Kcomb)
0200 stdUComb_26K=std(U_26Kcomb)
0201 
0202 temps=[12 16 18 26]
0203 stdValsQ=[stdQ_12K stdQ_16K stdQ_18K stdQ_26K]
0204 stdValsU=[stdU_12K stdU_16K stdU_18K stdU_26K]
0205 stdValsQ2=[stdQ_12K2 stdQ_16K2 stdQ_18K2 stdQ_26K2]
0206 stdValsU2=[stdU_12K2 stdU_16K2 stdU_18K2 stdU_26K2]
0207 stdValsQComb=[stdQComb_12K stdQComb_16K stdQComb_18K stdQComb_26K]
0208 stdValsUComb=[stdUComb_12K stdUComb_16K stdUComb_18K stdUComb_26K]
0209 
0210 
0211 temps=[12 16 18 ]
0212 stdValsQ=[stdQ_12K stdQ_16K stdQ_18K ]
0213 stdValsU=[stdU_12K stdU_16K stdU_18K ]
0214 stdValsQ2=[stdQ_12K2 stdQ_16K2 stdQ_18K2 ]
0215 stdValsU2=[stdU_12K2 stdU_16K2 stdU_18K2 ]
0216 
0217 
0218 
0219 plot(temps,stdValsQ)
0220 hold all
0221 plot(temps,stdValsU)
0222 plot(temps,stdValsQ2)
0223 plot(temps,stdValsU2)
0224 plot(temps,stdValsQComb)
0225 plot(temps,stdValsUComb)
0226 
0227 deltaFreq=(61-7)./64.*(500*1e6)
0228 tsysExpected=12+15; %trx+tresistor
0229 tau=0.01
0230 sigmaExpected = tsysExpected./(sqrt(deltaFreq*tau))
0231 
0232 sigmaExpected2 = tsysExpected./(sqrt(deltaFreq*tau*2))
0233