Blackman window magnitude compensation

I am trying to compute the fundamental phasor using sliding window DFT. I have employed a Blackman window in conjunction i.e. where is the time-domain signal sampled at a rate which is an integral multiple of the fundamental frequency i.e., where N is the number of samples per cycle of the fundamental, and represents the Blackman window.

I am referring a paper which is giving me an equation that should compensate for the magnitude attenuation due to this window whenever the time-domain signal frequency deviates from the nominal. The equation is given as follows,

The author also mentions that , where is the length of a simple moving average filter used in the pre-DFT computing process and also with being the computed time-domain signal frequency .

Using values such as

The list of values obtained for for a list of values of are given below:

1) 50 Hz=0.0639, 2) 51 Hz=0.0618, 3) 52 Hz=0.0594, 4) 53 Hz=0.0569, 5) 54 Hz=0.0543, 6) 55 Hz=0.0516.

I have a few questions regarding this, 1) Why is involved in computation. 2) Why are the values of magnitude gain of the blackman window not accurately compensating for the magnitude attenuation caused by the same, i.e. for 50 Hz the gain should have been unity exactly since the DFT is centered around it, and then close to unity as the frequency deviates from the nominal. Any help related to this will be really appreciated.