Bayesian Analysis of Femtosecond Pump-Probe Photoelectron-Photoion Coincidence Spectra with Fluctuating Laser Intensities.

This paper employs Bayesian probability theory for analyzing data generated in femtosecond pump-probe photoelectron-photoion coincidence (PEPICO) experiments. These experiments allow investigating ultrafast dynamical processes in photoexcited molecules. Bayesian probability theory is consistently applied to data analysis problems occurring in these types of experiments such as background subtraction and false coincidences. We previously demonstrated that the Bayesian formalism has many advantages, amongst which are compensation of false coincidences, no overestimation of pump-only contributions, significantly increased signal-to-noise ratio, and applicability to any experimental situation and noise statistics. Most importantly, by accounting for false coincidences, our approach allows running experiments at higher ionization rates, resulting in an appreciable reduction of data acquisition times. In addition to our previous paper, we include fluctuating laser intensities, of which the straightforward implementation highlights yet another advantage of the Bayesian formalism. Our method is thoroughly scrutinized by challenging mock data, where we find a minor impact of laser fluctuations on false coincidences, yet a noteworthy influence on background subtraction. We apply our algorithm to data obtained in experiments and discuss the impact of laser fluctuations on the data analysis.

1. Introduction

Coincidence measurements are a widely-used and powerful experimental technique in physics and chemistry. Photoelectron-photoion coincidence (PEPICO) spectroscopy utilizes not only information obtained from the detection of electrons and ions, but also the fact that they stem from the very same ionization event [1,2,3,4,5,6]. Frequently used in photoionization studies of gas phase molecules or clusters, this technique allows for conclusions about the ionization process such as the disentanglement of competing intramolecular relaxation channels [4,7,8,9] or multiple species [10], a depth of insight that cannot be achieved without assigning electrons to the ions from which they originate. Thus, the success of PEPICO is based on these recordings of pairs being unambiguous, energy-resolved for electrons, and mass-resolved for the cations. Yet, the correct pairwise assignment (true coincidence) may be affected by certain experimental conditions: If a laser pulse triggers a number of simultaneous ionization events arising from different neutral molecules, the assignment of correlated electron-ion pairs is impaired and causes so-called false coincidences [11]. The three possible events are: (1) Not exactly one electron and one ion are detected; in this case the event is rejected. (2) One electron and one ion are detected, which can originate from the same molecule (true coincidence) or (3) an electron and ion can originate from different molecules (false coincidence). Let this be illustrated using the example of exactly two ionization events. If two electrons and/or two ions are detected, the measurement would simply be discarded since no unambiguous assignments could be made. Yet, with imperfect detectors, there is a non-negligible probability of the following event: The electron from Molecule 1 is detected, the electron from Molecule 2 is not detected; Cation 1 is not detected; Cation 2 is detected. Hence, the experimentalist sees a false coincidence, where Electron 1 is wrongly assigned to Cation 2. Obviously, false coincidences only arise if both detectors are not perfect, i.e., the detection probabilities are less than unity, and are thus to some extent present in any such experiment. An easy way out is to work with low ionization rates, for the price of either a bad signal-to-noise ratio or time-consuming measurements. In principle, detector noise or ionization events not caused by the laser pulse might also lead to false coincidences, but, for the situation at hand, are sufficiently low to be neglected.

Time-resolved studies are typically carried out as pump-probe experiments [10,12] as depicted in Figure 1. Excitation by a laser pulse, commonly referred to as the pump pulse, triggers dynamical processes in the molecule after which a time-delayed second laser pulse, commonly referred to as the probe pulse, ionizes the molecule. The transient change of photo-electron and -ion signals associated with the excited states, as a function of the time-delay, provides insight into the underlying processes. Unfortunately, pump and/or probe pulses on their own can ionize the molecule as well, leading to signals that are referred to as pump-only and probe-only further on. This background signal is superimposed on the excited state signal, and in many experimental situations, the pump-only and/or the probe-only signals significantly contribute to the pump-probe signal, e.g., if multiphoton transitions are applied for pumping or probing, or if high photon energies are used for probing [13]. In order to extract the excited state transients, the pump-only and/or the probe-only signals are measured separately and usually subtracted from the pump-probe signals, obviously resulting in increased noise if the background and pump-probe spectra energetically overlap with each other. The application of the Bayesian formalism to background subtraction alone was already presented for astrophysical applications [14,15] and for photo-induced X-ray emission spectroscopy (PIXE) [16,17,18,19,20,21]. Yet, here, it has to be considered that the pump-only, the probe-only, and the pump-probe measurements have different ionization rates. Since the statistics of the coincidences depend on the ionization rates, the statistics of sole pump-only and probe-only measurements differ from those in the pump-probe measurement, and simple subtraction turns out not to be an unbiased estimator any longer [22].

