|
|
|||
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 */
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated by the 2.3.7 LXR engine. The LXR team |
|