FIMsqgxDipole.h

// -*- C++ -*-
//
// FIMsqgxDipole.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 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.
//
#ifndef HERWIG_FIMsqgxDipole_H
#define HERWIG_FIMsqgxDipole_H
//
// This is the declaration of the FIMsqgxDipole class.
//
 
#include "Herwig/MatrixElement/Matchbox/Dipoles/SubtractionDipole.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class FIMsqgxDipole: public SubtractionDipole {
 
public:
 
  FIMsqgxDipole();
 
public:
 
  virtual bool canHandleEmitter(const cPDVector& partons, int emitter) const {
    return 
      emitter > 1 &&
      (abs(abs(partons[emitter]->id()) - 1000000) < 7 || abs(abs(partons[emitter]->id()) - 2000000) < 7);
  }
 
  virtual bool canHandleSplitting(const cPDVector& partons, int emitter, int emission) const {
    return 
      canHandleEmitter(partons,emitter) &&
      partons[emission]->id() == ParticleID::g;
  }
 
  virtual bool canHandleSpectator(const cPDVector& partons, int spectator) const {
    return spectator < 2 && partons[spectator]->coloured();
  }
 
  virtual bool canHandle(const cPDVector& partons,
                         int emitter, int emission, int spectator) const;
 
  virtual int samplingZ() const {return 3;}
  
  virtual double me2() const;
 
  virtual double me2Avg(double ccme2) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
 
private:
 
  FIMsqgxDipole & operator=(const FIMsqgxDipole &) = delete;
 
};
 
}
 
#endif /* HERWIG_FIMsqgxDipole_H */