#NEWLINE# def midplane_cumulative_density( self , r ) : #NEWLINE# #INDENT# '\n        Find the cumulative column density as a function of radius.\n\n        The cumulative density is measured outwards from the origin, and in\n        the midplane.\n\n        Parameters\n        ----------\n        r : np.ndarray\n            Array of values of the radius up to which to tabulate the\n            cumulative density.\n\n        Returns\n        -------\n        rho : np.ndarray\n            Array of values of the cumulative density.\n        ' #NEWLINE# self . _check_all_set ( ) #NEWLINE# if ( self . rmax <= self . rmin ) : #NEWLINE# #INDENT# logger . warn ( 'Ignoring power-law envelope, since rmax < rmin' ) #NEWLINE# return np . zeros ( r . shape ) #NEWLINE# #UNINDENT# return ( ( self . rho_0 * integrate_powerlaw ( self . rmin , r . clip ( self . rmin , self . rmax ) , r . power ) ) / ( self . r_0 ** self . power ) )