hotColdLoadTests

script to generate Tsys with frequency

0001 %script to generate Tsys with frequency
0002 d_hot=read_arcSouth('20-Apr-2013:21:00:00','20-Apr-2013:21:01:00');% eccosorb in place
0003 d_cold=read_arcSouth('20-Apr-2013:21:13:30','20-Apr-2013:21:14:30');% eccosorb off
0004 
0005 d_hot=read_arcSouth('22-Apr-2013:14:40:00','22-Apr-2013:14:40:30');% eccosorb in place
0006 d_cold=read_arcSouth('22-Apr-2013:14:41:20','22-Apr-2013:14:41:50');% eccosorb off
0007 
0008 %check that no noise diode are in the data
0009 plot(d_hot.antenna0.roach1.LL)
0010 plot(d_cold.antenna0.roach1.LL)
0011 thot=290
0012 tcold=6
0013 tsys=tsysCalculation(d_hot,d_cold,thot,tcold) % calculate the tsys
0014 
0015 thot=repmat(thot,1,2000);
0016 tcold=repmat(tcold,1,2000);
0017 
0018 plot(tsys(1,:),tsys(2,:),'Linewidth',2)
0019 hold all
0020 plot(tsys(1,:),tsys(3,:),'Linewidth',2)
0021 set(gca,'FontSize',16)
0022 xlabel('Frequency [GHz]','Fontsize',16)
0023 ylabel('T_{rx} [K]','Fontsize',16)
0024 title('Digital Backend T_{rx}','Fontsize',18)
0025 legend('RCP','LCP')
0026 set(gca,'FontSize',16)
0027 grid on
0028 set(gca,'XTick',[4.5:0.1:5.5]) %
0029 
0030 print('-depsc','backendTrx')
0031 
0032 
0033 
0034 figure
0035 subplot(411)
0036 
0037 plot(tsys(1,:),tsys(8,:))
0038 hold all
0039 plot(tsys(1,:),tsys(9,:))
0040 plot(tsys(1,:),tsys(10,:))
0041 plot(tsys(1,:),tsys(11,:))
0042 title('Spectra')
0043 legend('RHot','RCold','LHot','LCold','Location','EastOutside')
0044 xlabel('Frequency [GHz]','Fontsize',14)
0045 ylabel('Spectra [B.U','Fontsize',14)
0046 grid on
0047 
0048 hold all
0049 plot(tsys(1,:),tsys(3,:))
0050 title('Tsys Eccosorb 22April2013','Fontsize',16)
0051 xlabel('Frequency [GHz]','Fontsize',14)
0052 ylabel('Tsys [K]','Fontsize',14)
0053 legend('RCP','LCP','Location','EastOutside')
0054 grid on
0055 subplot(412)
0056 plot(thot)
0057 hold all
0058 plot(tcold)
0059 title('Assumed Temperature')
0060 legend('Hot','Cold','Location','EastOutside')
0061 xlabel('Sample')
0062 ylabel('Temperature [K]')
0063 grid on
0064 
0065 subplot(413)
0066 plot(tsys(1,:),tsys(4,:))
0067 hold all
0068 plot(tsys(1,:),tsys(5,:))
0069 plot(tsys(1,:),tsys(6,:))
0070 plot(tsys(1,:),tsys(7,:))
0071 title('Spectra')
0072 legend('RHot','RCold','LHot','LCold','Location','EastOutside')
0073 xlabel('Frequency [GHz]','Fontsize',14)
0074 ylabel('Spectra [B.U','Fontsize',14)
0075 grid on
0076 
0077 
0078 
0079 
0080 
0081 lengthvec = length(d_hot.antenna0.roach1.RR)
0082 
0083 rfactorLL1 = mean(d_hot.antenna0.roach1.LL((1:lengthvec-0),4:60),1)./mean(d_hot.antenna0.roach1.load2((1:lengthvec-0),4:60),1)
0084 rfactorRR1 = mean(d_hot.antenna0.roach1.RR((1:lengthvec-0),4:60),1)./mean(d_hot.antenna0.roach1.load1((1:lengthvec-0),4:60),1)
0085 rfactorLL2 = mean(d_hot.antenna0.roach2.LL((1:lengthvec-0),4:60),1)./mean(d_hot.antenna0.roach2.load2((1:lengthvec-0),4:60),1)
0086 rfactorRR2 = mean(d_hot.antenna0.roach2.RR((1:lengthvec-0),4:60),1)./mean(d_hot.antenna0.roach2.load1((1:lengthvec-0),4:60),1)
0087 
0088 
0089 
0090 plot(smooth(rfactorLL1,3))
0091 hold all
0092 plot(smooth(rfactorLL2,3))
0093 plot(smooth(rfactorRR1,3))
0094 plot(smooth(rfactorRR2,3))
0095 
0096 subplot(414)
0097 plot(mean(d_hot.antenna0.roach1.LL((1:lengthvec-0),4:60),2));
0098 hold all
0099 plot(mean(d_hot.antenna0.roach1.RR((1:lengthvec-0),4:60),2));
0100 plot(mean(d_hot.antenna0.roach2.LL((1:lengthvec-0),4:60),2));
0101 plot(mean(d_hot.antenna0.roach2.RR((1:lengthvec-0),4:60),2));
0102 lengthvec = length(d_cold.antenna0.roach1.RR)
0103 plot(mean(d_cold.antenna0.roach1.LL((1:lengthvec-0),4:60),2));
0104 plot(mean(d_cold.antenna0.roach1.RR((1:lengthvec-0),4:60),2));
0105 plot(mean(d_cold.antenna0.roach2.LL((1:lengthvec-0),4:60),2));
0106 plot(mean(d_cold.antenna0.roach2.RR((1:lengthvec-0),4:60),2));
0107 grid on
0108 legend('Hot#1 LL','Hot#1RR','Hot#2LL','Hot#2RR','Cold#1 LL','Cold#1RR','Cold#2LL','Cold#2RR','Location','EastOutside')
0109 title('Time Series')
0110 
0111 xlabel('Sample','Fontsize',14)
0112 ylabel('B.U','Fontsize',14)
0113 
0114 print('-depsc','afterBrokenOMTCable.eps')
0115 
0116 %% check how much of the load leaks
0117 plot(tsys(1,:),tsys(8,:))
0118 hold all
0119 plot(tsys(1,:),tsys(9,:))
0120 plot(tsys(1,:),tsys(10,:))
0121 plot(tsys(1,:),tsys(11,:))
0122 title('Spectra')
0123 legend('RHot','RCold','LHot','LCold','Location','EastOutside')
0124 xlabel('Frequency [GHz]','Fontsize',14)
0125 ylabel('Spectra [B.U','Fontsize',14)
0126 grid on
0127 
0128