| 46 |
|
|
| 47 |
|
namespace oopse { |
| 48 |
|
|
| 49 |
< |
GofXyz::GofXyz(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2) |
| 50 |
< |
: RadialDistrFunc(info, filename, sele1, sele2){ |
| 49 |
> |
GofXyz::GofXyz(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2, double len, int nrbins) |
| 50 |
> |
: RadialDistrFunc(info, filename, sele1, sele2), len_(len), nRBins_(nrbins) { |
| 51 |
|
setOutputName(getPrefix(filename) + ".gxyz"); |
| 52 |
|
|
| 53 |
+ |
deltaR_ = len_ / nRBins_; |
| 54 |
+ |
|
| 55 |
+ |
histogram_.resize(nRBins_); |
| 56 |
+ |
for (int i = 0 ; i < nRBins_; ++i) { |
| 57 |
+ |
histogram_[i].resize(nRBins_); |
| 58 |
+ |
for(int j = 0; j < nRBins_; ++j) { |
| 59 |
+ |
histogram_[i][j].resize(nRBins_); |
| 60 |
+ |
} |
| 61 |
+ |
} |
| 62 |
|
} |
| 63 |
|
|
| 64 |
|
|
| 147 |
|
|
| 148 |
|
} else { |
| 149 |
|
|
| 150 |
< |
|
| 150 |
> |
sprintf(painCave.errMsg, "GofXyz: unable to open %s\n", outputFilename_.c_str()); |
| 151 |
> |
painCave.isFatal = 1; |
| 152 |
> |
simError(); |
| 153 |
|
} |
| 154 |
|
|
| 155 |
|
rdfStream.close(); |