Models/UED/UEDF1F0W1Vertex.h

0001 // -*- C++ -*-
0002 //
0003 // UEDF1F0W1Vertex.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_UEDF1F0W1Vertex_H
0010 #define HERWIG_UEDF1F0W1Vertex_H
0011 //
0012 // This is the declaration of the UEDF1F0W1Vertex class.
0013 //
0014
0015 #include "ThePEG/Helicity/Vertex/Vector/FFVVertex.h"
0016 #include "UEDBase.h"
0017
0018 namespace Herwig {
0019 using namespace ThePEG;
0020
0021 /**
0022  * This is the \f$ \bar{f}^{(1)} f^{(0)} W^{(1)}\f$ coupling, where
0023  * \f$W^{(1)}\f$ is any level 1 EW gauge boson. The class
0024  * inherits from FFVVertex and implements the setCoupling virtual member.
0025  *
0026  * @see \ref UEDF1F0W1VertexInterfaces "The interfaces"
0027  * defined for UEDF1F0W1Vertex.
0028  */
0029 class UEDF1F0W1Vertex: public FFVVertex {
0030
0031 public:
0032
0033   /**
0034    * The default constructor.
0035    */
0036   UEDF1F0W1Vertex();
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   /** Calculate the coupling
0065    *@param q2 The scale at which to evaluate the coupling
0066    *@param part1 The first interacting particle
0067    *@param part2 The second interacting particle
0068    *@param part3 The third interacting particle
0069    */
0070   virtual void setCoupling(Energy2 q2, tcPDPtr part1, tcPDPtr part2,
0071                tcPDPtr part3);
0072
0073 protected:
0074
0075   /** @name Clone Methods. */
0076   //@{
0077   /**
0078    * Make a simple clone of this object.
0079    * @return a pointer to the new object.
0080    */
0081   virtual IBPtr clone() const {return new_ptr(*this);}
0082
0083   /** Make a clone of this object, possibly modifying the cloned object
0084    * to make it sane.
0085    * @return a pointer to the new object.
0086    */
0087   virtual IBPtr fullclone() const {return new_ptr(*this);}
0088   //@}
0089
0090 protected:
0091
0092   /** @name Standard Interfaced functions. */
0093   //@{
0094   /**
0095    * Initialize this object after the setup phase before saving an
0096    * EventGenerator to disk.
0097    * @throws InitException if object could not be initialized properly.
0098    */
0099   virtual void doinit();
0100   //@}
0101
0102 private:
0103
0104   /**
0105    * The assignment operator is private and must never be called.
0106    * In fact, it should not even be implemented.
0107    */
0108   UEDF1F0W1Vertex & operator=(const UEDF1F0W1Vertex &) = delete;
0109
0110 private:
0111
0112   /**
0113    * The value of \f$\sin\theta_W\f$
0114    */
0115   double theSinW;
0116
0117   /**
0118    * The value of \f$\cos\theta_W\f$
0119    */
0120   double theCosW;
0121
0122   /**
0123    * The value of \f$\sin\theta_1\f$
0124    */
0125   double theSinOne;
0126
0127   /**
0128    * The value of \f$\cos\theta_1\f$
0129    */
0130   double theCosOne;
0131
0132   /**
0133    * The value of \f$\sin(\theta_W - \theta_1)\f$
0134    */
0135   double theSinWmO;
0136
0137   /**
0138    * The value of \f$\cos(\theta_W - \theta_1)\f$
0139    */
0140   double theCosWmO;
0141
0142   /**
0143    * Store the CKM matrix for the \f$W\f$ interactions
0144    */
0145   vector< vector<Complex> > theCKM;
0146
0147   /**
0148    * The scale at which the coupling was last evaluated.
0149    */
0150   Energy2 theq2last;
0151
0152   /**
0153    * The value of \f$\sqrt{4\pi\alpha}\f$ when it was last evaluated.
0154    */
0155   Complex theCouplast;
0156
0157   /**
0158    * The value of the left coupling when it was last evaluated.
0159    */
0160   Complex theLlast;
0161
0162   /**
0163    * The value of the right coupling when it was last evaluated.
0164    */
0165   Complex theRlast;
0166
0167   /**
0168    * The last gauge boson in the vertex
0169    */
0170   long theGBlast;
0171
0172   /**
0173    * The last KK-fermion in the vertex
0174    */
0175   long theKKlast;
0176
0177   /**
0178    * The last SM fermion in the vertex
0179    */
0180   long theSMlast;
0181 };
0182 }
0183
0184
0185 #endif /* HERWIG_UEDF1F0W1Vertex_H */