TwoKaonCzyzCurrent.h

// -*- C++ -*-
#ifndef Herwig_TwoKaonCzyzCurrent_H
#define Herwig_TwoKaonCzyzCurrent_H
//
// This is the declaration of the TwoKaonCzyzCurrent class.
//
 
#include "WeakCurrent.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class TwoKaonCzyzCurrent: public WeakCurrent {
 
public:
 
  TwoKaonCzyzCurrent();
 
 
  virtual bool createMode(int icharge, tcPDPtr resonance,FlavourInfo flavour,
                          unsigned int imode,PhaseSpaceModePtr mode,
                          unsigned int iloc,int ires,
                          PhaseSpaceChannel phase, Energy upp );
 
  virtual tPDVector particles(int icharge, unsigned int imode, int iq, int ia);
 
  virtual vector<LorentzPolarizationVectorE> 
  current(tcPDPtr resonance, FlavourInfo flavour,
          const int imode, const int ichan,Energy & scale,
          const tPDVector & outgoing,
          const vector<Lorentz5Momentum> & momenta,
          DecayIntegrator::MEOption meopt) const;
 
  virtual bool accept(vector<int> id);
 
  virtual unsigned int decayMode(vector<int> id);
 
  virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
  
  Complex Fkaon(Energy2 q2,const int imode, const int ichan,
                IsoSpin::IsoSpin Itotal,
                Strangeness::Strange strange,
                tcPDPtr resonance,
                Energy ma, Energy mb) const;
  
  Complex FkaonRemainderI1(Energy2 q2) const;
  
  Complex FkaonRemainderI0(Energy2 q2,Energy ma, Energy mb) 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 void doinit();
 
  void constructInterpolators() const;
 
private:
 
  TwoKaonCzyzCurrent & operator=(const TwoKaonCzyzCurrent &) = delete;
 
private:
 
  vector<Complex> rhoWgt_;
 
  vector<double> rhoMag_;
 
  vector<double> rhoPhase_;
 
  vector<Complex> omegaWgt_;
 
  vector<double> omegaMag_;
 
  vector<double> omegaPhase_;
 
  vector<Complex> phiWgt_;
 
  vector<double> phiMag_;
 
  vector<double> phiPhase_;
  
  vector<Energy> rhoMasses_;
 
  vector<Energy> rhoWidths_;
 
  vector<Energy> omegaMasses_;
 
  vector<Energy> omegaWidths_;
  
  vector<Energy> phiMasses_;
 
  vector<Energy> phiWidths_;
 
  double betaRho_;
 
  double betaOmega_;
 
  double betaPhi_;
  
  unsigned int nMax_;
 
  double etaPhi_;
 
  double gammaOmega_;
 
  double gammaPhi_;
  
  vector<vector<Energy> > mass_;
 
  vector<vector<Energy> > width_;
 
  vector<vector<Complex> > coup_;
 
  vector<double> dh_;
 
  vector<Energy2> hres_;
 
  vector<Energy2> h0_;
 
  Energy mpi_;
 
  Energy eMax_;
 
  mutable Interpolator<double,Energy2>::Ptr fKI0Re_, fKI0Im_,fKI1Re_, fKI1Im_;
  
};
 
}
 
#endif /* Herwig_TwoKaonCzyzCurrent_H */