/**
* Adds the zero-offset of the chopper to the slit times
*
* @param chopper :: PoldiAbstractChopper with slit times, not corrected with zero-offset
* @return vector with zero-offset-corrected chopper slit times
*/
std::vector<double> PoldiSpectrumDomainFunction::getChopperSlitOffsets(const PoldiAbstractChopper_sptr &chopper)
{
const std::vector<double> &chopperSlitTimes = chopper->slitTimes();
std::vector<double> offsets;
offsets.reserve(chopperSlitTimes.size());
for(std::vector<double>::const_iterator time = chopperSlitTimes.begin(); time != chopperSlitTimes.end(); ++time) {
offsets.push_back(*time + chopper->zeroOffset());
}
return offsets;
}
/**
* Adds the zero-offset of the chopper to the slit times
*
* @param chopper :: PoldiAbstractChopper with slit times, not corrected with
* zero-offset
* @return vector with zero-offset-corrected chopper slit times
*/
std::vector<double> PoldiSpectrumDomainFunction::getChopperSlitOffsets(
const PoldiAbstractChopper_sptr &chopper) {
const std::vector<double> &chopperSlitTimes = chopper->slitTimes();
std::vector<double> offsets;
offsets.reserve(chopperSlitTimes.size());
for (double chopperSlitTime : chopperSlitTimes) {
offsets.push_back(chopperSlitTime + chopper->zeroOffset());
}
return offsets;
}
void PoldiAutoCorrelation5::logConfigurationInformation(
boost::shared_ptr<PoldiDeadWireDecorator> cleanDetector,
PoldiAbstractChopper_sptr chopper) {
if (cleanDetector && chopper) {
g_log.information()
<< "____________________________________________________ " << std::endl;
g_log.information()
<< "_Poldi chopper conf ------------------------------ " << std::endl;
g_log.information() << "_Poldi - Chopper speed: "
<< chopper->rotationSpeed() << " rpm" << std::endl;
g_log.information() << "_Poldi - Number of slits: "
<< chopper->slitPositions().size() << std::endl;
g_log.information() << "_Poldi - Cycle time: "
<< chopper->cycleTime() << " µs" << std::endl;
g_log.information() << "_Poldi - Zero offset: "
<< chopper->zeroOffset() << " µs" << std::endl;
g_log.information() << "_Poldi - Distance: "
<< chopper->distanceFromSample() << " mm" << std::endl;
if (g_log.is(Poco::Message::PRIO_DEBUG)) {
for (size_t i = 0; i < chopper->slitPositions().size(); ++i) {
g_log.information() << "_Poldi - Slits: " << i
<< ": Position = " << chopper->slitPositions()[i]
<< "\t Time = " << chopper->slitTimes()[i] << " µs"
<< std::endl;
}
}
g_log.information()
<< "_Poldi detector conf ------------------------------ "
<< std::endl;
g_log.information() << "_Poldi - Element count: "
<< cleanDetector->elementCount() << std::endl;
g_log.information() << "_Poldi - Central element: "
<< cleanDetector->centralElement() << std::endl;
g_log.information() << "_Poldi - 2Theta(central): "
<< cleanDetector->twoTheta(199) / M_PI * 180.0 << "°"
<< std::endl;
g_log.information() << "_Poldi - Distance(central): "
<< cleanDetector->distanceFromSample(199) << " mm"
<< std::endl;
std::set<int> deadWires = cleanDetector->deadWires();
g_log.information() << "_Poldi - Number of dead wires: "
<< deadWires.size() << std::endl;
g_log.information() << "_Poldi - Wire indices: ";
for (std::set<int>::const_iterator dw = deadWires.begin();
dw != deadWires.end(); ++dw) {
g_log.information() << *dw << " ";
}
g_log.information() << std::endl;
}
}