Builtin/Amplitudes/MatchboxAmplitudelnuqqbargg.h

0001 // -*- C++ -*-
0002 //
0003 // MatchboxAmplitudelnuqqbargg.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_MatchboxAmplitudelnuqqbargg_H
0010 #define Herwig_MatchboxAmplitudelnuqqbargg_H
0011 //
0012 // This is the declaration of the MatchboxAmplitudelnuqqbargg class.
0013 //
0014
0015 #include "Herwig/MatrixElement/Matchbox/Base/MatchboxAmplitude.h"
0016 #include "Herwig/MatrixElement/Matchbox/Builtin/Amplitudes/MatchboxCurrents.h"
0017 #include "ThePEG/StandardModel/StandardModelBase.h"
0018 #include "ThePEG/Interface/ClassDocumentation.h"
0019 #include "ThePEG/EventRecord/Particle.h"
0020 #include "ThePEG/Repository/UseRandom.h"
0021 #include "ThePEG/Repository/EventGenerator.h"
0022 #include "ThePEG/Utilities/DescribeClass.h"
0023 #include "ThePEG/Persistency/PersistentOStream.h"
0024 #include "ThePEG/Persistency/PersistentIStream.h"
0025 #include "ThePEG/PDT/EnumParticles.h"
0026 #include "Herwig/MatrixElement/Matchbox/Utility/SU2Helper.h"
0027
0028 namespace Herwig {
0029
0030 using namespace ThePEG;
0031
0032 /**
0033  * \ingroup Matchbox
0034  * \author Thomas Schuh
0035  *
0036  * \brief MatchboxAmplitudelnuqqbargg
0037  */
0038 class MatchboxAmplitudelnuqqbargg: public MatchboxAmplitude, public MatchboxCurrents {
0039
0040 public:
0041
0042   /**
0043    * The default constructor.
0044    */
0045   MatchboxAmplitudelnuqqbargg();
0046
0047   /**
0048    * Return the (tree-level) order in \f$g_S\f$ in which this matrix
0049    * element is given.
0050    */
0051   virtual unsigned int orderInGs() const { return 2; }
0052
0053   /**
0054    * Return the (tree-level) order in \f$g_{EM}\f$ in which this matrix
0055    * element is given.
0056    */
0057   virtual unsigned int orderInGem() const { return 2; }
0058
0059   /**
0060    * Return true, if this amplitude is capable of calculating one-loop
0061    * (QCD) corrections.
0062    */
0063   virtual bool haveOneLoop() const { return false; }
0064
0065   /**
0066    * Return true, if this amplitude can handle the given process.
0067    */
0068   virtual bool canHandle(const PDVector&) const;
0069
0070   /**
0071    * Calculate the tree level amplitudes for the phasespace point
0072    * stored in lastXComb.
0073    */
0074   virtual void prepareAmplitudes(Ptr<MatchboxMEBase>::tcptr);
0075
0076   /**
0077    * Evaluate the amplitude for the given colour tensor id and
0078    * helicity assignment
0079    */
0080   virtual Complex evaluate(size_t, const vector<int>&, Complex&);
0081
0082   /**
0083    * Return true, if one loop corrections are given in the conventions
0084    * of BDK.
0085    */
0086   virtual bool isBDK() const { return true; }
0087
0088   /**
0089    * Return the value of the dimensional regularization
0090    * parameter. Note that renormalization scale dependence is fully
0091    * restored in DipoleIOperator.
0092    */
0093   virtual Energy2 mu2() const { return lastSHat(); }
0094
0095   /**
0096    * Flush all cashes.
0097    */
0098   virtual void flushCaches() {
0099     MatchboxCurrents::reset();
0100     MatchboxAmplitude::flushCaches();
0101   }
0102
0103   /**
0104    * Function used to write out object persistently.
0105    * @param os the persistent output stream written to.
0106    */
0107   void persistentOutput(PersistentOStream & os) const;
0108
0109   /**
0110    * Function used to read in object persistently.
0111    * @param is the persistent input stream read from.
0112    * @param version the version number of the object when written.
0113    */
0114   void persistentInput(PersistentIStream & is, int version);
0115
0116   /**
0117    * The standard Init function used to initialize the interfaces.
0118    * Called exactly once for each class by the class description system
0119    * before the main function starts or
0120    * when this class is dynamically loaded.
0121    */
0122   static void Init();
0123
0124 protected:
0125
0126   /**
0127    * Make a simple clone of this object.
0128    * @return a pointer to the new object.
0129    */
0130   virtual IBPtr clone() const;
0131
0132   /** Make a clone of this object, possibly modifying the cloned object
0133    * to make it sane.
0134    * @return a pointer to the new object.
0135    */
0136   virtual IBPtr fullclone() const;
0137
0138   /**
0139    * Initialize this object after the setup phase before saving an
0140    * EventGenerator to disk.
0141    * @throws InitException if object could not be initialized properly.
0142    */
0143   virtual void doinit();
0144
0145   /**
0146    * Initialize this object. Called in the run phase just before
0147    * a run begins.
0148    */
0149   virtual void doinitrun();
0150
0151 private:
0152
0153   /**
0154    * True, if a diagonal CKM matrix should be assumed. This ignores
0155    * the CKM object of the StandardModel.
0156    */
0157   bool theDiagonal;
0158
0159   /**
0160    * The ckm.
0161    */
0162   vector< vector<Complex> > theCKM;
0163
0164   /**
0165    * The assignment operator is private and must never be called.
0166    * In fact, it should not even be implemented.
0167    */
0168   MatchboxAmplitudelnuqqbargg & operator=(const MatchboxAmplitudelnuqqbargg &) = delete;
0169
0170 };
0171
0172 }
0173
0174 #endif /* Herwig_MatchboxAmplitudelnuqqbargg_H */