MAP

synapse_power=IHCpreSynapseParams.power;

% tauCa vector is established across channels to allow vectorization

tauCa=repmat(ANtauCas, nBFs,1);

tauCa=reshape(tauCa, nANchannels, 1);

    %% synapse -----------------------------

    % Compute the vesicle release rate for each fiber type at each BF

    Vsynapse=repmat(IHC_RP, nANfiberTypes,1);

    % look-up table of target fraction channels open for a given IHC_RP

    mICaINF=    1./( 1 + exp(-gamma  * Vsynapse)  /beta);

    for idx=1:segmentLength

        % mICaINF is the current 'target' value of mICa

        mICaCurrent=mICaCurrent+(mICaINF(:,idx)-mICaCurrent)*dt./tauM;

        mICa(:,idx)=mICaCurrent;

    end

    ICa=   (GmaxCa* mICa.^3) .* (Vsynapse- ECa);

    for idx=1:segmentLength

        CaCurrent=CaCurrent +  ICa(:,idx)*dt - CaCurrent*dt./tauCa;

        synapticCa(:,idx)=CaCurrent;

    %     end

    switch AN_spikesOrProbability

        case 'probability'

            for t = 1:segmentLength;

                M_Pq=PAN_M-Pavailable;

                M_Pq(M_Pq<0)=0;

            %% Apply refractory effect

               % the probability of a spike's occurring in the preceding

               %    pFired= 1 - II(1-p(t)),

               % we need a running account of cumProb=II(1-p(t))

               %   cumProb(t)= cumProb(t-1)*(1-p(t))/(1-p(t-refracPeriod))

               % not the relative refractory period. To approximate this it

               % is necessary to extend the refractory period by 50%

                for t = segmentStartPTR:segmentEndPTR;

                    if t>1

                    ANprobRateOutput(:,t)= ANprobRateOutput(:,t)...

%                 figure(88), plot(cumANnotFireProb'), title('cumNotFire')

%                 figure(89), plot(ANprobRateOutput'), title('ANprobRateOutput')

            %  acoustic reflex

            [r c]=size(ANrate);

            if r>nBFs % Only if LSR fibers are computed