hotColdLoadTests_jl

script to generate Tsys with frequency
d_hot=read_arcSouth('20-Apr-2013:21:00:00','20-Apr-2013:21:01:00');% eccosorb in place
d_cold=read_arcSouth('20-Apr-2013:21:13:30','20-Apr-2013:21:14:30');% eccosorb off

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('23-Apr-2013:17:03:10','23-Apr-2013:17:12:56');% eccosorb in place
0006 d_cold=read_arcSouth('23-Apr-2013:17:12:56','23-Apr-2013:17:20:42');% eccosorb off
0007 
0008 %check that no noise diode are in the data
0009 figure;
0010 plot(d_hot.antenna0.roach1.LL)
0011 hold all;
0012 plot(d_cold.antenna0.roach1.LL)
0013 thot=290
0014 tcold=6
0015 tload=11
0016   tsys=tsysCalculation_trad_jl(d_hot,d_cold,thot,tcold,tload); % calculate the tsys
0017 
0018 thot=repmat(thot,1,2000);
0019 tcold=repmat(tcold,1,2000);
0020 
0021 
0022 figure
0023 subplot(411)
0024 plot(tsys(1,:),tsys(2,:))
0025 hold all
0026 plot(tsys(1,:),tsys(3,:))
0027 plot(tsys(1,:),tsys(4,:))
0028 plot(tsys(1,:),tsys(5,:))
0029 filename ='be_tsys_23apr.txt';
0030 trec_mat=horzcat(tsys(1,:)',tsys(2,:)',tsys(3,:)',tsys(4,:)',tsys(5,:)');
0031 dlmwrite(filename,trec_mat,' ');
0032 
0033 
0034 title('Tsys Eccosorb 23April2013','Fontsize',16)
0035 xlabel('Frequency [GHz]','Fontsize',14)
0036 ylabel('Tsys [K]','Fontsize',14)
0037 legend('ch1','ch2','ch3','ch4','Location','EastOutside')
0038 ylim([0 100]);
0039 grid on
0040 subplot(412)
0041 plot(thot)
0042 hold all
0043 plot(tcold)
0044 title('Assumed Temperature')
0045 legend('Hot','Cold','Location','EastOutside')
0046 xlabel('Sample')
0047 ylabel('Temperature [K]')
0048 grid on
0049 
0050 subplot(413)
0051 plot(tsys(1,:),tsys(6,:))
0052 hold all
0053 plot(tsys(1,:),tsys(7,:))
0054 plot(tsys(1,:),tsys(8,:))
0055 plot(tsys(1,:),tsys(9,:))
0056 title('Spectra')
0057 legend('RHot','RCold','LHot','LCold','Location','EastOutside')
0058 xlabel('Frequency [GHz]','Fontsize',14)
0059 ylabel('Spectra [B.U','Fontsize',14)
0060 grid on
0061 
0062 
0063 
0064 
0065 
0066 lengthvec = length(d_hot.antenna0.roach1.RR);
0067 
0068 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);
0069 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);
0070 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);
0071 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);
0072 
0073 
0074 
0075 plot(smooth(rfactorLL1,3));
0076 hold all
0077 plot(smooth(rfactorLL2,3));
0078 plot(smooth(rfactorRR1,3));
0079 plot(smooth(rfactorRR2,3));
0080 
0081 subplot(414)
0082 plot(mean(d_hot.antenna0.roach1.LL((1:lengthvec-0),4:60),2));
0083 hold all
0084 plot(mean(d_hot.antenna0.roach1.RR((1:lengthvec-0),4:60),2));
0085 plot(mean(d_hot.antenna0.roach2.LL((1:lengthvec-0),4:60),2));
0086 plot(mean(d_hot.antenna0.roach2.RR((1:lengthvec-0),4:60),2));
0087 lengthvec = length(d_cold.antenna0.roach1.RR);
0088 plot(mean(d_cold.antenna0.roach1.LL((1:lengthvec-0),4:60),2));
0089 plot(mean(d_cold.antenna0.roach1.RR((1:lengthvec-0),4:60),2));
0090 plot(mean(d_cold.antenna0.roach2.LL((1:lengthvec-0),4:60),2));
0091 plot(mean(d_cold.antenna0.roach2.RR((1:lengthvec-0),4:60),2));
0092 grid on
0093 legend('Hot#1 LL','Hot#1RR','Hot#2LL','Hot#2RR','Cold#1 LL','Cold#1RR','Cold#2LL','Cold#2RR','Location','EastOutside')
0094 title('Time Series')
0095 
0096 xlabel('Sample','Fontsize',14)
0097 ylabel('B.U','Fontsize',14)
0098 
0099 print('-depsc','afterBrokenOMTCable.eps')
0100 
0101 %% check how much of the load leaks
0102 plot(tsys(1,:),tsys(10,:))
0103 hold all
0104 plot(tsys(1,:),tsys(11,:))
0105 plot(tsys(1,:),tsys(12,:))
0106 plot(tsys(1,:),tsys(13,:))
0107 title('Spectra')
0108 legend('RHot','RCold','LHot','LCold','Location','EastOutside')
0109 xlabel('Frequency [GHz]','Fontsize',14)
0110 ylabel('Spectra [B.U','Fontsize',14)
0111 grid on
0112 saveas(gcf,'tsys_plot_jl.png');
0113