| 93 |  | Vector3d pos1 = sd1->getPos(); | 
| 94 |  | Vector3d pos2 = sd2->getPos(); | 
| 95 |  | Vector3d r12 = pos1 - pos2; | 
| 96 | < | currentSnapshot_->wrapVector(r12); | 
| 96 | > | if (usePeriodicBoundaryConditions_) | 
| 97 | > | currentSnapshot_->wrapVector(r12); | 
| 98 |  | Vector3d dipole1 = sd1->getElectroFrame().getColumn(2); | 
| 99 |  | Vector3d dipole2 = sd2->getElectroFrame().getColumn(2); | 
| 100 |  |  | 
| 108 |  |  | 
| 109 |  | RealType halfBin = (nAngleBins_ - 1) * 0.5; | 
| 110 |  | int angleBin1 = halfBin * (cosAngle1 + 1.0); | 
| 111 | < | int angleBin2 = halfBin * (cosAngle1 + 1.0); | 
| 111 | > | int angleBin2 = halfBin * (cosAngle2 + 1.0); | 
| 112 |  |  | 
| 113 | < | ++histogram_[angleBin1][angleBin1]; | 
| 113 | > | ++histogram_[angleBin1][angleBin2]; | 
| 114 |  | ++npairs_; | 
| 115 |  | } | 
| 116 |  |  |