| 46 | 
  | 
#include "primitives/Molecule.hpp" | 
| 47 | 
  | 
namespace oopse { | 
| 48 | 
  | 
 | 
| 49 | 
< | 
  RadialDistrFunc::        RadialDistrFunc(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2) | 
| 50 | 
< | 
    : info_(info), currentSnapshot_(NULL), dumpFilename_(filename), step_(1),  | 
| 49 | 
> | 
  RadialDistrFunc::RadialDistrFunc(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2) | 
| 50 | 
> | 
    : StaticAnalyser(info, filename),  | 
| 51 | 
  | 
      selectionScript1_(sele1), selectionScript2_(sele2), evaluator1_(info), evaluator2_(info),  | 
| 52 | 
  | 
      seleMan1_(info), seleMan2_(info), common_(info), sele1_minus_common_(info), sele2_minus_common_(info){ | 
| 53 | 
  | 
           | 
| 131 | 
  | 
      //stage 3 are completely overlapping | 
| 132 | 
  | 
 | 
| 133 | 
  | 
      if (evaluator1_.isDynamic() || evaluator2_.isDynamic()) { | 
| 134 | 
– | 
 | 
| 134 | 
  | 
        common_ = seleMan1_ & seleMan2_; | 
| 135 | 
  | 
        sele1_minus_common_ = seleMan1_ - common_; | 
| 136 | 
  | 
        sele2_minus_common_ = seleMan2_ - common_;             | 
| 138 | 
  | 
        int nSelected2 = seleMan2_.getSelectionCount(); | 
| 139 | 
  | 
        int nIntersect = common_.getSelectionCount(); | 
| 140 | 
  | 
             | 
| 141 | 
< | 
        nPairs_ = nSelected1 * nSelected2 - (nIntersect +1) * nIntersect/2;                       | 
| 141 | 
> | 
        nPairs_ = nSelected1 * nSelected2 - (nIntersect +1) * nIntersect/2; | 
| 142 | 
  | 
      } | 
| 144 | 
– | 
 | 
| 143 | 
  | 
      processNonOverlapping(sele1_minus_common_, seleMan2_); | 
| 144 | 
  | 
      processNonOverlapping(common_, sele2_minus_common_);         | 
| 145 | 
  | 
      processOverlapping(common_); | 
| 146 | 
< | 
         | 
| 149 | 
< | 
         | 
| 146 | 
> | 
       | 
| 147 | 
  | 
      processHistogram(); | 
| 148 | 
  | 
         | 
| 149 | 
  | 
    } |