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

added tensor fit to pbar p benchmark channel

parent 95771d5e
No related branches found
No related tags found
No related merge requests found
......@@ -144,17 +144,21 @@ void IsobarTensorDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fs
motherTensor4Vec.Boost(all4Vec);
Vector4<double> daughter2HelMother(0.,0.,0.,0.);
if(_hasMotherPart){
if(fabs(mother_4Vec==all4Vec)){
daughter2HelMother=daughter2_4Vec;
daughter2HelMother.Boost(daughter2HelMother); //is this correct????
}
else daughter2HelMother=helicityVec(prodParticle4Vec, mother_4Vec, daughter2_4Vec);
}
else{
daughter2HelMother=daughter2_4Vec;
// if(_hasMotherPart){
// if(mother_4Vec==prodParticle4Vec){
// daughter2HelMother=daughter2_4Vec;
// daughter2HelMother.Boost(daughter2HelMother); //is this correct????
// }
// else{
// daughter2HelMother=helicityVec(prodParticle4Vec, mother_4Vec, daughter2_4Vec);
// }
// }
// else{
daughter2HelMother=daughter2_4Vec; //only needed for prod amps
daughter2HelMother.Boost(mother_4Vec);
}
// }
//rotate everything into the cms flight direction
// Vector4<double> allRot4Vec=all4Vec;
......
......@@ -148,11 +148,12 @@ complex<double> TensorDecAmps::lsLoop(AbsXdecAmp* grandmaAmp, Spin lamX, EvtData
tmpResult+=amp;
}
}
if (_absDyn->isLdependent()) tmpResult*=_cachedDynLSMap.at(std::this_thread::get_id()).at((*it)->L);
if (_absDyn->isLdependent()) tmpResult*=_cachedDynLSMap.at(std::this_thread::get_id()).at((*it)->L);
result+=tmpResult;
}
if (!_absDyn->isLdependent()) result *=_cachedDynMap.at(std::this_thread::get_id()).at(_absDyn->grandMaKey(grandmaAmp));
result*=_isospinCG;
return result;
}
......@@ -214,7 +215,7 @@ void TensorDecAmps::fillDefaultParams(std::shared_ptr<AbsPawianParameters> fitP
}
void TensorDecAmps::fillParamNameList(){
_paramNameList.clear();
std::vector< std::shared_ptr<const LScomb> >::const_iterator itLS;
for(itLS=_LSs.begin(); itLS!=_LSs.end(); ++itLS){
std::string magName=(*itLS)->name()+_key+"Mag";
......
......@@ -67,6 +67,11 @@ double pbarpTensorLh::calcEvtIntensity(EvtData* theData, std::shared_ptr<AbsPawi
double result=0.;
std::vector< std::shared_ptr<AbsXdecAmp> >::iterator itDecAll;
for (itDecAll=_decAmps.begin(); itDecAll!=_decAmps.end(); ++itDecAll){
(*itDecAll)->calcDynamics(theData);
}
std::map <std::shared_ptr<const JPCLS>, std::vector< std::shared_ptr<AbsXdecAmp> >, pawian::Collection::SharedPtrLess >::iterator it;
Spin lamSteps=1;
......
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