/*
 * Amplitude.h
 *
 *  Created on: Jan 14, 2012
 *      Author: Anastasia Karavdina
 */

#ifndef AMPLITUDE_H_
#define AMPLITUDE_H_

/*
 * Calculation of amplitude in helicity formalism
 */

#include "Chain.h"
#include <complex>
#include <vector>
#include "TLorentzVector.h"
class Amplitude {
private:
	Chain _ch;
	//vector<int> usedpos;//positions in the list [lfpcode] of used f.p

	complex<double> jtos1s2(double phi, double theta,
					 double j, double m,
					 double s1,double lya1,
					 double s2,double lya2, vector<double> l, double s,
					 double mass,
					 double mass1, double mass2, double mass0, double width);

	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);
	double momentum(double m, double m1, double m2, double m0);
	complex<double> pwa(double l, double mass, double mass1, double mass2,  double mass0, double width);
	double Fl(const int n, const double p);/// Calculation of barrier function

public:
	Amplitude(Chain ch);
	Amplitude(const Amplitude&);
	complex<double> Calc(list<TLorentzVector> lfp, int eJm, int elya);
	void Build(vector<int> lfpcode);
	vector<int> BusyParticle();
	Chain GetChain(){return _ch;};
	virtual ~Amplitude();

};

#endif /* AMPLITUDE_H_ */