y@0: function [binNb, width]  = refineScore(mag,binNb)
y@0: 
y@0: toneFactor = 2^(1/12);
y@0: 
y@0: firstBin = floor((binNb+1)/toneFactor)-2;
y@0: lastBin = ceil((binNb+1)*toneFactor)+2;
y@0: firstBin = max(firstBin,1);
y@0: width =0;
y@0: if (firstBin<lastBin-1)
y@0:     %pitch refinment
y@0:     [maxB,indMax ]= max(mag(firstBin:lastBin));
y@0:     maxIdx = firstBin+indMax-1;
y@0:     if (maxB <0.9)
y@0:     binNb =0;
y@0:     elseif (maxIdx ~= lastBin && maxIdx ~= firstBin)
y@0:         a = log(mag(maxIdx-1));
y@0:         b = log(mag(maxIdx));
y@0:         c = log(mag(maxIdx+1));
y@0:         binNb = (maxIdx + (a - c)/(2*(a-2*b+c)))-1;
y@0:         minIdx = floor(binNb);
y@0:         maxIdx = ceil(binNb);
y@0:         while(minIdx>0&&mag(minIdx)>0.1*maxB)
y@0:             minIdx = minIdx-1;
y@0:         end
y@0:         while(mag(maxIdx)>0.1*maxB)
y@0:             maxIdx = maxIdx+1;
y@0:         end
y@0:         width = 2*min( maxIdx -binNb+1,binNb- minIdx+1);
y@0:         width = max(width,3);
y@0:     else
y@0: %         plot(mag(1:maxIdx+10))
y@0: %         close all;
y@0:     end
y@0: end
y@0: