From 9a3761ff09a18335d224d0e9962630e798e6f1f1 Mon Sep 17 00:00:00 2001 From: Volodymyr Vovchenko Date: Sat, 13 Jul 2019 21:54:20 +0200 Subject: [PATCH] Proxy electric charge susceptibility back to electric charge from net-pion, needed to preserve correct cpc4 result --- src/library/HRGBase/ThermalModelBase.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/library/HRGBase/ThermalModelBase.cpp b/src/library/HRGBase/ThermalModelBase.cpp index 5ddee0d..6160552 100644 --- a/src/library/HRGBase/ThermalModelBase.cpp +++ b/src/library/HRGBase/ThermalModelBase.cpp @@ -946,8 +946,8 @@ 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) { @@ -955,8 +955,8 @@ namespace thermalfist { 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];