predefined_materials.hpp

// This file is part of the ACTS project.
//
// Copyright (C) 2016 CERN for the benefit of the ACTS project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
#pragma once
 
// Project include(s)
#include "detray/definitions/algebra.hpp"
#include "detray/definitions/units.hpp"
#include "detray/material/material.hpp"
 
// System include(s)
#include <limits>
 
namespace detray {

// Vacuum
DETRAY_DECLARE_MATERIAL(vacuum, std::numeric_limits<scalar_t>::max(),
                        std::numeric_limits<scalar_t>::max(), 0.f, 0.f, 0.f,
                        material_state::e_unknown);
 
// H₂ (1): Hydrogen Gas
DETRAY_DECLARE_MATERIAL(hydrogen_gas, 7.526E3f * unit<scalar_t>::m,
                        6.209E3f * unit<scalar_t>::m, 2.f * 1.008f, 2.f * 1.f,
                        static_cast<scalar_t>(8.376E-5 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// H₂ (1): Hydrogen Liquid
DETRAY_DECLARE_MATERIAL(hydrogen_liquid, 8.904f * unit<scalar_t>::m,
                        7.346f * unit<scalar_t>::m, 2.f * 1.008f, 2.f * 1.f,
                        static_cast<scalar_t>(0.07080f * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_liquid);
 
// He (2): Helium Gas
DETRAY_DECLARE_MATERIAL(helium_gas, 5.671E3f * unit<scalar_t>::m,
                        4.269E3f * unit<scalar_t>::m, 4.003f, 2.f,
                        static_cast<scalar_t>(1.663E-4 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// Be (4)
DETRAY_DECLARE_MATERIAL(beryllium, 352.8f * unit<scalar_t>::mm,
                        421.0f * unit<scalar_t>::mm, 9.012f, 4.f,
                        static_cast<scalar_t>(1.848 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// C (6): Carbon (amorphous)
DETRAY_DECLARE_MATERIAL(carbon_gas, 213.5f * unit<scalar_t>::mm,
                        429.0f * unit<scalar_t>::mm, 12.01f, 6.f,
                        static_cast<scalar_t>(2.0 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// N₂ (7): Nitrogen Gas
DETRAY_DECLARE_MATERIAL(nitrogen_gas, 3.260E+02f * unit<scalar_t>::m,
                        7.696E+02f * unit<scalar_t>::m, 2.f * 14.007f,
                        2.f * 7.f,
                        static_cast<scalar_t>(1.165E-03 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// O₂ (8): Oxygen Gas
DETRAY_DECLARE_MATERIAL(oxygen_gas, 2.571E+02f * unit<scalar_t>::m,
                        6.772E+02f * unit<scalar_t>::m, 2.f * 15.999f,
                        2.f * 8.f,
                        static_cast<scalar_t>(1.332E-3 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// O₂ (8): Oxygen liquid
DETRAY_DECLARE_MATERIAL(oxygen_liquid, 300.1f * unit<scalar_t>::mm,
                        790.3f * unit<scalar_t>::mm, 2.f * 15.999f, 2.f * 8.f,
                        static_cast<scalar_t>(1.141 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_liquid);
 
// Al (13)
DETRAY_DECLARE_MATERIAL(aluminium, 88.97f * unit<scalar_t>::mm,
                        397.0f * unit<scalar_t>::mm, 26.98f, 13.f,
                        static_cast<scalar_t>(2.699 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Si (14)
DETRAY_DECLARE_MATERIAL(silicon, 93.7f * unit<scalar_t>::mm,
                        465.2f * unit<scalar_t>::mm, 28.0855f, 14.f,
                        static_cast<scalar_t>(2.329 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Si (14) with density effect data
DETRAY_DECLARE_MATERIAL_WITH_DED(silicon_with_ded, 93.7f * unit<scalar_t>::mm,
                                 465.2f * unit<scalar_t>::mm, 28.0855f, 14.f,
                                 static_cast<scalar_t>(2.329 * unit<double>::g /
                                                       unit<double>::cm3),
                                 material_state::e_solid, 0.1492f, 3.2546f,
                                 0.2015f, 2.8716f, 173.0f, 4.4355f, 0.14f);
 
// Ar (18): Argon gas
DETRAY_DECLARE_MATERIAL(argon_gas, 1.176E+02f * unit<scalar_t>::m,
                        7.204E+02f * unit<scalar_t>::m, 39.948f, 18.f,
                        static_cast<scalar_t>(1.662E-03 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// Ar (18): Argon liquid
DETRAY_DECLARE_MATERIAL(argon_liquid, 14.f * unit<scalar_t>::cm,
                        85.77f * unit<scalar_t>::cm, 39.948f, 18.f,
                        static_cast<scalar_t>(1.396 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_liquid);
 
// Fe (26)
DETRAY_DECLARE_MATERIAL(iron, 1.757f * unit<scalar_t>::cm,
                        16.77f * unit<scalar_t>::cm, 55.845f, 26.f,
                        static_cast<scalar_t>(7.874 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Fe (26) with density effect data
DETRAY_DECLARE_MATERIAL_WITH_DED(iron_with_ded, 1.757f * unit<scalar_t>::cm,
                                 16.77f * unit<scalar_t>::cm, 55.845f, 26.f,
                                 static_cast<scalar_t>(7.874 * unit<double>::g /
                                                       unit<double>::cm3),
                                 material_state::e_solid, 0.14680f, 2.9632f,
                                 -0.0012f, 3.1531f, 286.0f, 4.2911f, 0.12f);
 
// Copper (29)
DETRAY_DECLARE_MATERIAL(copper, 1.436f * unit<scalar_t>::cm,
                        15.32f * unit<scalar_t>::cm, 63.546f, 29.f,
                        static_cast<scalar_t>(8.960 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Copper (29) with density effect data
DETRAY_DECLARE_MATERIAL_WITH_DED(copper_with_ded, 1.436f * unit<scalar_t>::cm,
                                 15.32f * unit<scalar_t>::cm, 63.546f, 29.f,
                                 static_cast<scalar_t>(8.960 * unit<double>::g /
                                                       unit<double>::cm3),
                                 material_state::e_solid, 0.14339f, 2.9044f,
                                 -0.0254f, 3.2792f, 322.0f, 4.4190f, 0.08f);
 
// W (74)
DETRAY_DECLARE_MATERIAL(tungsten, 3.504f * unit<scalar_t>::mm,
                        99.46f * unit<scalar_t>::mm, 183.84f, 74.f,
                        static_cast<scalar_t>(19.3 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Au (79)
DETRAY_DECLARE_MATERIAL(gold, 3.344f * unit<scalar_t>::mm,
                        101.6f * unit<scalar_t>::mm, 196.97f, 79.f,
                        static_cast<scalar_t>(19.32 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);

 
// Be (4)
DETRAY_DECLARE_MATERIAL(beryllium_tml, 352.8f * unit<scalar_t>::mm,
                        407.f * unit<scalar_t>::mm, 9.012f, 4.f,
                        static_cast<scalar_t>(1.848 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);
 
// Si (14)
DETRAY_DECLARE_MATERIAL(silicon_tml, 95.7f * unit<scalar_t>::mm,
                        465.2f * unit<scalar_t>::mm, 28.03f, 14.f,
                        static_cast<scalar_t>(2.32 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);

 
// Air (dry, 1 atm)
// @note:
// https://pdg.lbl.gov/2020/AtomicNuclearProperties/HTML/air_dry_1_atm.html
// @note: Ar from Wikipedia (https://en.wikipedia.org/wiki/Molar_mass)
DETRAY_DECLARE_MATERIAL(air, 3.039E+02f * unit<scalar_t>::m,
                        7.477E+02f * unit<scalar_t>::m, 28.97f, 14.46f,
                        static_cast<scalar_t>(1.205E-03 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// (CH3)2CHCH3 Gas
// @note: (X0, L0, mass_rho) from https://pdg.lbl.gov/2005/reviews/atomicrpp.pdf
// @note: Ar from Wikipedia (https://en.wikipedia.org/wiki/Isobutane)
// @note: Z was calculated by simply summing the number of atoms. Surprisingly
// it seems the right value because Z/A is 0.58496, which is the same with <Z/A>
// in the pdg reference
DETRAY_DECLARE_MATERIAL(isobutane, 1693E+02f * unit<scalar_t>::mm,
                        288.3f * unit<scalar_t>::mm, 58.124f, 34.f,
                        static_cast<scalar_t>(2.67 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// C3H8 Gas
// @note: (X0, L0, mass_rho) from
// https://pdg.lbl.gov/2020/AtomicNuclearProperties/HTML/propane.html
// @note: Ar from Wikipedia (https://en.wikipedia.org/wiki/Propane)
DETRAY_DECLARE_MATERIAL(propane, 2.429E+02f * unit<scalar_t>::m,
                        4.106E+02f * unit<scalar_t>::m, 44.097f, 26.f,
                        static_cast<scalar_t>(1.868E-03 * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_gas);
 
// Cesium Iodide (CsI)
// https://pdg.lbl.gov/2023/AtomicNuclearProperties/HTML/cesium_iodide_CsI.html
DETRAY_DECLARE_MATERIAL(cesium_iodide, 1.86f * unit<scalar_t>::cm,
                        38.04f * unit<scalar_t>::cm, 259.81f, 108.f,
                        static_cast<scalar_t>(4.510f * unit<double>::g /
                                              unit<double>::cm3),
                        material_state::e_solid);

DETRAY_DECLARE_MATERIAL_WITH_DED(
    cesium_iodide_with_ded, 1.86f * unit<scalar_t>::cm,
    38.04f * unit<scalar_t>::cm, 259.81f, 108.f,
    static_cast<scalar_t>(4.510f * unit<double>::g / unit<double>::cm3),
    material_state::e_solid, 0.25381f, 2.6657f, 0.0395f, 3.3353f, 553.1f,
    6.2807f, 0.00f);
 
}  // namespace detray