make_plots_I

 Read in the data:

0001 % Read in the data:
0002 function make_plots_I(filestring)
0003 
0004 dirname = '/home/heikoh/cbass/New_Codes/Full_Sky/Working_With_Stephen/Matlab';
0005 txt = sprintf('ls %s/%s*.mat', dirname, filestring);
0006   [s w] =unix(txt); % reads in the file list
0007 
0008   index = 1;
0009   while(~isempty(w)) % goes through file list and assigns each file an index
0010     [filename{index}, w] = strtok(w);
0011     index = index+1;
0012   end
0013   
0014   dall = [];
0015   for m=1:length(filename)-1
0016       eval(sprintf(' load %s drfi', filename{m}));% loads each file
0017       if(exist('drfi', 'var'))
0018       if(isempty(dall))
0019           dall = drfi; % for first value, assigns dall = drfi
0020       else
0021           dall = [dall, drfi]; % concantenates each drfi
0022       end
0023       end
0024 
0025   end
0026   
0027 %size(dall);
0028 %size(dall{1});
0029 %s = load(filenames);
0030 %size(s.drfi);
0031 %{
0032 for m=1:size(s.drfi,2)
0033     
0034     LL1 = s.drfi{m}.antenna0.roach1.LLfreq;
0035     LL2 = s.drfi{m}.antenna0.roach2.LLfreq;
0036     RR1 = s.drfi{m}.antenna0.roach1.RRfreq;
0037     RR2 = s.drfi{m}.antenna0.roach2.RRfreq;
0038     Intensity{m} = (LL1 + LL2 + RR1 + RR2)/4;
0039     
0040     LL1 = s.drfi{m}.antenna0.roach1.LL;
0041     LL2 = s.drfi{m}.antenna0.roach2.LL;
0042     RR1 = s.drfi{m}.antenna0.roach1.RR;
0043     RR2 = s.drfi{m}.antenna0.roach2.RR;
0044     Intensity{m} = (LL1 + LL2 + RR1 + RR2)/4
0045     
0046 end
0047  %}
0048 % check to see sizes are correct:
0049 %{
0050 size(s.drfi{2}.antenna0.roach1.LL)
0051 size(Intensity{2})
0052 size(Intensity)
0053 size(s.drfi)
0054 %}
0055 plotfigs_I(dall)
0056 
0057 end % read_data_HMH
0058 
0059 
0060 %%%  plotting figures
0061 function plotfigs_I(drfi)
0062 clc
0063 close all
0064 %cbass_startup
0065 
0066 % first, let's just find the peak and see what its spectrum looks like
0067 figure()
0068 allspectra = [];
0069 Intensity_freq1 = [];
0070 Intensity1 = [];
0071 Intensity_freq2 = [];
0072 Intensity2 = [];
0073 
0074 for m=1:size(drfi,2)
0075 
0076     LL1freq = drfi{m}.antenna0.roach1.LLfreq;
0077     LL2freq = drfi{m}.antenna0.roach2.LLfreq;
0078     RR1freq = drfi{m}.antenna0.roach1.RRfreq;
0079     RR2freq = drfi{m}.antenna0.roach2.RRfreq;
0080     Intensity_freq1{m} = (LL1freq + RR1freq)/2;
0081     Intensity_freq2{m} = (LL2freq + RR2freq)/2;
0082     
0083     LL1 = drfi{m}.antenna0.roach1.LL;
0084     LL2 = drfi{m}.antenna0.roach2.LL;
0085     RR1 = drfi{m}.antenna0.roach1.RR;
0086     RR2 = drfi{m}.antenna0.roach2.RR;
0087     Intensity1{m} = (LL1 + RR1)/2;
0088     Intensity2{m} = (LL2 + RR2)/2;
0089     
0090     aa1 = find(Intensity_freq1{m} == max(Intensity_freq1{m}));
0091     aa2 = find(Intensity_freq2{m} == max(Intensity_freq2{m}));
0092     
0093     
0094     % to visualize them (First step)
0095     subplot(2,1,1); plot(Intensity1{m}(aa1,:));hold on
0096     title('Peak Intensity Values Across All Channels: Roach1')
0097     xlabel('Channels')
0098     ylabel('Intensity')
0099     subplot(2,1,2); plot(Intensity2{m}(aa2,:)); hold on
0100     title('Peak Intensity Values Across All Channels: Roach2')
0101     xlabel('Channels')
0102     ylabel('Intensity')
0103     allspectra(m,:) = [Intensity1{m}(aa1,:), Intensity2{m}(aa2,:)];
0104     %allspectra(m,:) = [drfi{m}.Intensity(aa,:)];
0105     %hold on
0106 end
0107 hold off
0108 %{
0109 keyboard;
0110 
0111 figure()
0112 for m=1:length(drfi)
0113   disp('New Event');
0114   f = find(drfi{m}.antenna0.roach1.LLfreq == max(drfi{m}.antenna0.roach1.LLfreq));
0115   minInt = f-200;
0116   maxInt = f+200;
0117   if(minInt < 1)
0118     minInt = 1;
0119   end
0120   if(maxInt > length(drfi{m}.antenna0.roach1.LLfreq))
0121     maxInt = length(drfi{m}.antenna0.roach1.LLfreq);
0122   end
0123   
0124   maxval = max(drfi{m}.antenna0.roach1.LL(:));
0125   cc = [drfi{m}.antenna0.roach1.LL(f,:) , drfi{m}.antenna0.roach2.LL(f,:)];
0126   for mm=minInt:maxInt
0127      bb = [drfi{m}.antenna0.roach1.LL(mm,:) , drfi{m}.antenna0.roach2.LL(mm,:)];
0128      plot(bb);
0129      ylim([0,maxval]);
0130      hold on
0131      plot(cc, 'r-');
0132      txt = sprintf('Event # %d', m);
0133      title(txt);
0134      pause(0.025);
0135      hold off     
0136   end
0137   pause(1)
0138   
0139 end
0140 
0141 keyboard;
0142 %}
0143 figure()
0144 plot(allspectra');
0145 title('Each RFI Event Across All Channels')
0146 ylabel('Intensity')
0147 xlabel('Channels')
0148 
0149 
0150 %{
0151 
0152 % later, we can take the "rfi event", and compare them to each other (once
0153 % normalized)
0154 figure()
0155 for m=1:size(drfi,2)
0156     bb = drfi{m}.antenna0.roach1.LL;
0157     bb = bb./repmat(bb(1,:), [size(bb,1), 1]);
0158     bb = bb-1;
0159     plot(sum(bb,1));
0160     hold on;
0161 end
0162 title('Search For Peak RFI Across All Channels')
0163 xlabel('Channels')
0164 ylabel('LL Intensity')
0165 hold off
0166 %}
0167 end