Susy/RPV/RPVWSSVertex.h

0001 // -*- C++ -*-
0002 #ifndef Herwig_RPVWSSVertex_H
0003 #define Herwig_RPVWSSVertex_H
0004 //
0005 // This is the declaration of the RPVWSSVertex class.
0006 //
0007
0008 #include "ThePEG/Helicity/Vertex/Scalar/VSSVertex.h"
0009 #include "Herwig/Models/Susy/MixingMatrix.h"
0010
0011 namespace Herwig {
0012
0013 using namespace ThePEG;
0014
0015 /**
0016  * Here is the documentation of the RPVWSSVertex class.
0017  *
0018  * @see \ref RPVWSSVertexInterfaces "The interfaces"
0019  * defined for RPVWSSVertex.
0020  */
0021 class RPVWSSVertex: public Helicity::VSSVertex {
0022
0023 public:
0024
0025   /**
0026    * The default constructor.
0027    */
0028   RPVWSSVertex();
0029
0030  /**
0031   * Calculate the couplings.
0032   * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
0033   * @param part1 The ParticleData pointer for the first  particle.
0034   * @param part2 The ParticleData pointer for the second particle.
0035   * @param part3 The ParticleData pointer for the third  particle.
0036   */
0037   virtual void setCoupling(Energy2 q2,tcPDPtr part1,
0038                            tcPDPtr part2,tcPDPtr part3);
0039
0040 public:
0041
0042   /** @name Functions used by the persistent I/O system. */
0043   //@{
0044   /**
0045    * Function used to write out object persistently.
0046    * @param os the persistent output stream written to.
0047    */
0048   void persistentOutput(PersistentOStream & os) const;
0049
0050   /**
0051    * Function used to read in object persistently.
0052    * @param is the persistent input stream read from.
0053    * @param version the version number of the object when written.
0054    */
0055   void persistentInput(PersistentIStream & is, int version);
0056   //@}
0057
0058   /**
0059    * The standard Init function used to initialize the interfaces.
0060    * Called exactly once for each class by the class description system
0061    * before the main function starts or
0062    * when this class is dynamically loaded.
0063    */
0064   static void Init();
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 an
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   RPVWSSVertex & operator=(const RPVWSSVertex &) = delete;
0102
0103 private:
0104
0105   /**
0106    * Value of \f$sin(\theta_w)\f$
0107    */
0108   double _sw;
0109
0110   /**
0111    * Value of \f$cos(\theta_w)\f$
0112    */
0113   double _cw;
0114
0115   /**
0116    * Stau mixing matrix
0117    */
0118   tMixingMatrixPtr _stau;
0119
0120   /**
0121    * Stop mixing matrix
0122    */
0123   tMixingMatrixPtr _stop;
0124
0125   /**
0126    * Sbottom mixing matrix
0127    */
0128   tMixingMatrixPtr _sbottom;
0129
0130   /**
0131    * The value of \f$\sin2\theta_W\f$
0132    */
0133   double _s2w;
0134
0135   /**
0136    * The value of \f$\cos2\theta_W\f$
0137    */
0138   double _c2w;
0139
0140   /**
0141    *  Coupling of Z to scalar and pseudoscalar
0142    */
0143   vector<vector<Complex> > Cijeo_;
0144
0145   /**
0146    *  Coupling of W to scalar charged
0147    */
0148   vector<vector<Complex> > Cijec_;
0149
0150   /**
0151    *  Coupling of W to pseudoscalar charged
0152    */
0153   vector<vector<Complex> > Cijco_;
0154
0155   /**
0156    *  Coupling of Z to charged Higgs
0157    */
0158   vector<vector<Complex> > Cijc_;
0159
0160   /**
0161    *  Which interactions to include
0162    */
0163   unsigned int _interactions;
0164
0165   /**
0166    * Scale at which the coupling was last evaluated
0167    */
0168   Energy2 _q2last;
0169
0170    /**
0171     * The up type sfermion present when the vertex was evaluated.
0172     */
0173    long _ulast;
0174
0175    /**
0176     * The down type sfermion present when the vertex was evaluated.
0177     */
0178    long _dlast;
0179
0180    /**
0181     * The gauge boson present when the vertex was last evaluated.
0182     */
0183    long _gblast;
0184
0185   /**
0186    * The value of the mixing matrix dependent part when the vertex was
0187    * last evaluated
0188    */
0189   Complex _factlast;
0190
0191   /**
0192    * Value of coupling when last evaluated
0193    */
0194   Complex _couplast;
0195 };
0196
0197 }
0198
0199 #endif /* Herwig_RPVWSSVertex_H */