Let X be a ball quasi-Banach function space on R n . In this article, assuming that the powered Hardy-Littlewood maximal operator satisfies some Fefferman-Stein vector-valued maximal inequality on X as well as it is bounded on both the weak ball quasi-Banach function space WX and the associated space, the authors establish various Littlewood-Paley function characterizations of WH X (R n ) under some weak assumptions on the Littlewood-Paley functions. The authors also prove that the real interpolation intermediate space (H X (R n ), L ∞ (R n )) θ,∞ , between the Hardy space associated with X, H X (R n ), and the Lebesgue space L ∞ (R n ), is WH X 1/(1−θ) (R n ), where θ ∈ (0, 1). All these results are of wide applications. Particularly, when X := M p q (R n ) (the Morrey space), X := L p (R n ) (the mixed-norm Lebesgue space) and X := (E q Φ ) t (R n ) (the Orlicz-slice space), all these results are even new; when X := L Φ ω (R n ) (the weighted Orlicz space), the result on the real interpolation is new and, when X := L p(·) (R n ) (the variable Lebesgue space) and X := L Φ ω (R n ), the Littlewood-Paley function characterizations of WH X (R n ) obtained in this article improves the existing results via weakening the assumptions on the Littlewood-Paley functions.