UEDBase.h

// -*- C++ -*-
//
// UEDBase.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_UEDBase_H
#define HERWIG_UEDBase_H
//
// This is the declaration of the UEDBase class.
//
 
#include "Herwig/Models/General/BSMModel.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVVVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVSSVertex.h"
 
#include "UEDBase.fh"
 
namespace Herwig {
using namespace ThePEG;
 
class UEDBase: public BSMModel {
 
public:
 
  typedef pair<long, Energy> IDMassPair;
 
  typedef map<unsigned int, double> WAMap;
 
public:
 
  UEDBase();
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
public:
 
  InvEnergy compactRadius() const {
    return 1./theInvRadius;
  }
 
  double sinThetaN(const unsigned int n) const;
 
  double sinThetaOne() const {
    return theSinThetaOne;
  }
 
protected:
 
  void addMassElement(IDMassPair elem) {
    theMasses.push_back(elem);
  }
 
  void addMixingAngle(const unsigned int n, 
                      const double val) {
    theMixingAngles.insert(make_pair(n, val));
  }
  
private:
 
  void calculateKKMasses(const unsigned int n);
 
  void bosonMasses(const unsigned int n);
  
  void fermionMasses(const unsigned int n);
 
  void resetMass(long id, Energy value);
 
  double calculateMixingAngle(const unsigned int n);
  
  void writeSpectrum();
 
  static bool lowerMass(const pair<long, Energy> & p1, 
                        const pair<long, Energy> & p2) {
    return p1.second < p2.second;
  }
  
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
 
protected:
 
  virtual void doinit();
 
private:
 
  UEDBase & operator=(const UEDBase &) = delete;
 
private:
  
  bool theRadCorr;
  
  Energy theInvRadius;
  
  double theLambdaR;
 
  Energy theMbarH;
 
  WAMap theMixingAngles;
 
  double theSinThetaOne;
  
  vector<IDMassPair> theMasses;
  
  Energy theVeV;
 
  bool includeSMMass_;
 
  bool fixedCouplings_;
 
  bool includeGaugeMixing_;
  
  AbstractFFVVertexPtr theF1F1Z0Vertex;
 
  AbstractFFVVertexPtr theF1F1G0Vertex;
 
  AbstractFFVVertexPtr theF1F0G1Vertex;
 
  AbstractVVVVertexPtr theG1G1G0Vertex;
 
  AbstractVVVVVertexPtr theG0G0G1G1Vertex;
 
  AbstractFFVVertexPtr theF1F1P0Vertex;
 
  AbstractFFVVertexPtr theF1F1W0Vertex;
 
  AbstractFFVVertexPtr theF1F0W1Vertex;
 
  AbstractFFSVertexPtr theF1F0H1Vertex;
 
  AbstractVSSVertexPtr theP0H1H1Vertex;
 
  AbstractVSSVertexPtr theZ0H1H1Vertex;
 
  AbstractVSSVertexPtr theW0A1H1Vertex;
 
  AbstractVSSVertexPtr theZ0A1h1Vertex;
  
  AbstractVVVVertexPtr theW0W1W1Vertex;
};
 
 
}
 
#endif /* HERWIG_UEDBase_H */