Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2026-08-06 09:24:19

0001 // -*- C++ -*-
0002 //
0003 // SSGNGVertex.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_SSGNGVertex_H
0010 #define HERWIG_SSGNGVertex_H
0011 //
0012 // This is the declaration of the SSGNGVertex class.
0013 //
0014 
0015 #include "ThePEG/Helicity/Vertex/Vector/GeneralFFVVertex.h"
0016 #include "Herwig/Models/Susy/MSSM.h"
0017 #include "Herwig/Models/Susy/MixingMatrix.fh"
0018 
0019 namespace Herwig {
0020 using namespace ThePEG;
0021 
0022 /**
0023  * This class implements the coupling of a gluon to a gluino and a neutralino
0024  * via loop diagrams
0025  * It inherits from GeneralFFVertex and implements the setCoupling method.
0026  *
0027  * @see \ref SSGNGVertexInterfaces "The interfaces"
0028  * defined for SSGNGVertex.
0029  * @see FFVertex
0030  */
0031 class SSGNGVertex: public GeneralFFVVertex {
0032 
0033 public:
0034 
0035   /**
0036    * The default constructor.
0037    */
0038   SSGNGVertex();
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   /**
0067    * Calculate the couplings.
0068    * @param q2 The scale \f$q^2\f$ for the coupling at the vertex.
0069    * @param part1 The ParticleData pointer for the first  particle.
0070    * @param part2 The ParticleData pointer for the second particle.
0071    * @param part3 The ParticleData pointer for the third  particle.
0072    */
0073   virtual void setCoupling(Energy2 q2, tcPDPtr part1,
0074                            tcPDPtr part2, tcPDPtr part3);
0075 
0076 protected:
0077 
0078   /** @name Clone Methods. */
0079   //@{
0080   /**
0081    * Make a simple clone of this object.
0082    * @return a pointer to the new object.
0083    */
0084   virtual IBPtr clone() const {return new_ptr(*this);}
0085 
0086   /** Make a clone of this object, possibly modifying the cloned object
0087    * to make it sane.
0088    * @return a pointer to the new object.
0089    */
0090   virtual IBPtr fullclone() const {return new_ptr(*this);}
0091   //@}
0092 
0093 protected:
0094 
0095   /**
0096    *  Evaluate the loop integrals
0097    */
0098   void loopIntegrals(Energy Mi, Energy Mj, Energy M, Energy m,
0099              complex<InvEnergy2> & I, complex<InvEnergy2> & J,
0100              complex<InvEnergy2> & K, complex<InvEnergy2> & I2);
0101 
0102 protected:
0103 
0104   /** @name Standard Interfaced functions. */
0105   //@{
0106   /**
0107    * Initialize this object after the setup phase before saving an
0108    * EventGenerator to disk.
0109    * @throws InitException if object could not be initialized properly.
0110    */
0111   virtual void doinit();
0112 
0113   /**
0114    * Initialize this object. Called in the run phase just before
0115    * a run begins.
0116    */
0117   virtual void doinitrun();
0118 
0119   /**
0120    * Finalize this object. Called in the run phase just after a
0121    * run has ended. Used eg. to write out statistics.
0122    */
0123   virtual void dofinish();
0124   //@}
0125 
0126 private:
0127 
0128   /**
0129    * The assignment operator is private and must never be called.
0130    * In fact, it should not even be implemented.
0131    */
0132   SSGNGVertex & operator=(const SSGNGVertex &) = delete;
0133 
0134 private:
0135 
0136   /**
0137    *  Whether of not to include on-shell intermediate states
0138    */
0139   bool _includeOnShell;
0140 
0141   /**
0142    *  Only include the real part of the integral
0143    */
0144   bool _realIntegral;
0145 
0146   /**
0147    *  Option to omit light quark yukawas 
0148    */
0149   bool _omitLightQuarkYukawas;
0150 
0151   /**
0152    * Pointer to the stop mixing matrix
0153    */
0154   tMixingMatrixPtr _stop;
0155 
0156   /**
0157    * Pointer to the sbottom mixing matrix
0158    */
0159   tMixingMatrixPtr _sbot;
0160 
0161   /**
0162    * The value of \f$sin(\theta_W)\f$
0163    */
0164   double _sw;
0165   
0166   /**
0167    * The value of \f$cos(\theta_W)\f$
0168    */
0169   double _cw;
0170   
0171   /**
0172    * Mass of the W
0173    */
0174   Energy _mw;
0175 
0176   /**
0177    * \f$\sin(\beta)\f$
0178    */
0179   double _sb;
0180   
0181   /**
0182    * \f$\cos(\beta)\f$
0183    */
0184   double _cb;
0185 
0186   /**
0187    * Store the neutralino mixing matrix
0188    */
0189   tMixingMatrixPtr _theN;
0190 
0191   /**
0192    * Store the id of the neutralino when the coupling was last evaluated
0193    */
0194   long _idlast;
0195 
0196   /**
0197    * Store the value at which the coupling when it was last evaluated
0198    */
0199   Energy2 _q2last;
0200 
0201   /**
0202    * Store the value of the coupling when it was last evaluated
0203    */
0204   Complex _couplast;
0205 
0206   /**
0207    * Store the value of the left-coupling when it was last evaluated
0208    */
0209   complex<InvEnergy> _leftlast;
0210 
0211   /**
0212    * Store the value of the right-coupling when it was last evaluated
0213    */
0214   complex<InvEnergy> _rightlast;
0215 
0216   /**
0217    *  Whether or not initialised
0218    */
0219   bool _initLoops;
0220 };
0221 }
0222 
0223 #endif /* HERWIG_SSGNGVertex_H */