Cov4f.h

Go to the documentation of this file.

// AUTOMATICALLY GENERATED FILE - DO NOT EDIT
 
#ifndef EICD_Cov4f_H
#define EICD_Cov4f_H
 
#include <ostream>
 
#ifdef PODIO_JSON_OUTPUT
#include "nlohmann/json.hpp"
#endif
 
namespace eicd {
 
 
class Cov4f {
public:
  float xx{};
  float yy{};
  float zz{};
  float tt{};
  float xy{};
  float xz{};
  float xt{};
  float yz{};
  float yt{};
  float zt{};
 
 Cov4f() : xx{0}, yy{0}, zz{0}, tt{0}, xy{0}, xz{0}, xt{0}, yz{0}, yt{0}, zt{0} {}
 Cov4f(double vx, double vy, double vz, double vt,
 double vxy = 0, double vxz = 0, double vxt = 0,
 double vyz = 0, double vyt = 0, double vzt = 0)
 : xx{static_cast<float>(vx)}, yy{static_cast<float>(vy)}, zz{static_cast<float>(vz)}, tt{static_cast<float>(vt)},
 xy{static_cast<float>(vxy)}, xz{static_cast<float>(vxz)}, xt{static_cast<float>(vxt)},
 yz{static_cast<float>(vyz)}, yt{static_cast<float>(vyt)}, zt{static_cast<float>(vzt)} {}
 float operator()(unsigned i, unsigned j) const {
 // diagonal
 if (i == j) {
 return *(&xx + i);
 // off-diagonal, can probably be done with less if statements 
 } else {
 if (i > j) { 
 std::swap(i,j); 
 } 
 if (i == 0) { 
 return *(&xy + j - 1); 
 } else if (i == 1) { 
 return *(&yz + j - 2); 
 } else { 
 return zt; 
 } 
 } 
 }
 
 
};
 
inline std::ostream& operator<<(std::ostream& o, const eicd::Cov4f& value) {
  o << value.xx << " ";
  o << value.yy << " ";
  o << value.zz << " ";
  o << value.tt << " ";
  o << value.xy << " ";
  o << value.xz << " ";
  o << value.xt << " ";
  o << value.yz << " ";
  o << value.yt << " ";
  o << value.zt << " ";
 
  return o;
}
 
#ifdef PODIO_JSON_OUTPUT
inline void to_json(nlohmann::json& j, const Cov4f& value) {
  j = nlohmann::json{
  {"xx", value.xx}
  ,{"yy", value.yy}
  ,{"zz", value.zz}
  ,{"tt", value.tt}
  ,{"xy", value.xy}
  ,{"xz", value.xz}
  ,{"xt", value.xt}
  ,{"yz", value.yz}
  ,{"yt", value.yt}
  ,{"zt", value.zt}
  };
}
#endif
 
} // namespace eicd
 
 
#endif