[1]:
import clustertools as cts
import numpy as np
Setup¶
As discussed in the documentation, a wrapper has been written around the LIMEPY code (Gieles, M. & Zocchi, A. 2015, MNRAS, 454, 576 & Claydon, I., Gieles, M., Varri, A.L., Heggie, D.C., Zocchi, A. 2019, MNRAS, 487, 147) to set up a StarCluster that represents a Galactic Globular Cluster, where the structural information is taken from either de Boer, T. J. L., Gieles, M., Balbinot, E., Hénault-Brunet, V., Sollima, A., Watkins, L. L., Claydon, I. 2019, MNRAS, 485, 4906 (default) or Harris, W.E.
1996 (2010 Edition), AJ, 112, 1487. Orbital information is taken from Vasiliev E., 2019, MNRAS, 484,2832 via galpy’s Orbit.from_name() feature Bovy J., 2015, ApJS, 216, 29.
To setup NGC6101, for example:
[2]:
cluster=cts.setup_cluster(ctype='limepy',gcname='NGC6101')
print(cluster.ntot,cluster.rm,cluster.mtot)
print(cluster.xgc,cluster.ygc,cluster.zgc,cluster.vxgc,cluster.vygc,cluster.vzgc)
LOAD_CLUSTER MADE USE OF:
Gieles, M. & Zocchi, A. 2015, MNRAS, 454, 576
Vasiliev E., 2019, MNRAS, 484,2832
de Boer, T. J. L., Gieles, M., Balbinot, E., Hénault-Brunet, V., Sollima, A., Watkins, L. L., Claydon, I. 2019, MNRAS, 485, 4906
1936 11.8971760895 1936.332
-2005.21009948 -9348.18148437 -3945.46817625 -296.181213344 82.7743019402 -190.8475368
Unless otherwise specified, the cluster is setup to be in pckms units in clustercentric coordinates. The number of stars is set using mbar variable which has a default of 1.0 solar masses.
Clusters can easily be viewed as they would be in the sky using the skyplot command:
[3]:
cts.skyplot(cluster)
[3]:
[<matplotlib.lines.Line2D at 0x140654ac0>]
