Matchbox/Dipoles/IFqqxDipole.h

0001 // -*- C++ -*-
0002 //
0003 // IFqqxDipole.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_IFqqxDipole_H
0010 #define HERWIG_IFqqxDipole_H
0011 //
0012 // This is the declaration of the IFqqxDipole class.
0013 //
0014
0015 #include "Herwig/MatrixElement/Matchbox/Dipoles/SubtractionDipole.h"
0016
0017 namespace Herwig {
0018
0019 using namespace ThePEG;
0020
0021 /**
0022  * \ingroup Matchbox
0023  * \author Simon Platzer
0024  *
0025  * \brief IFqqxDipole implements the D^{q,qbar}_k subtraction dipole.
0026  *
0027  */
0028 class IFqqxDipole: public SubtractionDipole {
0029
0030 public:
0031
0032   /**
0033    * The default constructor.
0034    */
0035   IFqqxDipole();
0036
0037 public:
0038
0039   /**
0040    * Return true, if this dipole can possibly handle the indicated
0041    * emitter.
0042    */
0043   virtual bool canHandleEmitter(const cPDVector& partons, int emitter) const {
0044     return emitter < 2 && abs(partons[emitter]->id()) < 7;
0045   }
0046
0047   /**
0048    * Return true, if this dipole can possibly handle the indicated
0049    * splitting.
0050    */
0051   virtual bool canHandleSplitting(const cPDVector& partons, int emitter, int emission) const {
0052     return
0053       canHandleEmitter(partons,emitter) &&
0054       partons[emitter]->id() == partons[emission]->id();
0055   }
0056
0057   /**
0058    * Return true, if this dipole can possibly handle the indicated
0059    * spectator.
0060    */
0061   virtual bool canHandleSpectator(const cPDVector& partons, int spectator) const {
0062     return spectator > 1 && partons[spectator]->coloured();
0063   }
0064
0065   /**
0066    * Return true, if this dipole applies to the selected
0067    * configuration.
0068    */
0069   virtual bool canHandle(const cPDVector& partons,
0070              int emitter, int emission, int spectator) const;
0071
0072   /**
0073    *  How to sample the z-distribution.
0074    *  FlatZ = 1
0075    *  OneOverZ = 2
0076    *  OneOverOneMinusZ = 3
0077    *  OneOverZOneMinusZ = 4
0078    */
0079
0080   virtual int samplingZ() const {return 2;}
0081
0082   /**
0083    * Return the matrix element for the kinematical configuation
0084    * previously provided by the last call to setKinematics(), suitably
0085    * scaled by sHat() to give a dimension-less number.
0086    */
0087   virtual double me2() const;
0088
0089   /**
0090    * Return the matrix element averaged over spin correlations.
0091    */
0092   virtual double me2Avg(double ccme2) const;
0093
0094 public:
0095
0096   /** @name Functions used by the persistent I/O system. */
0097   //@{
0098   /**
0099    * Function used to write out object persistently.
0100    * @param os the persistent output stream written to.
0101    */
0102   void persistentOutput(PersistentOStream & os) const;
0103
0104   /**
0105    * Function used to read in object persistently.
0106    * @param is the persistent input stream read from.
0107    * @param version the version number of the object when written.
0108    */
0109   void persistentInput(PersistentIStream & is, int version);
0110   //@}
0111
0112   /**
0113    * The standard Init function used to initialize the interfaces.
0114    * Called exactly once for each class by the class description system
0115    * before the main function starts or
0116    * when this class is dynamically loaded.
0117    */
0118   static void Init();
0119
0120 protected:
0121
0122   /** @name Clone Methods. */
0123   //@{
0124   /**
0125    * Make a simple clone of this object.
0126    * @return a pointer to the new object.
0127    */
0128   virtual IBPtr clone() const;
0129
0130   /** Make a clone of this object, possibly modifying the cloned object
0131    * to make it sane.
0132    * @return a pointer to the new object.
0133    */
0134   virtual IBPtr fullclone() const;
0135   //@}
0136
0137
0138 // If needed, insert declarations of virtual function defined in the
0139 // InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
0140
0141
0142 private:
0143
0144   /**
0145    * The assignment operator is private and must never be called.
0146    * In fact, it should not even be implemented.
0147    */
0148   IFqqxDipole & operator=(const IFqqxDipole &) = delete;
0149
0150 };
0151
0152 }
0153
0154 #endif /* HERWIG_IFqqxDipole_H */