Small tuning for PCE-HRG fits

vlvovch · Aug 25, 2022 · 2ffde41 · 2ffde41
1 parent 974838c
commit 2ffde41
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 2 deletions.
diff --git a/src/gui/QtThermalFIST/fittoexperimenttab.cpp b/src/gui/QtThermalFIST/fittoexperimenttab.cpp
@@ -908,7 +908,7 @@ void FitToExperimentTab::finalize() {
   if (model->UsePartialChemicalEquilibrium()) {
     tableParameters->setItem(cRow, 0, new QTableWidgetItem("Tkin (MeV)"));
     tableParameters->setItem(cRow, 1, new QTableWidgetItem(QString::number(result.Tkin.value * 1.e3)));
-    if (result.T.toFit == true)
+    if (result.Tkin.toFit == true)
       tableParameters->setItem(cRow, 2, new QTableWidgetItem(QString::number(result.Tkin.error * 1.e3)));
     else
       tableParameters->setItem(cRow, 2, new QTableWidgetItem("--"));

diff --git a/src/library/HRGFit/ThermalModelFitParameters.cpp b/src/library/HRGFit/ThermalModelFitParameters.cpp
@@ -49,7 +49,7 @@ namespace thermalfist {
     S = params.S;
     C = params.C;
 
-    Tkin = FitParameter("Tkin", false, params.T, 0.05, 0.02, 0.300);
+    Tkin = FitParameter("Tkin", false, params.T, 0.05, 0.050, 0.200);
 
     FillParameterList();
   }

diff --git a/src/library/HRGPCE/ThermalModelPCE.cpp b/src/library/HRGPCE/ThermalModelPCE.cpp
@@ -336,12 +336,16 @@ namespace thermalfist {
   {
     for (int ipart = 0; ipart < m_EffectiveCharges.size(); ++ipart) {
       ThermalParticle& part = m_model->TPS()->Particle(ipart);
+      part.CalculateAndSetDynamicalThreshold();
       if (part.ZeroWidthEnforced() || part.Statistics() != -1)
         continue;
 
       double totmu = 0.;
+      int nonzerocharges = 0;
       for (int ifeed = 0; ifeed < m_EffectiveCharges[ipart].size(); ++ifeed) {
         totmu += m_EffectiveCharges[ipart][ifeed] * m_model->TPS()->Particle(m_StableMapTo[ifeed]).Mass();
+        if (m_EffectiveCharges[ipart][ifeed] != 0.0)
+          nonzerocharges++;
       }
 
       if (totmu > part.DecayThresholdMassDynamical()) {