#include "Resonance.h"
//ClassImp(Intensity);
#include <iostream> 
using namespace std;

Resonance::Resonance(){
  //  std::cout<<"new object Intensity is created!"<<std::endl;
}
Resonance::Resonance(double mass, double width, double J, int P){
  _J = J;
  _mass = mass;
  _width = width;
  for(int i=0;i<(2*J+1);i++) _Mj.push_back(i-J);
  // const int sizel = l.size();
  // for(int il=0;il<sizel;il++)
  //   _l.push_back(l[il]);
  _p = P;
}

double Resonance::M(){return _mass;}

double Resonance::W(){return _width;}

double Resonance::J(){return _J;}

int Resonance::P(){return _p;}

std::vector<double> Resonance::Mj(){
  // std::vector<double> Mjc;
  // for(int i=0;i<=(2*_J+1);i++) Mjc.push_back(_Mj[i]);
  // return Mjc;
  return _Mj;
}
std::vector<double> Resonance::l(){
  return _l;
}

double Resonance::Mj(int k){
  if(k<(2*_J+1)) return _Mj[k];
  else return -1000;
}

void Resonance::Print(){
  std::cout<<"Resonance with:"<<std::endl;
  std::cout<<" _mass = "<<_mass<<", _width = "<<_width<<std::endl;
  std::cout<<"J(P)="<<_J<<"("<<_p<<")"<<std::endl;
}