We observed several young stellar objects (YSOs) using the Photodetector Array Camera and Spectrometer (PACS) aboard the Herschel Space Observatory. CO, OH, $H_2O$, [O I], and [C II] lines were detected. CO rotational diagrams show two distinct gas components of ~400 K and ~1000 K with a break around 1500 K of the CO excitation energy, indicative of two different heating mechanisms: PDR and outflow shocks. OH and $H_2O$ line fluxes can be fitted with temperatures different from what are derived from the CO rotational diagrams. In order to understand the physical environment of line formation, the sources were modeled with the 3-D radiative transfer code, LIME. We present the results of observations, simple analysis, and modeling of Herschel-PACS spectra of the YSOs.