function dprime = eurich2022_decision(processed,mpar)
%eurich2022_decision  Decision stage of the Eurich et al. (2022) model 
%   Usage: dprime = eurich2022_decision(processed,mpar)
%
%   Input parameters:
%     processed : Matrix containing the processed stimulus from eurich2022_processing
%
%     mpar      : Structure with the following model parameters:
%
%                 - fs:                  sampling frequency of model
%                 - Filters_per_ERB_aud: spacing of peripheral filter central frequencies in ERB
%                 - flow:                lower bound of gammatone filterbank
%                 - fhigh:               higher bound of gammatone filterbank
%                 - gtorder:             filter order of gammatone filters (hohmann2002)
%                 - IPDnoise:            Internal noise affecting the IPD extraction
%                 - Dnoise:              Internal noise affecting the Decision stage
%                 - GaussSigma:          std of Gaussian across-channel smoothing
%                 - iKernelThresh:       treshold above which a value of the Gaussian filter window is used (below := 0)
%                 - n_hcbins:            number of bins to be used in IPD distribution computation via histcounts
%   
%   
%   Output parameters:
%     dprime                   : Sensitivity index :math:d' = frac{leftvert mu_a - mu_b rightvert}{sigma}
%
%
%   See also: eurich2022 exp_eurich2022
%
%   References:
%     B. Eurich, J. Encke, S. D. Ewert, and M. Dietz. Lower interaural
%     coherence in off-signal bands impairs binaural detection. The Journal
%     of the Acoustical Society of America, 151(6):3927--3936, 06 2022.
%     [1]arXiv | [2]http ]
%     
%     References
%     
%     1. http://arxiv.org/abs/https://pubs.aip.org/asa/jasa/article-pdf/151/6/3927/16528275/3927\_1\_online.pdf
%     2. https://doi.org/10.1121/10.0011673
%     
%
%   Url: http://amtoolbox.org/amt-1.4.0/doc/modelstages/eurich2022_decision.php


%   #Author: Bernhard Eurich (2022): original implementation
%   #Author: Piotr Majdak (2023): adaptation to AMT 1.4

% This file is licensed unter the GNU General Public License (GPL) either 
% version 3 of the license, or any later version as published by the Free Software 
% Foundation. Details of the GPLv3 can be found in the AMT directory "licences" and 
% at <https://www.gnu.org/licenses/gpl-3.0.html>. 
% You can redistribute this file and/or modify it under the terms of the GPLv3. 
% This file is distributed without any warranty; without even the implied warranty 
% of merchantability or fitness for a particular purpose. 

binaural_feature_ref  = processed(1,:);
binaural_feature_test = processed(3,:);

monaural_feature_ref  = processed(2,:);
monaural_feature_test = processed(4,:);

% binaural decision
dprime_binaural_multichannel = abs(binaural_feature_ref - binaural_feature_test) / mpar.binaural_internal_noise_sigma;

% monaural decision
delta_P = abs(monaural_feature_test - monaural_feature_ref);
dprime_monaural_multichannel = delta_P ./ mean([monaural_feature_test; monaural_feature_ref],1) / mpar.monaural_internal_noise_sigma;


if find(mpar.idx_of_target_channel)
    dprime_binaural = dprime_binaural_multichannel(mpar.idx_of_target_channel);
    dprime_monaural  = dprime_monaural_multichannel(mpar.idx_of_target_channel);
else
    dprime_binaural = max(dprime_binaural_multichannel);
    dprime_monaural = max(dprime_monaural_multichannel);

   % dprime_binaural = sqrt(sum(dprime_binaural_multichannel.^2));
    %dprime_monaural = sqrt(sum(dprime_monaural_multichannel.^2));
end


dprime = sqrt(dprime_binaural^2 + dprime_monaural^2);


end