Proxy electric charge susceptibility back to electric charge from net…

…-pion, needed to preserve correct cpc4 result
vlvovch · Jul 13, 2019 · 9a3761f · 9a3761f
1 parent 941bf69
commit 9a3761f
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/src/library/HRGBase/ThermalModelBase.cpp b/src/library/HRGBase/ThermalModelBase.cpp
@@ -946,17 +946,17 @@ namespace thermalfist {
             int c1 = 0;
             //if (i == 0) c1 = m_TPS->Particles()[k].BaryonCharge();
             if (i == 0) c1 = 1 * (m_TPS->Particles()[k].PdgId() == 2212) - 1 * (m_TPS->Particles()[k].PdgId() == -2212);
-            //if (i == 1) c1 = m_TPS->Particles()[k].ElectricCharge();
-            if (i == 1) c1 = 1 * (m_TPS->Particles()[k].PdgId() == 211) - 1 * (m_TPS->Particles()[k].PdgId() == -211);
+            if (i == 1) c1 = m_TPS->Particles()[k].ElectricCharge();
+            //if (i == 1) c1 = 1 * (m_TPS->Particles()[k].PdgId() == 211) - 1 * (m_TPS->Particles()[k].PdgId() == -211);
             if (i == 2) c1 = 1 * (m_TPS->Particles()[k].PdgId() == 321) - 1 * (m_TPS->Particles()[k].PdgId() == -321);
             if (i == 3) c1 = m_TPS->Particles()[k].Charm();
             for (size_t kp = 0; kp < m_TotalCorrel.size(); ++kp) {
               if (m_TPS->Particles()[kp].IsStable()) {
                 int c2 = 0;
                 //if (j == 0) c2 = m_TPS->Particles()[kp].BaryonCharge();
                 if (j == 0) c2 = 1 * (m_TPS->Particles()[kp].PdgId() == 2212) - 1 * (m_TPS->Particles()[kp].PdgId() == -2212);
-                //if (j == 1) c2 = m_TPS->Particles()[kp].ElectricCharge();
-                if (j == 1) c2 = 1 * (m_TPS->Particles()[kp].PdgId() == 211) - 1 * (m_TPS->Particles()[kp].PdgId() == -211);
+                if (j == 1) c2 = m_TPS->Particles()[kp].ElectricCharge();
+                //if (j == 1) c2 = 1 * (m_TPS->Particles()[kp].PdgId() == 211) - 1 * (m_TPS->Particles()[kp].PdgId() == -211);
                 if (j == 2) c2 = 1 * (m_TPS->Particles()[kp].PdgId() == 321) - 1 * (m_TPS->Particles()[kp].PdgId() == -321);
                 if (j == 3) c2 = m_TPS->Particles()[kp].Charm();
                 m_ProxySusc[i][j] += c1 * c2 * m_TotalCorrel[k][kp];