We calculate the chiral string amplitude in pure spinor formalism and take four point amplitude as an example. The method could be easily generalized to $N$ point amplitude by complicated calculations. By doing the usual calculations of string theory first and using a special singular gauge limit, we produce the amplitude with the integral over Dirac $delta$-functions. The Bosonic part of the amplitude matches the CHY amplitude and the Fermionic part gives us the supersymmetric generalization of CHY amplitude. Finally, we also check the dependence on boundary condition for heterotic chiral string amplitudes.