We give an interpretation of the t=1 specialization of the modified Macdonald
polynomial as a generating function of the energy statistics defined on the set
of paths arising in the context of Box-Ball Systems (BBS-paths for short). We
also introduce one parameter generalization of the energy statistics on the set
of BBS-paths which all, conjecturally, have the same distribution.