package affectors;

import main.Node;
import main.Globals;

/** This is one half of an affector that dephosphorylates a single node.  Use this
affector in the product species.  Use the A version in the phosphorylated species.
For either membrane-bound or cytoplasmic phosphatases, as long as target is cytoplasmic. 

<p> <b>Formula</b>
<br>dPRODUCT/dt = Pmax * SUBSTRATE_P * (PHOSPHO^nu / (K_PHOSPHO^nu + PHOSPHO^nu))

<p> <b>Parameters</b>
<table width="500" border="1">
<tr><td width="150">Substrate [SUBSTRATE_P]</td><td> The unphosphorylated Node </td></tr>
<tr><td> Phosphotase [PHOSPHO]</td><td> The phosphatase Node.</td></tr>
<tr><td> Half-max activation [K_PHOSPHO_SUBSTRATE_P]</td><td> The concentration of phosphatase at which it de-phosphorylates at half its maximum rate.</td></tr>
<tr><td> Cooperativity [nu_PHOSPHO_SUBSTRATE_P]</td><td> The non-linearity in phosphatase activity as phosphatase concentration increases.</td></tr>
<tr><td> Max reaction rate [P_PHOSPHO]</td><td> The maximum reaction rate.</td></tr>
</table>

<p> <b>Usage</b>
<code><font size = "3"><br>&SUBSTRATE_P
<blockquote>    &DephosphorylationB_ICAff  SUBSTRATE_P	PHOSPHO	K_PHOSPHO_SUBSTRATE_P	nu_PHOSPHO_SUBSTRATE_P	P_PHOSPHO
</blockquote>
&endSUBSTRATE_P</code></font>

<p>@see DephosphorylationA_ICAff
<br>@see PhosphorylationA_ICAff

*/
public class DephosphorylationB_ICAff extends Affector {	
	/** The maximal reaction rate. */
	int		kappaParam, nuParam, dMaxParam;
	
	static final String		desc = "Intracellular dephosphorylation of a substrate; version B for product";
	static final String	[]	nodeDescriptions = {"Phosphorylated substrate", "Phosphatase"};
	static final String	[]	paramDescriptions = {"Kappa: half-maximal activity of the phosphatase",
										"nu: cooperativity of reaction",
										"Dmax: maximal reaction rate"};
	
	public DephosphorylationB_ICAff()	{}

	protected void setLabelsAndTypes()	{
		setDescriptions(this.desc, nodeDescriptions, paramDescriptions);

		this.Type[GUI_CAPABLE] = 1;
		this.Type[CERTIFICATION] = Affector.RETURNS_DERIV;
		this.Type[MATHTYPE] = Affector.FF;
		this.Type[TERMTYPE] = Affector.CONVERSION;
	}
	
	public void setParameterNumbers(int [] param_nums)
	{
		kappaParam = param_nums[0];
		nuParam = param_nums[1];
		dMaxParam = param_nums[2];
	}

	public float getValue(Node which_node)	{
		float total_phosphatase = Nodes[1].getIntegrationValue();
		float source = Nodes[0].getIntegrationValue();
		float theTerm = params[dMaxParam] * source * 
					Phi(total_phosphatase, params[kappaParam], params[nuParam]);

		return Globals.characteristicTime * theTerm;
	}	
}