Figure 1

Utterly simplified sketch of a time-resolved photoionization study carried out with a pump-probe setup and a time-of-flight spectrometer. A commercial Ti:sapphire laser system delivers pulses of nm in center wavelength and fs in temporal length at a repetition rate of kHz. The delay stage is used to control the length of the optical path, and hence the time delay. The energy level diagram shows how the electron kinetic energy, given the energy of the states and the photons, identifies the state the system was in at the moment of ionization. A detailed description of the setup can be found in our previous publications [7,28].

Although it is feasible to distinguish between true and false coincidences by pure experimental finesse, such as in cold target recoil ion momentum spectroscopy (COLTRIMS) [23,24], it demands quite a technical and financial effort and is entirely impossible for time-of-flight detection, as used in the presented experiment. Covariance mapping, which is based on the calculation of the covariance for the photoelectron and mass spectra measured with each laser shot [25,26,27], does not guarantee that the reconstructed spectrum is positive, is restricted to Poisson processes, and leads to systematic deviations in other scenarios [26,27]. Further limitations are outlined in [25].

We recently presented a Bayesian approach to PEPICO, which treats both coincidences and background subtraction on the same footing [22,29]. In this work, we extend our theory to the prevalent experimental situation of unstable laser intensity, i.e., ionization rates fluctuating from pulse to pulse. We provide our software, including introductory examples, at https://github.com/fslab-tugraz/PEPICOBayes/. The experiment described in Section 2 is treated with the Bayesian formalism developed in Section 3. It will be tested by some challenging mock data in Section 4 and applied to real experimental data in Section 5.

2. Experiment

The analyzed experiment is of the type depicted in Figure 1 and described in detail in previous publications [7,28]. To apply our method, we choose our specimen and excitation-ionization-scheme according to Figure 2, since we expect the effects described above to be of particular importance in this scenario. Acetone molecules are excited by a three-photon transition to high-lying Rydberg states and ionized in the extraction region of a time-of-flight spectrometer, which measures both the electron kinetic energies and the ion masses [4,7]. The electron and ion flight times are then analyzed by a coincidence algorithm to produce photoelectron spectra corresponding to either an intact parent acetone ion or a fragmented acetyl ion. The concept and the necessity of time-resolved PEPICO (TR-PEPICO) in this context are further elucidated in Animations A1 and A2 in the Supplementary Material.

Figure 2

Pump-probe ionization scheme to investigate excited state dynamics in molecules.

The addressed excited state lies energetically close to the ionization continuum, resulting in a certain probability of four-photon ionization from the ground state caused by the pump pulse alone (measurement and Channel 1), leading to a background signal. In a separate pump-probe measurement (measurement ), the pump process is the same as in the pump-only case. The population is generated in the excited states, which in turn is ionized by a time-delayed probe pulse. Due to the low laser intensity of the probe pulse, ground state molecules are not ionized by the probe pulse alone, i.e., there is no probe-only background. Consequently, the measured pump-probe spectrum consists of pump-only ionization events (Channel 1) and pump-probe ionization events (Channel 2). Cations, produced in both channels, dependent on the ionization path, can either be stable and detected as parent ions or undergo fragmentation into neutral and ionic fragments. Coincidence detection of electrons and ions allows obtaining separate electron spectra for each ion, i.e., parent and fragment. The excited state of the molecule at the moment of ionization is identified by the measured electron kinetic energy in combination with the energy of the ionizing photon and knowledge of the vertical ionization energy of the excited state. In addition to the information of species and electronic state that is ionized, the related ion mass of the PEPICO spectrum provides insight into the fragmentation behavior. For example, the assignment of the photoelectron kinetic energy to an excited electronic state of the unfragmented molecule and coincidence detection of an ion fragment show that the molecule was intact at the moment of ionization and that fragmentation must have occurred afterwards. Moreover, the population in the excited state can decay to energetically lower states quite quickly, e.g., on a femtosecond timescale [4,7,28]. It is due to this decay that the Channel 2 signal can become significantly smaller than the Channel 1 background, in particular for long delay times, causing a poor signal-to-noise ratio.

3. Bayesian Data Analysis

3.1. Preliminary Considerations

