prevent double counting hits

interface/SimCluster.h

@@ -44,6 +44,8 @@ public:
       hits_.resize(nHits);
       fractions_.resize(nHits);
       for (unsigned int i = 0; i < nHits; ++i) {
+        if (std::find(hits_.begin(), hits_.end(), hits[i]) != hits_.end())
+          continue;
         hits_[i] = hits[i];
         fractions_[i] = fractions[i];
       }