Ignore:
Timestamp:
Apr 24, 2018 3:55:39 PM (5 years ago)
Author:
srkline
Message:

added procedures to compare files to see if they are from the same configuration, same wavelength, etc. so they can be properly chosen for transmission files, scattering files, and properly identified for the different resolution conditions.

Re-worked the logic of dispatching averaging, plotting, and saving in Execute_Protocol. Hopefully this will alow for easier dispatching for future conditions, including getting the correct resolution calculation.

File:
1 edited

Legend:

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

    r1095 r1097  
    785785// 
    786786// folderStr = WORK folder, type = the binning type (may include multiple detectors) 
    787 Function VC_fDoBinning_QxQy2D(folderStr,type) 
    788         String folderStr,type 
     787Function VC_fDoBinning_QxQy2D(folderStr,type,collimationStr) 
     788        String folderStr,type,collimationStr 
    789789         
    790790        Variable nSets = 0 
     
    14801480// 
    14811481 
    1482  
    1483         ii=0 
    1484         do 
    1485                 V_getResolution(qBin_qxqy[ii],lambda,lambdaWidth,DDet,apOff,S1,S2,L1,L2,BS,del_r,usingLenses,ret1,ret2,ret3) 
    1486                 sigmaq[ii] = ret1        
    1487                 qbar[ii] = ret2  
    1488                 fsubs[ii] = ret3         
    1489                 ii+=1 
    1490         while(ii<nq) 
    1491          
    1492          
    1493          
     1482// possible values are: 
     1483// 
     1484// pinhole 
     1485// pinhole_whiteBeam 
     1486// narrowSlit 
     1487// narrowSlit_whiteBeam 
     1488// convergingPinholes 
     1489 
     1490        if(cmpstr(collimationStr,"pinhole") == 0) 
     1491 
     1492                ii=0 
     1493                do 
     1494                        V_getResolution(qBin_qxqy[ii],lambda,lambdaWidth,DDet,apOff,S1,S2,L1,L2,BS,del_r,usingLenses,ret1,ret2,ret3) 
     1495                        sigmaq[ii] = ret1        
     1496                        qbar[ii] = ret2  
     1497                        fsubs[ii] = ret3         
     1498                        ii+=1 
     1499                while(ii<nq) 
     1500         
     1501        endif 
     1502         
     1503 
     1504        if(cmpstr(collimationStr,"pinhole_whiteBeam") == 0) 
     1505 
     1506//              set lambdaWidth == 0 so that the gaussian resolution calculates only the geometry contribution. 
     1507// the white beam distribution will need to be flagged some other way 
     1508// 
     1509                lambdaWidth = 0 
     1510                 
     1511                ii=0 
     1512                do 
     1513                        V_getResolution(qBin_qxqy[ii],lambda,lambdaWidth,DDet,apOff,S1,S2,L1,L2,BS,del_r,usingLenses,ret1,ret2,ret3) 
     1514                        sigmaq[ii] = ret1        
     1515                        qbar[ii] = ret2  
     1516                        fsubs[ii] = ret3         
     1517                        ii+=1 
     1518                while(ii<nq) 
     1519         
     1520        endif 
     1521 
     1522        if(cmpstr(collimationStr,"convergingPinholes") == 0) 
     1523 
     1524//              set usingLenses == 1 so that the Gaussian resolution calculation will be for a focus condition 
     1525// 
     1526                usingLenses = 1 
     1527                 
     1528                ii=0 
     1529                do 
     1530                        V_getResolution(qBin_qxqy[ii],lambda,lambdaWidth,DDet,apOff,S1,S2,L1,L2,BS,del_r,usingLenses,ret1,ret2,ret3) 
     1531                        sigmaq[ii] = ret1        
     1532                        qbar[ii] = ret2  
     1533                        fsubs[ii] = ret3         
     1534                        ii+=1 
     1535                while(ii<nq) 
     1536         
     1537        endif 
     1538                 
    14941539        SetDataFolder root: 
    14951540         
Note: See TracChangeset for help on using the changeset viewer.