lukesReadArcSouthCompare

% SOUTHERN SYSTEM
%

0001 %% SOUTHERN SYSTEM
0002 %%
0003 
0004 %% Read in time as fun of time read in
0005 t0 = cputime
0006 d=read_arcSouth('15-Jan-2014:19:00:00','15-Jan-2014:19:05:00')
0007 t1 = cputime
0008 d=read_arcSouth('15-Jan-2014:19:10:00','15-Jan-2014:19:20:00')
0009 t2 = cputime
0010 d=read_arcSouth('15-Jan-2014:19:30:00','15-Jan-2014:19:50:00')
0011 t3 = cputime
0012 d=read_arcSouth('15-Jan-2014:20:00:00','15-Jan-2014:20:40:00')
0013 t4 = cputime
0014 d=read_arcSouth('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0015 t5 = cputime
0016 d=read_arcSouth('15-Jan-2014:22:00:00','16-Jan-2014:00:00:00')
0017 t6 = cputime
0018 d=read_arcSouth('16-Jan-2014:01:00:00','16-Jan-2014:04:00:00')
0019 t7 = cputime
0020 d=read_arcSouth('15-Jan-2014:19:00:00','15-Jan-2014:23:00:00')
0021 t8 = cputime
0022 
0023 
0024 ta = [t1-t0,t2-t1,t3-t2,t4-t3,t5-t4,t6-t5,t7-t6,t8-t7];
0025 minutes = [5,10,20,40,60,120,180,240];
0026 
0027 figure
0028 plot( minutes, ta,'o' )
0029 xlabel( 'Minutes read in' )
0030 ylabel( 'cputime to run read_arcSouth' )
0031 title( 'Time to read in different amounts of southern data' )
0032 
0033 %% Reading in less data
0034 t17 = cputime
0035 d1 = read_arcSouthLL('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0036 t18 = cputime
0037 d1 = read_arcSouthLLRR('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0038 t19 = cputime
0039 d1 = read_arcSouthLLRRload1('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0040 t20 = cputime
0041 d1 = read_arcSouthLLRRload12('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0042 t21 = cputime
0043 d1 = read_arcSouthLLRRload12Q('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0044 t22 = cputime
0045 d1 = read_arcSouthLLRRload12QU('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0046 t23 = cputime
0047 d1 = read_arcSouthLL_2_RRload12QU('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0048 t24 = cputime
0049 d1 = read_arcSouthLLRR_2_load12QU('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0050 t25 = cputime
0051 d1 = read_arcSouthLLRRload1_2_2QU('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0052 t26 = cputime
0053 d1 = read_arcSouthLLRRload12_2_QU('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0054 t27 = cputime
0055 d1 = read_arcSouthLLRRload12Q_2_U('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0056 t28 = cputime
0057 d1 = read_arcSouthLLRRload12QU_2('15-Jan-2014:22:00:00','15-Jan-2014:23:00:00')
0058 t29 = cputime
0059 
0060 tb = [t5-t4,t18-t17,t19-t18,t20-t19,t21-t20,t22-t21,t23-t22,t24-t23,t25-t24,t26-t25,t27-t26,t28-t27,t29-t28]
0061 channels = [0,1,2,3,4,5,6,7,8,9,10,11,12]*64
0062 
0063 figure
0064 plot( channels, tb, 'o' )
0065 xlabel( 'No channels excluded' )
0066 ylabel( 'cputime to run read_arcSouth...' )
0067 title( 'Excluding LL, RR, load1, load2, Q and U from read_arcSouth (each with 64 channels) (roach1 then roach2)' )
0068 
0069 %% ROACH 1 LL vs LLfreq
0070 LLmean = mean(d.antenna0.roach1.LL' );
0071 %figure
0072 %plot(LLmean,'g')
0073 %hold all
0074 %plot(d.antenna0.roach1.LLfreq,'r')
0075 %ylabel( 'Raw units' )
0076 %xlabel( 'Index' )
0077 
0078 LLm = mean( reshape(LLmean, 100, []) );
0079 LLy = LLm( ones(100,1), :);
0080 LLy = LLy(:)';
0081 
0082 figure
0083 plot( LLy,'g' )
0084 hold
0085 plot(d.antenna0.roach1.LLfreq,'r')
0086 ylabel( 'Raw units' )
0087 xlabel( 'Index' )
0088 title( 'LLfreq (red) vs calculation of LLfreq from LL (green) (roach1)' )
0089 
0090 %% ROACH 1 RR vs RRfreq
0091 RRmean = mean(d.antenna0.roach1.RR' );
0092 %figure
0093 %plot(RRmean,'g')
0094 %hold all
0095 %plot(d.antenna0.roach1.RRfreq,'r')
0096 %ylabel( 'Raw units' )
0097 %xlabel( 'Index' )
0098 
0099 RRm = mean( reshape(RRmean, 100, []) );
0100 RRy = RRm( ones(100,1), :);
0101 RRy = RRy(:)';
0102 
0103 figure
0104 plot( RRy,'g' )
0105 hold
0106 plot(d.antenna0.roach1.RRfreq,'r')
0107 ylabel( 'Raw units' )
0108 xlabel( 'Index' )
0109 title( 'RRfreq (red) vs calculation of RRfreq from RR (green) (roach1)' )
0110 
0111 %% ROACH 2 LL vs LLfreq
0112 LLmean = mean(d.antenna0.roach2.LL' );
0113 %figure
0114 %plot(LLmean,'g')
0115 %hold all
0116 %plot(d.antenna0.roach1.LLfreq,'r')
0117 %ylabel( 'Raw units' )
0118 %xlabel( 'Index' )
0119 
0120 LLm = mean( reshape(LLmean, 100, []) );
0121 LLy = LLm( ones(100,1), :);
0122 LLy = LLy(:)';
0123 
0124 figure
0125 plot( LLy,'g' )
0126 hold
0127 plot(d.antenna0.roach2.LLfreq,'r')
0128 ylabel( 'Raw units' )
0129 xlabel( 'Index' )
0130 title( 'LLfreq (red) vs calculation of LLfreq from LL (green) (roach2)' )
0131 
0132 %% ROACH 2 RR vs RRfreq
0133 RRmean = mean(d.antenna0.roach2.RR' );
0134 %figure
0135 %plot(RRmean,'g')
0136 %hold all
0137 %plot(d.antenna0.roach1.RRfreq,'r')
0138 %ylabel( 'Raw units' )
0139 %xlabel( 'Index' )
0140 
0141 RRm = mean( reshape(RRmean, 100, []) );
0142 RRy = RRm( ones(100,1), :);
0143 RRy = RRy(:)';
0144 
0145 figure
0146 plot( RRy,'g' )
0147 hold
0148 plot(d.antenna0.roach2.RRfreq,'r')
0149 ylabel( 'Raw units' )
0150 xlabel( 'Index' )
0151 title( 'RRfreq (red) vs calculation of RRfreq from RR (green) (roach2)' )
0152 
0153 %% NORTHERN SYSTEM
0154 %%
0155 
0156 %% Reading in different lengths of data
0157 t40 = cputime
0158 d=read_arc('15-Jan-2014:19:00:00','15-Jan-2014:19:05:00')
0159 t41 = cputime
0160 d=read_arc('15-Jan-2014:19:10:00','15-Jan-2014:19:20:00')
0161 t42 = cputime
0162 d=read_arc('15-Jan-2014:19:30:00','15-Jan-2014:19:50:00')
0163 t43 = cputime
0164 d=read_arc('15-Jan-2014:20:00:00','15-Jan-2014:20:40:00')
0165 t44 = cputime
0166 d=read_arc('15-Jan-2014:21:00:00','15-Jan-2014:22:00:00')
0167 t45 = cputime
0168 d=read_arc('15-Jan-2014:22:00:00','16-Jan-2014:00:00:00')
0169 t46 = cputime
0170 d=read_arc('16-Jan-2014:01:00:00','16-Jan-2014:04:00:00')
0171 t47 = cputime
0172 d=read_arc('15-Jan-2014:19:00:00','15-Jan-2014:23:00:00')
0173 t48 = cputime
0174 
0175 
0176 ta = [t41-t40,t42-t41,t43-t42,t44-t43,t45-t44,t46-t45,t47-t46,t48-t47];
0177 minutes = [5,10,20,40,60,120,180,240];
0178 
0179 figure
0180 plot( minutes, ta,'o' )
0181 xlabel( 'Minutes read in' )
0182 ylabel( 'cputime to run read_arcSouth' )
0183 title( 'Time to read in different amounts of northern data' )
0184