Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2026-08-06 09:38:27

0001 // -*- C++ -*-
0002 //
0003 // MENCDIS.h is a part of ThePEG - Toolkit for HEP Event Generation
0004 // Copyright (C) 1999-2019 Leif Lonnblad
0005 //
0006 // ThePEG 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 ThePEG_MENCDIS_H
0010 #define ThePEG_MENCDIS_H
0011 // This is the declaration of the MENCDIS class.
0012 
0013 #include "ThePEG/MatrixElement/ME2to2QCD.h"
0014 // #include "MENCDIS.fh"
0015 // #include "MENCDIS.xh"
0016 
0017 namespace ThePEG {
0018 
0019 /**
0020  * The MENCDIS class implements the \f$e^\pm q\rightarrow e^\pm q\f$
0021  * matrix element. Both the gamma and \f$Z^0\f$ terms as well
0022  * as the interference term is included. Although not a strict QCD
0023  * matrix element the class inherits from ME2to2QCD, mainly to inherit
0024  * the parameter for the number of active quark flavours.
0025  *
0026  * @see \ref MENCDISInterfaces "The interfaces"
0027  * defined for MENCDIS.
0028  * @see ME2to2QCD
0029  */
0030 class MENCDIS: public ME2to2QCD {
0031 
0032 public:
0033 
0034   /** @name Standard constructors and destructors. */
0035   //@{
0036   /**
0037    * Default constructor.
0038    */
0039   MENCDIS();
0040 
0041   /**
0042    * Copy-constructor.
0043    */
0044   MENCDIS(const MENCDIS &);
0045   //@}
0046 
0047 public:
0048 
0049   /** @name Virtual functions required by the MEBase class. */
0050   //@{
0051   /**
0052    * Return the order in \f$\alpha_S\f$ in which this matrix element
0053    * is given. Returns .
0054    */
0055   virtual unsigned int orderInAlphaS() const;
0056 
0057   /**
0058    * Return the order in \f$\alpha_{EM}\f$ in which this matrix
0059    * element is given. Returns 2.
0060    */
0061   virtual unsigned int orderInAlphaEW() const;
0062 
0063   /**
0064    * The matrix element for the kinematical configuration
0065    * previously provided by the last call to setKinematics(), suitably
0066    * scaled by sHat() to give a dimension-less number.
0067    * @return the matrix element scaled with sHat() to give a
0068    * dimensionless number.
0069    */
0070   virtual double me2() const;
0071 
0072   /**
0073    * Add all possible diagrams with the add() function.
0074    */
0075   virtual void getDiagrams() const;
0076 
0077   /**
0078    * Return a Selector with possible colour geometries for the selected
0079    * diagram weighted by their relative probabilities.
0080    * @param diag the diagram chosen.
0081    * @return the possible colour geometries weighted by their
0082    * relative probabilities.
0083    */
0084   virtual Selector<const ColourLines *>
0085   colourGeometries(tcDiagPtr diag) const;
0086 
0087   /**
0088    * Get diagram selector. With the information previously supplied with the
0089    * setKinematics method, a derived class may optionally
0090    * override this method to weight the given diagrams with their
0091    * (although certainly not physical) relative probabilities.
0092    * @param dv the diagrams to be weighted.
0093    * @return a Selector relating the given diagrams to their weights.
0094    */
0095   virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
0096 
0097   /**
0098    * Return the scale associated with the last set phase space point.
0099    */
0100   virtual Energy2 scale() const;
0101   //@}
0102 
0103 public:
0104 
0105   /** @name Functions used by the persistent I/O system. */
0106   //@{
0107   /**
0108    * Function used to write out object persistently.
0109    * @param os the persistent output stream written to.
0110    */
0111   void persistentOutput(PersistentOStream & os) const;
0112 
0113   /**
0114    * Function used to read in object persistently.
0115    * @param is the persistent input stream read from.
0116    * @param version the version number of the object when written.
0117    */
0118   void persistentInput(PersistentIStream & is, int version);
0119   //@}
0120 
0121   /**
0122    * Standard Init function used to initialize the interface.
0123    */
0124   static void Init();
0125 
0126 protected:
0127 
0128   /** @name Clone Methods. */
0129   //@{
0130   /**
0131    * Make a simple clone of this object.
0132    * @return a pointer to the new object.
0133    */
0134   virtual IBPtr clone() const;
0135 
0136   /** Make a clone of this object, possibly modifying the cloned object
0137    * to make it sane.
0138    * @return a pointer to the new object.
0139    */
0140   virtual IBPtr fullclone() const;
0141   //@}
0142 
0143   /** @name Standard Interfaced functions. */
0144   //@{
0145   /**
0146    * Initialize this object after the setup phase before saving an
0147    * EventGenerator to disk.
0148    * @throws InitException if object could not be initialized properly.
0149    */
0150   virtual void doinit();
0151   //@}
0152 
0153 protected:
0154 
0155   /**
0156    * The squared mass of the Z0.
0157    */
0158   Energy2 mZ2;
0159 
0160   /**
0161    * The last gamma term to be used to select primary diagram.
0162    */
0163   mutable double lastG;
0164 
0165   /**
0166    * The last Z0 term to be used to select primary diagram.
0167    */
0168   mutable double lastZ;
0169 
0170 private:
0171 
0172   /**
0173    * Describe a concrete class with persistent data.
0174    */
0175   static ClassDescription<MENCDIS> initMENCDIS;
0176 
0177   /**
0178    *  Private and non-existent assignment operator.
0179    */
0180   MENCDIS & operator=(const MENCDIS &) = delete;
0181 
0182 };
0183 
0184 /** @cond TRAITSPECIALIZATIONS */
0185 
0186 /** This template specialization informs ThePEG about the
0187  *  base classes of MENCDIS. */
0188 template <>
0189 struct BaseClassTrait<MENCDIS,1>: public ClassTraitsType {
0190   /** Typedef of the first base class of MENCDIS. */
0191   typedef ME2to2QCD NthBase;
0192 };
0193 
0194 /** This template specialization informs ThePEG about the name of the
0195  *  MENCDIS class and the shared object where it is defined. */
0196 template <>
0197 struct ClassTraits<MENCDIS>: public ClassTraitsBase<MENCDIS> {
0198   /** Return a platform-independent class name */
0199   static string className() { return "ThePEG::MENCDIS"; }
0200   /** Return the name of the shared library be loaded to get
0201    *  access to the MENCDIS class and every other class it uses
0202    *  (except the base class). */
0203   static string library() { return "MENCDIS.so"; }
0204 };
0205 
0206 /** @endcond */
0207 
0208 }
0209 
0210 #endif /* ThePEG_MENCDIS_H */