timingCheckCJC

0001 d=read_arcSouth('15-Sep-2013:05:40:00','15-Sep-2013:06:40:00')
0002 
0003 plot(86400*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),mean(d.antenna0.roach1.LL-d.antenna0.roach1.load2,2))
0004 hold all;
0005 plot(86400*(d.antenna0.roach2.utc-d.antenna0.roach2.utc(1)),mean(d.antenna0.roach2.LL-d.antenna0.roach2.load2,2))
0006 
0007 legend('Roach1','Roach2')
0008 
0009 
0010 
0011 %%%check alignment with the servo time
0012 d=read_arcSouth('15-Sep-2013:05:42:00','15-Sep-2013:05:48:00') %%time 1
0013 
0014 t1='15-Sep-2013:05:42:00',t2='15-Sep-2013:05:48:00';
0015 t1='15-Sep-2013:06:02:00';t2='15-Sep-2013:06:08:00';
0016 
0017 d=read_arcSouth(t1,t2) %%time 1
0018 
0019 figure
0020 plot(86400*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),mean(d.antenna0.roach1.LL-d.antenna0.roach1.load2,2))
0021 hold all;
0022 plot(86400*(d.antenna0.roach2.utc-d.antenna0.roach2.utc(1)),mean(d.antenna0.roach2.LL-d.antenna0.roach2.load2,2))
0023 title(strcat('Roachs---',t1))
0024 legend('Roach 1','Roach2')
0025 ylim([10000 30000])
0026 grid on
0027 print( gcf, '-depsc', 'Time1' )
0028 
0029 
0030 
0031 
0032 
0033 
0034 figure
0035 h(1)=subplot(211)
0036 plot(86400*(d.antenna0.servo.utcFast-d.antenna0.servo.utcFast(1)),d.antenna0.servo.apparent(:,1))
0037 title('Servo position')
0038 
0039 h(2)=subplot(212)
0040 plot(86400*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)),mean(d.antenna0.roach1.LL-d.antenna0.roach1.load2,2))
0041  linkaxes(h,'x');
0042 
0043  
0044  %%% check roach vs servo position
0045  plot(86400*(d.antenna0.servo.utcFast-d.antenna0.servo.utcFast(1)),'.')
0046  hold all
0047  plot(86400*(d.antenna0.roach1.utc-d.antenna0.roach1.utc(1)))
0048  
0049  figure
0050  subplot(211)
0051   plot(d.antenna0.servo.apparent(:,1),mean(d.antenna0.roach1.LL-d.antenna0.roach1.load2,2))
0052  hold all
0053  subplot(212)
0054   plot(d.antenna0.servo.apparent(:,1),mean(d.antenna0.roach2.LL-d.antenna0.roach2.load2,2))
0055  
0056   
0057   t1='15-Sep-2013:05:42:00',t2='15-Sep-2013:07:38:00';
0058 t1='15-Sep-2013:06:02:00';t2='15-Sep-2013:06:08:00';
0059 
0060 d=read_arcSouth(t1,t2) %%time 1
0061 
0062   
0063   
0064    %%%%%%%%%NOW CHECK THE FIRST ROACH TIMING OFFSET
0065    close all
0066   roachdata=mean(d.antenna0.roach1.LL-d.antenna0.roach1.load2,2);
0067   
0068   
0069   
0070   indexaz=[1000:50000];
0071   posScan=find(diff(d.antenna0.servo.apparent(indexaz,1))>0);
0072   negScan=find(diff(d.antenna0.servo.apparent(indexaz,1))<0);
0073   
0074   figure
0075    offset=0;
0076   plot(d.antenna0.servo.apparent(posScan,1),roachdata(posScan+offset),'.','Markersize',5);
0077   hold all
0078   plot(d.antenna0.servo.apparent(negScan,1),roachdata(negScan+offset),'.r','Markersize',5);
0079   title(num2str(offset));
0080   ylim([16000 25000])
0081   %xlim([148 153])
0082   grid on
0083   
0084   hold all
0085   offset=218; 
0086   subplot(312)
0087   h(2)=plot(d.antenna0.servo.apparent(indexaz,1),roachdata(indexaz+offset),'.','Markersize',5);
0088   legend('+2.18s')
0089   ylim([16000 25000])
0090   xlim([148 153])
0091   grid on
0092   
0093   offset=216; 
0094   subplot(313)
0095   h(3)=plot(d.antenna0.servo.apparent(indexaz,1),roachdata(indexaz+offset),'.','Markersize',5);
0096   legend('+2.16s')
0097   ylim([16000 25000])
0098     xlim([148 153])
0099 
0100   grid on
0101   
0102   linkaxes(h,'x');
0103     
0104     
0105     
0106     
0107   %%%%%%%%%NOW CHECK THE SECOND ROACH TIMING OFFSET
0108     roachdata=mean(d.antenna0.roach2.LL-d.antenna0.roach2.load2,2)
0109       indexaz=[1000:30000]
0110   offset=150;
0111   subplot(311)
0112     h(1)=plot(d.antenna0.servo.apparent(indexaz,1),roachdata(indexaz+offset))
0113     legend('+2.5s')
0114     hold all
0115      indexaz=[1000:30000]
0116   offset=140; 
0117   subplot(312)
0118     h(2)=plot(d.antenna0.servo.apparent(indexaz,1),roachdata(indexaz+offset))
0119     legend('+2.4s')
0120            indexaz=[1000:30000]
0121   offset=130; 
0122   subplot(313)
0123     h(3)=plot(d.antenna0.servo.apparent(indexaz,1),roachdata(indexaz+offset))
0124     legend('+2.3s')
0125   
0126     linkaxes(h,'x');
0127     
0128     plot(d.antenna0.servo.apparent(:,2))
0129