| 328 | 
  | 
 | 
| 329 | 
  | 
  void SlabStatistics::processFrame(int istep) { | 
| 330 | 
  | 
    RealType z; | 
| 331 | 
+ | 
 | 
| 332 | 
  | 
    hmat_ = currentSnapshot_->getHmat(); | 
| 333 | 
  | 
    for (int i = 0; i < nBins_; i++) { | 
| 334 | 
  | 
      z = (((RealType)i + 0.5) / (RealType)nBins_) * hmat_(2,2); | 
| 354 | 
  | 
    SpatialStatistics(info, filename, sele, nbins), coordinateOrigin_(V3Zero) { | 
| 355 | 
  | 
     | 
| 356 | 
  | 
    binWidth_ = 1.0; | 
| 357 | 
+ | 
 | 
| 358 | 
+ | 
    Globals* simParams = info->getSimParams(); | 
| 359 | 
+ | 
    RNEMDParameters* rnemdParams = simParams->getRNEMDParameters(); | 
| 360 | 
+ | 
    bool hasCoordinateOrigin = rnemdParams->haveCoordinateOrigin(); | 
| 361 | 
  | 
     | 
| 362 | 
+ | 
    if (hasCoordinateOrigin) { | 
| 363 | 
+ | 
      coordinateOrigin_ = rnemdParams->getCoordinateOrigin(); | 
| 364 | 
+ | 
    } else { | 
| 365 | 
+ | 
      coordinateOrigin_ = V3Zero; | 
| 366 | 
+ | 
    } | 
| 367 | 
+ | 
     | 
| 368 | 
  | 
    r_ = new OutputData; | 
| 369 | 
  | 
    r_->units =  "Angstroms"; | 
| 370 | 
  | 
    r_->title =  "R"; | 
| 384 | 
  | 
  ShellStatistics::~ShellStatistics() { | 
| 385 | 
  | 
  } | 
| 386 | 
  | 
 | 
| 387 | 
< | 
  int ShellStatistics::getBin(Vector3d pos) {     | 
| 387 | 
> | 
  int ShellStatistics::getBin(Vector3d pos) {  | 
| 388 | 
  | 
    Vector3d rPos = pos - coordinateOrigin_; | 
| 389 | 
  | 
    return int(rPos.length() / binWidth_); | 
| 390 | 
  | 
  } |