ThreePionCzyzCurrent.h

// -*- C++ -*-//
// ThreePionCzyzCurrent.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_ThreePionCzyzCurrent_H
#define Herwig_ThreePionCzyzCurrent_H
//
// This is the declaration of the ThreePionCzyzCurrent class.
//
 
#include "WeakCurrent.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
 
class ThreePionCzyzCurrent: public WeakCurrent {
 
public:
 
  ThreePionCzyzCurrent();
 
  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;
  
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();
 
private:
 
  ThreePionCzyzCurrent & operator=(const ThreePionCzyzCurrent &) = delete;
 
private:
 
  vector<Energy> rhoMasses_;
  vector<Energy> rhoWidths_;
  
  vector<Energy> omegaMasses_;
  vector<Energy> omegaWidths_;
  
  Energy phiMass_;
  Energy phiWidth_;
 
  vector<InvEnergy3> coup_I0_;
 
  vector<Energy> rhoMasses_I1_;
  vector<Energy> rhoWidths_I1_;
  
  Energy omegaMass_I1_;
  Energy omegaWidth_I1_;
  
  double sigma_;
 
  InvEnergy GW_pre_;
  
  Energy GW_;
 
  double g_omega_pi_pi_;
 
  Energy mpip_, mpi0_;
 
};
 
}
 
#endif /* Herwig_ThreePionCzyzCurrent_H */