I fitted the dominant spectral peak in each power column with a Lorentzian plus a linear background. Its fitted intensity scales with power as \(I \propto P^b\), with \(b = 0.9210 \pm 0.0120\). The fitted peak positions and intensities are in [lorentzian_peak_parameters.csv](<BENCH>/raw/commercial/luna/bench/t-pl/run-25/lorentzian_peak_parameters.csv); the requested plots are [intensity vs power](<BENCH>/raw/commercial/luna/bench/t-pl/run-25/peak_intensity_vs_power.png) and [peak position vs power](<BENCH>/raw/commercial/luna/bench/t-pl/run-25/peak_position_vs_power.png).

RESULT: b = 0.9210 ± 0.0120