› Forums › Variable Stars › CG Dra: a VSS campaign
- This topic has 197 replies, 9 voices, and was last updated 4 months, 1 week ago by Jeremy Shears.
-
AuthorPosts
-
3 October 2022 at 8:34 am #612800Maxim UsatovParticipant
It did went into bright outburst mode. I’ve been only able to collect 58 measurements this night due to clouds.
Max
Attachments:
8 October 2022 at 11:09 am #612950Maxim UsatovParticipantHi All,
Strange things going on… The timing of minima keeps sliding. I’ve plotted the minimum a few days ago and it was right on the observed minimum. Tonight it is about 0.2 phase off. There is very little data, as it was shot through the hole in the clouds, but the eclipse was captured. Either this is a bug in my code or the period keeps changing. But the code here is rather straightforward – I get the epoch time and simply add periods to it to plot lines. How can a period drift by 0.2 phase in such a short period of time?
Code:
if self.plot_epoch > 0 and self.plot_period > 0:
extrema = []
current_extrema = self.plot_epoch
light_curve_start = light_curve[“DATE”].min()
light_curve_end = light_curve[“DATE”].max()
while current_extrema < light_curve_end:
if light_curve_end > current_extrema > light_curve_start:
extrema.append(current_extrema)
print(“Ephemeris extrema:”, current_extrema)
ax0.axvline(current_extrema – 2400000.50, color=”green”, linestyle=’–‘)
current_extrema = current_extrema + self.plot_periodMax
- This reply was modified 2 years, 1 month ago by Maxim Usatov.
Attachments:
8 October 2022 at 1:34 pm #612953Jeremy ShearsParticipantHmmm…that is very odd Max. Could it be the recent dip is a random fluctuation rather than an eclipse?
8 October 2022 at 3:59 pm #612954Maxim UsatovParticipantI don’t think so. Hopefully more data will follow.
Max14 October 2022 at 10:45 am #613045Maxim UsatovParticipantEclipses No(s): 94
State: QuiescenceThe eclipse minima definitely has shifted or there is an error in how I am calculating minima. I am puzzled. Can O-C analysis be done in Peranso or some other software?
Max
Attachments:
14 October 2022 at 12:48 pm #613050Jeremy ShearsParticipantPeriod analysis can be performed using one of Peranso spectrum analysis tools. But these analyse the data globally, resulting in significant periods being identified, rather than calculating the ToMin that would be required for O-C analysis.
14 October 2022 at 4:50 pm #613054Maxim UsatovParticipantJust double-checked my minima calculation: MetroPSF provides the same ephemeris as Peranso, so looks like there is something going on with periods here. I just wonder if there’s an automated way to identify observed minima of all ~ 90 eclipses we have, or is manual the only way? Anyone brave enough to do O-C analysis? 🙂
Max
15 October 2022 at 10:39 am #613057Maxim UsatovParticipantEclipses No(s): 95
State: QuiescenceTypical U/N/A eclipse for this state.
Max
Attachments:
16 October 2022 at 12:18 am #613059Maxim UsatovParticipantThinking how to automate ToM detection for the data collected so far, I have tried fitting a Gaussian function and finding its minimum – see screenshots attached for the last eclipse recorded. Unfortunately, this method is very sensitive to the data boundaries if the data is noisy. If I include more points to the left or right, I get a different curve fit. Even the Kwee-van Woerden method displays the same sensitivity in Peranso. Attached are a few examples. I was thinking to try to find ToM automatically using some predefined phase shift, say, delta phi = 0.2, from calculated minima to set search boundaries, but the choice of delta phi produces different ToM results with both Gaussian fits and the KvW.
Max
Attachments:
16 October 2022 at 2:22 pm #613064Maxim UsatovParticipantOK, I have managed to do Gaussian fits. First, I define a delta phi search region, then identify the faintest magnitude within the search region, smooth out the noise, determine eclipse well and then fit a Gaussian curve to the well. Presenting the first O-C table. Looks like it’s linear. ANOVA period determined incorrectly?
Max
Attachments:
16 October 2022 at 5:58 pm #613067Maxim UsatovParticipantEclipses No(s): 96
State: Quiescence
U/H?/PEHAttachments:
20 October 2022 at 8:29 am #613200Maxim UsatovParticipantEclipses No(s): 97
State: Rising
U/N/A-type eclipse.Attachments:
24 October 2022 at 11:26 pm #613249Maxim UsatovParticipantFlat O-C (ignore a few outliers) with P=0.188644 d.
Attachments:
25 October 2022 at 5:32 pm #613270Maxim UsatovParticipantEclipses No(s): 98
State: Quiescent.
U/N/A-type eclipse.Attachments:
25 October 2022 at 7:05 pm #613273Jeremy ShearsParticipantLooking good Max. Any ideas about outliers? I’ve sometimes struggled to assign proper ToMs in specific eclipses, especially small ones, when the underlying light curve has been busy with bumps and humps.
- This reply was modified 2 years ago by Jeremy Shears.
26 October 2022 at 9:57 am #613279Maxim UsatovParticipantThese are exactly the cases where the algorithm has failed to fit a Gaussian due to noise or partial data. I also need to resample the data to a finer grid, as described by KvW, to refine the result.
26 October 2022 at 11:51 pm #613284Jeremy ShearsParticipantThanks Max. I suspected this was the case. But ignoring these, the O-C is looking good 👍🏻
28 October 2022 at 10:29 am #613323Maxim UsatovParticipantEclipses No(s): 100 (didn’t show a pic for 99th, partial data)
State: Outburst
U/L/S-type eclipse.Attachments:
28 October 2022 at 11:48 am #613325Jeremy ShearsParticipantCongratulations on the century, Max! Very nice work, indeed.
28 October 2022 at 2:54 pm #613326Maxim UsatovParticipantThank you very much for organizing this, Jeremy!
-
AuthorPosts
- You must be logged in to reply to this topic.