--- trunk/src/applications/staticProps/GofRAngle.cpp 2005/02/13 20:08:29 327 +++ trunk/src/applications/staticProps/GofRAngle.cpp 2005/02/13 20:36:24 328 @@ -100,13 +100,14 @@ void GofRAngle::collectHistogram(StuntDouble* sd1, Stu double distance = r12.length(); int whichRBin = distance / deltaR_; - - double cosAngle = evaluateAngle(sd1, sd2); - double halfBin = (nAngleBins_ - 1) * 0.5; - int whichThetaBin = halfBin * (cosAngle + 1.0); - ++histogram_[whichRBin][whichThetaBin]; - - ++npairs_; + if (distance <= len_) { + double cosAngle = evaluateAngle(sd1, sd2); + double halfBin = (nAngleBins_ - 1) * 0.5; + int whichThetaBin = halfBin * (cosAngle + 1.0); + ++histogram_[whichRBin][whichThetaBin]; + + ++npairs_; + } } void GofRAngle::writeRdf() {