BaryonThreeQuarkModelFormFactor.h

// -*- C++ -*-
#ifndef HERWIG_BaryonThreeQuarkModelFormFactor_H
#define HERWIG_BaryonThreeQuarkModelFormFactor_H
//
// This is the declaration of the BaryonThreeQuarkModelFormFactor class.
//
#include "BaryonFormFactor.h"
#include "ThePEG/PDT/ParticleData.h"
 
namespace Herwig {
using namespace ThePEG;
 
class BaryonThreeQuarkModelFormFactor: public BaryonFormFactor {
  
public:
 
  BaryonThreeQuarkModelFormFactor();
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
public:
 
  virtual void SpinHalfSpinHalfFormFactor(Energy2 q2,int iloc, int id0, int id1,
                                          Energy m0, Energy m1,
                                          Complex & f1v,Complex & f2v,Complex & f3v,
                                          Complex & f1a,Complex & f2a,Complex & f3a,
                                          FlavourInfo flavour,
                                          Virtuality virt=SpaceLike);
 
  virtual void SpinHalfSpinThreeHalfFormFactor(Energy2 q2,int iloc, int id0, int id1,
                                               Energy m0, Energy m1,
                                               Complex & g1v,Complex & g2v,Complex & g3v,
                                               Complex & g4v,Complex & g1a,Complex & g2a,
                                               Complex & g3a,Complex & g4a,
                                               FlavourInfo flavour,
                                               Virtuality virt=SpaceLike);
 
  virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
 
  double operator ()(double x) const;
  typedef double ArgType;
  typedef double ValType;
 
  Energy widthIntegrand(double omega,Energy m0, Energy m1, int type, int imass,
                        int id0,int id1);
 
protected:
 
  vector<double> phiFunction(double);
 
  void SN(double y,int N,double & SNm2,double & SN) const;
 
protected:
 
 
  virtual void doinit();
 
private:
 
  BaryonThreeQuarkModelFormFactor & operator=(const BaryonThreeQuarkModelFormFactor &) = delete;
 
private:
 
  bool _initialize;
 
  unsigned int _order;
 
  Energy _mlight;
 
  Energy _mstrange;
 
  Energy _LambdaQ;
 
  Energy _Lambdaqq;
 
  Energy _Lambdasq;
 
 
  Energy _Lambdass;
 
  vector<double> _C0;
 
  vector<double> _C1;
 
  vector<double> _C2;
 
  double _a;
 
  double _b;
 
  double _mu2;
 
  int _N;
};
 
}
 
#endif /* HERWIG_BaryonThreeQuarkModelFormFactor_H */