tsysCalculation

function to calculate the Tsys of the system using the "traditional" Tsys formula 
 (i.e. averaging sky and load and warm and load).
inputs dhot - data sequence holding the hot data
dcold - data sequence with the cold data
thot - nominal hot temperature
tcold - nominal cold temperature

0001 %function to calculate the Tsys of the system using the "traditional" Tsys formula
0002 % (i.e. averaging sky and load and warm and load).
0003 %inputs dhot - data sequence holding the hot data
0004 %dcold - data sequence with the cold data
0005 %thot - nominal hot temperature
0006 %tcold - nominal cold temperature
0007 
0008 % example tsys=tsysCalculation(d_30k,d_12k,d_30k.antenna0.thermal.ccTemperatureLoad,d_12k.antenna0.thermal.ccTemperatureLoad)
0009 %plot(tsys(1,:),tsys(2,:))
0010 %hold all
0011 %plot(tsys(1,:),tsys(3,:))
0012 %title('Tsys with Freq','Fontsize',16)
0013 %xlabel('Frequency [GHz]','Fontsize',14)
0014 %ylabel('Tsys [K]','Fontsize',14)
0015 %print('-depsc','Tsys.eps')
0016 
0017   function [Tsys] = tsysCalculation_trad_jl(dhot,dcold,thot,tcold,tload)
0018 close all;
0019 timeroachhot=(dhot.antenna0.roach1.utc-dhot.antenna0.roach1.utc(1)).*24.*3600;
0020 timetempshot=dhot.antenna0.thermal.utc-dhot.antenna0.thermal.utc(1);
0021 timelsTempshot= timetempshot([1:5:length(timetempshot)]).*24;
0022 lengthvechot = length(dhot.antenna0.roach1.RR)
0023 timeroachshorthot=timeroachhot(1:lengthvechot-0);
0024 
0025 lengthvec=lengthvechot;
0026 
0027 ll1hot=mean(dhot.antenna0.roach1.LL((1:lengthvec-0),1:64),1);
0028 loadl1hot =mean(dhot.antenna0.roach1.load2((1:lengthvec-0),1:64),1);
0029 rr1hot=mean(dhot.antenna0.roach1.RR((1:lengthvec-0),1:64),1);
0030 loadr1hot =mean(dhot.antenna0.roach1.load1((1:lengthvec-0),1:64),1);
0031 Q1hot=mean(dhot.antenna0.roach1.Q((1:lengthvec-0),1:64),1);
0032 U1hot=mean(dhot.antenna0.roach1.U((1:lengthvec-0),1:64),1);
0033 
0034 ll2hot=mean(dhot.antenna0.roach2.LL((1:lengthvec-0),1:64),1);
0035 loadl2hot =mean(dhot.antenna0.roach2.load2((1:lengthvec-0),1:64),1);
0036 rr2hot=mean(dhot.antenna0.roach2.RR((1:lengthvec-0),1:64),1);
0037 loadr2hot =mean(dhot.antenna0.roach2.load1((1:lengthvec-0),1:64),1);
0038 Q2hot=mean(dhot.antenna0.roach2.Q((1:lengthvec-0),1:64),1);
0039 U2hot=mean(dhot.antenna0.roach2.U((1:lengthvec-0),1:64),1);
0040 
0041 
0042 timeroachcold=(dcold.antenna0.roach1.utc-dcold.antenna0.roach1.utc(1)).*24.*3600;
0043 timetempscold=dcold.antenna0.thermal.utc-dcold.antenna0.thermal.utc(1);
0044 timelsTempscold= timetempscold([1:5:length(timetempscold)]).*24;
0045 lengthveccold = length(dcold.antenna0.roach1.RR)
0046 timeroachshortcold=timeroachcold(1:lengthveccold-0);
0047 
0048 lengthvec=lengthveccold;
0049 ll1cold=mean(dcold.antenna0.roach1.LL((1:lengthvec-0),1:64),1);
0050 loadl1cold =mean(dcold.antenna0.roach1.load2((1:lengthvec-0),1:64),1);
0051 rr1cold=mean(dcold.antenna0.roach1.RR((1:lengthvec-0),1:64),1);
0052 loadr1cold =mean(dcold.antenna0.roach1.load1((1:lengthvec-0),1:64),1);
0053 Q1cold=mean(dcold.antenna0.roach1.Q((1:lengthvec-0),1:64),1);
0054 U1cold=mean(dcold.antenna0.roach1.U((1:lengthvec-0),1:64),1);
0055 
0056 ll2cold=mean(dcold.antenna0.roach2.LL((1:lengthvec-0),1:64),1);
0057 loadl2cold =mean(dcold.antenna0.roach2.load2((1:lengthvec-0),1:64),1);
0058 rr2cold=mean(dcold.antenna0.roach2.RR((1:lengthvec-0),1:64),1);
0059 loadr2cold =mean(dcold.antenna0.roach2.load1((1:lengthvec-0),1:64),1);
0060 Q2cold=mean(dcold.antenna0.roach2.Q((1:lengthvec-0),1:64),1);
0061 U2cold=mean(dcold.antenna0.roach2.U((1:lengthvec-0),1:64),1);
0062 
0063 
0064 %need to flip the upper band as this is actually inverted frequency wise
0065 ll2hot=flipdim(ll2hot,2);
0066 rr2hot=flipdim(rr2hot,2);
0067 ll2cold=flipdim(ll2cold,2);
0068 rr2cold=flipdim(rr2cold,2);
0069 
0070 loadl2hot=flipdim(loadl2hot,2);
0071 loadr2hot=flipdim(loadr2hot,2);
0072 loadl2cold=flipdim(loadl2cold,2);
0073 loadr2cold=flipdim(loadr2cold,2);
0074 
0075 
0076 llhot=[ll1hot ll2hot];
0077 rrhot = [rr1hot rr2hot];
0078 llcold=[ll1cold ll2cold];
0079 rrcold=[rr1cold rr2cold];
0080 loadlhot=[loadl1hot loadl2hot];
0081 loadrhot=[loadr1hot loadr2hot];
0082 loadlcold=[loadl1cold loadl2cold];
0083 loadrcold=[loadr1cold loadr2cold];
0084 
0085 Yll=llhot./llcold;
0086 Yrr=rrhot./rrcold;
0087 Yloadr=loadrhot./loadrcold;
0088 Yloadl=loadlhot./loadlcold;
0089 
0090 thotrep=mean(thot);
0091 tcoldrep=mean(tcold);
0092 tloadrep=mean(tload);
0093 
0094 thotrepmat=repmat(thotrep,1,128);
0095 tcoldrepmat=repmat(tcoldrep,1,128);
0096 tloadrepmat=repmat(tloadrep,1,128);
0097 
0098 Tsysll=((thotrepmat+tloadrepmat)/2 - Yll.*((tcoldrepmat+tloadrepmat)/2))./(Yll-1); % ch1
0099 Tsysrr=((thotrepmat+tloadrepmat)/2 - Yrr.*((tcoldrepmat+tloadrepmat)/2))./(Yrr-1); % ch2
0100 Tsysloadr=((thotrepmat+tloadrepmat)/2 - Yloadr.*((tcoldrepmat+tloadrepmat)/2))./(Yloadr-1); % ch3
0101 Tsysloadl=((thotrepmat+tloadrepmat)/2 - Yloadl.*((tcoldrepmat+tloadrepmat)/2))./(Yloadl-1); % ch4
0102 
0103 freq=[1:1000./128:1000]+4500;
0104 Tsys=[freq./1000;Tsysrr;Tsysll;Tsysloadr;Tsysloadl;rrhot;rrcold;llhot;llcold;loadrhot;loadrcold;loadlhot;loadlcold];
0105 
0106 
0107 
0108