Models/Susy/SSWGSSVertex.h

0001 // -*- C++ -*-
0002 //
0003 // SSWGSSVertex.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_SSWGSSVertex_H
0010 #define HERWIG_SSWGSSVertex_H
0011 //
0012 // This is the declaration of the SSWGSSVertex class.
0013 //
0014
0015 #include "ThePEG/Helicity/Vertex/Scalar/VVSSVertex.h"
0016 //#include "SusyBase.h"
0017 #include "MSSM.h"
0018
0019 namespace Herwig {
0020 using namespace ThePEG;
0021
0022 /**
0023  * This is the implementation of the 4 point gluon-gluon-squark-squark
0024  * vertex. It inherits from VVSSVertex and implements the setCouopling
0025  * method.
0026  *
0027  * @see VVSSVertex
0028  */
0029 class SSWGSSVertex: public VVSSVertex {
0030
0031 public:
0032   /** @name Functions used by the persistent I/O system. */
0033   //@{
0034   /**
0035    * Function used to write out object persistently.
0036    * @param os the persistent output stream written to.
0037    */
0038   void persistentOutput(PersistentOStream & os) const;
0039
0040   /**
0041    * Function used to read in object persistently.
0042    * @param is the persistent input stream read from.
0043    * @param version the version number of the object when written.
0044    */
0045   void persistentInput(PersistentIStream & is, int version);
0046   //@}
0047
0048   /**
0049    * The default constructor.
0050    */
0051   SSWGSSVertex();
0052
0053   /**
0054    * Calculate the couplings.
0055    * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
0056    * @param part1 The ParticleData pointer for the first  particle.
0057    * @param part2 The ParticleData pointer for the second particle.
0058    * @param part3 The ParticleData pointer for the third  particle.
0059    * @param part4 The ParticleData pointer for the fourth particle.
0060   */
0061   virtual void setCoupling(Energy2 q2,tcPDPtr part1,
0062                tcPDPtr part2,tcPDPtr part3,tcPDPtr part4);
0063
0064 public:
0065
0066   /**
0067    * The standard Init function used to initialize the interfaces.
0068    * Called exactly once for each class by the class description system
0069    * before the main function starts or
0070    * when this class is dynamically loaded.
0071    */
0072   static void Init();
0073
0074 protected:
0075
0076   /** @name Clone Methods. */
0077   //@{
0078   /**
0079    * Make a simple clone of this object.
0080    * @return a pointer to the new object.
0081    */
0082   virtual IBPtr clone() const {return new_ptr(*this);}
0083
0084   /** Make a clone of this object, possibly modifying the cloned object
0085    * to make it sane.
0086    * @return a pointer to the new object.
0087    */
0088   virtual IBPtr fullclone() const {return new_ptr(*this);}
0089   //@}
0090
0091 protected:
0092
0093   /** @name Standard Interfaced functions. */
0094   //@{
0095   /**
0096    * Initialize this object after the setup phase before saving an
0097    * EventGenerator to disk.
0098    * @throws InitException if object could not be initialized properly.
0099    */
0100   virtual void doinit();
0101   //@}
0102
0103 private:
0104
0105   /**
0106    * The assignment operator is private and must never be called.
0107    * In fact, it should not even be implemented.
0108    */
0109   SSWGSSVertex & operator=(const SSWGSSVertex &) = delete;
0110
0111 private:
0112
0113   /**
0114    * Value of \f$sin(\theta_w)\f$
0115    */
0116   double _sw;
0117
0118   /**
0119    * Value of \f$cos(\theta_w)\f$
0120    */
0121   double _cw;
0122
0123   /**
0124    * Scale at which the coupling was last evaluated
0125    */
0126   Energy2 _q2last;
0127
0128   /**
0129    * Value of em coupling when last evaluated
0130    */
0131   Complex _emcouplast;
0132
0133   /**
0134    * Value of strong coupling when last evaluated
0135    */
0136   Complex _scouplast;
0137
0138   /**
0139    * Stop mixing matrix
0140    */
0141   tMixingMatrixPtr _stop;
0142
0143   /**
0144    * Sbottom mixing matrix
0145    */
0146   tMixingMatrixPtr _sbottom;
0147
0148   /**
0149    * The up type squark present when the vertex was evaluated.
0150    */
0151   long _ulast;
0152
0153   /**
0154    * The down type squark present when the vertex was evaluated.
0155    */
0156   long _dlast;
0157
0158   /**
0159    * The gauge boson present when the vertex was last evaluated.
0160    */
0161   long _gblast;
0162
0163   /**
0164    * The value of the mixing matrix dependent part when the vertex was
0165    * last evaluated
0166    */
0167   Complex _factlast;
0168
0169 };
0170 }
0171
0172 #endif /* HERWIG_SSWGSSVertex_H */