KaonThreeMesonCurrent.h

// -*- C++ -*-
//
// KaonThreeMesonCurrent.h is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_KaonThreeMesonCurrent_H
#define HERWIG_KaonThreeMesonCurrent_H
//
// This is the declaration of the KaonThreeMesonCurrent class.
//

#include "ThreeMesonCurrentBase.h"

namespace Herwig {

using namespace ThePEG;

class KaonThreeMesonCurrent: public ThreeMesonCurrentBase {

public:

  KaonThreeMesonCurrent();

  virtual bool createMode(int icharge,unsigned int imode,DecayPhaseSpaceModePtr mode,
                          unsigned int iloc,unsigned int ires,
                          DecayPhaseSpaceChannelPtr phase,Energy upp);

  virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
  
  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;

protected:

  virtual bool acceptMode(int) const;

  virtual FormFactors calculateFormFactors(const int ichan,const int imode,
                                           Energy2 q2,Energy2 s1,Energy2 s2,Energy2 s3) 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();

  virtual void doinitrun();

  virtual void doupdate();

private:

  KaonThreeMesonCurrent & operator=(const KaonThreeMesonCurrent &);

private:

  Complex Trho1(Energy2 q2,int ires) const;

  Complex Trho2(Energy2 q2,int ires) const;

  Complex TKstar1(Energy2 q2,int ires) const;

  Complex TKstar2(Energy2 q2,int ires) const;

  Complex BWrho1(Energy2 q2, unsigned int ires) const;

  Complex BWrho2(Energy2 q2, unsigned int ires) const;

  Complex BWKstar1(Energy2 q2, unsigned int ires) const;

  Complex BWKstar2(Energy2 q2, unsigned int ires) const;

  Complex a1BreitWigner(Energy2 q2) const;

  Complex TK1(Energy2 q2,unsigned int iopt,int ires) const;

  Complex K1BreitWigner(Energy2 q2,unsigned int ires) const;

  Energy a1Width(Energy2 q2) const;

  double g(Energy2 q2) const;

  void inita1Width(int iopt);

  Complex Tomega(Energy2 q2, int ires) const;

  Complex OmegaPhiBreitWigner(Energy2 q2, unsigned int ires) const;

  Complex TOmegaKStar(Energy2 s1,Energy2 s2,int ires) const;

private:

  vector<double> _rho1wgts;

  vector<Energy> _rho1mass;

  vector<Energy> _rho1width;

  vector<double> _rho2wgts;

  vector<Energy> _rho2mass;

  vector<Energy> _rho2width;

  vector<double> _kstar1wgts;

  vector<Energy> _kstar1mass;

  vector<Energy> _kstar1width;

  vector<double> _kstar2wgts;

  vector<Energy> _kstar2mass;

  vector<Energy> _kstar2width;

  Energy _a1mass;
  
  Energy _a1width;

  vector<Energy> _k1mass;

  vector<Energy> _k1width;

  vector<double> _k1wgta;

  vector<double> _k1wgtb;

  vector<Energy>  _a1runwidth;

  vector<Energy2> _a1runq2;

  Interpolator<Energy,Energy2>::Ptr _a1runinter;

  double _epsomega;

  Energy _omegamass;

  Energy _omegawidth;

  Energy _phimass;

  Energy _phiwidth;

  double _omegaKstarwgt;

  Energy _fpi;

  Energy _mpi;

  Energy _mK;

  bool _initializea1;
  
  bool _rhoparameters;

  bool _kstarparameters;

  bool _a1parameters;

  bool _k1parameters;

  bool _a1opt;

  bool _omegaopt;

  Energy _maxmass;

  Energy _maxcalc;
};

}

#endif /* HERWIG_KaonThreeMesonCurrent_H */