Few samples of model C¹⁸O (J = 2–1) spectra toward G333.60 region obtained from Gausspy+ module. First, third, and fifth rows show the observed and model C¹⁸O (J = 2–1) spectra in G333.60 region. Black solid line denotes the observed spectra, whereas, red sold lines represent the decomposed and resultant model spectra in G333.60 region. In each of the figures two red dotted horizontal lines denote the ±1σ value of the noise. In addition one solid horizontal red line represents the +3σ noise level of the spectra. In the header of each of the figures idx represents the spectra number that is obtained from the Gausspy+ module. X and Y denote the pixel number in the image cube, N_comp is the number of decomposed components and ${\chi }_{\text{red\hspace{0.17em}}}^2$ is the reduced chi-squared value after fitting. Second, fourth, and sixth rows represent the residuals of the corresponding spectra after being fitted with model components. Here also two red dotted horizontal lines in each figure denotes the ±1σ value of the noise.