DalitzResonance.h

// -*- C++ -*-
#ifndef Herwig_DalitzResonance_H
#define Herwig_DalitzResonance_H
//
// This is the declaration of the DalitzResonance class.
//
 
#include "ThePEG/Config/ThePEG.h"
#include "DalitzResonance.fh"
#include <cassert>
 
namespace Herwig {
 
using namespace ThePEG;
 
namespace ResonanceType {
 
enum Type {NonResonant=0,
           Spin0=1,Spin1=3,Spin2=5,
           Spin0E691=11,Spin1E691=13,Spin2E691=15,
           BABARf0=21, Spin0Gauss=31, Flattef0=41, Spin0Complex=51,
           Flattea0=61, FlatteKstar0=71,
           Sigma=81,
           Spin1GS=23,
           Spin0NonResonant=101, Spin1NonResonant=103, Spin2NonResonant=105,
           Spin0MIPWA=-1, PiPiI2=-11, KMatrix=-21, LASS=-31};
}
 
class DalitzResonance: public Base {
 
public:
 
  DalitzResonance() : id(0), type(ResonanceType::NonResonant),mass(ZERO),width(ZERO),
                      daughter1(0),daughter2(0),spectator(0),amp(0.),R(ZERO)
  {}
 
  DalitzResonance(long pid, ResonanceType::Type rtype, Energy m, Energy w,
                  unsigned int d1, unsigned int d2, unsigned int s,
                  double mag, double phi, InvEnergy rr)
    : id(pid), type(rtype), mass(m),width(w),
      daughter1(d1),daughter2(d2),spectator(s),
      amp(mag*exp(Complex(0,phi))), R(rr)
  {}
 
public:
 
  virtual Complex BreitWigner(const Energy & mAB, const Energy & mA, const Energy & mB) const;
 
  virtual void dataBaseOutput(ofstream & output);
 
  static DalitzResonancePtr readResonance(string arg, string & error);
  
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
private:
 
  DalitzResonance & operator=(const DalitzResonance &) = delete;
 
public:
 
  long id;
 
  ResonanceType::Type type;
 
  Energy mass;
 
  Energy width;
 
  unsigned int daughter1,daughter2;
 
  unsigned int spectator;
 
  Complex amp;
 
  InvEnergy R;
  
};
 
}
 
#endif /* Herwig_DalitzResonance_H */