Mask/src/DecaySystem.cpp

#include "DecaySystem.h"
#include "RandomGenerator.h"

namespace Mask {

	DecaySystem::DecaySystem() :
		ReactionSystem()
	{
		nuclei.resize(3);
	}
	
	DecaySystem::DecaySystem(std::vector<int>& z, std::vector<int>& a) :
		ReactionSystem()
	{
		nuclei.resize(3);
		SetNuclei(z, a);
	}
	
	DecaySystem::~DecaySystem() {}
	
	bool DecaySystem::SetNuclei(std::vector<int>& z, std::vector<int>& a) {
		if(z.size() != a.size() || z.size() != 2) {
			return false;
		}

		step1.SetNuclei(z[0], a[0], 0, 0, z[1], a[1]);
		SetSystemEquation();
		return true;
	}

	const std::vector<Nucleus>& DecaySystem::GetNuclei()
	{
		nuclei[0] = step1.GetTarget();
		nuclei[1] = step1.GetEjectile();
		nuclei[2] = step1.GetResidual();

		return nuclei;
	}
	
	void DecaySystem::LinkTarget() {
		step1.SetLayeredTarget(&target);
	
		rxnLayer = target.FindLayerContaining(step1.GetTarget().GetZ(), step1.GetTarget().GetA());
		if(rxnLayer != -1) {
			step1.SetRxnLayer(rxnLayer);
			target_set_flag = true;
		} else {
			throw ReactionLayerException();
		}
	}
	
	void DecaySystem::SetSystemEquation() {
		m_sys_equation = step1.GetTarget().GetIsotopicSymbol();
		m_sys_equation += "-> ";
		m_sys_equation += step1.GetEjectile().GetIsotopicSymbol();
		m_sys_equation += "+";
		m_sys_equation += step1.GetResidual().GetIsotopicSymbol();
	}
	
	void DecaySystem::RunSystem() {
		//Link up the target if it hasn't been done yet
		if(!target_set_flag) {
			LinkTarget();
		}
	
		double rxnTheta = std::acos(decay1dist.GetRandomCosTheta());
		double rxnPhi = (*m_phi1Range)(RandomGenerator::GetInstance().GetGenerator());
		step1.SetPolarRxnAngle(rxnTheta);
		step1.SetAzimRxnAngle(rxnPhi);
	
		step1.TurnOnResidualEloss();
		step1.Calculate();

	}

}