We now introduce our notation and develop the Bayesian algorithm for analyzing the data generated in the experiment described in Section 1 and Section 2. We consider the following standard setup consisting of two experiments on the same target: pump-only and pump-probe, denoted by and , respectively. Each experiment consists of measurements. A measurement of the experiment is performed with exactly one pump pulse, while a measurement of the experiment comprises exactly one pump pulse and one probe pulse. During one measurement, two types of elementary coincidence events are detected, either a molecule is ionized from its ground state (referred to as Channel 1) or from its excited state (Channel 2). The latter is only possible in a pump-probe measurement (). We assume that the number of ionization events in channel is Poisson distributed with some mean ionization rate . For the sake of readability, we suppress the index j in the following considerations. Furthermore, we assume that in each experiment, characterized by a defined delay time between pump and probe pulse, is independent of the occupation of the states, which means that we neglect population depletion effects. We additionally presume the laser intensity, and thus the ionization rate , to fluctuate. Mikosch and Patchkovskii [26,27] proposed to describe with a Gaussian Probability Density Function (PDF). We rather resort to a -distribution of the latter instead, since: (1) the assignment of a Gaussian PDF for is inconsistent with the fact that , while the -distribution includes this constraint naturally; (2) it turns out to be quite convenient mathematically, while the effect of this choice on the result is deemed negligible, which will become apparent later on. We parametrize the -distribution, , with their expectation value and variance . All results match the findings of our previous paper [22] in the limit .

In one elementary event, the involved molecule can have mass and the emitted electron energy . For brevity, we will refer to this particular event as . The ion masses and the electron energies are discretized, , due to the finite resolution of the time-of-flight spectrometer. We will also use the symbol , if we refer to the measurements/experiments or , the symbol j for Channel 1 or 2, and x for the combination of both sets, i.e., . Given that an elementary event happens in channel , the probability that it corresponds to is denoted by: where denotes the conditional complex. The probabilities are properly normalized: In the pump-only measurement (), all molecules are in their respective ground state; therefore, only Channel 1 is allowed and:

We now introduce the spectrum: with a subtle distinction of in the explicit propositions. is conditioned on being a true single coincidence, whereas is conditioned on being a single coincidence, either true or false, and hence also on . The latter is what is actually observed in a pump-only PEPICO measurement, while the former is what is desired. It will become apparent in the following how this additional condition distorts our statistics. We summarize all unknown parameters in the variable . , , , and describe the fluctuations of and according to -distributions, and and are the detection probabilities of ions and electrons, respectively. The probability of detecting an electron with energy and an ion with mass in a single coincidence measurement is: with the marginals and and the detection failure probabilities and . In the second line, we use that was already derived in Appendix 1 of [22]. The appearing integrals are of the type: and solved in Appendix A. Note that the PDF describing the laser fluctuations enters the whole algorithm only within integrals of the above type. Thus, the theory can easily be adapted to different descriptions of the fluctuations by merely exchanging those integrals. The spectrum including false coincidences, , is then given by: turns out to be . The false coincidences are represented by the term , and therefore is a measure of the amount of false coincidences.

We now define analogously to and derive its dependence on and , and , and and , abbreviated in the appendix as q, , and . Therefore, we have to reassign the single coincidence event to the electron stemming from channel and the ion from channel , respectively. Distinguishing the case where the electron and ion are measured in coincidence and stem from the same () or a different (, or ) channel, we obtain in Appendix B: In summary, we have: with the parameters: where and: Up to now, we have defined all variables and dependencies we need for the derivation of the desired posterior distribution presented in the following sections.

3.2. The Posterior PDF

We use Bayes’ theorem to determine the posterior probability of the parameters we want to estimate, , and the spectra and . where capital P shall denote discrete and lower case p continuous distributions. and are two datasets, and contains the count rates and . counts how often the pair was detected as a single coincidence event during the experiment . Dataset contains and , which counts how many measurements lead to the detection of electrons and ions during the experiments and , respectively. In this case, it is expedient to use all detected events, not just single coincidences.

3.3. The Prior PDF

