We present spectroscopy of stars in the immediate vicinity of the dwarf nova (DN) KZ Gem to confirm its identification, which had been ambiguous in the literature. Analysis of 73 radial velocities spanning from 2014 to 2019 provides a high-precision orbital period of 0.2224628(2)\,d ($\sim5.34$\,hr) and shows KZ\,Gem to be a double-lined DN. Time series photometry taken from 2016 to 2018 shows a variable double-hump modulation with a full amplitude of $\sim0.3$\,mag, along with five Gaussian-like transient events lasting $\sim30$\,min or more. Using the light curve code XRBinary and nonlinear fitting code NMfit, we obtain an optimized binary model of the dwarf nova (DN) KZ Gem, from time series photometry, consisting of a Roche-lobe-filling K type dwarf with a mass transfer rate of $2.7\,-\,7.9\times10^{-10}\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ to a large, cool and thick disk surrounding a white dwarf, in an orbit with an inclination of $51^{\circ}.6(\pm1^{\circ}.4)$. Two hotspots on the disk are demonstrated to cause the observed variations in the ellipsoidal modulations from the secondary star. This physical model is compatible with the Gaia distance of KZ\,Gem.

We used the light curve code XRBinary to model the quiescent K2 light curves of three low-inclination cataclysmic variables (CVs): 1RXS\,J0632+2536 (J0632+2536), RZ\,Leo, TW\,Vir and the pre-CV WD\,1144+011. Optimized light curve models were obtained using a nonlinear fitting code NMfit and visualized by Phoebe 2.0. The disk model of J0632+2536 shows that one hotspot at the edge of the disk is enough to describe its light curve, while the other two dwarf nova (DN): RZ\,Leo and TW\,Vir require two hotspots. A typical pre-CV model with a weak irradiation effect for WD\,1144+011 can explain its single-hump modulation, and the newly observed spectrum confirms its previous classification. The synthetic analyses for the DN clearly indicate that phase zero of the double-hump modulations occurs around the secondary minimum and the primary hump is mainly caused by the hotspot at the edge of the disk. The quiescent disk has a flat temperature distribution with a power index of $\sim0.11$. The disk model of RZ\,Leo implies a truncated disk, supporting its previously speculated classification as an intermediate polar (IP). Except for the IP model of RZ\,Leo, which lacks a component related to the inferred accretion curtain, the models of J0632+2536, TW\,Vir and WD\,1144+011 are consistent with results from the Gaia mission. The derived masses and radii of the secondaries of the three DN are consistent with the semi-empirical relations for CV donor stars, while their effective temperatures are higher than the predictions. Irradiation of the donor stars is investigated to explain this discrepancy.