//#include "TComplex.h"
//#include "ClebschGordan.h"
#include <complex>
#include <vector>

#ifndef _FULLAMP_
#define _FULLAMP_

class FullAmp{
 private:
  /* //helicity amplitude */
  /* std::complex<double> NFj(double j, std::vector<double> l, std::vector<double> s, double s1, double lya1, double s2, double lya2); */

  /// one term of helicity amplitude for l,s
  std::complex<double> nfls(double j, double l, double s, double s1, double lya1, double s2, double lya2, 
			     double mass, double mass1, double mass2, double mass0, double width);
  /// for nonresonans case
  std::complex<double> nfls(double j, double l, double s, double s1, double lya1, double s2, double lya2);


  ///Partial wave amplitude a_ls, responcble for dynamics
  std::complex<double> pwa(double l,  double mass, double mass1, double mass2,  double mass0, double width);

  /// Calculation of momentum for decay mass m doubleo particles with masses m1, m2
  /// through resonance with mass m0
  double momentum(double m, double m1, double m2, double m0);

 public:
  ///Constructor & Dectructor
  FullAmp();
  virtual ~FullAmp(){};
  ///-------------------------
  
  /// j to s1 and s2
  std::complex<double> jtos1s2(double phi, double theta,
			       double j, double m,
			       double s1, double lya1,
			       double s2, double lya2, 
			       std::vector<double> l, double s,
			       double mass, double mass1, double mass2, double mass0, double width);

  ///j to s1 and s2 (noresonance decay)
  std::complex<double> jtos1s2(double phi, double theta, 
			       double j, double m,
			       double s1,double lya1,
			       double s2,double lya2, std::vector<double> l, double s);

  ///calculation of barrier function
  double Fl(const int n, const double p);

  //  ClassDef(FullAmp, 1); // Calculation of full amplitude

};
#endif