plots_for_raster_memo

0001 load('tau_newcal_apr14/06-Feb-2014:01:24:02_taua_I_pipelined.mat')
0002 
0003 % Extract raster scan data ONLY
0004 d.index.radio_point_scan.fast = d.index.beammap.fast;
0005 d.index.radio_point_scan.slow = d.index.beammap.slow;
0006 df = framecut(d,d.index.radio_point_scan.fast);
0007 
0008 tdate=(mjd2date_v2(df.antenna0.receiver.utc));% get in date form
0009 aa  = datenum(tdate.year,tdate.month,tdate.day,tdate.hour,tdate.minute,tdate.second);
0010 
0011 %Plot the raw data
0012 % .data is currently [I1,Q1,U1,Q2,U2,Q3,U3,I2]
0013 % We will plot I1,Q1,U1,I2,Q2,U2]
0014 
0015 plot_labels={'I1','Q1','U1','Q2','U2','Q3','U3','I2'}
0016 j=1
0017 for i=[1 2 3 8 4 5];
0018     subplot(2,3,j);
0019     plot(aa,df.antenna0.receiver.data(:,i));
0020     %xlabel(plot_labels(i));
0021     j=j+1;
0022     
0023     xlabel('Time,hours');
0024     ylabel('Astro units, K');
0025     title(plot_labels(i));
0026     xlim([min(aa) max(aa)]);
0027     datetick('x','hh','keeplimits','keepticks');
0028 end
0029 
0030 %% Now get detrended data
0031 clear
0032 close all
0033 
0034 load('tau_newcal_apr14/06-Feb-2014:01:24:02_taua_I_dgood.mat')
0035 d=dgood;
0036 clear dgood
0037 df = framecut(d,d.index.radio_point_scan.fast);
0038 clear d;
0039 
0040 % Flag spike in data
0041 abs_time=(df.antenna0.receiver.utc(:)-df.antenna0.receiver.utc(1));
0042 flagged=(abs_time >0.0327& abs_time<0.0331);
0043 
0044 
0045 tdate=(mjd2date_v2(df.antenna0.receiver.utc));% get in date form
0046 aa  = datenum(tdate.year,tdate.month,tdate.day,tdate.hour,tdate.minute,tdate.second);
0047 j=1
0048 plot_labels={'I1','Q1','U1','Q2','U2','Q3','U3','I2'}
0049 for i=[1 2 3 8 4 5];
0050     subplot(2,3,j);
0051     plot(aa(~flagged),df.antenna0.receiver.data(~flagged,i));
0052     %xlabel(plot_labels(i));
0053     j=j+1;
0054     
0055     xlabel('Time,hours');
0056     ylabel('Astro units, K');
0057     title(plot_labels(i));
0058     xlim([min(aa) max(aa)]);
0059     datetick('x','hh','keeplimits','keepticks');
0060 end
0061 
0062 %% Save dgood (detrended) as .fits file
0063 
0064 clear
0065 close all
0066 
0067 load('tau_newcal_apr14/06-Feb-2014:01:24:02_taua_I_dgood.mat')
0068 d = dgood;
0069 d.antenna0.receiver.data(:,9)=[];
0070 datadir_fits='n_raster_elliptical_newcal_apr14'
0071 dred = reduceData(d,'redScript_newcal_raster_2014_fits.red',0,datadir_fits)
0072