Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F19244127
ThreeMesonDefaultCurrent.icc
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Award Token
Flag For Later
Size
8 KB
Referenced Files
None
Subscribers
None
ThreeMesonDefaultCurrent.icc
View Options
// -*- C++ -*-
//
// This is the implementation of the inlined member functions of
// the ThreeMesonDefaultCurrent class.
//
namespace Herwig {
using namespace ThePEG;
inline ThreeMesonDefaultCurrent::
ThreeMesonDefaultCurrent(const ThreeMesonDefaultCurrent & x)
: ThreeMesonCurrentBase(x),_rhoF123wgts(x._rhoF123wgts),
_KstarF123wgts(x._KstarF123wgts),_rhoF5wgts(x._rhoF5wgts),
_KstarF5wgts(x._KstarF5wgts), _rhoKstarwgt(x._rhoKstarwgt),
_a1runwidth(x._a1runwidth),_a1runq2(x._a1runq2),_a1runinter(x._a1runinter),
_initializea1(x._initializea1), _a1mass(x._a1mass),_a1width(x._a1width),
_K1mass(x._K1mass),_K1width(x._K1width), _fpi(x._fpi),_mpi(x._mpi),_mK(x._mK),
_rhoparameters(x._rhoparameters), _rhoF123masses(x._rhoF123masses),
_rhoF5masses(x._rhoF5masses), _rhoF123widths(x._rhoF123widths),
_rhoF5widths(x._rhoF5widths), _Kstarparameters(x._Kstarparameters),
_KstarF123masses(x._KstarF123masses),_KstarF5masses(x._KstarF5masses),
_KstarF123widths(x._KstarF123widths),_KstarF5widths(x._KstarF5widths),
_a1parameters(x._a1parameters), _K1parameters(x._K1parameters) {}
inline IBPtr ThreeMesonDefaultCurrent::clone() const {
return new_ptr(*this);
}
inline IBPtr ThreeMesonDefaultCurrent::fullclone() const {
return new_ptr(*this);
}
inline void ThreeMesonDefaultCurrent::doupdate() throw(UpdateException) {
ThreeMesonCurrentBase::doupdate();
// First update base class.
bool redo = touched();
// redo if touched.
// UpdateChecker::check(aDependentMember, redo);
// Update referenced objects on which this depends redo is set to true
// if the dependent object is touched.
// for_each(ContainerOfDependencies, UpdateChecker(redo));
// Update a container of references.
// for_each(MapOfDependencies, UpdateMapChecker(redo));
// Update a map of references.
if ( !redo ) return;
// return if nothing has been touched. Otherwise do the actual update.
// touch()
// Touch if anything has changed.
}
inline void ThreeMesonDefaultCurrent::dofinish() {
ThreeMesonCurrentBase::dofinish();
}
inline void ThreeMesonDefaultCurrent::doinitrun() {
// set up the running a_1 width
inita1width(0);
ThreeMesonCurrentBase::doinitrun();
}
inline void ThreeMesonDefaultCurrent::rebind(const TranslationMap & trans)
throw(RebindException) {
// dummy = trans.translate(dummy);
ThreeMesonCurrentBase::rebind(trans);
}
inline IVector ThreeMesonDefaultCurrent::getReferences() {
IVector ret = ThreeMesonCurrentBase::getReferences();
// ret.push_back(dummy);
return ret;
}
// functions to return the Breit-Wigners
inline Complex ThreeMesonDefaultCurrent::BrhoF123(Energy2 q2,int ires) const
{
Complex output(0.),norm(0.);
for(unsigned int ix=0,N=min(3,int(_rhoF123wgts.size()));ix<N;++ix)
{norm+=_rhoF123wgts[ix];}
if(ires<0)
{for(unsigned int ix=0,N=min(3,int(_rhoF123wgts.size()));ix<N;++ix)
{output+=_rhoF123wgts[ix]*rhoKBreitWigner(q2,0,ix);}}
else
{
unsigned int temp(ires);
if(temp<_rhoF123wgts.size()&&temp<3)
{output=_rhoF123wgts[temp]*rhoKBreitWigner(q2,0,temp);}
else{output=0.;}
}
return output/norm;
}
inline Complex ThreeMesonDefaultCurrent::BrhoF5(Energy2 q2,int ires) const
{
Complex output(0.),norm(0.);
for(unsigned int ix=0,N=min(3,int(_rhoF5wgts.size()));ix<N;++ix)
{norm+=_rhoF5wgts[ix];}
if(ires<0)
{for(unsigned int ix=0,N=min(3,int(_rhoF5wgts.size()));ix<N;++ix)
{output+=_rhoF5wgts[ix]*rhoKBreitWigner(q2,1,ix);}}
else
{
unsigned int temp(ires);
if(temp<_rhoF5wgts.size()&&temp<3)
{output=_rhoF5wgts[temp]*rhoKBreitWigner(q2,1,temp);}
}
return output/norm;
}
inline Complex ThreeMesonDefaultCurrent::BKstarF123(Energy2 q2,int ires) const
{
Complex output(0.),norm(0.);
for(unsigned int ix=0,N=min(3,int(_KstarF123wgts.size()));ix<N;++ix)
{norm+=_KstarF123wgts[ix];}
if(ires<0)
{for(unsigned int ix=0,N=min(3,int(_KstarF123wgts.size()));ix<N;++ix)
{output+=_KstarF123wgts[ix]*rhoKBreitWigner(q2,2,ix);}}
else
{
unsigned int temp(ires);
if(temp<_KstarF123wgts.size()&&temp<3)
{output=_KstarF123wgts[temp]*rhoKBreitWigner(q2,2,temp);}
}
return output/norm;
}
inline Complex ThreeMesonDefaultCurrent::BKstarF5(Energy2 q2,int ires) const
{
Complex output(0.),norm(0.);
for(unsigned int ix=0,N=min(3,int(_KstarF5wgts.size()));ix<N;++ix)
{norm+=_KstarF5wgts[ix];}
if(ires<0)
{for(unsigned int ix=0,N=min(3,int(_KstarF5wgts.size()));ix<N;++ix)
{output+=_KstarF5wgts[ix]*rhoKBreitWigner(q2,3,ix);}}
else
{
unsigned int temp(ires);
if(temp<_KstarF5wgts.size()&&temp<3)
{output=_KstarF5wgts[ires]*rhoKBreitWigner(q2,3,temp);}
}
return output/norm;
}
// matrix element for the running a_1 width
inline double ThreeMesonDefaultCurrent::a1MatrixElement(Energy2 q2, Energy2 s3,
Energy2 s2,Energy2 s1,
Energy m1,Energy m2,Energy m3)
{
Energy2 mpi2(getParticleData(ParticleID::piplus)->mass());mpi2*=mpi2;
Complex propb(BrhoF123(s1,-1)),propa(BrhoF123(s2,-1));
// the matrix element
double output(0.);
// first resonance
output+= ((s1-4.*mpi2)+0.25*(s3-s2)*(s3-s2)/q2)*real(propb*conj(propb));
// second resonance
output+= ((s2-4.*mpi2)+0.25*(s3-s1)*(s3-s1)/q2)*real(propa*conj(propa));
// the interference term
output+= (0.5*q2-s3-0.5*mpi2+0.25*(s3-s2)*(s3-s1)/q2)*real(propa*conj(propb)+
propb*conj(propa));
return output;
}
// initialisation of the a_1 width
// (iopt=-1 initialises, iopt=0 starts the interpolation)
inline void ThreeMesonDefaultCurrent::inita1width(int iopt)
{
if(iopt==-1)
{
// parameters for the table of values
Energy mtau(getParticleData(ParticleID::tauplus)->mass());
Energy2 mtau2(mtau*mtau),step(mtau*mtau/199.);
// function to be integrated to give the matrix element
Genfun::AbsFunction *inme = new Defaulta1MatrixElement(this);
// integrator to perform the integral
vector<double> inweights;inweights.push_back(0.5);inweights.push_back(0.5);
vector<int> intype;intype.push_back(2);intype.push_back(3);
Energy mrho(getParticleData(ParticleID::rhoplus)->mass()),
wrho(getParticleData(ParticleID::rhoplus)->width());
vector<Energy> inmass(2,mrho),inwidth(2,wrho);
Energy m1(getParticleData(ParticleID::piplus)->mass());
ThreeBodyIntegrator *widthgen = new ThreeBodyIntegrator(inweights,intype,
inmass,inwidth,
inme,m1,m1,m1);
// nornmalisation constant to give physical width if on shell
double a1const(_a1width/(widthgen->width(_a1mass*_a1mass)));
// loop to give the values
Energy moff2(0.);
_a1runq2.resize(0);_a1runwidth.resize(0);
for(;moff2<=mtau2;moff2+=step)
{
_a1runwidth.push_back(widthgen->width(moff2)*a1const);
_a1runq2.push_back(moff2);
}
delete inme; delete widthgen;
}
// set up the interpolator
else if(iopt==0)
{_a1runinter = new Interpolator(_a1runwidth,_a1runq2,3);}
}
// running width of the a_1
inline Energy ThreeMesonDefaultCurrent::a1width(Energy2 q2) const
{return (*_a1runinter)(q2);}
// a_1 Breit Wigner
inline Complex ThreeMesonDefaultCurrent::a1BreitWigner(Energy2 q2) const
{
Complex ii(0.,1.);
Energy2 m2(_a1mass*_a1mass);
Energy q(sqrt(q2));
return m2/(m2-q2-ii*q*a1width(q2));
}
// K_1 Breit Wigner
inline Complex ThreeMesonDefaultCurrent::K1BreitWigner(Energy2 q2) const
{
Energy m2(_K1mass*_K1mass);
Complex ii(0.,1.);
complex<Energy2> fact(m2-ii*_K1mass*_K1width);
return fact/(fact-q2);
}
// mixed Breit Wigner
inline Complex ThreeMesonDefaultCurrent::FKrho(Energy2 si, Energy2 sj,int ires) const
{
Complex output;
if(ires<0){output = _rhoKstarwgt*BKstarF123(si,-1)+BrhoF123(sj,-1);}
else if(ires%2==0){output= _rhoKstarwgt*BKstarF123(si,ires/2);}
else if(ires%2==1){output=BrhoF123(sj,ires/2);}
output /=(1.+_rhoKstarwgt);
return output;
}
inline Complex ThreeMesonDefaultCurrent::rhoKBreitWigner(Energy2 q2,unsigned int itype,
unsigned int ires) const
{
Energy q(sqrt(q2)),mass,width,mout[2];
// get the mass and width of the requested resonance
if(itype==0){mass=_rhoF123masses[ires];width=_rhoF123widths[ires];}
else if(itype==1){mass=_rhoF5masses[ires];width=_rhoF5widths[ires];}
else if(itype==2){mass=_KstarF123masses[ires];width=_KstarF123widths[ires];}
else if(itype==3){mass=_KstarF5masses[ires];width=_KstarF5widths[ires];}
// calculate the momenta for the running widths
if(itype<=1){mout[0]=_mpi;mout[1]=_mpi;}
else{mout[0]=_mK;mout[1]=_mpi;}
Energy pcm0(Kinematics::pstarTwoBodyDecay(mass,mout[0],mout[1])),pcm(0.);
if(mout[0]+mout[1]<q){pcm=Kinematics::pstarTwoBodyDecay(q,mout[0],mout[1]);}
double ratio(pcm/pcm0);ratio*=ratio*ratio;
Energy gamrun(width*mass*ratio/q);
Complex ii(0.,1.);
complex<Energy2> denom(q2-mass*mass+ii*mass*gamrun),numer(-mass*mass);
return numer/denom;
}
}
File Metadata
Details
Attached
Mime Type
text/plain
Expires
Tue, Sep 30, 4:39 AM (15 h, 33 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6564733
Default Alt Text
ThreeMesonDefaultCurrent.icc (8 KB)
Attached To
Mode
rHERWIGHG herwighg
Attached
Detach File
Event Timeline
Log In to Comment