LocalCoordinatesSite

class LocalCoordinatesSite : public OpenMM::VirtualSite

This is a VirtualSite that uses the locations of several other particles to compute a local coordinate system, then places the virtual site at a fixed location in that coordinate system. The origin of the coordinate system and the directions of its x and y axes are each specified as a weighted sum of the locations of the other particles:

origin = wo1r1 + wo2r2 + wo3r3 + …

xdir = wx1r1 + wx2r2 + wx3r3 + …

ydir = wy1r1 + wy2r2 + wy3r3 + …

For the origin, the weights must add to one. For example if (wo1, wo2, wo3) = (1.0, 0.0, 0.0), the origin of the local coordinate system is at the location of particle 1. For xdir and ydir, the weights must add to zero. For example, if (wx1, wx2, wx3) = (-1.0, 0.5, 0.5), the x axis points from particle 1 toward the midpoint between particles 2 and 3.

The z direction is computed as zdir = xdir x ydir. To ensure the axes are all orthogonal, ydir is then recomputed as ydir = zdir x xdir. All three axis vectors are then normalized, and the virtual site location is set to

origin + x*xdir + y*ydir + z*zdir

Public Functions

LocalCoordinatesSite(const std::vector<int> &particles, const std::vector<double> &originWeights, const std::vector<double> &xWeights, const std::vector<double> &yWeights, const Vec3 &localPosition)

Create a new LocalCoordinatesSite virtual site.

Parameters
  • particles – the indices of the particle the site depends on

  • originWeights – the weight factors for the particles when computing the origin location

  • xWeights – the weight factors for the particles when computing xdir

  • yWeights – the weight factors for the particles when computing ydir

  • localPosition – the position of the virtual site in the local coordinate system

LocalCoordinatesSite(int particle1, int particle2, int particle3, const Vec3 &originWeights, const Vec3 &xWeights, const Vec3 &yWeights, const Vec3 &localPosition)

Create a new LocalCoordinatesSite virtual site. This constructor assumes the site depends on exactly three other particles.

Parameters
  • particle1 – the index of the first particle

  • particle2 – the index of the second particle

  • particle3 – the index of the third particle

  • originWeights – the weight factors for the three particles when computing the origin location

  • xWeights – the weight factors for the three particles when computing xdir

  • yWeights – the weight factors for the three particles when computing ydir

  • localPosition – the position of the virtual site in the local coordinate system

void getOriginWeights(std::vector<double> &weights) const

Get the weight factors for the particles when computing the origin location.

Vec3 getOriginWeights() const

Get the weight factors for the particles when computing the origin location. This version assumes the site depends on exactly three other particles, and throws an exception if that is not the case.

void getXWeights(std::vector<double> &weights) const

Get the weight factors for the particles when computing xdir.

Vec3 getXWeights() const

Get the weight factors for the particles when computing xdir. This version assumes the site depends on exactly three other particles, and throws an exception if that is not the case.

void getYWeights(std::vector<double> &weights) const

Get the weight factors for the particles when computing ydir.

Vec3 getYWeights() const

Get the weight factors for the particles when computing ydir. This version assumes the site depends on exactly three other particles, and throws an exception if that is not the case.

const Vec3 &getLocalPosition() const

Get the position of the virtual site in the local coordinate system.