public void CalculateAverageMagneticFieldStrength()
{
FireballParam param = CreateFireballParamForAverageFieldStrengths();
CollisionalElectromagneticField emf = new CollisionalElectromagneticField(param);
double result = emf.CalculateAverageMagneticFieldStrength(0.4, QGPConductivity);
AssertHelper.AssertApproximatelyEqual(0.16239, result, 5);
}
private SpatialVector[] CalculateMagneticFieldValues()
{
CollisionalElectromagneticField emf
= new CollisionalElectromagneticField(CreateFireballParam());
SpatialVector[] fieldValues = new SpatialVector[Positions.Length];
for (int i = 0; i < Positions.Length; i++)
{
fieldValues[i] = emf.CalculateMagneticField(
Time, Positions[i].X, Positions[i].Y, Positions[i].Z, QGPConductivity);
}
return(fieldValues);
}
/********************************************************************************************
* Public members, functions and properties
********************************************************************************************/
public double CalculateAverageSpinStateOverlap(
BottomiumState tripletState,
double properTime_fm
)
{
FireballParam param = CreateFireballParam();
LCFFieldFunction mixingCoefficientSquared = (x, y, rapidity) =>
{
CollisionalElectromagneticField emf = new CollisionalElectromagneticField(param);
double B_per_fm2 = emf.CalculateMagneticFieldInLCF(
properTime_fm, x, y, rapidity, QGPConductivity_MeV).Norm;
return(CalculateSpinStateOverlap(tripletState, B_per_fm2));
};
LCFFieldAverager avg = new LCFFieldAverager(param);
return(avg.AverageByBottomiumDistribution(mixingCoefficientSquared));
}