To define the prior PDF in Equation (9), we proceed with: For the parameters’ prior, , we use Jeffreys’ prior for scale variables for and and a flat prior for , , , and . The parameter ranges are , and . The spectra are normalized according to Equation (1) and positive, . Our choice for the prior of the spectra, are Dirichlet priors [30] in the spectra and , with and the normalization being the multivariate beta function. We can always choose the prior to be uninformative (flat) by setting all . The choice of these Dirichlet priors is without loss of generality and justified by the likelihood being multinomial in these variables, as we will see later. Previously, with Equation (7), we found to be a function of and , i.e., . With the usual transformation rules, the first term in Equation (10) becomes: The appearing transformation of the Dirac distribution, with , is derived in Appendix C, and the derivation of the Jacobian determinant, is shown in Appendix D. and denote the total number of bins for ions and electrons , respectively. Putting everything together, Equation (11) becomes: The second term in Equation (10), , is obtained by setting and replacing experiment by and Channel 2 by 1; hence, the prior is fully determined.

3.4. The Likelihood

In this paper, we only use the single coincidence events for estimating the spectra and given the parameters . This is justified by the fact that especially these events include relevant information about the spectrum, because the case of detecting more than one electron-ion pair does not allow one to link an electron to the ion it originates from, and the Bayesian approach would be different. The dataset will be used to determine the unknown parameters . Therefore, the likelihood splits into: The first term separates into: Marginalizing over and introduces the multinomial distribution: according to Appendix 2 of [22]. We now turn to the second term in Equation (14), With being the number of laser pulses; the terms are multinomial distributions, as well, Since the PDF is the same for both experiments, and , we suppress the index in the following considerations. The probability of detecting electrons and ions is: is already derived in Appendix 4 in [22] and has the general form: Inserting in Equation (15) produces: A list of for is given in Appendix E.

3.5. Remarks on the Posterior Sampling

In the previous sections, we fully determined the posterior . Ultimately, we are interested in the spectrum , or rather say its probability , which can easily be achieved by integrating out and in the posterior. For performing this integration and, more generally, computing expectation values, variances, and covariances of the posterior, we resort to numerical techniques, posterior sampling specifically. A suitable technique for sampling from this PDF is Markov Chain Monte Carlo (MCMC), which is based on a Markov chain that has the desired distribution as its equilibrium distribution. The technique is standard in Bayesian probability theory (see [30], especially Chapter 30 and the references therein). To generate the Markov chain, we used the Metropolis–Hastings algorithm with local updates in , , and . We chose the step size for the update of each parameter separately to achieve an acceptance probability of ca. 50%. We discarded the first 20% of the Markov chain as means of thermalization and tested the convergence to the equilibrium distribution with several initial states. Correlations were checked with binning.

4. Mock Data Analysis

In this section, we demonstrate the performance of our algorithm including -fluctuations. There are two disturbing influences in the reconstruction of the spectra: false coincidences and pump-only background. We study these influences separately. First, we investigate the influence of the -fluctuations on the false coincidences. In the second part, we scrutinize our algorithm in the presence of a challenging pump-only background and -fluctuations.

4.1. False Coincidences

In the simulation, we have two different ion masses , called parent (p) and fragment (f). We use the test spectra related to the ones reported in [22,26]. The electronic spectrum of the parent has a step-like form, and the electronic spectrum of the fragment consists of Gaussian peaks; see Figure 3. These challenging test spectra exhibit a strongly-varying parent/fragment-ratio as a function of the electron energy. Since we are firstly interested in the effect of -fluctuations on false coincidences, we suppress the background in our simulations ( and ). Before studying the simulation results in detail, we investigate the influence of -fluctuations on false coincidences by analyzing Equation (7). Since we have only signal contributions (), Equation (7) reduces to Equation (4) with the superscripts and 2 instead of and 1, respectively. The equation describes the connection between the measured spectrum , including false coincidences, and the true spectrum . gives the weight between false and true coincidences,

Figure 3

Simulation with mock data for studying the influence of -fluctuations on false coincidences. The black lines are the spectra used to generate the data; the green (blue) lines including error bands are the reconstructed spectra (not) including -fluctuations in the reconstruction. The parameters are and . For , differences between the algorithms are negligible even at relatively high -fluctuations with ; see spectra (a,b). When choosing (c,d), the algorithm not including -fluctuations produces small deviations, e.g., underestimation of the false coincidences at the first Gaussian in the fragment spectrum.

This equation shows that the influence of -fluctuations on false coincidences is negligible in two cases: first, if the relative fluctuations are small; and second, if the relation holds, the influence is negligible regardless of high relative fluctuations.

