VectorTo3PseudoScalarDalitz.h

// -*- C++ -*-
#ifndef Herwig_VectorTo3PseudoScalarDalitz_H
#define Herwig_VectorTo3PseudoScalarDalitz_H
//
// This is the declaration of the VectorTo3PseudoScalarDalitz class.
//
 
#include "DalitzBase.h"
#include "ThePEG/Helicity/LorentzPolarizationVector.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class VectorTo3PseudoScalarDalitz: public DalitzBase {
 
public:
 
  VectorTo3PseudoScalarDalitz() : coupling_(1./GeV)
  {}
 
  virtual double me2(const int ichan,const Particle & part,
                     const tPDVector & outgoing,
                     const vector<Lorentz5Momentum> & momenta,
                     MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
  virtual WidthCalculatorBasePtr threeBodyMEIntegrator(const DecayMode & dm) const;
 
  virtual double threeBodyMatrixElement(const int imode, const Energy2 q2,
                                        const  Energy2 s3, const Energy2 s2, const 
                                        Energy2 s1, const Energy m1,
                                        const Energy m2, const Energy m3) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
  
  virtual complex<InvEnergy2> amplitude(int ichan) const {
    complex<InvEnergy2> amp(ZERO);
    int iloc=-1;
    for(int ix=0;ix<int(resonances().size());++ix) {
      ++iloc;
      if(channel1()>=0) {
        if(ix!=channel1() && ix!=channel2()) continue;
      }
      if(ichan>=0&&ichan!=iloc) continue;
      amp += resAmp(ix);
    }
    return amp;
  }
 
  complex<InvEnergy2> resAmp(unsigned int i) const;
 
private:
 
  VectorTo3PseudoScalarDalitz & operator=(const VectorTo3PseudoScalarDalitz &) = delete;
 
private:
 
  InvEnergy coupling_;
 
  mutable Energy mD_;
 
  mutable Energy mOut_[3];
 
  mutable Energy m2_[3][3];
  
  mutable vector<Helicity::LorentzPolarizationVector> vectors_;
  
  mutable RhoDMatrix rho_;
 
};
 
}
 
#endif /* Herwig_VectorTo3PseudoScalarDalitz_H */