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)); }