RPV.h

// -*- C++ -*-
#ifndef HERWIG_RPV_H
#define HERWIG_RPV_H
//
// This is the declaration of the RPV class.
//
 
#include "Herwig/Models/Susy/MSSM.h"
#include "RPV.fh"
 
namespace Herwig {
 
using namespace ThePEG;
 
class RPV: public MSSM {
 
public:
 
  RPV() : lambdaLLE_(3,vector<vector<double> >(3,vector<double>(3,0.))),
          lambdaLQD_(3,vector<vector<double> >(3,vector<double>(3,0.))),
          lambdaUDD_(3,vector<vector<double> >(3,vector<double>(3,0.))),
          triLinearOnly_(false)
  {}
 
 
  const vector<vector<vector<double> > > & lambdaLLE() {return lambdaLLE_;}
 
  const vector<vector<vector<double> > > & lambdaLQD() {return lambdaLQD_;}
 
  const vector<vector<vector<double> > > & lambdaUDD() {return lambdaUDD_;}
 
  const vector<Energy> & sneutrinoVEVs() {return vnu_;}
 
  const vector<Energy> & epsilon() {return epsilon_;}
 
  const vector<Energy> & epsilonB() {return epsB_;}
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
  
protected:
 
  virtual void doinit();
 
protected:
 
  virtual void extractParameters(bool checkModel=true);
 
  virtual void createMixingMatrices();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  RPV & operator=(const RPV &) = delete;
 
private:
 
  //@
  vector<vector<vector<double> > > lambdaLLE_;
 
  vector<vector<vector<double> > > lambdaLQD_;
 
  vector<vector<vector<double> > > lambdaUDD_;
  //@
 
  vector<Energy> vnu_;
 
  vector<Energy> epsilon_;
 
  vector<Energy> epsB_;
 
  MixingMatrixPtr upSquarkMix_;
 
  MixingMatrixPtr downSquarkMix_;
 
  AbstractFFSVertexPtr LLEVertex_;
 
  AbstractFFSVertexPtr LQDVertex_;
 
  AbstractFFSVertexPtr UDDVertex_;
 
  bool triLinearOnly_;
};
 
}
 
#endif /* HERWIG_RPV_H */