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