function Ntrue = geiger(Nobs,tau,Ntest) % Da [Heslop2012], pagg.146--148
% `Nobs` esprime il conto osservato dell'attività
% (emissione di raggi gamma);
% `tau` esprime il tempo morto per ciascun impulso
% (circa 20e-6 secondi con strumenti vecchi);
% `Ntest` esprime un array di valori da cui
% va scelta la predizione che è più vicina al
% conto osservato `Nobs`.
Ncalc = Ntest .* exp( -Ntest * tau ); % `cps` predetti
Ndiff = abs(Ncalc-Nobs); % differenza rispetto a quanto osservato
figure, hold on,
plot([Ntest(1),Ntest(length(Ntest))],[Nobs,Nobs],'k');
plot(Ntest,Ncalc,'b'), xlabel('Campionamenti'), ylabel('Ticchettio'),
plot([Ntest(1),Ntest(length(Ntest))],[Ncalc(1),Ncalc(length(Ncalc))],'g'),
plot(Ntest,Ndiff,'m'),
[~,idx] = sort(Ndiff); % localizzare la minima distanza
Ntrue=Ntest(idx(1)); % scelta del valore che calza meglio di tutti