HardProcessConstructor.h

// -*- C++ -*-
#ifndef HERWIG_HardProcessConstructor_H
#define HERWIG_HardProcessConstructor_H
//
// This is the declaration of the HardProcessConstructor class.
//
 
#include "ThePEG/Interface/Interfaced.h"
#include "HPDiagram.h"
#include "Herwig/Models/StandardModel/StandardModel.h"
#include "ThePEG/Handlers/SubProcessHandler.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/Handlers/StandardEventHandler.h"
#include "Herwig/MatrixElement/General/GeneralHardME.h"
#include "HardProcessConstructor.fh"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HardProcessConstructor: public Interfaced {
 
public:
 
  typedef vector<HPDiagram> HPDVector;
 
public:
 
  HardProcessConstructor() : debug_(false) {}
 
  virtual void constructDiagrams() = 0;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  void fixFSOrder(HPDiagram & diag);
 
  void assignToCF(HPDiagram & diag);
 
  void sChannelCF(HPDiagram & diag);
 
  void tChannelCF(HPDiagram & diag);
 
  void uChannelCF(HPDiagram & diag);
 
  void fourPointCF(HPDiagram & diag);
 
  tHwSMPtr model() const {return model_;}
  
  tSubHdlPtr subProcess() const {return subProcess_;}
 
  bool debug() const {return debug_;}
 
  GeneralHardME::ColourStructure colourFlow(const tcPDVector & extpart) const;
 
  bool duplicate(const HPDiagram & diagram, 
                 const HPDVector & group) const;
 
  bool checkOrder(const HPDiagram & diag) const;
 
protected:
 
  virtual void doinit();
 
private:
 
  HardProcessConstructor & operator=(const HardProcessConstructor &) = delete;
 
private:
 
  tHwSMPtr model_;
  
   tSubHdlPtr subProcess_;
 
  bool debug_;
 
};
 
namespace HPC_helper {
  // Helper functor for find_if in duplicate function.
  class SameIncomingAs {
  public:
    SameIncomingAs(tPDPair in) : a(in.first->id()), b(in.second->id())  {}
    bool operator()(tPDPair ppair) const {
      long id1(ppair.first->id()), id2(ppair.second->id());
      return ( id1 == a && id2 == b ) || ( id1 == b && id2 == a );
    }
  private:
    long a, b;
  };
 
  inline bool duplicateIncoming(tPDPair ppair,const vector<tPDPair> &incPairs) {
    vector<tPDPair>::const_iterator it = 
      find_if( incPairs.begin(), incPairs.end(), SameIncomingAs(ppair) );
    return it != incPairs.end(); 
  }
}
 
}
 
#endif /* HERWIG_HardProcessConstructor_H */