Models/Susy/SSNFSVertex.h

0001 // -*- C++ -*-
0002 //
0003 // SSNFSVertex.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_SSNFSVertex_H
0010 #define HERWIG_SSNFSVertex_H
0011 //
0012 // This is the declaration of the SSNFSVertex class.
0013 //
0014
0015 #include "ThePEG/Helicity/Vertex/Scalar/FFSVertex.h"
0016 #include "ThePEG/StandardModel/StandardModelBase.h"
0017 #include "MSSM.h"
0018
0019 namespace Herwig {
0020 using namespace ThePEG;
0021
0022 /**
0023  * This is the implementation of the coupling of the neutralinos to
0024  * fermion-sfermions. It inherits from FFSVertex and implements the
0025  * virtual setCoupling() method.
0026  *
0027  * @see FFSVertex
0028  */
0029 class SSNFSVertex: public FFSVertex {
0030
0031 public:
0032
0033    /**
0034    * The default constructor.
0035    */
0036   SSNFSVertex();
0037
0038  public:
0039
0040   /** @name Functions used by the persistent I/O system. */
0041   //@{
0042   /**
0043    * Function used to write out object persistently.
0044    * @param os the persistent output stream written to.
0045    */
0046   void persistentOutput(PersistentOStream & os) const;
0047
0048   /**
0049    * Function used to read in object persistently.
0050    * @param is the persistent input stream read from.
0051    * @param version the version number of the object when written.
0052    */
0053   void persistentInput(PersistentIStream & is, int version);
0054   //@}
0055
0056   /**
0057    * The standard Init function used to initialize the interfaces.
0058    * Called exactly once for each class by the class description system
0059    * before the main function starts or
0060    * when this class is dynamically loaded.
0061    */
0062   static void Init();
0063
0064   /**
0065    * Calculate the couplings.
0066    * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
0067    * @param part1 The ParticleData pointer for the first  particle.
0068    * @param part2 The ParticleData pointer for the second particle.
0069    * @param part3 The ParticleData pointer for the third  particle.
0070    */
0071   virtual void setCoupling(Energy2 q2, tcPDPtr part1,
0072                            tcPDPtr part2, tcPDPtr part3);
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 and
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   SSNFSVertex & operator=(const SSNFSVertex &) = delete;
0110
0111   /**
0112    * Pointer to the stop mixing matrix
0113    */
0114   tMixingMatrixPtr _stop;
0115
0116   /**
0117    * Pointer to the sbottom mixing matrix
0118    */
0119   tMixingMatrixPtr _sbot;
0120
0121   /**
0122    * Pointer to the stau mixing matrix
0123    */
0124   tMixingMatrixPtr _stau;
0125
0126   /**
0127    * Pointer to the neutralino mixing matrix
0128    */
0129   tMixingMatrixPtr _nmix;
0130
0131   /**
0132    * Pointer to the Susy Model object
0133    */
0134   tMSSMPtr _theSS;
0135
0136   /**
0137    * \f$\sin(\theta_w)\f$
0138    */
0139   double _sw;
0140
0141   /**
0142    * \f$\cos(\theta_w)\f$
0143    */
0144   double _cw;
0145
0146   /**
0147    * Mass of the W
0148    */
0149   Energy _mw;
0150
0151   /**
0152    * \f$\sin(\beta)\f$
0153    */
0154   double _sb;
0155
0156   /**
0157    * \f$\cos(\beta)\f$
0158    */
0159   double _cb;
0160
0161   /**
0162    * The scale at which the coupling was last evaluated.
0163    */
0164   Energy2 _q2last;
0165
0166   /**
0167    * The value of the normalisation when it was evaluated at _q2last
0168    */
0169   Complex _couplast;
0170
0171   /**
0172    * Store the value of the left coupling when it was last evaluated
0173    */
0174   Complex _leftlast;
0175
0176   /**
0177    * Store the value of the right coupling when it was last evaluated
0178    */
0179   Complex _rightlast;
0180
0181   /**
0182    * Store the id of the last neutralino to be evaluate
0183    */
0184   long _id1last;
0185
0186   /**
0187    * Store the id of the last SM fermion to be evaluate
0188    */
0189   long _id2last;
0190
0191   /**
0192    * Store the id of the last scalar to be evaluate
0193    */
0194   long _id3last;
0195
0196   /**
0197    *  Include Yukawa's ?
0198    */
0199   unsigned int yukawa_;
0200 };
0201 }
0202
0203 #endif /* HERWIG_SSNFSVertex_H */