time   calls  line
                  1 function d = mainsWrapper(d, plotparams, parm, field);
                  2 
                  3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                  4 %
                  5 %  function d = mainsWrapper(d, plotparams, parm, field);
                  6 %
                  7 %    reduceData function that clears the data of the mains hum
                  8 %    This works on the switched data ie use before rfactor and alpha
                  9 %
                 10 %    sjcm
                 11 %    4-nov-2010 act : forced to work on switched data 
                 12 %
                 13 %    30-mar-2011 sjcm:  no plotting but can save plots
                 14 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                 15 
             1   16 useSwitch = 1; 
                 17 
             1   18 [d,q] = logcal(d, 'mains'); 
             1   19 if (q) 
                 20   disp('Mains correction has been applied');
                 21   disp(' ')
                 22 %  d=checkstatus(d, useSwitch);
                 23   return
                 24 end
                 25 
             1   26 if(useSwitch) 
             1   27   numChannels = 24; 
                 28 else
                 29   numChannels = 6;
                 30 end
                 31 
                 32 
                 33 
                 34 % check if we have the flagging parameters
             1   35 parFlag = checkpar(parm, 'mains'); 
             1   36 if(parFlag) 
             1   37   filterParams = reshape(parm.mains.tukey, [4 length(parm.mains.tukey)/4])'; 
                 38 else
                 39   display('Using Default mains filtering parameters about 40Hz (alias of 60Hz mains)');
                 40   filterParams = [38.9 39.9 40.1 41.1];
                 41 end
                 42 
                 43 % THIS METHOD TAKES TOO LONG.  LET'S PUT A FILTER ABOUT 60Hz.
                 44 % now we remove the mains
                 45 %[mains_factors, indices] = calculate_mains(d, useSwitch);
                 46 %d.correction.mains.factors = mains_factors;
                 47 %d.correction.mains.indices = indices;
                 48 
                 49 % apply the correction.
             1   50 dcorr = d; 
             1   51 for m=1:size(filterParams,1) 
             3   52   for nn=1:numChannels 
 45.56      72   53     dcorr = filter_data(dcorr, useSwitch, nn, filterParams(m,1), ... 
                 54 	filterParams(m,2), filterParams(m,3), filterParams(m,4), 0);
            72   55   end 
             3   56 end 
                 57 
                 58 
  0.03       1   59 setPlotDisplay(plotparams.plot); 
                 60 
                 61 % next we plot to see if there was any effect at all
                 62 
                 63 % only plot data that is on the blank stare
  0.01       1   64 if(isempty(find(dcorr.index.blank.fast))) 
                 65   display('No blank data to plot effect of filter');
             1   66 else 
             1   67   display('Plotting Power Spectrum'); 
             1   68   display('No flagging will be allowed'); 
             1   69   if(useSwitch) 
  1.91       1   70     d = packd(dcorr, [], 'none', 'mainsSwitch', plotparams, 'Main Plots', field, d); 
                 71       
                 72   else
                 73     d = packd(dcorr, [], 'none', 'mainsReg', plotparams, 'Main Plots', field, d);  
                 74   end
             1   75 end 
                 76 
                 77 
                 78 
                 79 % log that it has been applied
             1   80 d = dcorr; 
                 81 
             1   82 d = logcal(d, 'mains'); 
                 83 
                 84 
             1   85 return; 

Other subfunctions in this file are not included in this listing.