DalitzKMatrix.h

// -*- C++ -*-
#ifndef Herwig_DalitzKMatrix_H
#define Herwig_DalitzKMatrix_H
//
// This is the declaration of the DalitzKMatrix class.
//
 
#include "DalitzResonance.h"
#include "Herwig/Decay/FormFactors/KMatrix.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class DalitzKMatrix: public DalitzResonance {
 
public:
 
  DalitzKMatrix()
  {}
 
  DalitzKMatrix(long pid, ResonanceType::Type rtype, Energy m, Energy w,
                unsigned int d1, unsigned int d2, unsigned int s,
                double mag, double phi, InvEnergy rr,
                unsigned int imat, unsigned int chan,
                Energy2 sc, unsigned int itype,
                vector<pair<double,double> > beta,
                vector<pair<double,vector<double > > > coeffs)
    : DalitzResonance(pid,rtype,m,w,d1,d2,s,mag,phi,rr),
      imat_(imat), channel_(chan), sc_(sc), expType_(itype), coeffs_(coeffs) {
    beta_.clear();
    for(unsigned int ix=0;ix<beta.size();++ix) {
      beta_.push_back(beta[ix].first*exp(Complex(0.,beta[ix].second)));
    }
  }
 
public:
 
  virtual Complex BreitWigner(const Energy & mAB, const Energy & mA, const Energy & mB) const;
 
  virtual void dataBaseOutput(ofstream & output);
 
public:
 
  void setKMatrix(KMatrixPtr mat) {kMatrix_ = mat;}
 
  unsigned int imatrix() const {return imat_;}
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
private:
 
  DalitzKMatrix & operator=(const DalitzKMatrix &) = delete;
 
private:
 
  KMatrixPtr kMatrix_;
 
  unsigned int imat_;
 
  unsigned int channel_;
 
  Energy2 sc_;
 
  vector<Complex> beta_;
 
  unsigned int expType_;
 
  vector<pair<double,vector<double > > > coeffs_;
 
};
 
}
 
#endif /* Herwig_DalitzKMatrix_H */