Ignore:
Timestamp:
Jul 20, 2017 3:22:34 PM (5 years ago)
Author:
srkline
Message:

LOTS of changes to accommodate the beam center being reported in cm rather than pixels. Required a lot of changes to VCALC (to fill in simulated data), and in the reading and interpreting of data for display, and most importantly, the calculation of q.

There may still be a few residual bugs with this. I am still re-testing with new sample data sets.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/VC_DetectorBinning_Q.ipf

    r1038 r1055  
    177177        nPix_X = VCALC_get_nPix_X("FL") 
    178178        nPix_Y = VCALC_get_nPix_Y("FL") 
    179          
    180         xCtr = nPix_X+(F_LR_sep/2/pixSizeX)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
    181         yCtr = nPix_Y/2  
     179 
     180        if(kBCTR_CM) 
     181                xCtr = 0 
     182                yCtr = 0                        //values in cm 
     183        else     
     184                xCtr = nPix_X+(F_LR_sep/2/pixSizeX)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
     185                yCtr = nPix_Y/2  
     186        endif 
    182187        //put these  beam center values into the local folder 
    183188        V_putDet_beam_center_x("VCALC","FL",xCtr) 
     
    215220        nPix_X = VCALC_get_nPix_X("FR") 
    216221        nPix_Y = VCALC_get_nPix_Y("FR") 
    217          
    218         xCtr = -(F_LR_sep/2/pixSizeX)-1          
    219         yCtr = nPix_Y/2 
     222 
     223        if(kBCTR_CM) 
     224                xCtr = 0 
     225                yCtr = 0                        //values in cm 
     226        else     
     227                xCtr = -(F_LR_sep/2/pixSizeX)-1          
     228                yCtr = nPix_Y/2  
     229        endif 
     230 
    220231        //put these  beam center values into the local folder 
    221232        V_putDet_beam_center_x("VCALC","FR",xCtr) 
     
    248259        nPix_X = VCALC_get_nPix_X("FT") 
    249260        nPix_Y = VCALC_get_nPix_Y("FT") 
    250          
    251         xCtr = nPix_X/2 
    252         yCtr = -(F_TB_sep/2/pixSizeY)-1  
     261 
     262        if(kBCTR_CM) 
     263                xCtr = 0 
     264                yCtr = 0                        //values in cm 
     265        else     
     266                xCtr = nPix_X/2 
     267                yCtr = -(F_TB_sep/2/pixSizeY)-1  
     268        endif    
     269 
    253270                //put these  beam center values into the local folder 
    254271        V_putDet_beam_center_x("VCALC","FT",xCtr) 
     
    283300        nPix_X = VCALC_get_nPix_X("FB") 
    284301        nPix_Y = VCALC_get_nPix_Y("FB") 
    285                  
    286         xCtr = nPix_X/2 
    287         yCtr = nPix_Y+(F_TB_sep/2/pixSizeY)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
    288                 //put these  beam center values into the local folder 
     302 
     303        if(kBCTR_CM) 
     304                xCtr = 0 
     305                yCtr = 0                        //values in cm 
     306        else     
     307                xCtr = nPix_X/2 
     308                yCtr = nPix_Y+(F_TB_sep/2/pixSizeY)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
     309        endif    
     310                         
     311        //put these  beam center values into the local folder 
    289312        V_putDet_beam_center_x("VCALC","FB",xCtr) 
    290313        V_putDet_beam_center_y("VCALC","FB",yCtr) 
     
    596619        nPix_X = VCALC_get_nPix_X("ML") 
    597620        nPix_Y = VCALC_get_nPix_Y("ML") 
    598          
    599         xCtr = nPix_X+(M_LR_sep/2/pixSizeX)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
    600         yCtr = nPix_Y/2  
     621 
     622        if(kBCTR_CM) 
     623                xCtr = 0 
     624                yCtr = 0                        //values in cm 
     625        else     
     626                xCtr = nPix_X+(M_LR_sep/2/pixSizeX)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
     627                yCtr = nPix_Y/2  
     628        endif            
     629 
    601630                //put these  beam center values into the local folder 
    602631        V_putDet_beam_center_x("VCALC","ML",xCtr) 
     
    637666        nPix_X = VCALC_get_nPix_X("MR") 
    638667        nPix_Y = VCALC_get_nPix_Y("MR") 
    639          
    640         xCtr = -(M_LR_sep/2/pixSizeX)-1          
    641         yCtr = nPix_Y/2 
     668 
     669        if(kBCTR_CM) 
     670                xCtr = 0 
     671                yCtr = 0                        //values in cm 
     672        else     
     673                xCtr = -(M_LR_sep/2/pixSizeX)-1          
     674                yCtr = nPix_Y/2  
     675        endif    
     676                 
     677 
    642678                //put these  beam center values into the local folder 
    643679        V_putDet_beam_center_x("VCALC","MR",xCtr) 
     
    671707        nPix_X = VCALC_get_nPix_X("MT") 
    672708        nPix_Y = VCALC_get_nPix_Y("MT") 
    673          
    674         xCtr = nPix_X/2 
    675         yCtr = -(M_TB_sep/2/pixSizeY)-1  
     709 
     710        if(kBCTR_CM) 
     711                xCtr = 0 
     712                yCtr = 0                        //values in cm 
     713        else     
     714                xCtr = nPix_X/2 
     715                yCtr = -(M_TB_sep/2/pixSizeY)-1  
     716        endif 
     717                 
     718 
    676719                //put these  beam center values into the local folder 
    677720        V_putDet_beam_center_x("VCALC","MT",xCtr) 
     
    706749        nPix_X = VCALC_get_nPix_X("MB") 
    707750        nPix_Y = VCALC_get_nPix_Y("MB") 
    708                  
    709         xCtr = nPix_X/2 
    710         yCtr = nPix_Y+(M_TB_sep/2/pixSizeY)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)?? 
     751 
     752        if(kBCTR_CM) 
     753                xCtr = 0 
     754                yCtr = 0                        //values in cm 
     755        else     
     756                xCtr = nPix_X/2 
     757                yCtr = nPix_Y+(M_TB_sep/2/pixSizeY)             // TODO  -- check -- starting from 47 rather than 48 (but I'm in pixel units for centers)??  
     758        endif            
     759 
    711760        //put these  beam center values into the local folder 
    712761        V_putDet_beam_center_x("VCALC","MB",xCtr) 
Note: See TracChangeset for help on using the changeset viewer.