diff --git a/src/tools/NeighborList.cpp b/src/tools/NeighborList.cpp
index 00edab24fe8a83a47938d4e1b3110083e44967eb..55f2240a9c4b0eb21ab8074d83921e951ddd80fe 100644
--- a/src/tools/NeighborList.cpp
+++ b/src/tools/NeighborList.cpp
@@ -45,7 +45,7 @@ NeighborList::NeighborList(const vector<AtomNumber>& list0, const vector<AtomNum
   if(!do_pair) {
     nallpairs_=nlist0_*nlist1_;
   } else {
-    plumed_assert(nlist0_==nlist1_);
+    plumed_massert(nlist0_==nlist1_,"when using PAIR option, the two groups should have the same number of elements");
     nallpairs_=nlist0_;
   }
   initialize();