ThreeMesonCurrentBase.h

// -*- C++ -*-
//
// ThreeMesonCurrentBase.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_ThreeMesonCurrentBase_H
#define HERWIG_ThreeMesonCurrentBase_H
// This is the declaration of the ThreeMesonCurrentBase class.
 
#include "WeakDecayCurrent.h"
 
namespace Herwig {
using namespace ThePEG;
 
class ThreeMesonCurrentBase: public WeakDecayCurrent {
 
public:
 
  ThreeMesonCurrentBase();
 
  static void Init();
 
public:
 
 
  virtual vector<LorentzPolarizationVectorE> 
  current(const int imode,const int ichan,Energy & scale,
          const ParticleVector & decay,DecayIntegrator::MEOption meopt) const;
 
  virtual bool accept(vector<int> id);
 
  virtual unsigned int decayMode(vector<int> id);
 
  virtual tPDVector particles(int icharge, unsigned int imode, int iq, int ia);
 
  virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
 
protected:
 
  virtual bool acceptMode(int imode) const=0;
 
  struct FormFactors {
 
    complex<InvEnergy>  F1;
    
    complex<InvEnergy>  F2;
    
    complex<InvEnergy>  F3; 
    
    complex<InvEnergy>  F4;
    
    complex<InvEnergy3> F5;
 
    FormFactors(complex<InvEnergy>  f1 = InvEnergy(), 
                complex<InvEnergy>  f2 = InvEnergy(),
                complex<InvEnergy>  f3 = InvEnergy(),
                complex<InvEnergy>  f4 = InvEnergy(),
                complex<InvEnergy3> f5 = InvEnergy3())
      : F1(f1), F2(f2), F3(f3), F4(f4), F5(f5) {}
  };
 
  virtual FormFactors calculateFormFactors(const int ichan, const int imode,
                                           Energy2 q2,
                                           Energy2 s1, Energy2 s2, Energy2 s3) const = 0;
 
private:
 
  ThreeMesonCurrentBase & operator=(const ThreeMesonCurrentBase &) = delete;
 
};
 
}
 
#endif /* HERWIG_ThreeMesonCurrentBase_H */