Quick Start GuideΒΆ

[2]:
# Gaia DR2 source ids:
DSTucA = 6387058411482257536
DSTucB = 6387058411482257280
# Mass estimates, must be a tuple of (value,uncertainty)
# in solar masses:
massA = (0.97, 0.04)
massB = (0.87, 0.04)

# Import the Fitter and FitOrbit objects:
from lofti_gaia.lofti import Fitter, FitOrbit

# Initialize the fitter object:
fitterobject = Fitter(DSTucA,           # source id object 1
                      DSTucB,           # source id object 2
                      massA,            # mass object 1
                      massB,            # mass object 2
                      Norbits = 100     # number of desired accepted orbits for the posterior orbit sample
                     )

# run orbit fit:
orbits = FitOrbit(fitterobject)

# get fit results:
results = orbits.results

Created TAP+ (v1.2.1) - Connection:
        Host: gea.esac.esa.int
        Use HTTPS: True
        Port: 443
        SSL Port: 443
Created TAP+ (v1.2.1) - Connection:
        Host: geadata.esac.esa.int
        Use HTTPS: True
        Port: 443
        SSL Port: 443
Saving orbits in FitResults.2021.08.02.12.08.23.txt
inital chi min 56.70697569685297
100% (100 of 100): |####################|  Done...
Final Norbits: 100
[2]:
hists = results.PlotHists()
../_images/tutorials_QuickStart_2_0.png
[3]:
plotorbits = results.PlotOrbits()
../_images/tutorials_QuickStart_3_0.png
[4]:
plotseppa = results.PlotSepPA()
../_images/tutorials_QuickStart_4_0.png
[ ]: