/********************************************************************************************
* Public static members, functions and properties
********************************************************************************************/
/********************************************************************************************
* Constructors
********************************************************************************************/
public QQFreeState(
QQStateParam param
)
: base(param)
{
InitSolver();
AssertInputValid();
}
/********************************************************************************************
* Constructors
********************************************************************************************/
public QQState(
QQStateParam param
)
{
Param = param.Clone();
AssertValidParam();
SetCouplings();
SetDebyeMass();
SetPotential();
WaveFunction_fm = new Complex[Param.StepNumber + 1];
}
/********************************************************************************************
* Public static members, functions and properties
********************************************************************************************/
/********************************************************************************************
* Constructors
********************************************************************************************/
public QQBoundState(
QQStateParam param,
int quantumNumberN
)
: base(param)
{
InitSolver();
UseFixedAlpha = false;
QuantumNumberN = quantumNumberN;
Radius_fm = Solver.PositionValues;
WaveFunctionNorm = new double[Param.StepNumber + 1];
AssertInputValid();
}
private QQStateParam GetFreeStateParam(
double gluonEnergy_MeV,
int quantumNumberL
)
{
QQStateParam param = BoundParam.Clone();
param.AccuracyAlpha = 1e-6;
param.AccuracyWaveFunction = 1e-9;
param.AggressivenessAlpha = 0.5;
param.AggressivenessEnergy = quantumNumberL > 0 ? 0.02 : 40;
param.ColorState = ColorState.Octet;
param.Energy_MeV = gluonEnergy_MeV - MinEnergy_MeV;
param.GammaDamp_MeV = -GammaDamp_MeV / 8.0; /*GammaDamp_MeV = 0*/
param.MaxShootingTrials = 0;
param.PotentialType = PotentialType.Tzero_NoString;
param.QuantumNumberL = quantumNumberL;
param.Temperature_MeV = Temperature_MeV;
return(param);
}