Statistics of Non-Markovianity in the optimal control of an open quantum system described by hierarchical equations of motion