RFVVertex.h

// -*- C++ -*-
//
// RFVVertex.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 2003-2019 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.
//
#ifndef ThePEG_RFVVertex_H
#define ThePEG_RFVVertex_H
//
// This is the declaration of the RFVVertex class.

#include <ThePEG/Helicity/Vertex/AbstractRFVVertex.h>
#include <ThePEG/Helicity/WaveFunction/VectorWaveFunction.h>
#include <ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h>
#include <ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h>
#include "RFVVertex.fh"

namespace ThePEG {

namespace Helicity{

class RFVVertex: public AbstractRFVVertex {
  
public:
  
  static void Init();
  
public:

  Complex evaluate(Energy2 q2,const RSSpinorWaveFunction & sp1,
                   const SpinorBarWaveFunction & sbar2,
                   const VectorWaveFunction & vec3);

  Complex evaluate(Energy2 q2,const SpinorWaveFunction & sp1,
                   const RSSpinorBarWaveFunction & sbar2,
                   const VectorWaveFunction & vec3);

  SpinorBarWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                                 const RSSpinorBarWaveFunction & sbar2,
                                 const VectorWaveFunction & vec3,
                                 complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  RSSpinorBarWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                                   const SpinorBarWaveFunction & sbar2,
                                   const VectorWaveFunction & vec3,
                                   complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  VectorWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                              const RSSpinorWaveFunction & sp1,
                              const SpinorBarWaveFunction & sbar2,
                              complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  VectorWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                              const SpinorWaveFunction & sp1,
                              const RSSpinorBarWaveFunction & sbar2,
                              complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  RSSpinorWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                                const SpinorWaveFunction & sp1,
                                const VectorWaveFunction & vec3,
                                complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  SpinorWaveFunction evaluate(Energy2 q2,int iopt,tcPDPtr out,
                              const RSSpinorWaveFunction & sp1,
                              const VectorWaveFunction & vec3,
                              complex<Energy> mass=-GeV, complex<Energy> width=-GeV);

  virtual void setCoupling(Energy2 q2,tcPDPtr part1,
                           tcPDPtr part2,tcPDPtr part3)=0;

  virtual void setCoupling(Energy2,tcPDPtr,tcPDPtr,tcPDPtr,tcPDPtr) {
    assert(false);
  }

  const vector<Complex> &  left() const { return _left; }

  const vector<Complex> & right() const { return _right; }

protected:

  void left(const vector<Complex> & in) { _left = in; }

  void right(const vector<Complex> & in) { _right = in; }

private:
  
  RFVVertex & operator=(const RFVVertex &) = delete;
  
private:

  vector<Complex> _left;

  vector<Complex> _right;
  
};
}

}
#endif /* ThePEG_RFVVertex_H */