void StatisticsSampler::sampleKineticEnergy(System *system)
{
    m_kineticEnergy = 0;
    Particles *particles = system->particles();
    float constant = 0.5*system->settings()->mass*system->settings()->atomsPerParticle;

    for(unsigned int i=0; i<particles->numberOfParticles(); i++) {
        m_kineticEnergy += constant*particles->velocitySquared(i);
    }
}
float StatisticsSampler::averageVelocity(System *system)
{
    Particles *particles = system->particles();
    double averageVelocity = 0;
    for(unsigned int i=0; i<particles->numberOfParticles(); i++) {
        averageVelocity += sqrt(particles->velocitySquared(i));
    }

    averageVelocity /= particles->numberOfParticles();
    return averageVelocity;
}