Neighborhood_Criterion_Lattice.h

#ifndef QNANO_NEW_NEIGHBORHOOD_CRITERION_LATTICE_DEFINED_H
#define QNANO_NEW_NEIGHBORHOOD_CRITERION_LATTICE_DEFINED_H

#include "structure/Neighborhood_Criterion.h"
#include "structure/Lattice.h"

class Neighborhood_Criterion_Lattice: public Neighborhood_Criterion{
public: 
  Lattice lat;
  double tolerance;
  double upper_bound;

private:  
  void calculate_upper_bound();
 
public:
  Neighborhood_Criterion_Lattice(const Lattice &lat_, double tol_=0.20) : lat(lat_), tolerance(tol_){
    calculate_upper_bound();
  }

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

  virtual double upper_bound_radius() const{ return upper_bound;}

};


#endif