VectorMesonVectorScalarDecayer.cc
No OneTemporary
Actions

Size

6 KB

Referenced Files

None

Subscribers

None

VectorMesonVectorScalarDecayer.cc
View Options

	// -- C++ --
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the VectorMesonVectorScalarDecayer class.
	//

	#include "VectorMesonVectorScalarDecayer.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	#include "ThePEG/PDT/DecayMode.h"
	#include "ThePEG/Interface/ParVector.h"
	#include "ThePEG/Helicity/ScalarSpinInfo.h"
	#include "ThePEG/Helicity/VectorSpinInfo.h"
	#include "Herwig++/Helicity/WaveFunction/VectorWaveFunction.h"

	#ifdef ThePEG_TEMPLATES_IN_CC_FILE
	// #include "VectorMesonVectorScalarDecayer.tcc"
	#endif

	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"

	namespace Herwig {
	using namespace ThePEG;
	using namespace ThePEG::Helicity;
	using ThePEG::Helicity::ScalarSpinInfo;
	using ThePEG::Helicity::VectorSpinInfo;
	using Helicity::VectorWaveFunction;
	using Helicity::Direction;
	using Helicity::incoming;
	using Helicity::outgoing;

	VectorMesonVectorScalarDecayer::~VectorMesonVectorScalarDecayer() {}

	bool VectorMesonVectorScalarDecayer::accept(const DecayMode & dm) const {
	// is this mode allowed
	bool allowed=false;
	// must be two outgoing particles
	if(dm.products().size()!=2){return allowed;}
	// ids of the particles
	int id0=dm.parent()->id();
	ParticleMSet::const_iterator pit = dm.products().begin();
	int id1=(**pit).id();
	++pit;
	int id2=(**pit).id();
	unsigned int ix=0;
	do
	{
	if(id0==_incoming[ix])
	{
	if((id1==_outgoingV[ix]&&id2==_outgoingS[ix])\|\|
	(id2==_outgoingV[ix]&&id1==_outgoingS[ix])){allowed=true;}
	}
	++ix;
	}
	while(ix<_incoming.size()&&!allowed);
	return allowed;
	}

	ParticleVector VectorMesonVectorScalarDecayer::decay(const DecayMode & dm,
	const Particle & parent) const {
	int imode=-1;
	int id=parent.id();
	ParticleMSet::const_iterator pit = dm.products().begin();
	int id1=(**pit).id();
	++pit;
	int id2=(**pit).id();
	unsigned int ix=0;
	do
	{
	if(id==_incoming[ix])
	{
	if((id1==_outgoingV[ix]&&id2==_outgoingS[ix])\|\|
	(id2==_outgoingV[ix]&&id1==_outgoingS[ix])){imode=ix;}
	}
	++ix;
	}
	while(ix<_incoming.size()&&imode<0);
	// generate the decay
	bool cc=false;
	return generate(false,cc,imode,parent);
	}


	void VectorMesonVectorScalarDecayer::persistentOutput(PersistentOStream & os) const
	{os << _incoming << _outgoingV << _outgoingS << _maxweight << _coupling;}

	void VectorMesonVectorScalarDecayer::persistentInput(PersistentIStream & is, int)
	{is >> _incoming >> _outgoingV >> _outgoingS >> _maxweight >> _coupling;}

	ClassDescription<VectorMesonVectorScalarDecayer> VectorMesonVectorScalarDecayer::initVectorMesonVectorScalarDecayer;
	// Definition of the static class description member.

	void VectorMesonVectorScalarDecayer::Init() {

	static ClassDocumentation<VectorMesonVectorScalarDecayer> documentation
	("The\\classname{VectorMesonVectorScalarDecayer} class is designed for the "
	"decay of a vector meson to a vector meson, or the photon, and a "
	"scalar meson.");

	static ParVector<VectorMesonVectorScalarDecayer,int> interfaceIncoming
	("Incoming",
	"The PDG code for the incoming particle",
	&VectorMesonVectorScalarDecayer::_incoming,
	0, 0, 0, -10000, 10000, false, false, true);

	static ParVector<VectorMesonVectorScalarDecayer,int> interfaceOutcomingVector
	("OutgoingVector",
	"The PDG code for the outgoing spin-1 particle",
	&VectorMesonVectorScalarDecayer::_outgoingV,
	0, 0, 0, -10000, 10000, false, false, true);

	static ParVector<VectorMesonVectorScalarDecayer,int> interfaceOutcomingScalar
	("OutgoingScalar",
	"The PDG code for the outgoing spin-0 particle",
	&VectorMesonVectorScalarDecayer::_outgoingS,
	0, 0, 0, -10000, 10000, false, false, true);

	static ParVector<VectorMesonVectorScalarDecayer,InvEnergy> interfaceCoupling
	("Coupling",
	"The coupling for the decay mode",
	&VectorMesonVectorScalarDecayer::_coupling,
	0, 0, 0, -10000, 10000, false, false, true);

	static ParVector<VectorMesonVectorScalarDecayer,double> interfaceMaxWeight
	("MaxWeight",
	"The maximum weight for the decay mode",
	&VectorMesonVectorScalarDecayer::_maxweight,
	0, 0, 0, -10000, 10000, false, false, true);
	}

	// the hadronic currents
	vector<LorentzPolarizationVector>
	VectorMesonVectorScalarDecayer::decayCurrent(const bool vertex, const int,
	const Particle & inpart,
	const ParticleVector & decay) const
	{
	// storage for the current
	vector<LorentzPolarizationVector> temp;
	// work out which of the decay products is the vector and which is the scalar
	unsigned int ivec=1,isca=0;
	bool photon=_outgoingV[imode()]==ParticleID::gamma;
	if(decay[0]->id()==_outgoingV[imode()]){ivec=0;isca=1;}
	tVectorSpinPtr vecsp;
	// set up the spin information for the decay products
	if(vertex)
	{
	// scalar
	decay[isca]->spinInfo(new_ptr(ScalarSpinInfo(decay[isca]->momentum(),true)));
	// vector
	SpinPtr temp=new_ptr(VectorSpinInfo(decay[ivec]->momentum(),true));
	vecsp = dynamic_ptr_cast<tVectorSpinPtr>(temp);
	decay[ivec]->spinInfo(temp);
	}
	VectorWaveFunction vwave=VectorWaveFunction(decay[ivec]->momentum(),
	decay[ivec]->dataPtr(),outgoing);
	// calculate the currents
	Energy2 p0dotpv=inpart.momentum()*decay[ivec]->momentum();
	Complex epsdot=0.;
	LorentzPolarizationVector output;
	for(int ix=-1;ix<2;++ix)
	{
	if(ix==0&&photon)
	{
	temp.push_back(LorentzPolarizationVector());
	if(vertex){vecsp->setBasisState(ix,LorentzPolarizationVector());}
	}
	else
	{
	vwave.reset(ix);
	//if(vertex){vecsp->setBasisState(ix,vwave.Wave());}
	epsdot=vwave.Wave()*inpart.momentum();
	output = p0dotpvvwave.Wave()-epsdotdecay[ivec]->momentum();
	output *=_coupling[imode()]/inpart.mass();
	temp.push_back(output);
	}
	}
	return temp;
	}

	bool VectorMesonVectorScalarDecayer::twoBodyMEcode(const DecayMode & dm,
	int & mecode,
	double & coupling) const
	{
	// work out which mode it is for the coupling
	int imode=-1;
	int id=dm.parent()->id();
	ParticleMSet::const_iterator pit = dm.products().begin();
	int id1=(**pit).id();
	++pit;
	int id2=(**pit).id();
	unsigned int ix=0;
	do
	{
	if(id==_incoming[ix])
	{
	if((id1==_outgoingV[ix]&&id2==_outgoingS[ix])\|\|
	(id2==_outgoingV[ix]&&id1==_outgoingS[ix])){imode=ix;}
	}
	++ix;
	}
	while(ix<_incoming.size()&&imode<0);
	coupling = _coupling[imode]*dm.parent()->mass();
	mecode = 4;
	return id1==_outgoingV[imode]&&id2==_outgoingS[imode];
	}
	}

File Metadata

Mime Type: text/x-c
Expires: Tue, Sep 30, 4:48 AM (10 h, 13 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 6549572
Default Alt Text: VectorMesonVectorScalarDecayer.cc (6 KB)

VectorMesonVectorScalarDecayer.ccNo OneTemporaryActions

VectorMesonVectorScalarDecayer.ccView Options

File Metadata

Event Timeline

VectorMesonVectorScalarDecayer.cc
No OneTemporary
Actions

VectorMesonVectorScalarDecayer.cc
View Options