It is expected that the source redshift can be obtained from the identified EM counterparts from the NS-NS and BH-NS binary systems. The redshift distribution of GW sources is
taken from the form of [1]}, [2]}

P(z)4d2C(z)R(z)H(z)(1+z) ,


where \(H(z)\)  represents the Hubble parameter from the fiducial cosmological model, \(d_{\rm C}\)  is comoving distance, and \(R(z)\)  is the merger rate of binary systems with following expression [3]}, [1]},
\(R(z)={\left\lbrace \begin{array}{ll}1+2z, & (z\le 1) \\\frac{3}{4}(5-z), & (1<z<5) \\0, & (z\ge 5).\end{array}\right.}\) 
