A computer-implemented method for fitting scintillation pulse digitized signals, comprising the following steps: Step S1: Selecting a double-exponential model v(t) = a. e bt + c. e dt as a scintillation pulse prior model and determining said double-exponential model as a target function; Step S2: Selecting sampling points for classical pulses in a scintillation pulse database, said sampling points including eight voltage thresholds in sequence and eight time points in sequence, wherein said eight voltage thresholds in sequence are v1, v2, v3, v4, v5, v6, v7 and v8, and said eight time points in sequence are t1, t2, t3, t4, t5, t6, t7 and t8; said scintillation pulse database is obtained based on a MVT method; Step S3: storing said voltage thresholds in sequence as a voltage array v=[v1, v2, v3, v4, v5, v6, v7, v8], and carrying out time translation of a curve of said double-exponential model, said time translation comprising the following steps: setting a first time point t1 as 0, subtracting subsequent time points respectively by t1, and then storing obtained data as a time array t=[0, (t2-t1), (t3-t1), (t4-t1), (t5-t1), (t6-t1). (t7-t1), (t8-t1)]; Step S4: setting four initial default parameters a, b, c and d of said target function to be 1, and storing them as an initial parameter array para=[1, 1, 1, 1]; Step S5: inputting said voltage array v, said time array t, number n of sampling points for classical pulses, said initial parameter array para, and said target function V(t) into a Levenberg-Marquardt fitting function, and fitting to obtain four new fitting parameters a', b', c', d', and storing said four new fitting parameters as a second parameter array paranew=[a', b', c', d']; Step S6: as for each scintillation pulse in said scintillation pulse database, except replacing said initial parameter array para by said second parameter array paranew in Step S4, repeating said step S1- step S5 in sequence until all of scintillation pulses have been fitted, and recording a fitted function of a k-th scintillation pulse as V k (t) = a k · e b k t + c k · e d k t; Step S7: integrating said target function V k (t) of said scintillation pulse with said fitting parameters having been determined after fitting to get an energy value E k of a k-th pulse; Step S8: obtaining a zero-crossing point of a curve by solving an equation V k (t) = 0, obtaining an arrival time of said k-th pulse t k = t1 + t k0, and finishing extracting time and energy information of said pulse. Step S1: Selecting a double-exponential model v(t) = a. e bt + c. e dt as a scintillation pulse prior model and determining said double-exponential model as a target function; Step S2: Selecting sampling points for classical pulses in a scintillation pulse database, said sampling points including eight voltage thresholds in sequence and eight time points in sequence, wherein said eight voltage thresholds in sequence are v1, v2, v3, v4, v5, v6, v7 and v8, and said eight time points in sequence are t1, t2, t3, t4, t5, t6, t7 and t8; said scintillation pulse database is obtained based on a MVT method; Step S3: storing said voltage thresholds in sequence as a voltage array v=[v1, v2, v3, v4, v5, v6, v7, v8], and carrying out time translation of a curve of said double-exponential model, said time translation comprising the following steps: setting a first time point t1 as 0, subtracting subsequent time points respectively by t1, and then storing obtained data as a time array t=[0, (t2-t1), (t3-t1), (t4-t1), (t5-t1), (t6-t1). (t7-t1), (t8-t1)]; Step S4: setting four initial default parameters a, b, c and d of said target function to be 1, and storing them as an initial parameter array para=[1, 1, 1, 1]; Step S5: inputting said voltage array v, said time array t, number n of sampling points for classical pulses, said initial parameter array para, and said target function V(t) into a Levenberg-Marquardt fitting function, and fitting to obtain four new fitting parameters a', b', c', d', and storing said four new fitting parameters as a second parameter array paranew=[a', b', c', d']; Step S6: as for each scintillation pulse in said scintillation pulse database, except replacing said initial parameter array para by said second parameter array paranew in Step S4, repeating said step S1- step S5 in sequence until all of scintillation pulses have been fitted, and recording a fitted function of a k-th scintillation pulse as V k (t) = a k · e b k t + c k · e d k t; Step S7: integrating said target function V k (t) of said scintillation pulse with said fitting parameters having been determined after fitting to get an energy value E k of a k-th pulse; Step S8: obtaining a zero-crossing point of a curve by solving an equation V k (t) = 0, obtaining an arrival time of said k-th pulse t k = t1 + t k0, and finishing extracting time and energy information of said pulse.