Skip to content

Commit

Permalink
refs #5016. Fix input angles.
Browse files Browse the repository at this point in the history
  • Loading branch information
OwenArnold committed May 30, 2012
1 parent 15e9cc8 commit 28d55f4
Showing 1 changed file with 6 additions and 6 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -73,19 +73,19 @@ namespace Mantid
ws->splitBox();

auto spectraAxis = eventWs->getAxis(1);
const double two_pi = 6.28318531;
const double c_cos_theta_i = cos(m_incidentTheta);
const double c_sin_theta_i = sin(m_incidentTheta);

const double two_pi = 6.28318531;
const double to_radians_factor = 3.14159265/180;
const double c_cos_theta_i = cos(m_incidentTheta*to_radians_factor);
const double c_sin_theta_i = sin(m_incidentTheta*to_radians_factor);
for(size_t index = 0; index < eventWs->getNumberHistograms(); ++index)
{
auto counts = eventWs->readY(index);
auto wavelengths = eventWs->readX(index);
auto errors = eventWs->readE(index);
size_t nInputBins = eventWs->isHistogramData() ? wavelengths.size() -1 : wavelengths.size();
const double theta_final = spectraAxis->getValue(index)/2;
const double c_sin_theta_f = sin(theta_final);
const double c_cos_theta_f = cos(theta_final);
const double c_sin_theta_f = sin(theta_final*to_radians_factor);
const double c_cos_theta_f = cos(theta_final*to_radians_factor);
const double dirQx = (c_cos_theta_f - c_cos_theta_i);
const double dirQz = (c_sin_theta_f + c_sin_theta_i);
//Loop over all bins in spectra
Expand Down

0 comments on commit 28d55f4

Please sign in to comment.