public double CalculateAverageMagneticFieldStrength(
double properTime_fm,
double conductivity_MeV
)
{
LCFFieldFunction function = (x_fm, y_fm, rapidity) => CalculateMagneticFieldInLCF(
properTime_fm, x_fm, y_fm, rapidity, conductivity_MeV).Norm;
return(LCFFieldAverager.AverageByBottomiumDistribution(function));
}
public double CalculateAverageMagneticFieldStrength(
double properTime_fm,
double x_fm,
double y_fm,
double conductivity_MeV
)
{
Func <double, double> integrand = rapidity => CalculateMagneticFieldInLCF(
properTime_fm, x_fm, y_fm, rapidity, conductivity_MeV).Norm;
return(LCFFieldAverager.AverageRapidityDependence(integrand));
}
/********************************************************************************************
* Constructors
********************************************************************************************/
public CollisionalElectromagneticField(
FireballParam param
)
{
NucleusPositionA = param.NucleusPositionA;
NucleusPositionB = param.NucleusPositionB;
Nucleus.CreateNucleusPair(param, out Nucleus nucleusA, out Nucleus nucleusB);
NucleusEMFA = new NucleusElectromagneticField(
param.EMFCalculationMethod,
param.PartonPeakRapidity,
nucleusA,
param.EMFQuadratureOrder);
NucleusEMFB = new NucleusElectromagneticField(
param.EMFCalculationMethod,
-param.PartonPeakRapidity,
nucleusB,
param.EMFQuadratureOrder);
LCFFieldAverager = new LCFFieldAverager(param);
}
