diff --git a/Models/StandardModel/SMFFGVertex.cc b/Models/StandardModel/SMFFGVertex.cc --- a/Models/StandardModel/SMFFGVertex.cc +++ b/Models/StandardModel/SMFFGVertex.cc @@ -1,72 +1,71 @@ // -*- C++ -*- // // SMFFGVertex.cc is a part of Herwig - A multi-purpose Monte Carlo event generator // Copyright (C) 2002-2017 The Herwig Collaboration // // Herwig is licenced under version 3 of the GPL, see COPYING for details. // Please respect the MCnet academic guidelines, see GUIDELINES for details. // // // This is the implementation of the non-inlined, non-templated member // functions of the StandardModelFFGVertex class. // #include "SMFFGVertex.h" #include "ThePEG/Utilities/DescribeClass.h" #include "ThePEG/Interface/ClassDocumentation.h" #include "ThePEG/Persistency/PersistentOStream.h" #include "ThePEG/Persistency/PersistentIStream.h" #include "ThePEG/StandardModel/StandardModelBase.h" using namespace Herwig; using namespace ThePEG; // The following static variable is needed for the type // description system in ThePEG. DescribeNoPIOClass describeHerwigSMFFGVertex("Herwig::SMFFGVertex", "Herwig.so"); void SMFFGVertex::Init() { static ClassDocumentation documentation ("The SMFFGVertex class is the implementation of" "the coupling of the gluon to the Standard Model fermions"); } // coupling for FFG vertex #ifndef NDEBUG void SMFFGVertex::setCoupling(Energy2 q2,tcPDPtr aa,tcPDPtr bb,tcPDPtr) { #else void SMFFGVertex::setCoupling(Energy2 q2,tcPDPtr aa,tcPDPtr,tcPDPtr) { #endif // first the overall normalisation if(q2!=_q2last||_couplast==0.) { _couplast = -strongCoupling(q2); - cerr << "testing coupling " << sqr(_couplast)/4./Constants::pi << "\n"; _q2last=q2; } // the left and right couplings assert( abs(aa->id()) >= 1 && abs(aa->id()) <= 6 ); assert( aa->id()==-bb->id()); if(aa->id()<0) norm( _couplast); else norm(-_couplast); left(1.); right(1.); } SMFFGVertex::SMFFGVertex() : _couplast(0.), _q2last(ZERO) { orderInGs(1); orderInGem(0); colourStructure(ColourStructure::SU3TFUND); } void SMFFGVertex::doinit() { // PDG codes for the particles for(int ix=1;ix<7;++ix) { addToList(-ix,ix,21); } FFVVertex::doinit(); }