void IntegrateEllipsoids::calculateE1( const Geometry::DetectorInfo &detectorInfo) { for (size_t i = 0; i < detectorInfo.size(); ++i) { if (detectorInfo.isMonitor(i)) continue; // skip monitor if (!detectorInfo.isMasked(i)) continue; // edge is masked so don't check if not masked const auto &det = detectorInfo.detector(i); double tt1 = det.getTwoTheta(V3D(0, 0, 0), V3D(0, 0, 1)); // two theta double ph1 = det.getPhi(); // phi V3D E1 = V3D(-std::sin(tt1) * std::cos(ph1), -std::sin(tt1) * std::sin(ph1), 1. - std::cos(tt1)); // end of trajectory E1 = E1 * (1. / E1.norm()); // normalize E1Vec.push_back(E1); } }