Models/Susy/SSNNZVertex.h

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