Ignore:
Timestamp:
Mar 11, 2016 3:36:43 PM (7 years ago)
Author:
srkline
Message:

lots of changes to 1D averaging and the plotting routines, detector corrections, and basic reads

File:
1 edited

Legend:

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

    r983 r984  
    4040// 
    4141// TODO: -- is there an extra "entry" heading? Am I adding this by mistake by setting base_name="entry" for RAW data? 
    42 //                      -- as dumb as it is -- do I just leave it now, or break everything. ont the plus side, removing the extra entry 
     42//                      -- as dumb as it is -- do I just leave it now, or break everything. On the plus side, removing the extra "entry" 
    4343//          layer may catch a lot of the hard-wired junk that is present... 
    4444Function V_LoadHDF5Data(file,folder) 
     
    109109                                NonLinearCorrection(w,w_calib,tube_width,detStr,destPath) 
    110110                                 
     111                                 
     112                                //(2.4) Convert the beam center values from pixels to mm 
     113                                // TODO -- there needs to be a permanent location for these values?? 
     114                                // 
     115                                ConvertBeamCtr_to_mm(folder,detStr,destPath) 
     116                                 
    111117                                // (2.5) Calculate the q-values 
    112118                                // calculating q-values can't be done unless the non-linear corrections are calculated 
     
    122128                                 
    123129                        endfor 
     130                         
     131                        //"B" is separate 
     132                        NonLinearCorrection_B(folder,"B",destPath) 
     133                        ConvertBeamCtr_to_mmB(folder,"B",destPath) 
     134                        V_Detector_CalcQVals(folder,"B",destPath) 
     135                         
    124136                else 
    125137                        Print "Non-linear correction not done" 
     
    157169        det_B=tmpw 
    158170        det_B += 2 
     171        Wave distance=distance 
     172        distance = 2200 
     173 
    159174                         
    160175        Variable ctr=20,npix=128 
     
    170185        det_MT=tmpw 
    171186        det_MT += 2 
     187        Wave distance=distance 
     188        distance = 1030 
     189 
    172190         
    173191        SetDataFolder root:Packages:NIST:VSANS:RAW:entry:entry:instrument:detector_MB 
     
    182200        det_MB=tmpw 
    183201        det_MB += 2 
     202        Wave distance=distance 
     203        distance = 1030 
     204 
    184205         
    185206        ctr=30 
     
    195216        det_ML=tmpw 
    196217        det_ML += 2 
     218        Wave distance=distance 
     219        distance = 1000 
     220 
    197221                 
    198222        SetDataFolder root:Packages:NIST:VSANS:RAW:entry:entry:instrument:detector_MR 
     
    206230        det_MR=tmpw 
    207231        det_MR += 2 
     232        Wave distance=distance 
     233        distance = 1000 
     234         
    208235         
    209236        ctr=30 
     
    216243        wave tmpw=$"root:Packages:NIST:VSANS:VCALC:entry:entry:instrument:detector_FT:det_FT" 
    217244        det_FT=tmpw 
     245        Wave distance=distance 
     246        distance = 180 
     247 
    218248 
    219249        SetDataFolder root:Packages:NIST:VSANS:RAW:entry:entry:instrument:detector_FB 
     
    225255        wave tmpw=$"root:Packages:NIST:VSANS:VCALC:entry:entry:instrument:detector_FB:det_FB" 
    226256        det_FB=tmpw 
     257        Wave distance=distance 
     258        distance = 180 
     259 
    227260                         
    228261        SetDataFolder root:Packages:NIST:VSANS:RAW:entry:entry:instrument:detector_FL 
     
    234267        wave tmpw=$"root:Packages:NIST:VSANS:VCALC:entry:entry:instrument:detector_FL:det_FL" 
    235268        det_FL=tmpw 
     269        Wave distance=distance 
     270        distance = 150 
     271 
    236272         
    237273        SetDataFolder root:Packages:NIST:VSANS:RAW:entry:entry:instrument:detector_FR 
     
    243279        wave tmpw=$"root:Packages:NIST:VSANS:VCALC:entry:entry:instrument:detector_FR:det_FR" 
    244280        det_FR=tmpw 
     281        Wave distance=distance 
     282        distance = 150 
     283 
    245284         
    246285// get rid of zeros 
     
    250289        det_FB += 2 
    251290 
    252 V_RescaleToBeamCenter("RAW","MB",64,55) 
    253 V_RescaleToBeamCenter("RAW","MT",64,-8.7) 
    254 V_RescaleToBeamCenter("RAW","MR",-8.1,64) 
    255 V_RescaleToBeamCenter("RAW","ML",55,64) 
    256 V_RescaleToBeamCenter("RAW","FL",55,64) 
    257 V_RescaleToBeamCenter("RAW","FR",-8.1,64) 
    258 V_RescaleToBeamCenter("RAW","FT",64,-8.7) 
    259 V_RescaleToBeamCenter("RAW","FB",64,55) 
     291 
     292// fake beam center values 
     293        V_putDet_beam_center_x("RAW","B",160) 
     294        V_putDet_beam_center_y("RAW","B",160) 
     295 
     296 
     297        V_putDet_beam_center_x("RAW","MB",64) 
     298        V_putDet_beam_center_y("RAW","MB",55) 
     299        V_putDet_beam_center_x("RAW","MT",64) 
     300        V_putDet_beam_center_y("RAW","MT",-8.7) 
     301        V_putDet_beam_center_x("RAW","MR",-8.1) 
     302        V_putDet_beam_center_y("RAW","MR",64) 
     303        V_putDet_beam_center_x("RAW","ML",55) 
     304        V_putDet_beam_center_y("RAW","ML",64) 
     305 
     306        V_putDet_beam_center_x("RAW","FB",64) 
     307        V_putDet_beam_center_y("RAW","FB",55) 
     308        V_putDet_beam_center_x("RAW","FT",64) 
     309        V_putDet_beam_center_y("RAW","FT",-8.7) 
     310        V_putDet_beam_center_x("RAW","FR",-8.1) 
     311        V_putDet_beam_center_y("RAW","FR",64) 
     312        V_putDet_beam_center_x("RAW","FL",55) 
     313        V_putDet_beam_center_y("RAW","FL",64) 
     314 
     315 
     316 
     317 
     318 
     319 
     320        V_RescaleToBeamCenter("RAW","MB",64,55) 
     321        V_RescaleToBeamCenter("RAW","MT",64,-8.7) 
     322        V_RescaleToBeamCenter("RAW","MR",-8.1,64) 
     323        V_RescaleToBeamCenter("RAW","ML",55,64) 
     324        V_RescaleToBeamCenter("RAW","FL",55,64) 
     325        V_RescaleToBeamCenter("RAW","FR",-8.1,64) 
     326        V_RescaleToBeamCenter("RAW","FT",64,-8.7) 
     327        V_RescaleToBeamCenter("RAW","FB",64,55) 
    260328 
    261329 
     
    269337// this will load in the whole HDF file all at once. 
    270338// Attributes are NOT loaded at all. 
     339// 
     340// -- the Gateway function H5GW_ReadHDF5(parentFolder, fileName, [hdf5Path]) 
     341//    reads in the attributes too, but is very slow 
     342//   -- the H5GW function is called by: H_HDF5Gate_Read_Raw(file) 
    271343// 
    272344// TODO:  
     
    275347// -- be sure I'm using all of the correct flags in the HDF5LoadGroup operation 
    276348// -- settle on how the base_name is to be used. "entry" for the RAW, fileName for the "rawVSANS"? 
     349// x- error check for path existence 
    277350// 
    278351// passing in "" for base_name will take the name from the file name as selected 
     
    284357//              hdf5Path = "/" 
    285358//      endif 
     359 
     360        PathInfo/S catPathName 
     361        if(V_flag == 0) 
     362                V_PickPath() 
     363        endif 
    286364 
    287365        String hdf5path = "/"           //always read from the top 
Note: See TracChangeset for help on using the changeset viewer.