2,718
edits
Get rot from xparm: Difference between revisions
no edit summary
No edit summary |
No edit summary |
||
| Line 31: | Line 31: | ||
xyz=xyz/x | xyz=xyz/x | ||
x=dot_product(xyz,rotation) | x=dot_product(xyz,rotation) | ||
print*,axnam(i),' axis | print*,'angle between ',axnam(i),' axis and ROTATION_AXIS: ',acos(x)*r2d | ||
x=dot_product(xyz,detector_x) | x=dot_product(xyz,detector_x) | ||
print*,axnam(i),' axis | print*,'angle between ',axnam(i),' axis and DETECTOR_X-AXIS: ',acos(x)*r2d | ||
x=dot_product(xyz,detector_y) | x=dot_product(xyz,detector_y) | ||
print*,axnam(i),' axis | print*,'angle between ',axnam(i),' axis and DETECTOR_Y-AXIS: ',acos(x)*r2d | ||
x=dot_product(xyz,beam) | x=dot_product(xyz,beam) | ||
print*,axnam(i),' axis | print*,'angle between ',axnam(i),' axis and INCIDENT_BEAM DIRECTION:',acos(x)*r2d | ||
end do | end do | ||
end | end | ||