SimpleColourBasis2.h

// -*- C++ -*-
//
// SimpleColourBasis2.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_SimpleColourBasis2_H
#define Herwig_SimpleColourBasis2_H
//
// This is the declaration of the SimpleColourBasis2 class.
//
 
#include "Herwig/MatrixElement/Matchbox/Utility/ColourBasis.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class SimpleColourBasis2: public ColourBasis {
 
public:
 
  virtual size_t prepareBasis(const vector<PDT::Colour>&);
 
  virtual double scalarProduct(size_t a, size_t b,
                               const vector<PDT::Colour>& abBasis) const;
 
  virtual double tMatrixElement(size_t i, size_t a, size_t b,
                                const vector<PDT::Colour>& aBasis,
                                const vector<PDT::Colour>& bBasis,
                                size_t k, size_t l,
                                const map<size_t,size_t>& dict) const;
 
  /*
   * Temporary to make it compile
   */
  virtual double sMatrixElement(size_t, size_t, size_t,
                                const vector<PDT::Colour>&,
                                const vector<PDT::Colour>&,
                                size_t, size_t,
                                const map<size_t,size_t>&
                                ) const {
    assert( 0 == 1 );
    return 0;
  }
 
  virtual bool canSplitGluons() const {
    return false;
  }
 
  virtual bool colourConnected(const cPDVector&,
                               const vector<PDT::Colour>&,
                               const pair<int,bool>&, 
                               const pair<int,bool>&, 
                               size_t) const;
 
  virtual bool haveColourFlows() const { return true; }
 
  void makeIds() const;
 
  virtual map<size_t,vector<vector<size_t> > > basisList(const vector<PDT::Colour>&) 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:
 
  mutable vector<PDT::Colour> id88;
 
  mutable vector<PDT::Colour> id33bar;
 
  mutable vector<PDT::Colour> id888;
 
  mutable vector<PDT::Colour> id33bar8;
 
  mutable vector<PDT::Colour> id8888;
 
  mutable vector<PDT::Colour> id33bar88;
 
  mutable vector<PDT::Colour> id33bar33bar;
 
  mutable vector<PDT::Colour> id88888;
 
  mutable vector<PDT::Colour> id33bar888;
 
  mutable vector<PDT::Colour> id33bar33bar8;
 
private:
 
  SimpleColourBasis2 & operator=(const SimpleColourBasis2 &) = delete;
 
};
 
}
 
#endif /* Herwig_SimpleColourBasis2_H */