ObservableTypeId.h

Go to the documentation of this file.

#ifndef OBSERVABLE_TYPE_ID_H
#define OBSERVABLE_TYPE_ID_H
 
#include <stdexcept>
 
enum class PolarizationKind : long {
    Fermion = 1,
    LongitudinalVector = 2,
    TransverseVector = 3,
};
 
inline constexpr long kTransverseFractionTypeId = 932;
 
inline constexpr long kAlphaKTypeId = 933;
 
inline constexpr long kLegacyFermionPolarizationTypeId = 92015;
 
inline constexpr long kLegacyLongitudinalVectorPolarizationTypeId = 921423;
 
constexpr long make_polarization_type_id(
    PolarizationKind kind,
    unsigned int daughter_position
) {
    const long kind_digit = static_cast<long>(kind);
    if (kind_digit < 1 || kind_digit > 3) {
        throw std::invalid_argument("polarization kind must be encoded by i=1,2,3");
    }
    if (daughter_position < 1 || daughter_position > 9) {
        throw std::invalid_argument("polarized daughter position must be in [1,9]");
    }
    return 9200L + 10L * kind_digit + static_cast<long>(daughter_position);
}
 
static_assert(make_polarization_type_id(PolarizationKind::Fermion, 2) == 9212);
static_assert(make_polarization_type_id(PolarizationKind::LongitudinalVector, 1) == 9221);
static_assert(make_polarization_type_id(PolarizationKind::TransverseVector, 1) == 9231);
 
#endif // OBSERVABLE_TYPE_ID_H