TypeDefs.h

Go to the documentation of this file.

#pragma once
 
#include <optional>
#include <string>
#include <unordered_map>
 
namespace iguana {

  using vector_element_t = double;

  struct Momentum3 {
      vector_element_t px;
      vector_element_t py;
      vector_element_t pz;
  };

  struct Momentum4 {
      vector_element_t px;
      vector_element_t py;
      vector_element_t pz;
      vector_element_t E;
  };


  namespace particle {
    // clang-format off

    enum PDG {
      electron     = 11,
      photon       = 22,
      proton       = 2212,
      antiproton   = -2212,
      neutron      = 2112,
      antineutron  = -2112,
      pi_plus      = 211,
      pi_minus     = -211,
      kaon_plus    = 321,
      kaon_minus   = -321
    };

    const std::unordered_map<PDG, std::string> name{
      { electron, "electron" },
      { photon, "photon" },
      { proton, "proton" },
      { antiproton, "antiproton" },
      { neutron, "neutron" },
      { antineutron, "antineutron" },
      { pi_plus, "pi_plus" },
      { pi_minus, "pi_minus" },
      { kaon_plus, "kaon_plus" },
      { kaon_minus, "kaon_minus" }
    };

    const std::unordered_map<PDG, std::string> title{
      { electron, "e^{-}" },
      { photon, "#gamma" },
      { proton, "p" },
      { antiproton, "#bar{p}" },
      { neutron, "n" },
      { antineutron, "#bar{n}" },
      { pi_plus, "#pi^{+}" },
      { pi_minus, "#pi^{-}" },
      { kaon_plus, "K^{+}" },
      { kaon_minus, "K^{-}" }
    };

    const std::unordered_map<PDG, double> mass{
      { electron, 0.000511 },
      { photon, 0.0 },
      { proton, 0.938272 },
      { antiproton, 0.938272 },
      { neutron, 0.939565 },
      { antineutron, 0.939565 },
      { pi_plus, 0.139570 },
      { pi_minus, 0.139570 },
      { kaon_plus, 0.493677 },
      { kaon_minus, 0.493677 }
    };

    template <typename VALUE_TYPE>
    std::optional<VALUE_TYPE> const get(std::unordered_map<PDG,VALUE_TYPE> const& property, PDG const& pdg_code)
    {
      if(auto const& it = property.find(pdg_code); it != property.end())
        return it->second;
      return std::nullopt;
    }

    template <typename VALUE_TYPE>
    std::optional<VALUE_TYPE> const get(std::unordered_map<PDG,VALUE_TYPE> const& property, int const& pdg_code)
    {
      return get(property, static_cast<particle::PDG>(pdg_code));
    }
 
    // clang-format on
  }


  inline bool IsValidSector(int const& sec)
  {
    return sec >= 1 && sec <= 6;
  }

  enum DetectorType {
    UNDEFINED = 0,
    BMT       = 1,
    BST       = 2,
    CND       = 3,
    CTOF      = 4,
    CVT       = 5,
    DC        = 6,
    ECAL      = 7,
    FMT       = 8,
    FT        = 9,
    FTCAL     = 10,
    FTHODO    = 11,
    FTOF      = 12,
    FTTRK     = 13,
    HTCC      = 15,
    LTCC      = 16,
    RF        = 17,
    RICH      = 18,
    RTPC      = 19,
    HEL       = 20,
    BAND      = 21,
    RASTER    = 22,
    URWELL    = 23,
    AHDC      = 24,
    ATOF      = 25,
    RECOIL    = 26,
    TARGET    = 100,
    MAGNETS   = 101,
  };

  class DetectorLayer
  {
    public:
      static int const CND_INNER  = 1;
      static int const CND_MIDDLE = 2;
      static int const CND_OUTER  = 3;
 
      static int const PCAL_U = 1;
      static int const PCAL_V = 2;
      static int const PCAL_W = 3;
      static int const PCAL_Z = 9; // layer number used to define the longitudinal coordinate of the cluster
 
      static int const EC_INNER_U = 4;
      static int const EC_INNER_V = 5;
      static int const EC_INNER_W = 6;
      static int const EC_INNER_Z = 9; // layer number used to define the longitudinal coordinate of the cluster
 
      static int const EC_OUTER_U = 7;
      static int const EC_OUTER_V = 8;
      static int const EC_OUTER_W = 9;
      static int const EC_OUTER_Z = 9; // layer number used to define the longitudinal coordinate of the cluster
 
      static int const PCAL     = PCAL_U;
      static int const EC_INNER = EC_INNER_U;
      static int const EC_OUTER = EC_OUTER_U;
 
      static int const FTOF1A = 1;
      static int const FTOF1B = 2;
      static int const FTOF2  = 3;
 
      static int const TARGET_CENTER     = 1;
      static int const TARGET_DOWNSTREAM = 2;
      static int const TARGET_UPSTREAM   = 3;
 
      static int const FTTRK_MODULE1 = 1;
      static int const FTTRK_MODULE2 = 2;
      static int const FTTRK_LAYER1  = 1;
      static int const FTTRK_LAYER2  = 2;
      static int const FTTRK_LAYER3  = 3;
      static int const FTTRK_LAYER4  = 4;
 
      static int const RICH_MAPMT      = 1;
      static int const RICH_AEROGEL_B1 = 2;
      static int const RICH_AEROGEL_B2 = 3;
      static int const RICH_AEROGEL_L1 = 4;
  };
 
}