A multi-Band Monte Carlo approach to Coulomb interaction for device analysis