TwoBodyAllOnCalculator.h

// -*- C++ -*-
//
// TwoBodyAllOnCalculator.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_TwoBodyAllOnCalculator_H
#define HERWIG_TwoBodyAllOnCalculator_H
// This is the declaration of the TwoBodyAllOnCalculator class.
 
#include "WidthCalculatorBase.h"
#include "TwoBodyAllOnCalculator.fh"
#include "GenericWidthGenerator.h"
 
namespace Herwig {
using namespace ThePEG;
 
class TwoBodyAllOnCalculator: public WidthCalculatorBase {
 
public:
 
  friend class GenericWidthGenerator;
 
public:
 
  TwoBodyAllOnCalculator(tGenericWidthGeneratorPtr inwidth,int imode,
                         Energy m1,Energy m2)
    : _mode(imode),_mass1(m1),_mass2(m2),_widthgen(inwidth)
  {}
 
  Energy partialWidth(Energy2 scale) const;
 
  void resetMass(int imass,Energy mass) {
    if(imass==1)      _mass1=mass;
    else if(imass==2) _mass2=mass;
    else throw Exception() << "Unknown particle in " 
                           << "TwoBodyAllOnCalculator::resetMass()"
                           << Exception::runerror;
  }
 
  Energy getMass(const int imass) const {
    if(imass==1)      return _mass1;
    else if(imass==2) return _mass2;
    else throw Exception() << "Unknown particle in " 
                           << "TwoBodyAllOnCalculator::getMass()"
                           << Exception::runerror;
  }
 
  Energy otherMass(const int imass) const {
    if(imass==1)      return _mass2;
    else if(imass==2) return _mass1;
    else throw Exception() << "Unknown particle in " 
                           << "TwoBodyAllOnCalculator::otherMass()"
                           << Exception::runerror;
  }
 
private:
 
  TwoBodyAllOnCalculator & operator=(const TwoBodyAllOnCalculator &) = delete;
 
private:
 
  int _mode;
 
  Energy _mass1;
 
  Energy _mass2;
 
  GenericWidthGeneratorPtr _widthgen;
 
};
 
}
 
#endif /* HERWIG_TwoBodyAllOnCalculator_H */