def radialDistribution(this, bins, simulation = -1, states = "all", label = ''): print "Calculating radial distribution"; simPath = this.simulations[simulation]; length = this.stateNumbers[simulation]; mdf.combine(simPath, length, ); r, g = mda.calculateRadialDistribution(simPath, bins, length, states); mdplot.plotRadialDistribution(this.path, g, r, label = label); return r, g;