matteolucchini1

Hi @capy-on-caffeine ! One thing that might be worth trying if you want something more fancy is a more realistic case where your lightcurve is a sum of a QPO...

@matteobachetti Ah fair enough. I do still think it would be nice to show different ways of loading event files (e.g. here's what happens when you load an heasoft-compatible file,...

I really like that tutorial! I left a couple more minor comments otherwise it's getting close :)

Hi @matteobachetti , I would be very happy to be a co-author, thanks for including me! My current affiliation is "Anton Pannekoek Institute, University of Amsterdam, Science Park 904, 1098...

Sure! Following from the code above: ``` stingray_bins = nicer_rebin.emin[nicer_rebin.emin>0.5] stingray_bins = stingray_bins[stingray_bins

The value of plasma beta for velsw=15 is unphysical, it should be 0, or an appropriate small value (see figure 8 here https://arxiv.org/pdf/2108.12011 and the discussion in that section)

The output directory being empty is because of infosw being set to 0, it needs to be at least 1 to produce files (otherwise the code just returns the array...

@matteobachetti yea I agree. I think the easiest implementation to make this a bit more bulletproof would be to replace the calls with np.rint() to just int() - there's a...

I'm confused - the implementation in `avg_pds_from_timeseries` can only work if we are passing a time series with a known size (because n_bin is set explicitely from there). If instead...