Untitled

ionSpeed = Vect();
    speedRoll_final = sqrt(threashold_rolling * 2.0 /
                           (n_ioni * param_maille * Fevap)); // threshold in eV
    bar = bar - maxCharged->m_coord;
    normbar = bar.length();
    ionSpeed = (bar)*speedRoll_final / normbar;