#include <complex>
#include <vector>

#ifndef _RESONANCE_
#define _RESONANCE_

class Resonance{
 public:
  ///Constructor & Dectructor
  Resonance();
  Resonance(double mass, double width, double J, int P);
  virtual ~Resonance(){};

  ///Get parameters of resonance
  ///mass
  double M();
  ///width
  double W();
  ///momenta
  double J();
  ///momenta projection
  std::vector<double>  Mj();// obtain whole vector
  double Mj(int k);//obtain k-th element from vector;
  ///P parity
  int P();
  /* ///Spin */
  /* double S(); */
  std::vector<double>  l();// obtain whole vector 
  void setl( std::vector<double> l){ _l = l;};// set vector

  void Print();

 private:
  double _J;//total moment
  std::vector<double> _Mj;// moment projection
  double _mass; // mass of resonance
  double _width;// width of resonance
  std::vector<double> _l;   //possible angular momenta for this reaction
  int _p;
};
#endif