mne_get_current_comp.m
View Options

			function [comp] = mne_get_current_comp(info)
			%
			% [comp] = mne_get_current_comp(info)
			%
			% Get the current compensation in effect in the data
			%


			%
			% Author : Matti Hamalainen, MGH Martinos Center
			% License : BSD 3-clause
			%
			%
			% Revision 1.2 2006/04/23 15:29:40 msh
			% Added MGH to the copyright
			%
			% Revision 1.1 2006/04/12 10:51:19 msh
			% Added projection writing and compensation routines
			%
			%

			me='MNE:mne_get_current_comp';

			global FIFF;
			if isempty(FIFF)
			FIFF = fiff_define_constants();
			end

			comp = 0;
			first_comp = -1;
			for k = 1:info.nchan
			if info.chs(k).kind == FIFF.FIFFV_MEG_CH
			comp = bitshift(double(info.chs(k).coil_type),-16);
			if first_comp < 0
			first_comp = comp;
			elseif comp ~= first_comp
			error(me,'Compensation is not set equally on all MEG channels');
			end
			end
			end

			return;

Slepi/Code/Utilities/spm12/external/mne/mne_get_current_comp.m8a186cd82567master