########################################################
# Started Logging At: 2025-04-04 13:53:12
########################################################
########################################################
# # Started Logging At: 2025-04-04 13:53:12
########################################################
get_ipython().run_line_magic('pip', 'install daschlab')
import daschlab

# Open a session
session = daschlab.open_session(".")

# Define the coordinates for FU Orionis (J2000)
ra = "05h 45m 22.4s"
dec = "+09° 04′ 12″"

# Retrieve the light curve data
lightcurve = session.lightcurve(ra=ra, dec=dec, radius=5/60)  # radius in degrees
import daschlab

# Open a session
session = daschlab.open_session(".")

# Define the coordinates for FU Orionis (J2000)
ra = "05h 45m 22.4s"
dec = "+09° 04′ 12″"

# Retrieve the light curve data
lightcurve = session.lightcurve()#ra=ra, dec=dec, radius=5/60)  # radius in degrees
import daschlab
from astropy.coordinates import SkyCoord

# Open a session
session = daschlab.open_session(".")

# Define the coordinates for FU Orionis (J2000)
ra = "05h 45m 22.4s"
dec = "+09° 04′ 12″"

# Retrieve the light curve data
lightcurve = session.lightcurve(SkyCoord(ra=ra, dec=dec, frame='fk5'), radius=5/60)  # radius in degrees
import daschlab
from astropy.coordinates import SkyCoord

# Open a session
session = daschlab.open_session(".")

# Define the coordinates for FU Orionis (J2000)
ra = "05h 45m 22.4s"
dec = "+09° 04′ 12″"

# Retrieve the light curve data
lightcurve = session.lightcurve(SkyCoord(ra=ra, dec=dec, frame='fk5'),)  # radius in degrees
from daschlab import open_session

sess = open_session(".")
from daschlab import open_session

sess = open_session(".")import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")

import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")

import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

sess.select_target(fu_ori_coords)

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")

import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')


# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")

import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

lightcurve = session.lightcurve("FU Orionis")

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
from daschlab import open_session

sess = open_session(".")

import daschlab
from astropy.coordinates import SkyCoord
import astropy.units as u

# Open a session with DASCH
session = daschlab.Session(".")

# Define the coordinates for FU Orionis (J2000)
fu_ori_coords = SkyCoord(ra=86.343333*u.deg, dec=9.070000*u.deg, frame='icrs')

#lightcurve = session.lightcurve("FU Orionis")

# Query the reference catalog (e.g., APASS) around FU Orionis
refcat_results = session.refcat(
    refcat="apass",
    ra_deg=fu_ori_coords.ra.deg,
    dec_deg=fu_ori_coords.dec.deg,
    radius_arcsec=300  # Search radius of 5 arcminutes
)

# Display the reference catalog results
print(refcat_results)

# Assuming FU Orionis is the first entry in the results
fu_ori_entry = refcat_results[0]

# Retrieve the light curve using the reference number and GSC bin index
lightcurve = session.lightcurve(
    refcat="apass",
    ref_number=fu_ori_entry["ref_number"],
    gsc_bin_index=fu_ori_entry["gsc_bin_index"]
)
# Fill in the SIMBAD-resolvable name of your target source:
SOURCE = "source name goes here"

# Leave this unchanged unless you're sure that you want something else:
REFCAT = "apass"

# Fill in the SIMBAD-resolvable name of your target source:
SOURCE = "FU Orionis"

# Leave this unchanged unless you're sure that you want something else:
REFCAT = "apass"

# fu orionis

# Get the main module:
import daschlab

# Set up Bokeh plots:
from bokeh.io import output_notebook
output_notebook()

# Get some other imports that will be useful:
from astropy import units as u
from bokeh.plotting import figure, show
import numpy as np
from daschlab.photometry import AFlags, BFlags

