One of the unresolved mysteries in neutrino physics is the neutrino mass hierarchy. We present a new method to determine neutrino mass hierarchy by comparing the events of inverse beta decays (IBD), {bar{nu }}_e + prightarrow n + e^+, and neutral current (NC) interactions, nu ({overline{nu }}) + prightarrow nu ({overline{nu }}) + p, of supernova neutrinos from accretion and cooling phases in scintillation detectors. Supernova neutrino flavor conversions depend on the neutrino mass hierarchy. On account of Mikheyev–Smirnov–Wolfenstein effects, the full swap of the {bar{nu }}_e flux with the {bar{nu }}_x (x=mu ,~tau ) one occurs in the inverted hierarchy, while such a swap does not occur in the normal hierarchy. In consequence, the ratio of high energy IBD events to NC events for the inverted hierarchy is higher than in the normal hierarchy. Since the luminosity of {bar{nu }}_e is larger than that of nu _x in accretion phase while the luminosity of {bar{nu }}_e becomes smaller than that of nu _x in cooling phase, we calculate this ratio for both accretion and cooling phases. By analyzing the change of this event ratio from accretion phase to cooling phase, one can determine the neutrino mass hierarchy.