(79) where κ is the molecular line width and ⟨cos 2 θ⟩ (see Eq. 76) is the effective modification of polariton frequencies depending on orientation disorder. In Fig. , we present the numerical results obtained from our simulations (solid colors) when normalized to the peak height, as well as the analytic answer from Eq. 79 (black dashed lines). Overall, the analytic answer gives a reasonable estimate of the lineshape behavior across different angle disorders but slightly overestimates the linewidth. The narrower linewidth from the numerical results is likely due to the subaverage behavior (linewidth smaller than Γ/2 + κ/2 for the resonant case, see Eq. 34 in Ref. 50), which cannot be captured by the simple expression of the lineshape in Eq. 79 that predicts the linewidth to be Γ/2 + κ/2 under the resonance case.