RoughMultiLayerComputation Class Reference

Collaboration diagram for RoughMultiLayerComputation:

Public Member Functions
	RoughMultiLayerComputation (const ProcessedSample *p_sample)
void	compute (SimulationElement &elem) const

Private Member Functions
complex_t	get_refractive_term (size_t ilayer, double wavelength) const
complex_t	get_sum8terms (size_t ilayer, const SimulationElement &sim_element) const

Private Attributes
const ProcessedSample *	mp_sample

Detailed Description

Computes the diffuse reflection from the rough interfaces of a multilayer.

Used by DWBAComputation.

Definition at line 27 of file RoughMultiLayerComputation.h.

Constructor & Destructor Documentation

RoughMultiLayerComputation()

RoughMultiLayerComputation::RoughMultiLayerComputation

(

const ProcessedSample *

p_sample

)

Definition at line 43 of file RoughMultiLayerComputation.cpp.

    : mp_sample{p_sample}
{
}

Member Function Documentation

compute()

void RoughMultiLayerComputation::compute

(

SimulationElement &

elem

)

const

Definition at line 48 of file RoughMultiLayerComputation.cpp.

{
    if (elem.getAlphaMean() < 0.0)
        return;
    auto n_slices = mp_sample->numberOfSlices();
    kvector_t q = elem.getMeanQ();
    double wavelength = elem.getWavelength();
    double autocorr(0.0);
    complex_t crosscorr(0.0, 0.0);
 
    std::vector<complex_t> rterm(n_slices - 1);
    std::vector<complex_t> sterm(n_slices - 1);
 
    for (size_t i = 0; i + 1 < n_slices; i++) {
        rterm[i] = get_refractive_term(i, wavelength);
        sterm[i] = get_sum8terms(i, elem);
    }
    for (size_t i = 0; i + 1 < n_slices; i++) {
        const LayerRoughness* rough = mp_sample->bottomRoughness(i);
        if (rough)
            autocorr += std::norm(rterm[i]) * std::norm(sterm[i]) * rough->getSpectralFun(q);
    }
    // cross correlation between layers
    if (mp_sample->crossCorrelationLength() != 0.0) {
        for (size_t j = 0; j < n_slices - 1; j++) {
            for (size_t k = 0; k < n_slices - 1; k++) {
                if (j == k)
                    continue;
                crosscorr += rterm[j] * sterm[j] * mp_sample->crossCorrSpectralFun(q, j, k)
                             * std::conj(rterm[k]) * std::conj(sterm[k]);
            }
        }
    }
    //! @TODO clarify complex vs double
    elem.addIntensity((autocorr + crosscorr.real()) * M_PI / 4. / wavelength / wavelength);
}

References SimulationElement::addIntensity(), ProcessedSample::bottomRoughness(), ProcessedSample::crossCorrelationLength(), ProcessedSample::crossCorrSpectralFun(), get_refractive_term(), get_sum8terms(), SimulationElement::getAlphaMean(), SimulationElement::getMeanQ(), LayerRoughness::getSpectralFun(), SimulationElement::getWavelength(), M_PI, mp_sample, anonymous_namespace{SlicedCylindersBuilder.cpp}::n_slices(), ProcessedSample::numberOfSlices(), and anonymous_namespace{SlicedCylindersBuilder.cpp}::wavelength().

Here is the call graph for this function:

get_refractive_term()

complex_t RoughMultiLayerComputation::get_refractive_term ( size_t  ilayer, double  wavelength  ) const

private

Definition at line 85 of file RoughMultiLayerComputation.cpp.

{
    auto& slices = mp_sample->slices();
    return slices[ilayer].material().refractiveIndex2(wavelength)
           - slices[ilayer + 1].material().refractiveIndex2(wavelength);
}

References mp_sample, ProcessedSample::slices(), and anonymous_namespace{SlicedCylindersBuilder.cpp}::wavelength().

Referenced by compute().

Here is the call graph for this function:

get_sum8terms()

complex_t RoughMultiLayerComputation::get_sum8terms ( size_t  ilayer, const SimulationElement &  sim_element  ) const

private

Definition at line 92 of file RoughMultiLayerComputation.cpp.

{
    auto& slices = mp_sample->slices();
    auto p_fresnel_map = mp_sample->fresnelMap();
    const auto P_in_plus = p_fresnel_map->getInCoefficients(sim_element, ilayer);
    const auto P_out_plus = p_fresnel_map->getOutCoefficients(sim_element, ilayer);
 
    const auto P_in_minus = p_fresnel_map->getInCoefficients(sim_element, ilayer + 1);
    const auto P_out_minus = p_fresnel_map->getOutCoefficients(sim_element, ilayer + 1);
 
    complex_t kiz_plus = P_in_plus->getScalarKz();
    complex_t kfz_plus = P_out_plus->getScalarKz();
    complex_t qz1_plus = -kiz_plus - kfz_plus;
    complex_t qz2_plus = -kiz_plus + kfz_plus;
    complex_t qz3_plus = -qz2_plus;
    complex_t qz4_plus = -qz1_plus;
    double thickness = slices[ilayer].thickness();
    complex_t T_in_plus = P_in_plus->getScalarT() * exp_I(kiz_plus * thickness);
    complex_t R_in_plus = P_in_plus->getScalarR() * exp_I(-kiz_plus * thickness);
    complex_t T_out_plus = P_out_plus->getScalarT() * exp_I(kfz_plus * thickness);
    complex_t R_out_plus = P_out_plus->getScalarR() * exp_I(-kfz_plus * thickness);
 
    complex_t kiz_minus = P_in_minus->getScalarKz();
    complex_t kfz_minus = P_out_minus->getScalarKz();
    complex_t qz1_minus = -kiz_minus - kfz_minus;
    complex_t qz2_minus = -kiz_minus + kfz_minus;
    complex_t qz3_minus = -qz2_minus;
    complex_t qz4_minus = -qz1_minus;
 
    double sigma(0.0);
    if (const LayerRoughness* roughness = mp_sample->bottomRoughness(ilayer))
        sigma = roughness->getSigma();
    complex_t term1 = T_in_plus * T_out_plus * h_plus(qz1_plus * sigma);
    complex_t term2 = T_in_plus * R_out_plus * h_plus(qz2_plus * sigma);
    complex_t term3 = R_in_plus * T_out_plus * h_plus(qz3_plus * sigma);
    complex_t term4 = R_in_plus * R_out_plus * h_plus(qz4_plus * sigma);
    complex_t term5 =
        P_in_minus->getScalarT() * P_out_minus->getScalarT() * h_min(qz1_minus * sigma);
    complex_t term6 =
        P_in_minus->getScalarT() * P_out_minus->getScalarR() * h_min(qz2_minus * sigma);
    complex_t term7 =
        P_in_minus->getScalarR() * P_out_minus->getScalarT() * h_min(qz3_minus * sigma);
    complex_t term8 =
        P_in_minus->getScalarR() * P_out_minus->getScalarR() * h_min(qz4_minus * sigma);
 
    return term1 + term2 + term3 + term4 + term5 + term6 + term7 + term8;
}

References ProcessedSample::bottomRoughness(), exp_I(), ProcessedSample::fresnelMap(), IFresnelMap::getInCoefficients(), anonymous_namespace{RoughMultiLayerComputation.cpp}::h_min(), anonymous_namespace{RoughMultiLayerComputation.cpp}::h_plus(), mp_sample, and ProcessedSample::slices().

Referenced by compute().

Here is the call graph for this function:

Member Data Documentation

mp_sample

const ProcessedSample* RoughMultiLayerComputation::mp_sample

private

Definition at line 35 of file RoughMultiLayerComputation.h.

Referenced by compute(), get_refractive_term(), and get_sum8terms().

---

The documentation for this class was generated from the following files:

• /home/www/ba/Core/Computation/RoughMultiLayerComputation.h
• /home/www/ba/Core/Computation/RoughMultiLayerComputation.cpp