function [output] = hsfilter(theta,fs,input)
%This function processes a mono input signal applying the transfer function
%for head shadowing and delaying the output for compensating the interaural
%time differences. It outputs a mono signal so it has to be applied twice
%to a same input signal but using opposite values for theta to create
%output for the left and right channels of the stereo signal.
theta = theta + 90; %angle adjustment so that range is: -90