A method for determining a network topology (NT) of a hierarchical network , the method comprising the steps of: (a) providing (S1) an outage state matrix (X) representing an outage state of leaf nodes at the lowest hierarchical level of said hierarchical network ; (b) decomposing (S2) the state matrix (X) into a first probability matrix ( Ŵ ) indicating for each inner node of said hierarchical network the probability that the respective inner node forms the origin of an outage at the lowest hierarchical level of said hierarchical network and into a second probability matrix (H) indicating for each leaf node at the lowest hierarchical level of said hierarchical network the probability that an inner node forms a hierarchical superordinate node of the respective leaf node at the lowest hierarchical level of said hierarchical network ; and (c) evaluating (S3) the second probability matrix (H) to determine the network topology of said hierarchical network wherein the hierarchical network is a hierarchical power supply network (PSN), wherein the hierarchical power supply network (PSN) comprises at its lowest hierarchy level leaf nodes formed by smart meters generating event data streams (EDS) indicating a current outage state of the leaf nodes of said power supply network (PSN), wh erein the outage state matrix (X) derived from event data streams (EDS) is a two-dimensional matrix, wherein the first dimension of said outage state matrix (X) represents time, and wherein the second dimension of said outage state matrix (X) represents smart meters forming leaf nodes at the lowest hierarchical level of said power supply network (PSN), wherein each column of said two-dimensional outage state matrix (X) represents states of a respective smart meter at particular points in time, and wherein each row of said outage state matrix represents the states of all smart meters at a particular point in time. (a) providing (S1) an outage state matrix (X) representing an outage state of leaf nodes at the lowest hierarchical level of said hierarchical network ; (b) decomposing (S2) the state matrix (X) into a first probability matrix ( Ŵ ) indicating for each inner node of said hierarchical network the probability that the respective inner node forms the origin of an outage at the lowest hierarchical level of said hierarchical network and into a second probability matrix (H) indicating for each leaf node at the lowest hierarchical level of said hierarchical network the probability that an inner node forms a hierarchical superordinate node of the respective leaf node at the lowest hierarchical level of said hierarchical network ; and (c) evaluating (S3) the second probability matrix (H) to determine the network topology of said hierarchical network wherein the hierarchical network is a hierarchical power supply network (PSN), wherein the second dimension of said outage state matrix (X) represents smart meters forming leaf nodes at the lowest hierarchical level of said power supply network (PSN), wherein each column of said two-dimensional outage state matrix (X) represents states of a respective smart meter at particular points in time, and wherein each row of said outage state matrix represents the states of all smart meters at a particular point in time. <EOS>