In our mock data analysis, we choose the parameters and , , and . Now, we turn to the parameter estimates summarized in Table 1. The algorithm reported in [22], which does not include -fluctuations, provides estimates for , , and with about (for ) and (for ) relative deviation from the real values. The differences are caused by the assumption of having a constant without any statistical noise, which is simply wrong, since we generated our mock data by drawing from the -distribution with the parameters and . As expected, estimated by the algorithm ignoring -fluctuations lies within . The algorithm reported in this paper includes the -fluctuations and is therefore able to give reliable estimates for all parameters.

Table 1

Estimated parameters , , , and . In the lines showing the results of the algorithm presented in [22], is shown instead of . Each value denotes the mean and standard deviation of the parameter’s distribution.

	λ_2	σ_2	ξ_i	ξ_e
Parameters (Figure 3a,b)	1.5	0.5	0.5	0.5
Algorithm in [22]	1.4177±0.0008	-	0.5235±0.0003	0.5235±0.0003
Algorithm with λ-fluctuations	1.499±0.001	0.501±0.002	0.4998±0.0003	0.5000±0.0003
Parameters (Figure 3c,d)	0.5	0.5	0.5	0.5
Algorithm in [22]	0.4365±0.0003	-	0.5618±0.0004	0.5621±0.0004
Algorithm with λ-fluctuations	0.4992±0.0005	0.499±0.001	0.5000±0.0004	0.5004±0.0004

We will now show that even if the parameter estimation of the algorithm not including -fluctuations is slightly off, there is only a small difference in the reconstructed spectra. While at , -fluctuations do not affect the reconstruction (Figure 3a,b), at , the algorithm including -fluctuations (green line) is slightly closer to the true spectrum than the algorithm not including -fluctuations (blue line); see Figure 3c,d. In this parameter regime, neglecting -fluctuations leads to , while Equation (16) produces . Therefore, the algorithm not including -fluctuations underestimates the amount of false coincidences.

Altogether, the influence of -fluctuations on false coincidences in the reconstructed spectra is small in the experimentally-relevant parameter regime, which is in accordance with Mikosch et al. [26]. Still, including -fluctuations is important in the parameter estimation.

4.2. Background Subtraction

To analyze the influence of -fluctuations on the background subtraction, we use the same spectra as in Section 4.1. Now, the step-like spectrum is chosen as the signal and the spectrum consisting of Gaussian peaks as the background. We restrict ourselves to one ion mass to exclude the influence of false coincidences stemming from two ion masses. If we neglect false coincidences, the prefactor in defined in Equation (8) is: This factor represents the ratio between the background (1) and signal (2). The equation shows that the corrections cancel if and ; see Figure 4a,b. If (Figure 4c), the algorithm without -fluctuations (blue line) underestimates the background, leading to peaks in every step. In the opposite case, (Figure 4d), the algorithm without -fluctuations (blue line) overestimates the background, leading to notches in every step. The algorithm including -fluctuations (green line) is able to reconstruct the signal correctly in all scenarios.

Figure 4

Simulated test spectra for studying the influence of -fluctuations on the background subtraction. The parameters are , , and . and are different for every sub-figure. If (a) or (b), both algorithms (with (green line) and without (blue line) including -fluctuations) reconstruct the spectra correctly. and lead to an underestimation of the background when neglecting -fluctuations (c). Overestimation of the background happens in the case of and (d).

Similar to the mock data analysis in the previous section, not including -fluctuations leads to deviations in the parameter estimation in all of the simulated test cases; see Table 2. As expected, estimated by the algorithm ignoring -fluctuations lies within . Obviously, including -fluctuations in the parameter estimation becomes important especially for high or . In summary, the influence of -fluctuations on the background subtraction is important if the magnitude of the relative -fluctuations of one channel is large and does not compare to the other.

Table 2

Estimated parameters , , , , , and . The parameter regimes denoted by the identifications (a–d) are according to Figure 4. For each parameter set, the first line denotes the true value, while Line 2 (3) contains the parameter estimation performed with the algorithm without (with) -fluctuations, respectively. In the lines showing the results of the algorithm ignoring -fluctuations, is shown instead of . Each value denotes the mean and standard deviation of the parameter’s distribution.

	λ_1	λ_2	σ_1	σ_2	ξ_i	ξ_e
(a)	0.5	0.5	0.1	0.1	0.5	0.5
	0.4970±0.0003	0.4972±0.0005	-	-	0.5032±0.0002	0.5029±0.0002
	0.5000±0.0003	0.5007±0.0005	0.098±0.002	0.101±0.005	0.5003±0.0003	0.5000±0.0003
(b)	0.5	0.5	0.5	0.5	0.5	0.5
	0.4367±0.0003	0.4317±0.0004	-	-	0.5618±0.0002	0.5620±0.0002
	0.5002±0.0004	0.5009±0.0006	0.5026±0.0008	0.496±0.002	0.4996±0.0003	0.4998±0.0003
(c)	0.5	0.5	0.1	0.5	0.5	0.5
	0.4840±0.0003	0.4608±0.0005	-	-	0.5201±0.0002	0.5201±0.0002
	0.4991±0.0003	0.5013±0.0006	0.098±0.002	0.500±0.002	0.5001±0.0002	0.5001±0.0003
(d)	0.5	0.5	0.5	0.1	0.5	0.5
	0.4452±0.0003	0.4675±0.0004	-	-	0.5441±0.0002	0.5440±0.0002
	0.4998±0.0004	0.5005±0.0006	0.5002±0.0009	0.102±0.004	0.4998±0.0003	0.4997±0.0003

5. Application to Experimental Data

With our setup for time-resolved PEPICO studies on gas phase molecules (see Figure 2), we were not able to produce stable -fluctuations with a sufficiently large to demonstrate a systematic influence on the reconstructed spectra. The output of our commercial laser system (Coherent Vitara oscillator and Legend Elite Duo amplifier) turned out to be too stable for this purpose. We tried to simulate a less stable setup by operating the Optical Parametric Amplifier (OPA) at intensities below the saturation threshold, which induces -fluctuations. In this setting, however, the pulse energy also undergoes temporal drifts on time scales of the measurement, which are different from the statistical fluctuations and not accounted for in our analysis as a setup would not be operated in this regime.

Nevertheless, by using the algorithm including the -fluctuations, we find that the fluctuations were in the order of in the experimentally-relevant parameter regime. The parameters estimated by the different algorithms are depicted in Table 3. Finally, we note that less stable systems are expected to suffer from significant -fluctuations, which have to be accounted for in the reconstruction process of the spectrum. This will in particular be the case if multiple nonlinear optical processes like in the OPA are used in the pump path of a pump-probe setup. The fluctuations get even more prominent if multiphoton processes are used for the excitation process, because these processes depend nonlinearly on the pulse power.

Table 3

Estimated parameters , , , , , and . Line 1 (2) contains the parameter estimations performed with the algorithm without (with) -fluctuations, respectively. In the line showing the results of the algorithm ignoring -fluctuations, is shown instead of . Each value denotes the mean and standard deviation of the parameter’s distribution.

λ_1	λ_2	σ_1	σ_2	ξ_i	ξ_e
0.3328±0.0009	0.132±0.001	-	-	0.3247±0.0008	0.542±0.001
0.3328±0.0009	0.1335±0.0003	0.0004±0.0010	0.0004±0.0012	0.3245±0.0008	0.541±0.001

6. Conclusions

We used Bayesian probability theory to analyze data obtained from pump-probe photoionization experiments with photoelectron-photoion coincidence detection. We extended the algorithm developed previously in [22] by including fluctuations of the laser intensity as a random variable . Based on challenging mock data, we have demonstrated the reliability of the developed algorithm. In accordance with Mikosch et al. [26], the influence of -fluctuations on false coincidences was small in a certain parameter regime. We derived the condition for negligible influence of -fluctuations on false coincidences, even at high relative fluctuations. and denote the complement of the detection probabilities of electrons () and ions (), respectively. In the case of the pump-only background (Channel 1) and the signal (Channel 2) contained in the pump-probe measurements, neglecting -fluctuations underestimates the background signal in the case . The other way around, the background signal is overestimated. In both scenarios, including -fluctuations is important for estimating the parameters , , , , , and correctly. In our application to the experimental data, we find that the relative laser fluctuations in the experimental setup are fortunately too small to see effects on the reconstructed spectra. With the developed algorithm, we were able to determine the experimental -fluctuations to be in the order of . Compared to conventional subtraction of the pump-only spectrum from the pump-probe spectrum, the Bayesian approach provides several important advantages: (i) It results in a significant increase of the signal-to-noise ratio. (ii) It does not overestimate the pump-only contribution and never leads to negative spectra because the relative weight of the pump-only contribution is self-consistently determined. (iii) Spectral signatures based on false coincidences are eliminated, allowing for higher signal rates. (iv) It includes consistently all prior knowledge, such as positivity, and (v) a confidence interval is obtained for the estimated spectrum. (vi) It is applicable to any experimental situation and noise statistics, as demonstrated in this paper for the case of -fluctuations.