linse/sics
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

We're now saving two theta, not polar angle.

Browse Source 
r2230 | ffr | 2007-11-05 14:23:46 +1100 (Mon, 05 Nov 2007) | 2 lines
This commit is contained in:
Ferdi Franceschini
2007-11-05 14:23:46 +11:00
committed by Douglas Clowes
parent a8e95faacf
commit 0b53e28ef0
2 changed files with 0 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
site_ansto/instrument/hipd/config/nexus/nxscripts.tcl
View File

@@ -1,17 +1 @@
source $cfPath(nexus)/nxscripts_common_1.tcl
# The polar angle should be calculated on the
# horizontal spatial bin boundaries.
proc put_polar_angle {nxobj point dim0 dim1} {
global det_height
set det_angle_rad [SplitReply [detector_angle_rad]]
set det_active_ht_mm [SplitReply [detector_active_height_mm]]
set det_rot_rad [ expr [SplitReply [stth]]/[SplitReply [deg_per_rad]] ]
set det_zero_col [ SplitReply [detector_zero_col]]
set det_last_hor_pixel [SplitReply [detector_last_hor_pixel]]
set col_zero [expr ($dim1 -1.0)*$det_zero_col/$det_last_hor_pixel]
set angsep [expr $det_angle_rad / ($dim1-1.0)]
$nxobj putpolararray dtheta [list $point 0] [list 1 $dim1] $angsep $det_active_ht_mm $det_rot_rad $col_zero $dim0 $dim1
}

13
site_ansto/instrument/hrpd/config/nexus/nxscripts.tcl
View File

@@ -1,14 +1 @@
source $cfPath(nexus)/nxscripts_common_1.tcl
proc put_polar_angle {nxobj point dim0 dim1} {
global det_height;
set det_angle_rad [SplitReply [detector_angle_rad]]
set det_active_ht_mm [SplitReply [detector_active_height_mm]]
set det_rot_rad [ expr [SplitReply [stth]]/[SplitReply [deg_per_rad]] ]
set det_zero_col [ SplitReply [detector_zero_col]]
set det_last_hor_pixel [SplitReply [detector_last_hor_pixel]]
set col_zero [expr ($dim1 -1.0)*$det_zero_col/$det_last_hor_pixel]
set angsep [expr $det_angle_rad / ($dim1-1)]
$nxobj putpolararray dtheta [list $point 0] [list 1 $dim1] $angsep $det_active_ht_mm $det_rot_rad $col_zero $dim0 $dim1
}