From 61ab9719065cba722f8742601a57dfe48fcc9fbb Mon Sep 17 00:00:00 2001
From: Volodymyr Vovchenko <vlvovch@gmail.com>
Date: Thu, 18 Apr 2019 17:28:44 +0200
Subject: [PATCH] Fixed bug in CCE with Vc<V event generator

---
 CMakeLists.txt                                       | 2 +-
 src/library/HRGEventGenerator/EventGeneratorBase.cpp | 6 +-----
 2 files changed, 2 insertions(+), 6 deletions(-)

diff --git a/CMakeLists.txt b/CMakeLists.txt
index f005c9c..d3209d5 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -8,7 +8,7 @@ project (ThermalFIST)
 # The version number.
 set (ThermalFIST_VERSION_MAJOR 1)
 set (ThermalFIST_VERSION_MINOR 1)
-set (ThermalFIST_VERSION_DEVEL 4)
+set (ThermalFIST_VERSION_DEVEL 5)
 
 # configure a header file to pass some of the CMake settings
 # to the source code
diff --git a/src/library/HRGEventGenerator/EventGeneratorBase.cpp b/src/library/HRGEventGenerator/EventGeneratorBase.cpp
index 4547d21..5e01d47 100644
--- a/src/library/HRGEventGenerator/EventGeneratorBase.cpp
+++ b/src/library/HRGEventGenerator/EventGeneratorBase.cpp
@@ -277,10 +277,6 @@ namespace thermalfist {
         m_AntiChargeMesons.push_back(std::make_pair(densities[i], i));
         m_MeanACM += densities[i];
       }
-      else if (m_THM->TPS()->Particles()[i].BaryonCharge() == 0 && m_THM->TPS()->Particles()[i].Strangeness() == 0 && m_THM->TPS()->Particles()[i].ElectricCharge() == 1) {
-        m_ChargeMesons.push_back(std::make_pair(densities[i], i));
-        m_MeanCM += densities[i];
-      }
       else if (m_THM->TPS()->Particles()[i].BaryonCharge() == 0 && m_THM->TPS()->Particles()[i].Strangeness() == 0 && m_THM->TPS()->Particles()[i].ElectricCharge() == 0 && m_THM->TPS()->Particles()[i].Charm() == 1) {
         m_CharmMesons.push_back(std::make_pair(densities[i], i));
         m_MeanCHRMM += densities[i];
@@ -698,7 +694,7 @@ namespace thermalfist {
           std::vector<int> totalsaux2(m_THM->TPS()->Particles().size(), 0);
           int netC = 0;
           for (size_t i = 0; i < totalsaux.size(); ++i) {
-            if (m_THM->TPS()->Particles()[i].Strangeness() > 0) {
+            if (m_THM->TPS()->Particles()[i].Charm() > 0) {
               for (int j = 0; j < totalsaux[i]; ++j) {
                 if (RandomGenerators::randgenMT.rand() < fraction) {
                   totalsaux2[i]++;