MSSM right gaugino interaction fixed

marty/models/include/pmssm_lem.h

@@ -213,6 +213,10 @@ protected:
         void initInteractions164();
         void initInteractions165();
         void initInteractions166();
+        void initInteractions167();
+        void initInteractions168();
+        void initInteractions169();
+        void initInteractions170();
     void initSpectrum();
 
 protected:

marty/models/src/pmssm_content.cpp

@@ -923,6 +923,10 @@ void PMSSM_LEM::initInteractions()
     initInteractions164();
     initInteractions165();
     initInteractions166();
+    initInteractions167();
+    initInteractions168();
+    initInteractions169();
+    initInteractions170();
 }
 
 void PMSSM_LEM::initSpectrum()

marty/src/MSSM.cpp

@@ -400,11 +400,15 @@ void MSSM_Model::initGauginoInteractions(
                 break;
             }
 
+        bool isLeft = (fermion->getChirality() == Chirality::Left);
+        csl::Expr scal = (isLeft) ? cc(scalar(s_index)) : scalar(s_index);
+        csl::Expr ferm = (isLeft) ? fermion(f_index) : cc(fermion(f_index));
+        csl::Expr connexion = (isLeft) ? 
+            dirac4.C_matrix({a, b}) : -dirac4.getDelta()({a, b});
         addLagrangianTerm(
-                - csl::sqrt_s(2) * coupling 
-                * cc(scalar(s_index)) * T({A, i, j}) * fermion(f_index)
-                * dirac4.C_matrix({a, b})
-                * gaugino({A, b}),
+                csl::DeepRefreshed(csl::prod_s({
+                    -csl::sqrt_s(2), coupling, T({A, i, j}),
+                    scal, ferm, connexion, gaugino({A, b})})),
                 true
                 );
     }
@@ -434,11 +438,15 @@ void MSSM_Model::initU1GauginoInteractions(
         csl::Index b = a.rename();
         s_index.erase(s_index.end() - 1); // Dirac index
 
+        bool isLeft = (fermion->getChirality() == Chirality::Left);
+        csl::Expr scal = (isLeft) ? cc(scalar(s_index)) : scalar(s_index);
+        csl::Expr ferm = (isLeft) ? fermion(f_index) : cc(fermion(f_index));
+        csl::Expr connexion = (isLeft) ? 
+            dirac4.C_matrix({a, b}) : -dirac4.getDelta()({a, b});
         addLagrangianTerm(
-                -csl::sqrt_s(2) * coupling 
-                * charge * cc(scalar(s_index)) * fermion(f_index)
-                * dirac4.C_matrix({a, b})
-                * gaugino(b),
+                csl::DeepRefreshed(csl::prod_s({
+                    -csl::sqrt_s(2), coupling, charge,
+                    scal, ferm, connexion, gaugino(b)})),
                 true
                 );
     }

marty/src/PMSSM.cpp

@@ -35,13 +35,13 @@ PMSSM_Model::PMSSM_Model(
 {
     addAllowedMixing({"st_L", "st_R"});
     addAllowedMixing({"sb_L", "sb_R"});
-    addAllowedMixing({"sc_L", "sc_R"});
-    addAllowedMixing({"ss_L", "ss_R"});
-    addAllowedMixing({"su_L", "su_R"});
-    addAllowedMixing({"sd_L", "sd_R"});
+    // addAllowedMixing({"sc_L", "sc_R"});
+    // addAllowedMixing({"ss_L", "ss_R"});
+    // addAllowedMixing({"su_L", "su_R"});
+    // addAllowedMixing({"sd_L", "sd_R"});
     addAllowedMixing({"stau_L", "stau_R"});
-    addAllowedMixing({"smu_L", "smu_R"});
-    addAllowedMixing({"se_L", "se_R"});
+    // addAllowedMixing({"smu_L", "smu_R"});
+    // addAllowedMixing({"se_L", "se_R"});
     if (init) {
         std::cout << "Initializing gauge and particle content ..." << std::endl;
         initContent();

marty/src/interactionTerm.cpp

@@ -445,6 +445,8 @@ csl::Expr InteractionTerm::getFullExpression() const
     csl::Expr res = applyFactorAndSymmetriesOn(getFieldProduct());
     res = DeepRefreshed(res);
     if (csl::Abbrev::getFreeStructure(res->getIndexStructure()).size() > 0) {
+        std::cerr << "Term : " << getTerm() << '\n';
+        std::cerr << "Res  : " << res << '\n';
         CallHEPError(mty::error::IndexError, "")
     }
     return res;