Models/Susy/SSWWHVertex.h

0001 // -*- C++ -*-
0002 //
0003 // SSWWHVertex.h is a part of Herwig - A multi-purpose Monte Carlo event generator
0004 // Copyright (C) 2002-2019 The Herwig Collaboration
0005 //
0006 // Herwig is licenced under version 3 of the GPL, see COPYING for details.
0007 // Please respect the MCnet academic guidelines, see GUIDELINES for details.
0008 //
0009 #ifndef HERWIG_SSWWHVertex_H
0010 #define HERWIG_SSWWHVertex_H
0011 //
0012 // This is the declaration of the SSWWHVertex class.
0013 //
0014
0015 #include "ThePEG/Helicity/Vertex/Scalar/VVSVertex.h"
0016 #include "MSSM.h"
0017
0018 namespace Herwig {
0019 using namespace ThePEG;
0020
0021 /**
0022  * This class implements the coupling of a higgs in the MSSM
0023  * to a pair of SM gauge bosons.
0024  *
0025  * @see \ref SSWWHVertexInterfaces "The interfaces"
0026  * defined for SSWWHVertex.
0027  */
0028 class SSWWHVertex: public VVSVertex {
0029
0030 public:
0031
0032   /**
0033    * The default constructor.
0034    */
0035   SSWWHVertex();
0036
0037 public:
0038
0039   /** @name Functions used by the persistent I/O system. */
0040   //@{
0041   /**
0042    * Function used to write out object persistently.
0043    * @param os the persistent output stream written to.
0044    */
0045   void persistentOutput(PersistentOStream & os) const;
0046
0047   /**
0048    * Function used to read in object persistently.
0049    * @param is the persistent input stream read from.
0050    * @param version the version number of the object when written.
0051    */
0052   void persistentInput(PersistentIStream & is, int version);
0053   //@}
0054
0055   /**
0056    * The standard Init function used to initialize the interfaces.
0057    * Called exactly once for each class by the class description system
0058    * before the main function starts or
0059    * when this class is dynamically loaded.
0060    */
0061   static void Init();
0062
0063   /**
0064    * Calculate the coupling for the vertex
0065    * @param q2 The scale to at which evaluate the coupling.
0066    * @param particle1 The first particle in the vertex.
0067    * @param particle2 The second particle in the vertex.
0068    * @param particle3 The third particle in the vertex.
0069    */
0070   virtual void setCoupling(Energy2 q2, tcPDPtr particle1, tcPDPtr particle2,
0071                tcPDPtr particle3);
0072
0073 protected:
0074
0075   /** @name Clone Methods. */
0076   //@{
0077   /**
0078    * Make a simple clone of this object.
0079    * @return a pointer to the new object.
0080    */
0081   virtual IBPtr clone() const {return new_ptr(*this);}
0082
0083   /** Make a clone of this object, possibly modifying the cloned object
0084    * to make it sane.
0085    * @return a pointer to the new object.
0086    */
0087   virtual IBPtr fullclone() const {return new_ptr(*this);}
0088   //@}
0089
0090 protected:
0091
0092   /** @name Standard Interfaced functions. */
0093   //@{
0094   /**
0095    * Initialize this object after the setup phase before saving an
0096    * EventGenerator to disk.
0097    * @throws InitException if object could not be initialized properly.
0098    */
0099   virtual void doinit();
0100   //@}
0101
0102 private:
0103
0104   /**
0105    * The assignment operator is private and must never be called.
0106    * In fact, it should not even be implemented.
0107    */
0108   SSWWHVertex & operator=(const SSWWHVertex &) = delete;
0109
0110 private:
0111
0112
0113   /**
0114    * The value of the factor \f$\frac{m_W \sin(\beta-\alpha)}{\sin\theta_W}\f$
0115    */
0116   Energy theh0Wfact;
0117
0118   /**
0119    * The value of the factor \f$\frac{m_W \cos(\beta-\alpha)}{\sin\theta_W}\f$
0120    */
0121   Energy theH0Wfact;
0122
0123   /**
0124    * The value of the factor
0125    * \f$\frac{m_Z\sin(\beta-\alpha)}{\sin\theta_W\cos\theta_W}\f$
0126    */
0127   Energy theh0Zfact;
0128
0129   /**
0130    * The value of the factor
0131    * \f$\frac{m_Z\cos(\beta-\alpha)}{\sin\theta_W\cos\theta_W}\f$
0132    */
0133   Energy theH0Zfact;
0134
0135   /**
0136    * The value of the coupling when it was last evaluated
0137    */
0138   complex<Energy> theCoupLast;
0139
0140   /**
0141    * The value of \f$\sqrt{4\pi \alpha(q2)}\f$ when it was last evaluated.
0142    */
0143   Complex theElast;
0144
0145
0146   /**
0147    * The scale at which the coupling was last evaluated
0148    */
0149   Energy2 theq2last;
0150
0151   /**
0152    * The ID of the last higgs for which the vertex was evaluated
0153    */
0154   long theHlast;
0155
0156   /**
0157    * The ID of the last gauge boson for which the vertex was evaluated
0158    */
0159   long theGBlast;
0160 };
0161 }
0162
0163
0164 #endif /* HERWIG_SSWWHVertex_H */