FourPionNovosibirskCurrent.h

// -*- C++ -*-
//
// FourPionNovosibirskCurrent.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_FourPionNovosibirskCurrent_H
#define HERWIG_FourPionNovosibirskCurrent_H
//
// This is the declaration of the FourPionNovosibirskCurrent class.
//
#include "WeakCurrent.h"
#include "Herwig/Utilities/Interpolator.h"
#include "Herwig/Utilities/Kinematics.h"
 
namespace Herwig {
using namespace ThePEG;
 
class FourPionNovosibirskCurrent: public WeakCurrent {
 
  friend class FourPionDefaultMatrixElement;
 
public:
 
  FourPionNovosibirskCurrent();
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
public:
 
  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;
 
  double threeBodyMatrixElement(const int iopt, const Energy2 q2,
                                const Energy2 s3, const Energy2 s2, 
                                const Energy2 s1, const Energy  m1,
                                const Energy  m2, const Energy  m3) const;
  
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
  
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
  virtual void doupdate();
 
private:
 
  FourPionNovosibirskCurrent & operator=(const FourPionNovosibirskCurrent &) = delete;
      
protected:
  
 
  void inita1width(int iopt);
  
  double a1FormFactor(Energy2 q2) const {
    return sqr((1.+_a1massolam2)/(1.+q2*_onedlam2));   
  }
 
  Complex sigmaBreitWigner(Energy2 q2,unsigned int iopt) const;
 
  Complex a1BreitWigner(Energy2 q2) const;
 
  Complex omegaBreitWigner(Energy2 q2) const;
 
  Complex rhoBreitWigner(Energy2 q2) const;
 
  Energy a1width(Energy2 q2) const {return (*_a1runinter)(q2);}
 
  LorentzVector<complex<Energy5> > 
  t1(Lorentz5Momentum & q1,Lorentz5Momentum & q2,
     Lorentz5Momentum & q3,Lorentz5Momentum & q4) const;
 
  LorentzVector<complex<Energy5> > 
  t2(Lorentz5Momentum & q1,Lorentz5Momentum & q2,
     Lorentz5Momentum & q3,Lorentz5Momentum & q4,
     unsigned int iopt) const;
 
  LorentzVector<complex<Energy5> > 
  t3(Lorentz5Momentum & q1,Lorentz5Momentum & q2,
     Lorentz5Momentum & q3,Lorentz5Momentum & q4) const;
 
  InvEnergy6 gFunction(Energy2 q2, int ichan) const;
 
  Energy2 DParameter() const;
 
  double dhdq2Parameter() const;
 
  Energy2 hFunction(const Energy q) const;
 
private:
  
  Interpolator<double,Energy>::Ptr _Fomega;
 
  Interpolator<double,Energy>::Ptr _Fthreec;
 
  Interpolator<double,Energy>::Ptr _Fonec;
 
  Interpolator<double,Energy2>::Ptr _Fsigma;
 
  Energy _mpic;
 
  Energy _mpi0;
 
  Energy _rhomass;
 
  Energy _a1mass;
 
  Energy _omegamass;
 
  Energy _sigmamass;
 
  Energy _rhowidth;
 
  Energy _a1width;
 
  Energy _omegawidth;
 
  Energy _sigmawidth;
 
  Energy _intmass;
 
  Energy _intwidth; 
 
  Complex _zsigma;
 
  double _zmag;
 
  double _zphase;
 
  Energy2 _lambda2;
 
  InvEnergy2 _onedlam2;
 
  double _a1massolam2;
 
  vector<Energy> _psigma;
 
  Energy2 _mpic2;
 
  Energy2 _mpi02;
 
  Energy2 _hm2;
 
  Energy2 _rhoD;
 
  Energy _prho;
 
  double _dhdq2m2;
 
  InvEnergy _aomega;
 
  InvEnergy _athreec;
 
  InvEnergy _aonec;
 
  double _bomega;
 
  double _bthreec;
 
  double _bonec;
 
  double _comega;
 
  double _cthreec;
 
  double _conec;
 
  vector<double> _omegaparam;
 
  bool _initializea1;
 
  bool _localparameters;
 
  vector<Energy> _a1runwidth;
 
  vector<Energy2> _a1runq2;
 
  Interpolator<Energy,Energy2>::Ptr _a1runinter;
 
  Energy _maxmass;
 
  Energy _maxcalc;
};
 
}
 
#endif /* HERWIG_FourPionNovosibirskCurrent_H */