function J = J(tauM, w, n)

%function J = J(tauM, w, n, tauE, S, FULLMODELFREE)

% spectral density function 
% w in Hz, tau in ns

%tauAv = 1/(1/tauM+1/tauE);
%J = S^2*tauM*1e-9/(1+(n*w*2*pi*tauM*1e-9)^2) + (1-S^2)*tauAv*1e-9/(1+(n*w*2*pi*tauAv*1e-9)^2)*FULLMODELFREE ;

J = tauM*1e-9/(1+(n*w*2*pi*tauM*1e-9)^2)  ;