Neighborhood_Criterion_Vector.h

#ifndef QNANO_NEW_NEIGHBORHOOD_CRITERION_VECTOR_DEFINED_H
#define QNANO_NEW_NEIGHBORHOOD_CRITERION_VECTOR_DEFINED_H

#include "structure/Neighborhood_Criterion.h"

class Neighborhood_Criterion_Vector: public Neighborhood_Criterion{
public: 
  Vec3d vec;
  double length;
  double tolerance;

private:  
  void calculate_upper_bound(){  length=vec.norm();  }
 
public:
  Neighborhood_Criterion_Vector(const Vec3d &vec_, double tol_=0.20) : vec(vec_), tolerance(tol_){
    calculate_upper_bound();
  }

  virtual bool is_neighbor(const Atom &a1, const Atom &a2) const;

  virtual double upper_bound_radius() const{ return length*(1+tolerance);}

};


#endif