I treated each power column as one merged spectrum and fitted the resolved peak near 1.336 eV with a Lorentzian plus a linear baseline. The peak intensity is the fitted Lorentzian amplitude. The allometric fits use \(y=aP^b\); the reported uncertainties are standard errors from the fit covariance.

- Peak intensity: \(b = 1.208 \pm 0.032\)
- Peak position \(x_c\): \(b = 0.00356 \pm 0.00024\)

The high-energy rise reaches the scan boundary, so it was not treated as a resolved peak.

Files: [plots](<<BENCH>/raw/commercial/luna/bench/t-pl/run-04/power_scaling.png>) · [Lorentzian fit results](<<BENCH>/raw/commercial/luna/bench/t-pl/run-04/lorentzian_peak_fits.csv>) · [spectra near the fitted peak](<<BENCH>/raw/commercial/luna/bench/t-pl/run-04/spectra_peak_region.png>)

RESULT: b = 1.208 ± 0.032