added protection against calculating scattering properties (T-, S-matrix) below thresholds

Examples/benchmark/coupledChannel/global.cfg

@@ -8,7 +8,7 @@ minimumTolerance = 10
 
 pdgTableFile = /Particle/pdtNew.table
 
-serverPort = 38854
+serverPort = 20301
 serverAddress = localhost
 noOfClients = 4
 noOfThreads = 1

PwaUtils/TMatrixDynamics.cc

@@ -96,8 +96,14 @@ complex<double> TMatrixDynamics::eval(EvtData* theData, AbsXdecAmp* grandmaAmp,
     while( (deltaRel-phiData) > 90.) deltaRel-=180.;
     theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::PHIFIT_PIPISCAT_NAME))=deltaRel;
   }
-
-  theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::ETAFIT_PIPISCAT_NAME))=sqrt(norm(SijRel));
+  if(thePhpVecs[_decProjectionIndex]->factor(currentMass).real()<1.e-12){//protection against calculations below threshold
+    theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::ETAFIT_PIPISCAT_NAME))=0.;
+    theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::PHIFIT_PIPISCAT_NAME))=0.;
+  }
+  else{
+    double normSijRel=sqrt(norm(SijRel));
+    theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::ETAFIT_PIPISCAT_NAME))=sqrt(norm(SijRel));
+  }
   return (*_tMatr)(_prodProjectionIndex,_decProjectionIndex);
 }
 

pipiScatteringUtils/pipiScatteringBaseLh.cc

@@ -112,9 +112,11 @@ double pipiScatteringBaseLh::calcEvtIntensity(EvtData* theData, std::shared_ptr<
 
   double result=0.5*(chi2_phi+chi2_eta);
 // if(chi2_phi>1000.){
-//   InfoMsg << "chi2_phi: " << chi2_phi << "\tchi2_eta: " << chi2_eta << endmsg;
-//   InfoMsg << "phiData: " << phiData << "\tphiFit: " << phiFit << "\tphiErrData: " << phiErrData 
-//           << "\netaData: " << etaData << "\tetaFit: " << etaFit << endmsg;
+  // InfoMsg << "etaData: " << etaData << "\tetaFit: " << etaFit << endmsg;
+  // InfoMsg << "etaErrData: " << etaErrData << endmsg;
+  //  InfoMsg << "chi2_phi: " << chi2_phi << "\tchi2_eta: " << chi2_eta << endmsg;
+  //  InfoMsg << "phiData: " << phiData << "\tphiFit: " << phiFit << "\tphiErrData: " << phiErrData 
+  //          << "\netaData: " << etaData << "\tetaFit: " << etaFit << endmsg;
 // }
    return result;