| 46 |  |  | 
| 47 |  | namespace oopse { | 
| 48 |  |  | 
| 49 | < | GofAngle2::GofAngle2(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2) | 
| 50 | < | : RadialDistrFunc(info, filename, sele1, sele2){ | 
| 49 | > | GofAngle2::GofAngle2(SimInfo* info, const std::string& filename, const std::string& sele1, | 
| 50 | > | const std::string& sele2, int nangleBins) | 
| 51 | > | : RadialDistrFunc(info, filename, sele1, sele2), nAngleBins_(nangleBins) { | 
| 52 | > |  | 
| 53 |  | setOutputName(getPrefix(filename) + ".gto"); | 
| 54 |  |  | 
| 55 | + | deltaCosAngle_ = 2.0 / nAngleBins_; | 
| 56 | + |  | 
| 57 | + | histogram_.resize(nAngleBins_); | 
| 58 | + | avgGofr_.resize(nAngleBins_); | 
| 59 | + | for (int i = 0 ; i < nAngleBins_; ++i) { | 
| 60 | + | histogram_[i].resize(nAngleBins_); | 
| 61 | + | avgGofr_[i].resize(nAngleBins_); | 
| 62 | + | } | 
| 63 | + |  | 
| 64 |  | } | 
| 65 |  |  | 
| 66 |  |  | 
| 124 |  | double cosAngle1 = -1.0 + (i + 0.5)*deltaCosAngle_; | 
| 125 |  |  | 
| 126 |  | for(int j = 0; j < avgGofr_[i].size(); ++j) { | 
| 127 | < | double cosAngle2 = -1.0 + (i + 0.5)*deltaCosAngle_; | 
| 127 | > | double cosAngle2 = -1.0 + (j + 0.5)*deltaCosAngle_; | 
| 128 |  | rdfStream << cosAngle1 << "\t" << cosAngle2 << "\t" << avgGofr_[i][j]/nProcessed_ << "\n"; | 
| 129 |  | } | 
| 130 |  | } |