sess = daschlab.open_session(source=SOURCE)
sess.select_target(name=SOURCE)
sess.select_refcat(REFCAT)
#[Out]# <daschlab.Session at 0x1522474ead50>
await sess.connect_to_wwt()
sess.refcat()[:12]
#[Out]# <RefcatSources masked=True length=12>
#[Out]#        ref_text           ref_number   gsc_bin_index ... refcat local_id
#[Out]#                                                      ...                
#[Out]#         str22               uint64         uint32    ...  str5   int64  
#[Out]# ---------------------- --------------- ------------- ... ------ --------
#[Out]# APASS_J054522.4+090412 405452241090412      97795899 ...  apass        0
#[Out]#          N300323123094  11300323123094      97818650 ...  apass        1
#[Out]#          N300323122221  11300323122221      97773147 ...  apass        2
#[Out]#          N300323122224  11300323122224      97773148 ...  apass        3
#[Out]#          N300323122222  11300323122222      97818650 ...  apass        4
#[Out]#          N300323122219  11300323122219      97818649 ...  apass        5
#[Out]#          N300323123042  11300323123042      97773148 ...  apass        6
#[Out]#          N300323113884  11300323113884      97773148 ...  apass        7
#[Out]#          N300323113816  11300323113816      97773147 ...  apass        8
#[Out]#          N300323123112  11300323123112      97818650 ...  apass        9
#[Out]#          N300323123103  11300323123103      97818651 ...  apass       10
#[Out]#          N300323122235  11300323122235      97818650 ...  apass       11
sess.refcat().show()
TARGET_ID = 0
lc = sess.lightcurve(TARGET_ID)
lc.plot()
#[Out]# figure(id='p1001', ...)
await sess.connect_to_wwt()
# Fill in the SIMBAD-resolvable name of your target source:
SOURCE = "Westerlund 2"

# Leave this unchanged unless you're sure that you want something else:
REFCAT = "apass"

# fu orionis

# Get the main module:
import daschlab

# Set up Bokeh plots:
from bokeh.io import output_notebook
output_notebook()

# Get some other imports that will be useful:
from astropy import units as u
from bokeh.plotting import figure, show
import numpy as np
from daschlab.photometry import AFlags, BFlags

sess = daschlab.open_session(source=SOURCE)
sess.select_target(name=SOURCE)
sess.select_refcat(REFCAT)
#[Out]# <daschlab.Session at 0x15222ddb0e00>
await sess.connect_to_wwt()
sess.refcat()[:12]
#[Out]# <RefcatSources masked=True length=12>
#[Out]#        ref_text           ref_number   gsc_bin_index ... refcat local_id
#[Out]#                                                      ...                
#[Out]#         str22               uint64         uint32    ...  str5   int64  
#[Out]# ---------------------- --------------- ------------- ... ------ --------
#[Out]#          S111111026652  12111111026652      13025859 ...  apass        0
#[Out]# APASS_J102358.5-574545 410235852574545      13025859 ...  apass        1
#[Out]#          S111111026619  12111111026619      13025860 ...  apass        2
#[Out]#          S111111026825  12111111026825      13025859 ...  apass        3
#[Out]#          S111111026835  12111111026835      13025859 ...  apass        4
#[Out]#          S111111026827  12111111026827      13025860 ...  apass        5
#[Out]#          S111111037981  12111111037981      13013571 ...  apass        6
#[Out]#          S111111026828  12111111026828      13025860 ...  apass        7
#[Out]# APASS_J102356.7-574528 410235672574528      13025859 ...  apass        8
#[Out]#          S111111026743  12111111026743      13025859 ...  apass        9
#[Out]#          S111111037976  12111111037976      13013571 ...  apass       10
#[Out]# APASS_J102356.2-574530 410235622574530      13025859 ...  apass       11
sess.refcat().show()
# Fill in the SIMBAD-resolvable name of your target source:
SOURCE = "Westerlund 2"

# Leave this unchanged unless you're sure that you want something else:
REFCAT = "apass"

# fu orionis

# Get the main module:
import daschlab

# Set up Bokeh plots:
from bokeh.io import output_notebook
output_notebook()

# Get some other imports that will be useful:
from astropy import units as u
from bokeh.plotting import figure, show
import numpy as np
from daschlab.photometry import AFlags, BFlags

sess = daschlab.open_session(source=SOURCE)
sess.select_target(name=SOURCE)
sess.select_refcat(REFCAT)
#[Out]# <daschlab.Session at 0x15222ddb0e00>