Changeset acb37d9 in sasview


Ignore:
Timestamp:
Feb 12, 2009 12:05:11 PM (15 years ago)
Author:
Jae Cho <jhjcho@…>
Branches:
master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
Children:
31e3298
Parents:
ca10d8e
Message:

get_q_compo() for qx, and qy is added.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • DataLoader/manipulations.py

    rca10d8e racb37d9  
    3131    theta      = 0.5*math.atan(plane_dist/det_dist) 
    3232    return (4.0*math.pi/wavelength)*math.sin(theta) 
     33 
     34def get_q_compo(dx, dy, det_dist, wavelength,compo=None): 
     35    #This reduces tiny error at very large q. 
     36    #Implementation of this func is not started yet.<--ToDo 
     37    if dy==0: 
     38        if dx>=0: 
     39            angle_xy=0 
     40        else: 
     41            angle_xy=math.pi 
     42    else: 
     43        angle_xy=math.atan(dx/dy) 
     44         
     45    if compo=="x": 
     46        out=get_q(dx, dy, det_dist, wavelength)*cos(angle_xy) 
     47    elif compo=="y": 
     48        out=get_q(dx, dy, det_dist, wavelength)*sin(angle_xy) 
     49    else: 
     50        out=get_q(dx, dy, det_dist, wavelength) 
     51    return out 
    3352     
    3453class _Slab(object): 
Note: See TracChangeset for help on using the changeset viewer.