We perform a hierarchical Bayesian inference to investigate the population properties of the coalesc- ing compact binaries involving at least one neutron star (NS). With the current observation data, we can not rule out either of the Double Gaussian, Single Gaussian and Uniform NS mass distribution models, although the mass distribution of the Galactic NSs is slightly preferred by the gravitational wave (GW) observations. The mass distribution of black holes (BHs) in the neutron star-black hole (NSBH) population is found to be similar to that for the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be about 3 : 7. The spin properties of the binaries, though constrained relatively poor, play nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed.