GeoMaterialConverter.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/

// GeoMaterialConverter.cxx, (c) ATLAS Detector software
 
// Trk
#include "TrkDetDescrGeoModelCnv/GeoMaterialConverter.h"
// GeoModelKernel
#include "GeoModelKernel/GeoElement.h"
#include "GeoModelKernel/GeoMaterial.h"
#include "GeoModelKernel/Units.h"
namespace {
    constexpr double s_densityCnvFactor = 1. / GeoModelKernelUnits::gram;
}
Trk::Material Trk::GeoMaterialConverter::convert(const GeoMaterial* gm) {
    // get the obvious things
    float x0 = gm->getRadLength();
    float l0 = gm->getIntLength();
    float rho = gm->getDensity() * s_densityCnvFactor;
    // float dEdX = gm->getDeDxMin();
    float A = 0.;
    float Z = 0.;
    // get the number of elements that build the material
    int numberOfElements = gm->getNumElements();
    // loop over the elements
    for (int iEl = 0; iEl < numberOfElements; iEl++) {
        const GeoElement* geoEl = gm->getElement(iEl);
        float fraction = gm->getFraction(iEl);
        A += fraction * (geoEl->getA() / (GeoModelKernelUnits::gram));
        Z += fraction * (geoEl->getZ());
    }
    // return the result
    return {x0, l0, A, Z, rho};
}
 
bool Trk::GeoMaterialConverter::dummy_material(const GeoMaterial* mat) {
 
    bool dummyMat = false;
 
    const std::string& matName = mat->getName();
 
    if (matName == "special::Ether")
        dummyMat = true;
    if (matName == "WorldLog:Air")
        dummyMat = true;
    if (matName == "std::Air")
        dummyMat = true;
    if (matName == "Air")
        dummyMat = true;
 
    return dummyMat;
}