| 70 | 
  | 
 | 
| 71 | 
  | 
void GofR::processHistogram() { | 
| 72 | 
  | 
 | 
| 73 | 
+ | 
    int nRealPairs = getNRealPairs(); | 
| 74 | 
  | 
    double volume = info_->getSnapshotManager()->getCurrentSnapshot()->getVolume(); | 
| 75 | 
< | 
    double pairDensity = npairs_ /volume; | 
| 75 | 
> | 
    double pairDensity = nRealPairs /volume; | 
| 76 | 
  | 
    double pairConstant = ( 4.0 * NumericConstant::PI * pairDensity ) / 3.0; | 
| 77 | 
  | 
 | 
| 78 | 
  | 
    for(int i = 0 ; i < histogram_.size(); ++i){ | 
| 100 | 
  | 
 | 
| 101 | 
  | 
    double distance = r12.length(); | 
| 102 | 
  | 
 | 
| 103 | 
< | 
    int whichBin = distance / deltaR_; | 
| 104 | 
< | 
    ++histogram_[whichBin]; | 
| 105 | 
< | 
    ++npairs_; | 
| 103 | 
> | 
    if (distance < len_) { | 
| 104 | 
> | 
        int whichBin = distance / deltaR_; | 
| 105 | 
> | 
        ++histogram_[whichBin]; | 
| 106 | 
> | 
        ++npairs_; | 
| 107 | 
> | 
    } | 
| 108 | 
  | 
} | 
| 109 | 
  | 
 | 
| 110 | 
  | 
 |