MatchboxInsertionOperator.h

// -*- C++ -*-
//
// MatchboxInsertionOperator.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_MatchboxInsertionOperator_H
#define HERWIG_MatchboxInsertionOperator_H
//
// This is the declaration of the MatchboxInsertionOperator class.
//
 
#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/Handlers/StandardXComb.h"
#include "ThePEG/Handlers/LastXCombInfo.h"
#include "Herwig/MatrixElement/Matchbox/InsertionOperators/MatchboxInsertionOperator.fh"
#include "Herwig/MatrixElement/Matchbox/Utility/LastMatchboxXCombInfo.h"
#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.fh"
#include "Herwig/MatrixElement/Matchbox/MatchboxFactory.fh"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MatchboxInsertionOperator: 
    public HandlerBase, 
    public LastXCombInfo<StandardXComb>,
    public LastMatchboxXCombInfo {
 
public:
 
  MatchboxInsertionOperator();
 
  virtual ~MatchboxInsertionOperator();
 
public:
 
  Ptr<MatchboxFactory>::tptr factory() const;
 
 
  virtual bool apply(const cPDVector&) const = 0;
 
  Ptr<MatchboxMEBase>::tptr lastBorn() const ;
 
  virtual void setXComb(tStdXCombPtr xc) { 
    theLastXComb = xc;
    lastMatchboxXComb(xc);
  }
      
      
  virtual void setAlpha(double )const {assert(false);}
  
    
 
  virtual int nDimAdditional() const { return 0; }
 
 
 
 
      virtual bool isDRbar() const;
      virtual bool isDR() const;
      virtual bool isCS() const;
      virtual bool isBDK() const;
      virtual bool isExpanded() const;
  virtual double oneLoopDoublePole() const { return 0.; }
 
  virtual double oneLoopSinglePole() const { return 0.; }
 
 
 
  virtual double me2() const = 0;
 
  virtual CrossSection dSigHatDR() const;
      
      
  virtual CrossSection dSigHatDRAlphaDiff(double alpha) const{
    setAlpha(alpha);
    CrossSection res=dSigHatDR();
    setAlpha(1.);
    return res-dSigHatDR();
  }
 
 
 
  virtual void flushCaches() {}
 
  Ptr<MatchboxInsertionOperator>::ptr cloneMe() const {
    return dynamic_ptr_cast<Ptr<MatchboxInsertionOperator>::ptr>(clone());
  }
 
  virtual void cloneDependencies(const std::string& prefix = "");
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
private:
 
  MatchboxInsertionOperator & operator=(const MatchboxInsertionOperator &) = delete;
 
};
 
}
 
#endif /* HERWIG_MatchboxInsertionOperator_H */