ISGWFormFactor.h

// -*- C++ -*-
//
// ISGWFormFactor.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_ISGWFormFactor_H
#define HERWIG_ISGWFormFactor_H
//
// This is the declaration of the <ISGWFormFactor class.
//
#include "ScalarFormFactor.h"
#include "ThePEG/StandardModel/StandardModelBase.h"
#include "ThePEG/PDT/ParticleData.h"
#include "ThePEG/PDT/EnumParticles.h"
#include "ThePEG/PDT/EnumParticles.h"
#include "ThePEG/Repository/EventGenerator.h"
 
namespace Herwig {
using namespace ThePEG;
 
class ISGWFormFactor: public ScalarFormFactor {
 
public:
 
  ISGWFormFactor();
 
  virtual void ScalarScalarFormFactor(Energy2 q2,unsigned int iloc,int id0,int id1,Energy m0,
                                      Energy m1,Complex & f0,Complex & fp) const;
 
  virtual void ScalarVectorFormFactor(Energy2 q2, unsigned int iloc, int id0, int id1,
                                      Energy m0, Energy m1, Complex & V,
                                      Complex & A0,Complex & A1,Complex & A2) const;
 
  virtual void ScalarTensorFormFactor(Energy2 q2,unsigned int iloc,int id0,int id1,
                                      Energy m0,
                                      Energy m1, complex<InvEnergy2> & h,
                                      Complex & k, complex<InvEnergy2> & bp,
                                      complex<InvEnergy2> & bm) const;
 
  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:
 
  void formFactor(Energy2 q2,unsigned int iloc,int id0,int id1,Energy m0,
                  Energy m1,Complex & f1,Complex & f2,
                  Complex & f3,Complex & f4) const;
  // general member to calculate all the form-factors
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
private:
 
  ISGWFormFactor & operator=(const ISGWFormFactor &) = delete;
 
private:
 
  double _kappa;
 
  Energy _mdown;
 
  Energy _mup;
 
  Energy _mstrange;
 
  Energy _mcharm;
 
  Energy _mbottom;
 
  vector<Energy> _mquark;
 
 
  Energy _betaSud;
 
  Energy _betaSus;
 
  Energy _betaSuc;
 
  Energy _betaSub;
 
  vector<vector<Energy> > _betaS;
 
  Energy _betaPud;
 
  Energy _betaPus;
 
  Energy _betaPuc;
 
  vector<vector<Energy> > _betaP;
 
  double _thetaeta;
};
 
}
 
#endif /* HERWIG_ISGWFormFactor_H */