/**
*
* @param detectorInfo A reference to the workspace's detector info
* @param detIndex The index of the detector whose parameters should be updated
* @param pmap A reference to the ParameterMap instance to update
* @param l2 The new l2 value
* @param theta The new theta value
* @param phi The new phi value
* @param delay The new delay time
* @param pressure The new pressure value
* @param thickness The new thickness value
*/
void LoadDetectorInfo::updateParameterMap(Geometry::DetectorInfo &detectorInfo,
const size_t detIndex,
Geometry::ParameterMap &pmap,
const double l2, const double theta,
const double phi, const double delay,
const double pressure,
const double thickness) const {
const auto detCompID = detectorInfo.detector(detIndex).getComponentID();
// store detector params that are different to instrument level
if (fabs(delay - m_instDelta) > 1e-06)
pmap.addDouble(detCompID, DELAY_PARAM, delay);
if (fabs(pressure - m_instPressure) > 1e-06)
pmap.addDouble(detCompID, PRESSURE_PARAM, pressure);
if (fabs(thickness - m_instThickness) > 1e-06)
pmap.addDouble(detCompID, THICKNESS_PARAM, thickness);
// move
if (m_moveDets) {
V3D newPos;
newPos.spherical(l2, theta, phi);
// The sample position may not be at 0,0,0
newPos += m_samplePos;
detectorInfo.setPosition(detIndex, newPos);
}
}
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);
}
}