#include "Particle.h"
//ClassImp(Intensity);
#include <iostream> 

using namespace std;

Particle::Particle(){
  //  std::cout<<"new object Intensity is created!"<<std::endl;
}
Particle::Particle(double mass, double s, int p){
  _s = s;
  _mass = mass;
  _p = p;
  for(int i=0;i<=int(2*s);i++){
    if(_mass==0 && (_s==1) && (-_s+i)==0) continue;
    _lya.push_back(-_s+i);
  }
}

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

double Particle::S(){return _s;}

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

std::vector<double> Particle::lya(){
  std::vector<double> lyac;
  const int sizelya = _lya.size();
  for(int i=0;i<sizelya;i++) lyac.push_back(_lya[i]);
  return lyac;
}

double Particle::lya(int k){
  const int sizelya = _lya.size();
  if(k<sizelya) return _lya[k];
  else return -1000;
}

void Particle::Print(){
  std::cout<<"Particle with:"<<std::endl;
  std::cout<<" _mass = "<<_mass<<std::endl;
  std::cout<<"J(P)="<<_s<<"("<<_p<<")"<<std::endl;
}