SSHSFSFVertex.h

// -*- C++ -*-
//
// SSHSFSFVertex.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_SSHSFSFVertex_H
#define HERWIG_SSHSFSFVertex_H
//
// This is the declaration of the SSHSFSFVertex class.
//
 
#include "ThePEG/Helicity/Vertex/Scalar/SSSVertex.h"
#include "MSSM.h"
 
namespace Herwig {
using namespace ThePEG;
 
class SSHSFSFVertex: public SSSVertex {
 
  typedef vector<MixingMatrixPtr> MMPVector;
 
public:
 
  SSHSFSFVertex();
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
  
  void setCoupling(Energy2 q2, tcPDPtr particle1, tcPDPtr particle2, 
                   tcPDPtr particle3);
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
  
protected:
 
  virtual void doinit();
 
private:
 
  SSHSFSFVertex & operator=(const SSHSFSFVertex &) = delete;
 
private:
 
  void downSF(Energy2 q2, long higgs, long smID, unsigned int alpha, unsigned int beta);
  
  void upSF(Energy2 q2, long higgs, long smID, unsigned int alpha, unsigned int beta);
  
  void leptonSF(Energy2 q2, long higgs, long smID, unsigned int alpha, unsigned int beta);
  
  void chargedHiggs(Energy2 q2, long id1, long id2);
  
 
private:
  
  MMPVector theMix;
 
  vector<complex<Energy> > theTriC;
    
  double theSinA;
 
  double theCosA;
 
  double theSinB;
 
  double theCosB;
 
  double theTanB;
 
  double theSinAB;
 
  double theCosAB;
 
  Energy theMw;
 
  Energy theMz;
  
  Energy theMu;
 
  double theSw;
 
  double theCw;
 
  complex<Energy> theCoupLast;
  
  Energy2 theq2Last;
  
  double thegLast;
  
  long theHLast;
  
  long theSF1Last;
  
  long theSF2Last;
 
  tMSSMPtr theMSSM;
};
}
 
 
#endif /* HERWIG_SSHSFSFVertex_H */