| 47 |  | namespace oopse { | 
| 48 |  |  | 
| 49 |  | GofR::GofR(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2, double len) | 
| 50 | < | : RadialDistrFunc(info, filename, sele1, sele2, len){ | 
| 50 | > | : RadialDistrFunc(info, filename, sele1, sele2){ | 
| 51 |  |  | 
| 52 | < | histogram_.resize(nbins_); | 
| 53 | < | avgGofr_.resize(nbins_); | 
| 52 | > | deltaR_ = len_ /nRBins_; | 
| 53 | > |  | 
| 54 | > | histogram_.resize(nRBins_); | 
| 55 | > | avgGofr_.resize(nRBins_); | 
| 56 |  | } | 
| 57 |  |  | 
| 58 |  |  | 
| 59 |  | void GofR::preProcess() { | 
| 58 | – | avgGofr_.resize(nbins_); | 
| 60 |  | std::fill(avgGofr_.begin(), avgGofr_.end(), 0.0); | 
| 61 |  | } | 
| 62 |  |  | 
| 98 |  | double distance = r12.length(); | 
| 99 |  |  | 
| 100 |  | int whichBin = distance / deltaR_; | 
| 101 | < | histogram_[whichBin] ++; | 
| 102 | < | npairs_++; | 
| 101 | > | ++histogram_[whichBin]; | 
| 102 | > | ++npairs_; | 
| 103 |  | } | 
| 104 |  |  | 
| 105 |  |  |