IAP GITLAB

Skip to content

Trajectories crossing observation plane

Running hybrid_MC as at today (with cut lowered to 3 GeV), a sizeable fraction (> 20 %) of the showers fail because some trajectories exceed the observation plane and go below ground. gdb output:

#7  corsika::ShowerAxis::getProjectedX (this=this@entry=0x7fff07926b30, p=...) at /home/ka/ka_iap/ka_xxxxxx/corsika/corsika/detail/media/ShowerAxis.inl:109
109         return getX(projectedLength);
(gdb) p p
$3 = (const corsika::Point &) @0x7fff07926320: {<corsika::BaseVector<phys::units::dimensions<1, 0, 0, 0, 0, 0, 0, 0> >> = {quantityVector_ = {
      eigenVector_ = {<Eigen::PlainObjectBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >> = {<Eigen::MatrixBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >> = {<Eigen::DenseBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >> = {<Eigen::DenseCoeffsBase<Eigen::Matrix<double, 3, 1, 0, 3, 1>, 3>> = {<Eigen::DenseCoeffsBase<Eigen::Matrix<double, 3, 1, 0, 3, 1>, 1>> = {<Eigen::DenseCoeffsBase<Eigen::Matrix<double, 3, 1, 0, 3, 1>, 0>> = {<Eigen::EigenBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >> = {<No data fields>}, <No data fields>}, <No data fields>}, <No data fields>}, <No data fields>}, <No data fields>}, m_storage = {m_data = {array = {1116.5140376118152, 490.03453658239732, 6368384.0377964033}}}}, <No data fields>}},
    cs_ = {<std::__shared_ptr<corsika::CoordinateSystem const, (__gnu_cxx::_Lock_policy)2>> = {<std::__shared_ptr_access<corsika::CoordinateSystem const, (__gnu_cxx::_Lock_policy)2, false, false>> = {<No data fields>}, _M_ptr = 0x173b1f70, _M_refcount = {_M_pi = 0x173aa4b0}}, <No data fields>}}, <No data fields>}

Note that 6368384.0377964033 < mean Earth radius. The error seems to occur only for the magnetic tracking (default settings).

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information