Not sure whether modifying the pair.cpp was done correctly.
But whatever you are doing can done by introducing two types for single monomer and essentialy
each monomer will then have two types and you have define pair coeff accordingly. So if you have 100 monomers you will have 100 X2 types.
Was this helpful ?