In this paper, the second in a series, we continue to study the generalized Lamé equation with the Treibich-Verdier potentialfrom the monodromy aspect. We prove the existence of a pre-modular form Z n r,s (τ) of weight 1 2 ∑ n k (n k + 1) such that the monodromy data (r, s) is characterized by Z n r,s (τ) = 0. This generalizes the result in [16], where the Lamé case (i.e. n 1 = n 2 = n 3 = 0) was studied by Wang and the third author. As applications, we prove among other things that the following two mean field equations ∆u + e u = 16πδ 0 and ∆u + e u = 8π 3 ∑ k=1 δ ω k 2 on a flat torus E τ := C/(Z + Zτ) has the same number of even solutions. This result is quite surprising from the PDE point of view.