diff parameterStore/MAPparamsOHC.m @ 26:b03ef38fe497

refractory (prob) effect added

--- a/parameterStore/MAPparamsOHC.m	Fri Jun 17 14:30:28 2011 +0100
+++ b/parameterStore/MAPparamsOHC.m	Tue Jun 21 14:58:12 2011 +0100
@@ -1,5 +1,5 @@
 function method=MAPparamsOHC ...
-    (BFlist, sampleRate, showParams)
+    (BFlist, sampleRate, showParams, paramChanges)
 % MAPparams<> establishes a complete set of MAP parameters
 % Parameter file names must be of the form <MAPparams> <name>
 %
@@ -58,11 +58,11 @@
 % i.e. a minimum ratio of 0.056.
 % 'spikes' model: AR based on brainstem spiking activity (LSR)
 OMEParams.rateToAttenuationFactor=0.006;   % * N(all ICspikes)
-%     OMEParams.rateToAttenuationFactor=0;   % * N(all ICspikes)
+% OMEParams.rateToAttenuationFactor=0;   % * N(all ICspikes)
 
 % 'probability model': Ar based on AN firing probabilities (LSR)
 OMEParams.rateToAttenuationFactorProb=0.01;% * N(all ANrates)
-%     OMEParams.rateToAttenuationFactorProb=0;% * N(all ANrates)
+% OMEParams.rateToAttenuationFactorProb=0;% * N(all ANrates)
 
 % asymptote should be around 100-200 ms
 OMEParams.ARtau=.05; % AR smoothing function
@@ -104,9 +104,9 @@
 DRNLParams.rateToAttenuationFactor = .01;  % strength of MOC
 %      DRNLParams.rateToAttenuationFactor = 0;  % strength of MOC
 % 'probability' model: MOC based on AN spiking activity (HSR)
-DRNLParams.rateToAttenuationFactorProb = .005;  % strength of MOC
+DRNLParams.rateToAttenuationFactorProb = .0055;  % strength of MOC
 % DRNLParams.rateToAttenuationFactorProb = .0;  % strength of MOC
-DRNLParams.MOCrateThreshold =70;                % spikes/s probability only
+DRNLParams.MOCrateThresholdProb =70;                % spikes/s probability only
 
 DRNLParams.MOCtau =.1;                         % smoothing for MOC
 
@@ -263,6 +263,16 @@
 end
 
 
+%% now accept last minute parameter changes required by the calling program
+% paramChanges
+if nargin>3 && ~isempty(paramChanges)
+    nChanges=length(paramChanges);
+    for idx=1:nChanges
+        eval(paramChanges{idx})
+    end
+end
+
+
 %% write all parameters to the command window
 % showParams is currently set at the top of htis function
 if showParams
@@ -276,6 +286,14 @@
             eval(['UTIL_showStructureSummary(' nm{i} ', ''' nm{i} ''', 10)'])
         end
     end
+
+    % highlight parameter changes made locally
+    if nargin>3 && ~isempty(paramChanges)
+        fprintf('\n Local parameter changes:\n')
+        for i=1:length(paramChanges)
+            disp(paramChanges{i})
+        end
+    end
 end