massless after refresh corrected

marty/src/expander.cpp

@@ -308,7 +308,6 @@ Amplitude Expander::expandLagrangian()
             HEPAssert(csl::TestIndexSanity(ampl),
                     mty::error::RuntimeError,
                     "Bad index property for " + toString(ampl))
-            std::cout << ampl << std::endl;
             std::vector<Wilson> newWilsons 
                 = parseExpression(ampl, CSL_1, false);
             for (const auto &w : newWilsons) {

marty/src/modelBuilder.cpp

@@ -533,6 +533,9 @@ void ModelBuilder::rotateFields(
                 );
         for (int i = 0; i < nMassLessFields; ++i) {
             newFields[i]->setMass(CSL_0);
+            clearDependencies(L.mass, [&](Lagrangian::TermType const &t) {
+                return t->contains(newFields[i].get());
+            });
         }
     }
     applyUnitaryCondition(rotation);
@@ -586,6 +589,10 @@ void ModelBuilder::rotateFields(
     for (int i = 0; i < nMassLessFields; ++i) {
         newFields1[i]->setMass(CSL_0);
         newFields2[i]->setMass(CSL_0);
+        clearDependencies(L.mass, [&](Lagrangian::TermType const &t) {
+            return t->contains(newFields1[i].get())
+                || t->contains(newFields2[i].get());
+        });
     }
 }