afeinstein20 / eleanor

A tool for light curve extraction from the TESS FFIs.
MIT License
92 stars 39 forks source link

Getting different TPFs for the same star in different sectors #268

Open iamasciencer opened 1 year ago

iamasciencer commented 1 year ago

Hi,

I am having some issues when trying to perform series analysis for some TIC sources. To get the best out of the data, I am trying to extract all possible light curve measurements. However, when i plotted the light curves and periodograms, i found some weird behaviour. When i the looked at the TPFs and the apertures for a few stars, they varied.

I am using this code for TIC 28328233:

stars = eleanor.multi_sectors(tic=28328233, sectors="all")

    for star in stars:

        print(star)

        data = eleanor.TargetData(star, height=15, width=15, bkg_size=31, do_psf=False, do_pca=False)

        print(data)

        q = data.quality == 0
        eleanor.TargetData.custom_aperture(data, shape='circle', r=2, pos=[7,6], method='exact')
        eleanor.TargetData.get_lightcurve(data)

        fig, (ax1, ax2) = plt.subplots(ncols=2, figsize=(15,4), gridspec_kw={'width_ratios':[1,3]})
        ax1.imshow(data.tpf[0])
        ax1.imshow(data.all_apertures[0], cmap='Greys', alpha=0.5)
        ax1.set_title('Aperture over TPF')
        #plt.imshow(data.aperture)

        ax2.plot(data.time[q], data.raw_flux[q]/np.nanmedian(data.raw_flux[q]), 'k', label='Raw')
        ax2.plot(data.time[q], data.corr_flux[q]/np.nanmedian(data.corr_flux[q]) - 0.015, 'r', label='Corrected')
        ax2.legend()
        ax2.set_xlabel('Time [BJD - 2457000]')
        ax2.set_ylabel('Normalized Flux')
        plt.show()

so, basically the same procedure as in the tutorial at https://adina.feinste.in/eleanor/.

The output in the terminal looks like this:

2023-09-06 22:17:13.459783: W tensorflow/stream_executor/platform/default/dso_loader.cc:64] Could not load dynamic library 'libcudart.so.11.0'; dlerror: libcudart.so.11.0: cannot open shared object file: No such file or directory 2023-09-06 22:17:13.459814: I tensorflow/stream_executor/cuda/cudart_stub.cc:29] Ignore above cudart dlerror if you do not have a GPU set up on your machine. Found star in Sector(s) 5 32 44 INFO: Found cached file ./mastDownload/HLSP/hlsp_eleanor_tess_ffi_postcard-s0005-1-4-cal-0216-0592_tess_v2_pc/hlsp_eleanor_tess_ffi_postcard-s0005-1-4-cal-0216-0592_tess_v2_bkg.fits with expected size 73658880. [astroquery.query] INFO: Found cached file ./mastDownload/HLSP/hlsp_eleanor_tess_ffi_postcard-s0005-1-4-cal-0216-0592_tess_v2_pc/hlsp_eleanor_tess_ffi_postcard-s0005-1-4-cal-0216-0592_tess_v2_pc.fits with expected size 147421440. [astroquery.query] INFO: Found cached file ./mastDownload/HLSP/hlsp_eleanor_tess_ffi_postcard-s0005-1-4-cal-0216-0592_tess_v2_pc/hlsp_eleanor_tess_ffi_postcard-s0005-1-4_tess_v2_pm.txt with expected size 222713. [astroquery.query] No eleanor postcard has been made for your target (yet). Using TessCut instead. No eleanor postcard has been made for your target (yet). Using TessCut instead. <eleanor.source.Source object at 0x7f07e46216d0> <eleanor.targetdata.TargetData object at 0x7f07e4630640> <eleanor.source.Source object at 0x7f07e4630d90> <eleanor.targetdata.TargetData object at 0x7f07e1685c10> <eleanor.source.Source object at 0x7f07e4630cd0> <eleanor.targetdata.TargetData object at 0x7f07e05e7eb0>

Which is fine so far. however, looking at the images, they look vastly different from each other. For example, stars appear that are not present in the background when i look at the Aladin tool at Simbad (although that may be biased as well..)

Here are the plots:

eleanor_aperture_test_1 For sector 5

eleanor_aperture_test_2 for sector 32

eleanor_aperture_test_3 For sector 44

Could there be an issue with TessCut? or did I do something majorly wrong and haven't found the issue (I am not the most experienced programmer)?

As for my setup: I am using VSCode Version 1.81,1 (last updated about a month ago) on Ubuntu 20.04 on my laptop.

Thank you very much! Cheers