Page MenuHomeHEPForge
Files F10881712

No OneTemporary
Actions

View Options

diff --git a/Helicity/Vertex/AbstractVVVTVertex.h b/Helicity/Vertex/AbstractVVVTVertex.h
--- a/Helicity/Vertex/AbstractVVVTVertex.h
+++ b/Helicity/Vertex/AbstractVVVTVertex.h
@@ -1,113 +1,75 @@
// -*- C++ -*-
#ifndef HELICITY_AbstractVVVTVertex_H
#define HELICITY_AbstractVVVTVertex_H
//
// This is the declaration of the AbstractVVVTVertex class.
//
#include "VertexBase.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/TensorWaveFunction.h"
#include "AbstractVVVTVertex.fh"
namespace ThePEG {
namespace Helicity {
/**
* The AbstractVVVTVertex class is the base class for all vector-vector-vector-tensor
* interactions in ThePEG.
*/
class AbstractVVVTVertex: public VertexBase {
public:
/**
* Default constructor
*/
AbstractVVVTVertex() : VertexBase(VertexType::VVVT) {}
/**
* Members to calculate the helicity amplitude expressions for vertices
* and off-shell particles.
*/
//@{
/**
* Evalulate the vertex.
* @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
* @param vec1 The wavefunction for the first vector.
* @param vec2 The wavefunction for the second vector.
* @param vec3 The wavefunction for the third vector.
* @param ten4 The wavefunction for the tensor.
*/
virtual Complex evaluate(Energy2 q2,const VectorWaveFunction & vec1,
const VectorWaveFunction & vec2,
const VectorWaveFunction & vec3,
const TensorWaveFunction & ten4) = 0;
-
- /**
- * Evaluate the off-shell tensor coming from the vertex.
- * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
- * @param iopt Option of the shape of the Breit-Wigner for the off-shell tensor.
- * @param out The ParticleData pointer for the off-shell tensor.
- * @param vec1 The wavefunction for the first vector.
- * @param vec2 The wavefunction for the second vector.
- * @param vec3 The wavefunction for the third vector.
- * @param mass The mass of the off-shell particle if not taken from the ParticleData
- * object
- * @param width The width of the off-shell particle if not taken from the ParticleData
- * object
- */
- virtual TensorWaveFunction evaluate(Energy2 q2,int iopt, tcPDPtr out,
- const VectorWaveFunction & vec1,
- const VectorWaveFunction & vec2,
- const VectorWaveFunction & vec3,
- complex<Energy> mass=-GeV, complex<Energy> width=-GeV) = 0;
-
- /**
- * Evaluate the off-shell vector coming from the vertex.
- * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
- * @param iopt Option of the shape of the Breit-Wigner for the off-shell vector.
- * @param out The ParticleData pointer for the off-shell vector.
- * @param vec1 The wavefunction for the first vector.
- * @param vec2 The wavefunction for the second vector.
- * @param ten4 The wavefunction for the tensor.
- * @param mass The mass of the off-shell particle if not taken from the ParticleData
- * object
- * @param width The width of the off-shell particle if not taken from the ParticleData
- * object
- */
- virtual VectorWaveFunction evaluate(Energy2 q2,int iopt, tcPDPtr out,
- const VectorWaveFunction & vec1,
- const VectorWaveFunction & vec2,
- const TensorWaveFunction & ten4,
- complex<Energy> mass=-GeV, complex<Energy> width=-GeV) = 0;
//@}
public:
/**
* The standard Init function used to initialize the interfaces.
* Called exactly once for each class by the class description system
* before the main function starts or
* when this class is dynamically loaded.
*/
static void Init();
private:
/**
* The assignment operator is private and must never be called.
* In fact, it should not even be implemented.
*/
AbstractVVVTVertex & operator=(const AbstractVVVTVertex &);
};
}
}
namespace ThePEG {
}
#endif /* HELICITY_AbstractVVVTVertex_H */
diff --git a/Helicity/Vertex/Tensor/VVVTVertex.cc b/Helicity/Vertex/Tensor/VVVTVertex.cc
--- a/Helicity/Vertex/Tensor/VVVTVertex.cc
+++ b/Helicity/Vertex/Tensor/VVVTVertex.cc
@@ -1,153 +1,131 @@
// -*- C++ -*-
//
// VVVTVertex.cc is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 2003-2017 Peter Richardson, Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the VVVTVertex class.
//
#include "VVVTVertex.h"
#include "ThePEG/Utilities/DescribeClass.h"
#include "ThePEG/Interface/ClassDocumentation.h"
using namespace ThePEG;
using namespace Helicity;
// The following static variable is needed for the type
// description system in ThePEG.
DescribeAbstractNoPIOClass<VVVTVertex,AbstractVVVTVertex>
describeThePEGVVVTVertex("ThePEG::VVVTVertex", "libThePEG.so");
void VVVTVertex::Init() {
static ClassDocumentation<VVVTVertex> documentation
("The VVVTVertex class is the implementation f the"
" helicity amplitude calculation of the vector-vector-vector-tensor"
" vertex. All such vertices should inherit from it.");
}
// function to evaluate the vertex
Complex VVVTVertex::evaluate(Energy2 q2, const VectorWaveFunction & vec1,
const VectorWaveFunction & vec2,
const VectorWaveFunction & vec3,
const TensorWaveFunction & ten) {
// set the couplings
setCoupling(q2,vec1.particle(),vec2.particle(),
vec3.particle(),ten.particle());
Complex ii(0.,1.);
// dot products of the wavefunctions
Complex dotv1v2 = vec1.wave().dot(vec2.wave());
Complex dotv1v3 = vec1.wave().dot(vec3.wave());
Complex dotv2v3 = vec2.wave().dot(vec3.wave());
// dot product of wavefunctions and momenta
complex<Energy> dotv1k23 = vec1.wave().dot(vec2.momentum()-vec3.momentum());
complex<Energy> dotv2k31 = vec2.wave().dot(vec3.momentum()-vec1.momentum());
complex<Energy> dotv3k12 = vec3.wave().dot(vec1.momentum()-vec2.momentum());
// components of the tensor
Complex tentx = ten.tx()+ten.xt();
Complex tenty = ten.ty()+ten.yt();
Complex tentz = ten.tz()+ten.zt();
Complex tenxy = ten.xy()+ten.yx();
Complex tenxz = ten.xz()+ten.zx();
Complex tenyz = ten.yz()+ten.zy();
// dot product of wavefunctions and momenta with the tensor
Complex tenv1v2 =
2.*(+ten.tt()*vec1.t()*vec2.t()+ten.xx()*vec1.x()*vec2.x()
+ten.yy()*vec1.y()*vec2.y()+ten.zz()*vec1.z()*vec2.z())
-tentx*(vec1.t()*vec2.x()+vec1.x()*vec2.t())
-tenty*(vec1.t()*vec2.y()+vec1.y()*vec2.t())
-tentz*(vec1.t()*vec2.z()+vec1.z()*vec2.t())
+tenxy*(vec1.x()*vec2.y()+vec1.y()*vec2.x())
+tenxz*(vec1.x()*vec2.z()+vec1.z()*vec2.x())
+tenyz*(vec1.y()*vec2.z()+vec1.z()*vec2.y());
Complex tenv1v3 =
2.*(+ten.tt()*vec1.t()*vec3.t()+ten.xx()*vec1.x()*vec3.x()
+ten.yy()*vec1.y()*vec3.y()+ten.zz()*vec1.z()*vec3.z())
-tentx*(vec1.t()*vec3.x()+vec1.x()*vec3.t())
-tenty*(vec1.t()*vec3.y()+vec1.y()*vec3.t())
-tentz*(vec1.t()*vec3.z()+vec1.z()*vec3.t())
+tenxy*(vec1.x()*vec3.y()+vec1.y()*vec3.x())
+tenxz*(vec1.x()*vec3.z()+vec1.z()*vec3.x())
+tenyz*(vec1.y()*vec3.z()+vec1.z()*vec3.y());
Complex tenv2v3 =
2.*(+ten.tt()*vec2.t()*vec3.t()+ten.xx()*vec2.x()*vec3.x()
+ten.yy()*vec2.y()*vec3.y()+ten.zz()*vec2.z()*vec3.z())
-tentx*(vec2.t()*vec3.x()+vec2.x()*vec3.t())
-tenty*(vec2.t()*vec3.y()+vec2.y()*vec3.t())
-tentz*(vec2.t()*vec3.z()+vec2.z()*vec3.t())
+tenxy*(vec2.x()*vec3.y()+vec2.y()*vec3.x())
+tenxz*(vec2.x()*vec3.z()+vec2.z()*vec3.x())
+tenyz*(vec2.y()*vec3.z()+vec2.z()*vec3.y());
complex<Energy> tenv1k23 =
2.*(+ten.tt()*vec1.t()*(vec2.e() -vec3.e() )
+ten.xx()*vec1.x()*(vec2.px()-vec3.px())
+ten.yy()*vec1.y()*(vec2.py()-vec3.py())
+ten.zz()*vec1.z()*(vec2.pz()-vec3.pz()))
-tentx*(vec1.t()*(vec2.px()-vec3.px())+vec1.x()*(vec2.e() -vec3.e() ))
-tenty*(vec1.t()*(vec2.py()-vec3.py())+vec1.y()*(vec2.e() -vec3.e() ))
-tentz*(vec1.t()*(vec2.pz()-vec3.pz())+vec1.z()*(vec2.e() -vec3.e() ))
+tenxy*(vec1.x()*(vec2.py()-vec3.py())+vec1.y()*(vec2.px()-vec3.px()))
+tenxz*(vec1.x()*(vec2.pz()-vec3.pz())+vec1.z()*(vec2.px()-vec3.px()))
+tenyz*(vec1.y()*(vec2.pz()-vec3.pz())+vec1.z()*(vec2.py()-vec3.py()));
complex<Energy> tenv2k31 =
2.*(+ten.tt()*vec2.t()*(vec3.e() -vec1.e() )
+ten.xx()*vec2.x()*(vec3.px()-vec1.px())
+ten.yy()*vec2.y()*(vec3.py()-vec1.py())
+ten.zz()*vec2.z()*(vec3.pz()-vec1.pz()))
-tentx*(vec2.t()*(vec3.px()-vec1.px())+vec2.x()*(vec3.e() -vec1.e() ))
-tenty*(vec2.t()*(vec3.py()-vec1.py())+vec2.y()*(vec3.e() -vec1.e() ))
-tentz*(vec2.t()*(vec3.pz()-vec1.pz())+vec2.z()*(vec3.e() -vec1.e() ))
+tenxy*(vec2.x()*(vec3.py()-vec1.py())+vec2.y()*(vec3.px()-vec1.px()))
+tenxz*(vec2.x()*(vec3.pz()-vec1.pz())+vec2.z()*(vec3.px()-vec1.px()))
+tenyz*(vec2.y()*(vec3.pz()-vec1.pz())+vec2.z()*(vec3.py()-vec1.py()));
complex<Energy> tenv3k12 =
2.*(+ten.tt()*vec3.t()*(vec1.e() -vec2.e() )
+ten.xx()*vec3.x()*(vec1.px()-vec2.px())
+ten.yy()*vec3.y()*(vec1.py()-vec2.py())
+ten.zz()*vec3.z()*(vec1.pz()-vec2.pz()))
-tentx*(vec3.t()*(vec1.px()-vec2.px())+vec3.x()*(vec1.e() -vec2.e() ))
-tenty*(vec3.t()*(vec1.py()-vec2.py())+vec3.y()*(vec1.e() -vec2.e() ))
-tentz*(vec3.t()*(vec1.pz()-vec2.pz())+vec3.z()*(vec1.e() -vec2.e() ))
+tenxy*(vec3.x()*(vec1.py()-vec2.py())+vec3.y()*(vec1.px()-vec2.px()))
+tenxz*(vec3.x()*(vec1.pz()-vec2.pz())+vec3.z()*(vec1.px()-vec2.px()))
+tenyz*(vec3.y()*(vec1.pz()-vec2.pz())+vec3.z()*(vec1.py()-vec2.py()));
// trace of the tensor
Complex trace = ten.tt()-ten.xx()-ten.yy()-ten.zz();
// compute the vertex
Complex vertex= -0.5 * ii * norm() * UnitRemoval::InvE *
(
+dotv3k12*(tenv1v2-trace*dotv1v2)
+dotv2k31*(tenv1v3-trace*dotv1v3)
+dotv1k23*(tenv2v3-trace*dotv2v3)
+dotv2v3*tenv1k23+dotv1v3*tenv2k31+dotv1v2*tenv3k12
);
// return the answer
return vertex;
}
-
-TensorWaveFunction VVVTVertex::evaluate(Energy2 ,int , tcPDPtr ,
- const VectorWaveFunction & ,
- const VectorWaveFunction & ,
- const VectorWaveFunction &,
- complex<Energy>, complex<Energy> ) {
- throw Exception() << "VVVTVertex::evaluate() only implemented for the "
- << "member which returns the amplitude, "
- << "not the off-shell wavefunctions"
- << Exception::runerror;
-}
-
-VectorWaveFunction VVVTVertex::evaluate(Energy2 ,int , tcPDPtr ,
- const VectorWaveFunction & ,
- const VectorWaveFunction & ,
- const TensorWaveFunction &,
- complex<Energy>, complex<Energy> ) {
- throw Exception() << "VVVTVertex::evaluate() only implemented for the "
- << "member which returns the amplitude, "
- << "not the off-shell wavefunctions"
- << Exception::runerror;
-}

File Metadata

Mime Type
text/x-diff
Expires
Sat, May 3, 6:44 AM (1 d, 4 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4974567
Default Alt Text
(10 KB)

Event Timeline