Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Generalized lightcurves #754

Merged
merged 96 commits into from
Dec 15, 2023
Merged

Generalized lightcurves #754

merged 96 commits into from
Dec 15, 2023

Conversation

matteobachetti
Copy link
Member

@matteobachetti matteobachetti commented Sep 15, 2023

An attempt at more general time series objects. Not taking for granted counts or count rates, having multiple array attributes, etc.
It might be, for example, spectroscopic or polarimetric light curves with multiple measurements per time point

Documentation for these changes can be found at StingraySoftware/notebooks#73

(Depends on #774) Merged

@matteobachetti matteobachetti marked this pull request as draft September 15, 2023 10:10
@codecov
Copy link

codecov bot commented Sep 18, 2023

Codecov Report

All modified and coverable lines are covered by tests ✅

Comparison is base (424fa59) 96.18% compared to head (29cf3a0) 94.22%.

Additional details and impacted files
@@            Coverage Diff             @@
##             main     #754      +/-   ##
==========================================
- Coverage   96.18%   94.22%   -1.96%     
==========================================
  Files          43       43              
  Lines        8020     8473     +453     
==========================================
+ Hits         7714     7984     +270     
- Misses        306      489     +183     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@pep8speaks
Copy link

pep8speaks commented Sep 18, 2023

Hello @matteobachetti! Thanks for updating this PR. We checked the lines you've touched for PEP 8 issues, and found:

Line 92:101: E501 line too long (101 > 100 characters)
Line 908:101: E501 line too long (166 > 100 characters)
Line 947:101: E501 line too long (166 > 100 characters)
Line 964:101: E501 line too long (166 > 100 characters)

Line 1490:101: E501 line too long (103 > 100 characters)
Line 1504:101: E501 line too long (103 > 100 characters)
Line 1509:101: E501 line too long (116 > 100 characters)
Line 1519:101: E501 line too long (103 > 100 characters)

Comment last updated at 2023-12-15 13:24:42 UTC

stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
@matteobachetti matteobachetti force-pushed the generalized_lightcurves branch 10 times, most recently from f9f1421 to e29dcbb Compare September 26, 2023 16:20
@matteobachetti matteobachetti force-pushed the generalized_lightcurves branch 2 times, most recently from 23694cb to ea4bdb1 Compare September 29, 2023 08:47
@matteobachetti
Copy link
Member Author

Fixing things for making them work with the Generalized Light curves produced a monster PR. And there are more modifications to do. I will spin off the fixes to smaller PRs, and try to keep this PR on topic.

@matteobachetti matteobachetti force-pushed the generalized_lightcurves branch 2 times, most recently from 25eaa51 to faa4af7 Compare October 7, 2023 07:39
@matteobachetti matteobachetti requested a review from mgullik October 9, 2023 15:42
@matteobachetti matteobachetti marked this pull request as ready for review October 9, 2023 15:42
@matteobachetti matteobachetti force-pushed the generalized_lightcurves branch 3 times, most recently from 237b7af to 59ec903 Compare October 11, 2023 16:03
@matteobachetti matteobachetti force-pushed the generalized_lightcurves branch from 6b9407c to 2ddb7a6 Compare October 25, 2023 09:36
Copy link
Collaborator

@mgullik mgullik left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi @matteobachetti,

this is an important change in the structure of the code. I will probably divide the review into a few steps

_join_timeseries and join method:
Can this method join StingrayTimeseries with EventList objects?
If so, it might be useful to specify it (maybe I missed it)
Does this method take into account when the MJDREF is present only in one or a portion of the objects that are going to be joint?

In lightcurve.py
Is Lightcurve class a child class of StingrayTimeseries?
I don't see any change in the definition of the class Lightcurve, but I see a few methods that are deleted from here, I think this happened because Lightcurve is going to inherit all the methods of StingrayTimeseries, correct?

in crossspectrum.py and powerspectrum.py
a few comments directly in the text about names

I'm going to review the tests later

stingray/base.py Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/powerspectrum.py Outdated Show resolved Hide resolved
stingray/crossspectrum.py Outdated Show resolved Hide resolved
stingray/crossspectrum.py Show resolved Hide resolved
stingray/crossspectrum.py Outdated Show resolved Hide resolved
stingray/crossspectrum.py Outdated Show resolved Hide resolved
stingray/crossspectrum.py Outdated Show resolved Hide resolved
@matteobachetti
Copy link
Member Author

@mgullik thanks for the review! I fixed all the docstrings. Indeed, many of those methods were moved from EventList and Lightcurve, and the docstrings were left with the original class names.

In general, all child classes (EventList, Lightcurve) inherit the methods from StingrayTimeseries or re-write them if incompatible. join is partially re-written, but e.g. EventList.join does call _join_timeseries. However, note that there is a check that you cannot join different types. If you start from an event list, you cannot merge it with a light curve, even if both of them are children of StingrayTimeseries.

Changes of MJDREF are accounted for. Note that you cannot be a StingrayTimeseries without defining MJDREF, at least as 0 (It is a field automatically defined in the StingrayTimeseries class, see line 1085).

Copy link
Collaborator

@mgullik mgullik left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi @matteobachetti

this is the second part of the review.
You find every comment on the code.

I'm still not sure why StingrayTimeseries won't replace the EventList.
If the time array in StingrayTimeseries doesn't have to be well-sampled, it is possible to lead an event file in StingrayTimeseries and then use rebin (dt = float) to produce/transform the StingrayTimeseries into a light curve.
Of course, I don't want to suggest getting rid of EventList and Lightcurve, but I'd like to know if there is any functionality that EventList has and StingrayTimeseries hasn't

stingray/tests/test_base.py Outdated Show resolved Hide resolved
stingray/tests/test_base.py Outdated Show resolved Hide resolved
stingray/tests/test_base.py Outdated Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
Copy link
Collaborator

@mgullik mgullik left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi @matteobachetti,

thank you for addressing all my comments, I approve your Pull Request

Copy link
Member

@dhuppenkothen dhuppenkothen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Great work Matteo! I think this is a huge improvement! I have lots of minor nitpicks, a question about the relationship between EventList and TimeSeries, and a question for the future: with this added functionality, I'm expecting that people will want to somehow connect things like RMFs and ARFs to stingray. I know Matteo Lucchini is working on some code to do that. The question is: have you thought about how these would fit into the context of the current objects? Would they be separate objects? Would they get attached to a TimeSeries object as an attribute?

stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/base.py Outdated Show resolved Hide resolved
stingray/gti.py Show resolved Hide resolved
stingray/lightcurve.py Show resolved Hide resolved
stingray/powerspectrum.py Show resolved Hide resolved
stingray/tests/test_base.py Outdated Show resolved Hide resolved
stingray/tests/test_base.py Show resolved Hide resolved
@matteobachetti
Copy link
Member Author

@dhuppenkothen I should have addressed all your comments in one way or the other.
About responses: they could definitely be added as attributes, and then used for whatever the user wants. Or we could just create another Spectrum object, maybe compatible with the equivalent xspec or sherpa objects. I think feedback from people like @matteolucchini1 , actually doing these things, should drive the choice though.

@matteolucchini1
Copy link
Collaborator

matteolucchini1 commented Dec 15, 2023

@matteobachetti You can find the response stuff I've done so far here https://github.com/matteolucchini1/neXTsPec_prototype
Currently what it does is simply take an impulse response in units of count rate as a function of energy and time (although the time axis could easily be e.g. pulse phase or Fourier frequency), and convolve that with response files to give units of count rate as a function instrument channel and time/whatever else. In practice, what I actually chose to do is simply a matrix multiplication pretty much like x-psi, because it's much more straightforward than the xspec/isis/sherpa way of handling things and the computational cost is negligible.

I'm not too familiar with Stingray, but in principle I think it should be possible to adapt my code to have a more robust energy:channel mapping, and then replace the rough_calibration function.

Copy link
Member

@dhuppenkothen dhuppenkothen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this is good to go now! Two requests for separate PRs:
(1) add an rst page to the documentation that explains the overall hierarchy of stingray objects, including where TimeSeries fits in and how it relates to LightCurve etc
(2) I haven't looked at the notebooks yet, but it would be useful to make some of the additional functionality explicit there that might be confusing to the user. I'm thinking of the functionality to create binned TimeSeries from EventLists and the different syntax for creating PSDs out of light curves versus time series.

@matteobachetti
Copy link
Member Author

@matteolucchini1 thanks! We have a function that does something similar in HENDRICS: https://github.com/StingraySoftware/HENDRICS/blob/c0e2c68c965cd54404d7f1b859f1ad3bc96aa9f8/hendrics/calibrate.py#L133
It's probable that yours will be more advanced, but feel free to take a look and recycle/propose changes at will

@matteobachetti matteobachetti added this pull request to the merge queue Dec 15, 2023
Merged via the queue into main with commit 753eb2d Dec 15, 2023
14 of 15 checks passed
@matteobachetti matteobachetti deleted the generalized_lightcurves branch December 15, 2023 13:33
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

5 participants