sjaffa/q_plus

Tool for astronomical data analysis, quantifying the structure of young star clusters using a fractal dimension.

Q_plus

Q+: an algorithm for statistically quantifying sub-structure in star clusters

INPUTS:

x,y: two-dimensional position of stars in parsecs

OUTPUTS

filename (Possbile Q warning if centrally concentrated) PC1 and PC2 values (Possibe edge of parameter space warning) [D,D_err],[L,L_err],[C,C_err]

Where D,L,C areestimated fractal dimension, number of levels, density scaling exponent and D_err,L_err,C_err are errors on these.

DOCUMENTATION

This algorithm: -builds the complete graph and minimum spanning tree of the points, -derives various statistics of structure from these, -compares these statistics to a set with known structure, and -estimates the most likely structure parameters and their error. See paper for a full explanation of the algorithm.

Usage:

In a terminal run: python gridrecon.py

Where is the path to the text file you want to run and is a list of the indices if the x and y positions. Example:

python gridrecon.py /home/sjaffa/Dropbox/Data/IC348.txt [0,1]

LICENSE

This work is shared under and MIT license, please see the LICENSE file

q_plus's Issues

Is G equivalent to L in the code?

Hi Sarah,

It's Becky Arnold, we've met a couple of times at conferences, I did my PhD at Sheffield with Simon Goodwin, I'm at Keele now. I'm looking at applying Q+ to compare some simulations and real clusters and I'm confused because your paper talks about the L parameter, put the code outputs a parameter called G (in addition to D and C).

Is G just another name for L, or is it something different?

Thanks,

Becky