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()
[3]:
plotorbits = results.PlotOrbits()
[4]:
plotseppa = results.PlotSepPA()
[ ]: