@@ -16,12 +16,12 @@ Code in the [functions folder](https://ict.inaf.it/gitlab/heag-oar/scsearch/-/tr
powdist.m defines a function which briefly analyses a light curve to search for crosstalk signatures.
It plots the light curve together with a simple Leahy-normalized power spectrum (PS) of the data and the same PS rebinned onto a quasi-logarithmic grid. Then, assuming that the observation is mostly composed of Poissonian noise modified by crosstalk, it generates the expected distribution in three different crosstalk models, and plots the counts-per-bin distribution both in the data and the models.
By checking that the high-frequency PSs are flat (i.e., the noise tends to white noise), and that the alignment between at least one of the models and the data count distribution is good, users can safely just renormalize the powers in their PSs then by dividing them by half the average power at high frequencies.
By checking that the high-frequency PSs are flat (i.e., the noise tends to white noise), and that the alignment between at least one of the models and the data count distribution is good, users can safely just renormalize the powers in their PSs, even after correcting them or the underlying light curve for orbital motions: by dividing the PSs by half the average power at high frequencies (returned by the function as `HalfMeanPow`), they can be treated as normally done for the purpose of power significance and sensitivity estimation.
The function takes five arguments as input (NyqFreq, LC, Binned, PSname, Noisename) and returns three (extrVoverE, CrossTalkProb, HalfMeanPow).
#### Input arguments:
-`NyqFreq` is the Nyquist frequency in Hz
-`NyqFreq` is the Nyquist frequency in Hz (we suggest to use at least 1500 Hz)
-`LC` (light curve) is the vector of either binned counts, or photon times of arrival (ToAs)
- If `Binned` is true, we assume uniform binning and dt = 1/(2*NyqFreq), otherwise we assume the elements of LC are actually the ToAs, and the function bins the data with the aforementioned dt
-`PSname` and `Noisename` are the filenames (or complete paths) of graphs to be saved: the first figure will contain the power spectra (PSs) and light curve, and the second will compare the counts distribution with the one estimated through assumptions on the crosstalk distributions. Any of the [types supported](https://it.mathworks.com/help/matlab/ref/exportgraphics.html#mw_3db70293-26d3-4e3f-9abf-6b96488e4dde) by MATLAB's exportgraphics function should be fine, but formatting tests were carried out with pdfs
