BlobDiagram.h

// -*- C++ -*-
//
// BlobDiagram.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_BlobDiagram_H
#define ThePEG_BlobDiagram_H
// This is the declaration of the BlobDiagram class.

#include "ThePEG/MatrixElement/DiagramBase.h"
#include "ThePEG/MatrixElement/ColourLines.h"
#include "ThePEG/Handlers/StandardXComb.fh"

namespace ThePEG {

class BlobDiagram: public DiagramBase {

public:

  typedef cPDVector::size_type size_type;

public:

  BlobDiagram() 
    : DiagramBase() {}

  BlobDiagram(int id, tcPDPtr first, tcPDPtr second) 
    : DiagramBase() {
    addParton(first);
    addParton(second);
    diagramInfo(2,id);
  }

  ~BlobDiagram();

public:

  BlobDiagram& operator,(PDPtr pd) { addParton(pd); return *this; }

  BlobDiagram& operator,(cPDPtr pd) { addParton(pd); return *this; }

  BlobDiagram& operator,(tPDPtr pd) { addParton(pd); return *this; }

  BlobDiagram& operator,(tcPDPtr pd) { addParton(pd); return *this; }

  virtual tPVector construct(SubProPtr sb, const StandardXComb&,
                             const ColourLines&) const;

  tcPDPair incoming() const {
    return tcPDPair(partons()[0],partons()[1]); 
  }

  tcPDVector outgoing() const {
    return tcPDVector(partons().begin()+2,partons().end());
  }

  tcPDVector external() const {
    return tcPDVector(partons().begin(),partons().end());
  }

  size_type nOutgoing() const { return partons().size() - 2; }

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

private:

  static ClassDescription<BlobDiagram> initBlobDiagram;

  BlobDiagram & operator=(const BlobDiagram &) = delete;

};

}

namespace ThePEG {

template <>
struct BaseClassTrait<BlobDiagram,1>: public ClassTraitsType {
  typedef DiagramBase NthBase;
};

template <>
struct ClassTraits<BlobDiagram>: public ClassTraitsBase<BlobDiagram> {
  static string className() { return "ThePEG::BlobDiagram"; }
};

}

#endif /* ThePEG_BlobDiagram_H */