Models/Susy/SSNCTVertex.h

0001 // -*- C++ -*-
0002 #ifndef HERWIG_SSNCTVertex_H
0003 #define HERWIG_SSNCTVertex_H
0004 //
0005 // This is the declaration of the SSNCTVertex class.
0006 //
0007
0008 #include "ThePEG/Helicity/Vertex/Scalar/FFSVertex.h"
0009 #include "MSSM.h"
0010
0011 namespace Herwig {
0012
0013 using namespace ThePEG;
0014
0015 /**
0016  * Here is the documentation of the SSNCTVertex class.
0017  *
0018  * @see \ref SSNCTVertexInterfaces "The interfaces"
0019  * defined for SSNCTVertex.
0020  */
0021 class SSNCTVertex: public Helicity::FFSVertex {
0022
0023 public:
0024
0025   /**
0026    * The default constructor.
0027    */
0028   SSNCTVertex();
0029
0030 public:
0031
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 standard Init function used to initialize the interfaces.
0050    * Called exactly once for each class by the class description system
0051    * before the main function starts or
0052    * when this class is dynamically loaded.
0053    */
0054   static void Init();
0055
0056   /**
0057    * Calculate the couplings.
0058    * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
0059    * @param part1 The ParticleData pointer for the first  particle.
0060    * @param part2 The ParticleData pointer for the second particle.
0061    * @param part3 The ParticleData pointer for the third  particle.
0062    */
0063   virtual void setCoupling(Energy2 q2, tcPDPtr part1,
0064                            tcPDPtr part2, tcPDPtr part3);
0065
0066 protected:
0067
0068   /** @name Clone Methods. */
0069   //@{
0070   /**
0071    * Make a simple clone of this object.
0072    * @return a pointer to the new object.
0073    */
0074   virtual IBPtr clone() const;
0075
0076   /** Make a clone of this object, possibly modifying the cloned object
0077    * to make it sane.
0078    * @return a pointer to the new object.
0079    */
0080   virtual IBPtr fullclone() const;
0081   //@}
0082
0083 protected:
0084
0085   /** @name Standard Interfaced functions. */
0086   //@{
0087   /**
0088    * Initialize this object after the setup phase before saving and
0089    * EventGenerator to disk.
0090    * @throws InitException if object could not be initialized properly.
0091    */
0092   virtual void doinit();
0093   //@}
0094
0095 private:
0096
0097   /**
0098    * The assignment operator is private and must never be called.
0099    * In fact, it should not even be implemented.
0100    */
0101   SSNCTVertex & operator=(const SSNCTVertex &) = delete;
0102
0103 private:
0104
0105   /**
0106    *  High scale for the loops
0107    */
0108   Energy MX_;
0109
0110   /**
0111    * Pointer to the neutralino mixing matrix
0112    */
0113   tMixingMatrixPtr nmix_;
0114
0115   /**
0116    * \f$\sin(\theta_w)\f$
0117    */
0118   double sw_;
0119
0120   /**
0121    * \f$\cos(\theta_w)\f$
0122    */
0123   double cw_;
0124
0125   /**
0126    * Mass of the W
0127    */
0128   Energy mw_;
0129
0130   /**
0131    * \f$\sin(\beta)\f$
0132    */
0133   double sb_;
0134
0135   /**
0136    * \f$\cos(\beta)\f$
0137    */
0138   double cb_;
0139
0140   /**
0141    * The scale at which the coupling was last evaluated.
0142    */
0143   Energy2 q2last_;
0144
0145   /**
0146    * The value of the normalisation when it was evaluated at _q2last
0147    */
0148   Complex couplast_;
0149
0150   /**
0151    * Store the value of the left coupling when it was last evaluated
0152    */
0153   Complex leftlast_;
0154
0155   /**
0156    * Store the value of the right coupling when it was last evaluated
0157    */
0158   Complex rightlast_;
0159
0160   /**
0161    * Store the id of the last neutralino to be evaluate
0162    */
0163   long idlast_;
0164
0165   /**
0166    *  Mixing parameter
0167    */
0168   Complex epsilon_;
0169
0170 };
0171
0172 }
0173
0174 #endif /* HERWIG_SSNCTVertex_H */