ResonantProcessConstructor.h

// -*- C++ -*-
//
// ResonantProcessConstructor.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_ResonantProcessConstructor_H
#define HERWIG_ResonantProcessConstructor_H
//
// This is the declaration of the ResonantProcessConstructor class.
//
 
#include "HardProcessConstructor.h"
#include "ThePEG/Utilities/Exception.h"
#include "ResonantProcessConstructor.fh"
 
namespace Herwig {
using namespace ThePEG;
 
class ResonantProcessConstructor: public HardProcessConstructor {
 
public:
 
  typedef set<tPDPtr> tPDSet;
 
  typedef vector<vector<double> > CFMatrix;
 
  enum direction {incoming, outgoing};
 
public:
 
  ResonantProcessConstructor() :
    processOption_(0), scaleChoice_(1), scaleFactor_(1.),
    incoming_(0), intermediates_(0),
    outgoing_(0), diagrams_(0)
  {}
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
  void constructDiagrams() ;
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
private:
 
  void constructVertex2(IDPair in, VertexBasePtr vertex, 
                        PDPtr partc);
  
  void makeResonantDiagram(IDPair in, PDPtr offshell, long outa, 
                           long outb, VBPair vertices);
  
  tPDSet search(VertexBasePtr vertex, long part1, direction d1, 
               long part2, direction d2, direction d3 = outgoing);
 
  IDPair find(long part, const PDVector & out) const;
 
  void createMatrixElement(const HPDiagram & diag) const;
 
  string MEClassname(const tcPDVector & extpart, tcPDPtr inter,
                     string & objname) const;
 
private:
 
  ResonantProcessConstructor & operator=(const ResonantProcessConstructor &) = delete;
 
private:
 
  unsigned int processOption_;
 
  unsigned int scaleChoice_;
 
  double scaleFactor_;
  
  vector<PDPtr> incoming_;
 
  vector<PDPtr> intermediates_;
 
  vector<PDPtr> outgoing_; 
 
  vector<HPDiagram> diagrams_;
 
};
 
  class RPConstructorError : public Exception {};
 
}
 
 
#endif /* HERWIG_ResonantProcessConstructor_H */