Skip to content
Snippets Groups Projects
Commit e45ec176 authored by Bertram Kopf's avatar Bertram Kopf
Browse files

modified/fixed calculation in TMatrixDynamics

parent 7759c26d
No related branches found
No related tags found
No related merge requests found
......@@ -343,7 +343,7 @@ std::cout << "phpDescriptionVec.size(): " << phpDescriptionVec.size() << std::en
InfoMsg << "The k-Matrix input parameter are: " << endmsg;
_params->print(std::cout);
fillParams(_params);
if(_pathToFitParams != "") fillParams(_params);
_tMatr=std::shared_ptr<TMatrixRel>(new TMatrixRel(_kMatr));
......
......@@ -52,7 +52,7 @@ int main(int __argc,char *__argv[]){
int numStepsForSheetScan = 500;
std::string pathToConfigParser;
std::string pathToFitParams;
std::string pathToFitParams="";
std::vector<double> energyPlaneBorders;
energyPlaneBorders.resize(4);
......
......@@ -74,14 +74,15 @@ complex<double> TMatrixDynamics::eval(EvtData* theData, AbsXdecAmp* grandmaAmp,
_tMatr->evalMatrix(currentMass);
complex<double> currentTijRel=(*_tMatr)(_prodProjectionIndex,_decProjectionIndex);
// complex<double> SijRel=complex<double>(1.,0.)+2.*complex<double>(0.,1.)*sqrt(thePhpVecs[_prodProjectionIndex]->factor(currentMass)*thePhpVecs[_decProjectionIndex]->factor(currentMass))*currentTijRel;
complex<double> SijRel=complex<double>(1.,0.)+2.*complex<double>(0.,1.)*sqrt(thePhpVecs[_prodProjectionIndex]->factor(currentMass).real()*thePhpVecs[_decProjectionIndex]->factor(currentMass).real())*currentTijRel;
complex<double> SijRel=complex<double>(1.,0.)+2.*complex<double>(0.,1.)*sqrt(thePhpVecs[_prodProjectionIndex]->factor(currentMass)*thePhpVecs[_decProjectionIndex]->factor(currentMass))*currentTijRel;
// complex<double> SijRel=complex<double>(1.,0.)+2.*complex<double>(0.,1.)*sqrt(thePhpVecs[_prodProjectionIndex]->factor(currentMass).real()*thePhpVecs[_decProjectionIndex]->factor(currentMass).real())*currentTijRel;
theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::ETAFIT_PIPISCAT_NAME))=sqrt(norm(SijRel));
//phase
//note: this is a workaround
if(_prodProjectionIndex!=_decProjectionIndex){
SijRel=PawianConstants::i*2.*complex<double>(0.,1.)*sqrt(thePhpVecs[_prodProjectionIndex]->factor(currentMass)*thePhpVecs[_decProjectionIndex]->factor(currentMass))*currentTijRel;
theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::PHIFIT_PIPISCAT_NAME))=theData->DoubleId.at(IdStringMapRegistry::instance()->stringId(EvtDataScatteringList::PHI_PIPISCAT_NAME));
}
else{
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment