Cov4f.h

Go to the documentation of this file.

// AUTOMATICALLY GENERATED FILE - DO NOT EDIT
 
#ifndef EDM4EIC_Cov4f_H
#define EDM4EIC_Cov4f_H
 
#include <ostream>
 
#if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
#include "nlohmann/json_fwd.hpp"
#endif
 
namespace edm4eic {
 
/** @class Cov4f
 *  
 *  @author: 
 */
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; 
 } 
 } 
 }
 
 
};
 
std::ostream& operator<<(std::ostream& o, const edm4eic::Cov4f& value);
 
#if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
void to_json(nlohmann::json& j, const Cov4f& value);
#endif
 
} // namespace edm4eic
 
 
#endif