using System;

namespace Science.Physics.GeneralPhysics
{
	/// <summary>
	/// Decay
	/// </summary>
	public class Decay
	{
		public Decay()
		{
		}
		private double lambda, t;
		public double DecayConstant
		{
			get{return lambda;}
			set{lambda = value;t=Math.Log(2.0)/lambda;}
		}
		public double HalfLife
		{
			get{return t;}
			set{t = value;lambda=Math.Log(2.0)/t;}
		}
		public static Nucleus AlphaDecay(Nucleus X)
		{
			Nucleus Y = new Nucleus();
			Y.MassNumber = X.MassNumber - 4;
			Y.AtomicNumber = X.AtomicNumber - 2;
			Y.NeutronNumber = X.NeutronNumber - 2;
			return Y;
		}
		public static Nucleus BetaDecay(Nucleus X)
		{
			Nucleus Y = new Nucleus();
			Y.MassNumber = X.MassNumber;
			Y.AtomicNumber = X.AtomicNumber + 1;
			Y.NeutronNumber = X.NeutronNumber - 1;
			return Y;
		}
		public static Nucleus GammaDecay(Nucleus X)
		{
			Nucleus Y = new Nucleus();
			Y.MassNumber = X.MassNumber;
			Y.AtomicNumber = X.AtomicNumber;
			Y.NeutronNumber = X.NeutronNumber;
			return Y;
		}
	}
}