MSSM.h

// -*- C++ -*-
//
// MSSM.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_MSSM_H
#define HERWIG_MSSM_H
//
// This is the declaration of the MSSM class.
//
 
#include "SusyBase.h"
#include "MSSM.fh"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MSSM: public SusyBase {
 
public:
 
  MSSM() : createDiagonalMixing_(false) {}
 
public:
 
  double higgsMixingAngle() const {return theAlpha;}
 
  const complex<Energy> & topTrilinear() const {return theAtop;}
 
  const complex<Energy> & bottomTrilinear() const {return theAbottom;}
 
  const complex<Energy> & tauTrilinear() const {return theAtau;}
 
  const MixingMatrixPtr & stopMix() const {return theStopMix;}
 
  const MixingMatrixPtr & sbottomMix() const {return theSbotMix;}
 
  const MixingMatrixPtr & stauMix() const {return theStauMix;}
 
  const MixingMatrixPtr & CPevenHiggsMix() const {return theHiggsMix;}
 
  const MixingMatrixPtr & CPoddHiggsMix() const {
    return HiggsAMix_;
  }
 
  const MixingMatrixPtr & ChargedHiggsMix() const {
    return HiggsPMix_;
  }
 
public:
 
  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);}
 
protected:
 
  virtual void createMixingMatrices();
 
  virtual void extractParameters(bool checkModel=true);
 
  virtual void adjustMixingMatrix(long id);
 
  void CPevenHiggsMix(MixingMatrixPtr in) {theHiggsMix = in;}
 
  void CPoddHiggsMix(MixingMatrixPtr in) {HiggsAMix_ = in;}
 
  void ChargedHiggsMix(MixingMatrixPtr in) {HiggsPMix_ = in;}
 
  void stopMix(MixingMatrixPtr in) {theStopMix =in;}
 
  void sbottomMix(MixingMatrixPtr in) {theSbotMix = in;}
 
  void stauMix(MixingMatrixPtr in) {theStauMix = in;}
 
  void higgsMixingAngle(double in) {theAlpha = in;}
 
private:
 
  MSSM & operator=(const MSSM &) = delete;
 
private:
 
  MixingMatrixPtr theStopMix;
 
  MixingMatrixPtr theSbotMix;
 
  MixingMatrixPtr theStauMix;
 
  complex<Energy> theAtop;
 
  complex<Energy> theAbottom;
 
  complex<Energy> theAtau;
 
  double theAlpha;
 
  MixingMatrixPtr theHiggsMix;
  MixingMatrixPtr HiggsAMix_;
 
  MixingMatrixPtr HiggsPMix_;
 
  bool createDiagonalMixing_;
 
};
 
}
 
#endif /* HERWIG_MSSM_H */