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

added simple Blatt-Weisskopf factors w/o any mass shape for dynamics; qR value...

added simple Blatt-Weisskopf factors w/o any mass shape for dynamics; qR value as argument for BreitWignerBlattWRelDynamics
parent d1e47cfd
No related branches found
No related tags found
No related merge requests found
Showing
with 58 additions and 14 deletions
...@@ -22,7 +22,7 @@ ...@@ -22,7 +22,7 @@
//************************************************************************// //************************************************************************//
#include "PwaDynamics/BreitWignerFunction.hh" #include "PwaDynamics/BreitWignerFunction.hh"
#include "PwaDynamics/BarrierFactor.hh" //#include "PwaDynamics/BarrierFactor.hh"
complex<double> BreitWignerFunction::NonRel(double currentMass,double mass0, double width){ complex<double> BreitWignerFunction::NonRel(double currentMass,double mass0, double width){
complex<double> i(0.,1.); complex<double> i(0.,1.);
...@@ -35,14 +35,14 @@ complex<double> BreitWignerFunction::Rel(double currentMass,double mass0, double ...@@ -35,14 +35,14 @@ complex<double> BreitWignerFunction::Rel(double currentMass,double mass0, double
return mass0*width/(mass0*mass0 - currentMass*currentMass- i*rho0_m*mass0*width); return mass0*width/(mass0*mass0 - currentMass*currentMass- i*rho0_m*mass0*width);
} }
complex<double> BreitWignerFunction::BlattWRel(int orbMom, double currentMass,double mass0, double width, double massA, double massB){ complex<double> BreitWignerFunction::BlattWRel(int orbMom, double currentMass,double mass0, double width, double massA, double massB, double qR){
complex<double> i(0.,1.); complex<double> i(0.,1.);
complex<double> rho0=phaseSpaceFac(mass0, massA, massB); complex<double> rho0=phaseSpaceFac(mass0, massA, massB);
complex<double> rho=phaseSpaceFac(currentMass, massA, massB); complex<double> rho=phaseSpaceFac(currentMass, massA, massB);
complex<double> momQ0=breakupMomQ(mass0, massA, massB); complex<double> momQ0=breakupMomQ(mass0, massA, massB);
complex<double> momQ=breakupMomQ(currentMass, massA, massB); complex<double> momQ=breakupMomQ(currentMass, massA, massB);
std::complex<double> bwkFactor=BarrierFactor::BlattWeisskopfRatio(orbMom, momQ, momQ0, BarrierFactor::qRDefault); std::complex<double> bwkFactor=BarrierFactor::BlattWeisskopfRatio(orbMom, momQ, momQ0, qR);
return mass0*width*bwkFactor/(mass0*mass0 - currentMass*currentMass- i*(rho/rho0)*mass0*width*bwkFactor*bwkFactor); return mass0*width*bwkFactor/(mass0*mass0 - currentMass*currentMass- i*(rho/rho0)*mass0*width*bwkFactor*bwkFactor);
......
...@@ -35,6 +35,7 @@ ...@@ -35,6 +35,7 @@
#include <memory> #include <memory>
#include "qft++/relativistic-quantum-mechanics/Utils.hh" #include "qft++/relativistic-quantum-mechanics/Utils.hh"
#include "PwaDynamics/BarrierFactor.hh"
//_____________________________________________________________________________ //_____________________________________________________________________________
//_____________________________________________________________________________ //_____________________________________________________________________________
...@@ -42,7 +43,7 @@ ...@@ -42,7 +43,7 @@
namespace BreitWignerFunction { namespace BreitWignerFunction {
complex<double> NonRel(double currentMass,double mass0, double width); complex<double> NonRel(double currentMass,double mass0, double width);
complex<double> Rel(double currentMass,double mass0, double width, double massA, double massB); complex<double> Rel(double currentMass,double mass0, double width, double massA, double massB);
complex<double> BlattWRel(int orbMom, double currentMass,double mass0, double width, double massA, double massB); complex<double> BlattWRel(int orbMom, double currentMass,double mass0, double width, double massA, double massB, double qR=BarrierFactor::qRDefault);
}; // namespace BreitWignerFunction }; // namespace BreitWignerFunction
......
...@@ -26,6 +26,7 @@ ...@@ -26,6 +26,7 @@
#include <getopt.h> #include <getopt.h>
#include <fstream> #include <fstream>
#include <sstream>
#include <algorithm> #include <algorithm>
#include "PwaUtils/AbsDecay.hh" #include "PwaUtils/AbsDecay.hh"
...@@ -40,6 +41,7 @@ ...@@ -40,6 +41,7 @@
#include "PwaUtils/DynRegistry.hh" #include "PwaUtils/DynRegistry.hh"
#include "PwaUtils/AbsDynamics.hh" #include "PwaUtils/AbsDynamics.hh"
#include "PwaUtils/GlobalEnv.hh" #include "PwaUtils/GlobalEnv.hh"
#include "PwaDynamics/BarrierFactor.hh"
#include "ErrLogger/ErrLogger.hh" #include "ErrLogger/ErrLogger.hh"
AbsDecay::AbsDecay(Particle* mother, Particle* daughter1, Particle* daughter2, ChannelID channelId) : AbsDecay::AbsDecay(Particle* mother, Particle* daughter1, Particle* daughter2, ChannelID channelId) :
...@@ -55,6 +57,7 @@ AbsDecay::AbsDecay(Particle* mother, Particle* daughter1, Particle* daughter2, C ...@@ -55,6 +57,7 @@ AbsDecay::AbsDecay(Particle* mother, Particle* daughter1, Particle* daughter2, C
,_daughter1IGJPCPtr(getIGJPCPtr(daughter1)) ,_daughter1IGJPCPtr(getIGJPCPtr(daughter1))
,_daughter2IGJPCPtr(getIGJPCPtr(daughter2)) ,_daughter2IGJPCPtr(getIGJPCPtr(daughter2))
,_isospinClebschG(1.) ,_isospinClebschG(1.)
,_qR(BarrierFactor::qRDefault)
,_name(mother->name()+"To"+daughter1->name()+"_"+daughter2->name()) ,_name(mother->name()+"To"+daughter1->name()+"_"+daughter2->name())
,_fitParamSuffix(_name) ,_fitParamSuffix(_name)
,_massParamKey(_mother->name()) ,_massParamKey(_mother->name())
...@@ -235,6 +238,17 @@ void AbsDecay::enableDynamics(std::string& dynString, std::vector<std::string>& ...@@ -235,6 +238,17 @@ void AbsDecay::enableDynamics(std::string& dynString, std::vector<std::string>&
} }
_decPair2ndChannel=make_pair(firstParticle,secondParticle); _decPair2ndChannel=make_pair(firstParticle,secondParticle);
} }
else if(_dynType=="BlattWBarrier" || _dynType=="BreitWignerBlattWRel"){
if(additionalStringVec.size()>0){
std::istringstream currentqRIstr(additionalStringVec[0]);
_qR=stof(additionalStringVec[0]);
if ( _qR <1.e-5 || _qR > 20.){
Alert << "radius for barrier factor too high or too low qr=" << _qR << endmsg;
exit(0);
}
}
Info << "AmpName: " << name() << " radius for barrier factor qr= " << _qR << endmsg;
}
_absDynPtr=DynRegistry::instance()->getDynamics(shared_from_this()); _absDynPtr=DynRegistry::instance()->getDynamics(shared_from_this());
} }
...@@ -338,12 +352,6 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect ...@@ -338,12 +352,6 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect
// _absDecDaughter2->fillWignerDs(fsMap, daughter2_4Vec, evtData); // _absDecDaughter2->fillWignerDs(fsMap, daughter2_4Vec, evtData);
_absDecDaughter2->fillWignerDs(fsMap, mother4Vec, evtData, _refKey); _absDecDaughter2->fillWignerDs(fsMap, mother4Vec, evtData, _refKey);
} }
if (_useProdBarrier){
double qVal=daughter2HelMother.P();
double qValNorm=breakupMomQ(mother4Vec.M(), massSumFsParticlesDec1(), massSumFsParticlesDec2()).real();
evtData->DoubleString[_wignerDKey]=qVal;
evtData->DoubleString[_wignerDKey+"qNorm"] = qValNorm;
}
} }
else{ else{
if (!_daughter1IsStable){ if (!_daughter1IsStable){
...@@ -355,6 +363,19 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect ...@@ -355,6 +363,19 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect
_absDecDaughter2->fillWignerDs(fsMap, mother4Vec, evtData, _refKey); _absDecDaughter2->fillWignerDs(fsMap, mother4Vec, evtData, _refKey);
} }
} }
bool fillqVals=false;
if(_isProdAmp && _useProdBarrier) fillqVals=true;
else if(0!=_absDynPtr){
if(_absDynPtr->type()=="BlattWBarrierDynamics") fillqVals=true;
}
if(fillqVals){
double qVal=daughter2HelMother.P();
double qValNorm=breakupMomQ(mother4Vec.M(), massSumFsParticlesDec1(), massSumFsParticlesDec2()).real();
evtData->DoubleString[_wignerDKey]=qVal;
evtData->DoubleString[_wignerDKey+"qNorm"] = qValNorm;
}
// _alreadyFilledMap[evtNo]=true; // _alreadyFilledMap[evtNo]=true;
} }
......
...@@ -73,6 +73,7 @@ public: ...@@ -73,6 +73,7 @@ public:
virtual void fillWignerDs(std::map<std::string , Vector4<double> >& fsMap, Vector4<double>& prodParticle4Vec, EvtData* evtData, std::string& refKey); virtual void fillWignerDs(std::map<std::string , Vector4<double> >& fsMap, Vector4<double>& prodParticle4Vec, EvtData* evtData, std::string& refKey);
void enableDynamics(std::string& dynString, std::vector<std::string>& additionalStringVec); void enableDynamics(std::string& dynString, std::vector<std::string>& additionalStringVec);
std::shared_ptr<AbsDynamics> getDynamics(){return _absDynPtr;} std::shared_ptr<AbsDynamics> getDynamics(){return _absDynPtr;}
double barrierqR() {return _qR;}
virtual void print(std::ostream& os) const; virtual void print(std::ostream& os) const;
Particle* motherPart() {return _mother;} Particle* motherPart() {return _mother;}
...@@ -127,6 +128,7 @@ protected: ...@@ -127,6 +128,7 @@ protected:
std::shared_ptr<const IGJPC> _daughter2IGJPCPtr; std::shared_ptr<const IGJPC> _daughter2IGJPCPtr;
double _isospinClebschG; double _isospinClebschG;
double _qR;
std::string _name; std::string _name;
std::string _fitParamSuffix; std::string _fitParamSuffix;
......
...@@ -49,6 +49,7 @@ public: ...@@ -49,6 +49,7 @@ public:
virtual ~AbsDynamics(); virtual ~AbsDynamics();
virtual std::string name() {return _name;}; virtual std::string name() {return _name;};
virtual std::string type()=0;
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0)=0; virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0)=0;
virtual void cacheAmplitudes(); virtual void cacheAmplitudes();
virtual void fillMasses(EvtData* theData); virtual void fillMasses(EvtData* theData);
......
...@@ -34,8 +34,9 @@ ...@@ -34,8 +34,9 @@
#include "Particle/Particle.hh" #include "Particle/Particle.hh"
#include "PwaDynamics/BreitWignerFunction.hh" #include "PwaDynamics/BreitWignerFunction.hh"
BreitWignerBlattWRelDynamics::BreitWignerBlattWRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2) : BreitWignerBlattWRelDynamics::BreitWignerBlattWRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2, double qR) :
BreitWignerRelDynamics(name, fsParticles, mother, massSumDaughter1, massSumDaughter2) BreitWignerRelDynamics(name, fsParticles, mother, massSumDaughter1, massSumDaughter2)
,_qR(qR)
{ {
} }
...@@ -50,7 +51,7 @@ complex<double> BreitWignerBlattWRelDynamics::eval(EvtData* theData, AbsXdecAmp* ...@@ -50,7 +51,7 @@ complex<double> BreitWignerBlattWRelDynamics::eval(EvtData* theData, AbsXdecAmp*
return _cachedLMap.at(evtNo).at(orbMom); return _cachedLMap.at(evtNo).at(orbMom);
} }
complex<double> result=BreitWignerFunction::BlattWRel(orbMom, theData->DoubleString.at(_dynKey), _currentMass, _currentWidth, _fsp1Mass, _fsp2Mass); complex<double> result=BreitWignerFunction::BlattWRel(orbMom, theData->DoubleString.at(_dynKey), _currentMass, _currentWidth, _fsp1Mass, _fsp2Mass, _qR);
if ( _cacheAmps){ if ( _cacheAmps){
theMutex.lock(); theMutex.lock();
_cachedLMap[evtNo][orbMom]=result; _cachedLMap[evtNo][orbMom]=result;
......
...@@ -34,17 +34,20 @@ ...@@ -34,17 +34,20 @@
#include <memory> #include <memory>
#include "PwaUtils/BreitWignerRelDynamics.hh" #include "PwaUtils/BreitWignerRelDynamics.hh"
#include "PwaDynamics/BarrierFactor.hh"
class BreitWignerBlattWRelDynamics : public BreitWignerRelDynamics{ class BreitWignerBlattWRelDynamics : public BreitWignerRelDynamics{
public: public:
BreitWignerBlattWRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2); BreitWignerBlattWRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2, double qR=BarrierFactor::qRDefault);
virtual ~BreitWignerBlattWRelDynamics(); virtual ~BreitWignerBlattWRelDynamics();
virtual std::string type() {return "BreitWignerBlattWRelDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
protected: protected:
std::map<int, std::map<int, complex<float> > > _cachedLMap; std::map<int, std::map<int, complex<float> > > _cachedLMap;
double _qR;
private: private:
......
...@@ -41,6 +41,7 @@ public: ...@@ -41,6 +41,7 @@ public:
BreitWignerDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother); BreitWignerDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~BreitWignerDynamics(); virtual ~BreitWignerDynamics();
virtual std::string type() {return "BreitWignerDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -41,6 +41,7 @@ public: ...@@ -41,6 +41,7 @@ public:
BreitWignerRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2); BreitWignerRelDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, double massSumDaughter1, double massSumDaughter2);
virtual ~BreitWignerRelDynamics(); virtual ~BreitWignerRelDynamics();
virtual std::string type() {return "BreitWignerRelDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
......
...@@ -41,6 +41,7 @@ ...@@ -41,6 +41,7 @@
#include "PwaUtils/KMatrixDynamics.hh" #include "PwaUtils/KMatrixDynamics.hh"
#include "PwaUtils/VoigtDynamics.hh" #include "PwaUtils/VoigtDynamics.hh"
#include "PwaUtils/K0star1430LassDynamics.hh" #include "PwaUtils/K0star1430LassDynamics.hh"
#include "PwaUtils/BlattWBarrierDynamics.hh"
#include "PwaUtils/GlobalEnv.hh" #include "PwaUtils/GlobalEnv.hh"
#include "PwaUtils/WoDynamics.hh" #include "PwaUtils/WoDynamics.hh"
...@@ -90,7 +91,7 @@ std::shared_ptr<AbsDynamics> DynRegistry::getDynamics(std::shared_ptr<AbsDecay> ...@@ -90,7 +91,7 @@ std::shared_ptr<AbsDynamics> DynRegistry::getDynamics(std::shared_ptr<AbsDecay>
else if(theDec->dynType()=="BreitWignerRel") else if(theDec->dynType()=="BreitWignerRel")
result= std::shared_ptr<AbsDynamics>(new BreitWignerRelDynamics(theName, fsParticles, theDec->motherPart(), theDec->massSumFsParticlesDec1(), theDec->massSumFsParticlesDec2() )); result= std::shared_ptr<AbsDynamics>(new BreitWignerRelDynamics(theName, fsParticles, theDec->motherPart(), theDec->massSumFsParticlesDec1(), theDec->massSumFsParticlesDec2() ));
else if(theDec->dynType()=="BreitWignerBlattWRel") else if(theDec->dynType()=="BreitWignerBlattWRel")
result= std::shared_ptr<AbsDynamics>(new BreitWignerBlattWRelDynamics(theName, fsParticles, theDec->motherPart(), theDec->massSumFsParticlesDec1(), theDec->massSumFsParticlesDec2())); result= std::shared_ptr<AbsDynamics>(new BreitWignerBlattWRelDynamics(theName, fsParticles, theDec->motherPart(), theDec->massSumFsParticlesDec1(), theDec->massSumFsParticlesDec2(), theDec->barrierqR()));
else if(theDec->dynType()=="KMatrix"){ else if(theDec->dynType()=="KMatrix"){
std::string pathToConfigFile=theDec->pathToConfigParser(); std::string pathToConfigFile=theDec->pathToConfigParser();
result= std::shared_ptr<AbsDynamics>(new KMatrixDynamics(theName, fsParticles, theDec->motherPart(), pathToConfigFile)); result= std::shared_ptr<AbsDynamics>(new KMatrixDynamics(theName, fsParticles, theDec->motherPart(), pathToConfigFile));
...@@ -107,6 +108,10 @@ std::shared_ptr<AbsDynamics> DynRegistry::getDynamics(std::shared_ptr<AbsDecay> ...@@ -107,6 +108,10 @@ std::shared_ptr<AbsDynamics> DynRegistry::getDynamics(std::shared_ptr<AbsDecay>
result= std::shared_ptr<AbsDynamics>(new VoigtDynamics(theName, fsParticles, theDec->motherPart())); result= std::shared_ptr<AbsDynamics>(new VoigtDynamics(theName, fsParticles, theDec->motherPart()));
else if(theDec->dynType()=="K0star1430Lass") else if(theDec->dynType()=="K0star1430Lass")
result= std::shared_ptr<AbsDynamics>(new K0star1430LassDynamics(theName, fsParticles, theDec->motherPart())); result= std::shared_ptr<AbsDynamics>(new K0star1430LassDynamics(theName, fsParticles, theDec->motherPart()));
else if(theDec->dynType()=="BlattWBarrier")
result= std::shared_ptr<AbsDynamics>(new BlattWBarrierDynamics(theName, fsParticles, theDec->motherPart(), theDec->massSumFsParticlesDec1(), theDec->massSumFsParticlesDec2(), theDec->wignerDKey(), theDec->barrierqR()));
else if(theDec->dynType()=="WoDynamics") result= std::shared_ptr<AbsDynamics>(new WoDynamics(theName, fsParticles, theDec->motherPart())); else if(theDec->dynType()=="WoDynamics") result= std::shared_ptr<AbsDynamics>(new WoDynamics(theName, fsParticles, theDec->motherPart()));
else{ else{
Alert << "Dyn type:\t" << theDec->dynType() << "\tdoes not exist" << endmsg; Alert << "Dyn type:\t" << theDec->dynType() << "\tdoes not exist" << endmsg;
......
...@@ -43,6 +43,7 @@ public: ...@@ -43,6 +43,7 @@ public:
FlatteDynamics(std::string& key, std::vector<Particle*>& fsParticles, Particle* mother, std::pair<Particle*, Particle*>& decPair1stChannel, std::pair<Particle*, Particle*>& decPair2ndChannel); FlatteDynamics(std::string& key, std::vector<Particle*>& fsParticles, Particle* mother, std::pair<Particle*, Particle*>& decPair1stChannel, std::pair<Particle*, Particle*>& decPair2ndChannel);
virtual ~FlatteDynamics(); virtual ~FlatteDynamics();
virtual std::string type() {return "FlatteDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -43,6 +43,7 @@ public: ...@@ -43,6 +43,7 @@ public:
K0star1430LassDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother); K0star1430LassDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~K0star1430LassDynamics(); virtual ~K0star1430LassDynamics();
virtual std::string type() {return "K0star1430LassDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -50,6 +50,7 @@ public: ...@@ -50,6 +50,7 @@ public:
KMatrixDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, std::string& pathToConfigParser); KMatrixDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother, std::string& pathToConfigParser);
virtual ~KMatrixDynamics(); virtual ~KMatrixDynamics();
virtual std::string type() {return "KMatrixDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -46,6 +46,7 @@ public: ...@@ -46,6 +46,7 @@ public:
KPiSWaveIso12Dynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother); KPiSWaveIso12Dynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~KPiSWaveIso12Dynamics(); virtual ~KPiSWaveIso12Dynamics();
virtual std::string type() {return "KPiSWaveIso12Dynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -46,6 +46,7 @@ public: ...@@ -46,6 +46,7 @@ public:
KPiSWaveIso32Dynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother); KPiSWaveIso32Dynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~KPiSWaveIso32Dynamics(); virtual ~KPiSWaveIso32Dynamics();
virtual std::string type() {return "KPiSWaveIso32Dynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -49,6 +49,7 @@ public: ...@@ -49,6 +49,7 @@ public:
PiPiSWaveASDynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother, ParticleTable* thePdtTable); PiPiSWaveASDynamics(std::string& massKey, std::vector<Particle*>& fsParticles, Particle* mother, ParticleTable* thePdtTable);
virtual ~PiPiSWaveASDynamics(); virtual ~PiPiSWaveASDynamics();
virtual std::string type() {return "PiPiSWaveASDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -43,6 +43,7 @@ public: ...@@ -43,6 +43,7 @@ public:
VoigtDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother); VoigtDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~VoigtDynamics(); virtual ~VoigtDynamics();
virtual std::string type() {return "VoigtDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
...@@ -42,6 +42,7 @@ public: ...@@ -42,6 +42,7 @@ public:
WoDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother); WoDynamics(std::string& name, std::vector<Particle*>& fsParticles, Particle* mother);
virtual ~WoDynamics(); virtual ~WoDynamics();
virtual std::string type() {return "WoDynamics";}
virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0); virtual complex<double> eval(EvtData* theData, AbsXdecAmp* grandmaAmp, Spin OrbMom=0);
virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr); virtual void getDefaultParams(fitParams& fitVal, fitParams& fitErr);
......
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