Influence of x-ray pulse generated from gamma spectrometers should be eliminated in applications, which typically uses pulse shape techniques between gamma and x-ray pulses. In this study, we proposed and tested several algorithms aiming to eliminate this influence. The algorithms are based on curve fitting (CF), artificial neural network (ANN), system identification, peak shape, amplitude search with curve fitting and pulse tracking methods. Gamma pulses and X-ray pulses are detected by NaI(TI) scintillator detector and Silicon lithium Si(Li) detector, respectively. The developed algorithms are tested using 32,000 total instantaneous detector events of acquired gamma pulses and 65,536 total instantaneous detector events of x-ray source. An algorithm using the least square curve fitting method is applied for differentiation between gamma and x-ray pulses. ANN is employed as a classifier for identification of extracted spectrum and Bispectrum features of gamma and x-ray pulses. A comparison between identification results due to extracted spectrum and Bispectrum features is established. System identification algorithm is then built to determine the detection system response of each radiation pulse, which includes various models to attain best fitting. These models are Auto-regressive model with external input (ARX), the linear parametric model (IV) and process models (P1D). The peak shape algorithm is also tried, which depends on the individual classification of pulse width. The amplitude search with curve fitting algorithm is implemented. Moreover, the pulse tracking algorithm is investigated for PSD between gamma and x-ray pulses. The maximum peak of contaminated pulse is tracked using a suggested peak search method. Then, pulse position is estimated using matrix method. Comparison between these algorithms is conducted based on the evaluation of light of residuals, fitting error and processing time. The results confirm that peak shape algorithm is the best one from computational speed point of view, while ANN algorithm using Bispectrum feature extraction method is the most appropriate one that yields 100% accuracy over noisy environment with longer processing time. In addition, the system identification algorithm is the optimal algorithm that achieves zero fitting error under clean environment. These proposed algorithms for PSD between gamma and x-ray pulses lead to design efficient spectrometers with optimal applicability in various environments.