TwoBodyDecay.h

// -*- C++ -*-
//
// TwoBodyDecay.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_TwoBodyDecay_H
#define HERWIG_TwoBodyDecay_H
//
// This is the declaration of the TwoBodyDecay struct.
//
 
namespace Herwig {
using namespace ThePEG;
using Helicity::tVertexBasePtr;
 
struct TwoBodyDecay {
  
public:
  
  TwoBodyDecay(tPDPtr pa, tPDPtr pb, tPDPtr pc,
               tVertexBasePtr vertex) : parent_(pa), vertex_(vertex) {
    ParticleOrdering order;
    if( order(pb, pc) ) {
      children_.first = pb;
      children_.second = pc;
    }
    else {
      children_.first = pc;
      children_.second = pb;
    }
  }
  
  tPDPtr parent_;
  
  tPDPair children_;
  
  tVertexBasePtr vertex_;
  
private:
  
  TwoBodyDecay();
};
 
inline bool operator<(const TwoBodyDecay & x, const TwoBodyDecay & y) {
  if(x.parent_->id()!=y.parent_->id())
    return x.parent_->id()<y.parent_->id();
  if(x.children_.first->id()!=y.children_.first->id())
    return x.children_.first->id() < y.children_.first->id();
  if(x.children_.second->id()!=y.children_.second->id())
    return x.children_.second->id() < y.children_.second->id();
  return false;
}
inline bool operator==(const TwoBodyDecay & x, const TwoBodyDecay & y) {
  if(x.parent_->id()!=y.parent_->id())
    return false;
  if(x.children_.first->id()!=y.children_.first->id())
    return false;
  if(x.children_.second->id()!=y.children_.second->id())
    return false;
  return true;
}
 
}
#endif /* HERWIG_TwoBodyDecay_H */