/*
 * TwoParticleState.h
 *
 *  Created on: Jan 11, 2012
 *      Author: Anastasia Karavdina
 */

#ifndef TWOPARTICLESTATE_H_
#define TWOPARTICLESTATE_H_

/*
 * Two-Particle Plane-Wave Helicity States in the rest frame of resonance
 */
#include "Particle.h"
#include "ResonanceP.h"
#include "vector"
//#include "TLorentzVector.h"
#include "TVector3.h"
#include <iostream>

class TwoParticleState {
public:

	TwoParticleState(Particle p1, TVector3 mom1,
									   Particle p2, TVector3 mom2,
									   ResonanceP res);

	TwoParticleState(Particle p1, Particle p2, ResonanceP res);

	TwoParticleState(const TwoParticleState& st);
	//TwoParticleState();
	virtual ~TwoParticleState();

	//Recalculation of particle momentum in the "rest frame of resonance"
	void ResRestFrame(TVector3 mom1, TVector3 mom2);
	///Getting value:
	    vector<double> AngMom();
	    double Mom();
	    double Theta();
	    double Phi();
	    double S1();
	    double S2();
	    double J();
	    double mR();
	    double wR();
	    vector<double> Lya1();
	    vector<double> Lya2();
	    vector<double> Jm();
	    double mass1();
	    double mass2();
	    int CodeR();
	    int Code1();
	    int Code2();
	    vector<int> ListFP1(){return _lfp1;}
	    vector<int> ListFP2(){return _lfp2;}
	  /*  void setListFP1(vector<int> lfp){_lfp1 = lfp;}
	    void setListFP2(vector<int> lfp){_lfp2 = lfp;}*/
	    void setListFP1(vector<int> lfp);
	    void setListFP2(vector<int> lfp);
	    pair<Particle,Particle> GetFP();//get pair of final particles
	    ResonanceP GetIP();//get initial particles

private:
	double _p;//Momentum magnitude
	double _theta;//theta angle of 1st particle
	double _phi;//phi angle of 1st particle
	double _s1;//spin of 1st particle
	vector<double> _lya1;//helicity of 1st particle
	double _s2;//spin of 2nd particle
	vector<double> _lya2;//helicity of 2nd particle
	double _J;//Total moment of system
	vector<double> _Jm;//Total moment projection on z-axis
	vector<double> _l;//angular moment
	double _mR;//resonance mass
	double _wR;//resonance width
	double _m1;//mass of 1st particle
	double _m2;//mass of second particle
	//_codeR->_code1+_code2
	int _codeR;
	int _code1;
	int _code2;
	///for momentum of f.p calculation if one of final particles is resonance
	vector<int> _lfp1;//list of particles which built particle#1
	vector<int> _lfp2;//list of particles which built particle#2
	Particle _p1;
	Particle _p2;
	ResonanceP _res;
	void AngMomCalc(Particle p1, Particle p2, Particle res);//Calculation of angular moment

};

#endif /* TWOPARTICLESTATE_H_ */