I fitted the dominant peak in each spectrum with a Lorentzian plus a local linear background. The reported intensity is the fitted peak height. The allometric fits use \(y=aP^b\); the uncertainties below are standard errors on \(b\) from the log–log fits.

| Power (µW) | Peak intensity (counts) | Peak position xc (eV) |
|---:|---:|---:|
| 5 | 193.3 | 1.31326 |
| 10 | 241.2 | 1.32258 |
| 15 | 361.3 | 1.32409 |
| 20 | 444.0 | 1.32621 |
| 25 | 557.8 | 1.32719 |
| 30 | 676.9 | 1.32815 |
| 35 | 751.4 | 1.32883 |
| 40 | 881.2 | 1.33026 |
| 45 | 963.7 | 1.33035 |
| 50 | 1,063.8 | 1.33150 |
| 60 | 1,295.0 | 1.33194 |
| 70 | 1,520.0 | 1.33343 |
| 80 | 1,735.2 | 1.33345 |
| 90 | 1,994.2 | 1.33415 |
| 100 | 2,297.0 | 1.33472 |

- [Peak intensity vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-02/peak_intensity_vs_power.png>): \(b=0.8736\pm0.0330\)
- [Peak position vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-02/peak_position_vs_power.png>): \(b=0.00487\pm0.00027\)
- [Lorentzian fit results for each power](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-02/lorentzian_peak_fits.csv>)

RESULT: b = 0.8736 ± 0.0330