Holomorphic instantons in target torus

For computing instantons contributions from worldsheet torus to target torus, one can evaluate zero modes contribution of genus 1 partition function given by following expression:
$$Tr(-1)^FF_LF_Rq^{L_0}bar{q}^{bar{L_0}}= frac{t+bar t}{4pitau_2}sum_{m.n,r,s}expleft[-frac{t}{4tau_2rho_2}|(m+rrho)-bar{tau}(n+srho)|^2-frac{bar t}{4tau_2 rho_2}|(m+rrho)-{tau}(n+srho)|^2right]…….(1)$$
This is nothing but Narain lattice summed part of bosonic momenta contributing to the partition function $Tr(-1)^FF_LF_Rq^{L_0}bar{q}^{bar{L_0}}$ and can be obtained in a similar way one calculates threshold correction in heterotic string theory as done in Dixon et al."Moduli dependence of string loop corrections to gauge coupling constants."

Now the task at hand is to extract holomorphic and zero instanton contributions from the above expression which actually has following expression as per (BCOV1 eqn. 9):
$$ I_F/tau_2=Tr(-1)^FF_LF_Rq^{L_0}bar{q}^{bar{L_0}}simeqfrac{t+bar{t}}{4pitau_2}+sum_{Min
GL(2,Z)}frac{tau_2}{|detM|}e^{-|detM|t}delta(tau- M(rho))+ Order(e^{-ibar t})………..(2)$$
where $rho$ is complex structure of target torus and $tau$ is that of worldsheet torus. t is taken as Kaehler parameter here.

First term in the above expression is easy to see with M=0 but the second term is hard to derive though all the terms involved have very natural interpretation as $e^{-|detM|t}$ is instanton supression of path integral by holomorphic maps of degree |detM| and |detM| in the denominator corresponds to a multicovering contribution in target torus from worldsheet torus. However it does not seem to be possible to derive the second part just by taking $bar trightarrow infty$ limit of (1).

I looked it up in Claymath "Mirror Symmetry" where through a redefintion of Kaehler parmater in (1) can be rewritten as (eqn. (35.15) in Claymath),
$$I_F=frac{t-bar t}{2i}sum_M expleft[-2pi it|detM|-frac{pi A}{tau_2rho_2}left|(1 ,rho)Mbegin{pmatrix}
tau
1
end{pmatrix}right|^2right]$$
where M is
$$begin{pmatrix}
m_2 & k_2
m_1 & k_1
end{pmatrix}$$ and $A=(t-bar t)/2i.$

Now after taking the holomorphic limit of I_F, one would arrive at zero instanton part which is nothing but A and other part which is supposed to be,
$$I_{Fin M}= sum_{Min
GL(2,Z)}frac{tau_2^2}{|detM|}e^{2pi it|detM|}delta(tau- M(rho))….(3)$$ apart from an infinite constant but I am unable to derive this.

What I am doing is to send $m_i$‘s to $-m_i$‘s giving $e^{2pi it|detM|}$ part and other part becomes $$Asum_M expleft[-frac{pi A}{tau_2rho_2}left|(k_2+rho k_1)-tau(m_2+rho m_1)right|^2right]$$

Now this Gaussian in $bar{t}rightarrow-iinfty$ can be approximated by a delta function with $M=frac{k_2+rho k_1}{m_2+rho m_1}$ giving overall with factor $e^{2pi it|detM|}$,
$$Asum_Msqrt{frac{tau_2rho_2}{A|m_2+rho m_1|^2}}e^{2pi it|detM|}delta(tau- M(rho))$$ but this is not same as (3)!!

Anyone has any idea??