Changeset 947 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS
- Timestamp:
- Oct 16, 2014 10:59:26 AM (8 years ago)
- Location:
- sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/VCALC_PanelMockUp.ipf
r936 r947 567 567 NewDataFolder/O root:Packages:NIST:VSANS 568 568 NewDataFolder/O root:Packages:NIST:VSANS:VCALC 569 NewDataFolder/O root:Packages:NIST:VSANS:VCALC:Front 570 NewDataFolder/O root:Packages:NIST:VSANS:VCALC:Middle 571 NewDataFolder/O root:Packages:NIST:VSANS:VCALC:Back 569 572 570 573 … … 581 584 ///// FRONT DETECTOR BANKS 582 585 // dimensions for the detector banks (then get them in the drawing functions) 583 Variable/G front_LR_w = 384 //front bank, nominal LR panel width (mm)584 Variable/G front_LR_h = 1000585 Variable/G front_TB_w = 500586 Variable/G front_TB_h = 384586 Variable/G gFront_LR_w = 384 //front bank, nominal LR panel width (mm) 587 Variable/G gFront_LR_h = 1000 588 Variable/G gFront_TB_w = 500 589 Variable/G gFront_TB_h = 384 587 590 588 591 // SDD offset of T/B (decide on units??) 589 Variable/G front_SDDOffset = 300 // (mm)592 Variable/G gFront_SDDOffset = 300 // (mm) 590 593 591 594 // detector resolution (xy for each bank!) 592 Variable/G front_L_pixelX = 0.8 // (cm) these tubes are vertical593 Variable/G front_L_pixelY = 0.4 // (cm)594 Variable/G front_R_pixelX = 0.8 // (cm)595 Variable/G front_R_pixelY = 0.4 // (cm)596 597 Variable/G front_T_pixelX = 0.4 // (cm) these tubes are horizontal598 Variable/G front_T_pixelY = 0.8 // (cm)599 Variable/G front_B_pixelX = 0.4 // (cm)600 Variable/G front_B_pixelY = 0.8 // (cm)595 Variable/G gFront_L_pixelX = 0.8 // (cm) these tubes are vertical 596 Variable/G gFront_L_pixelY = 0.4 // (cm) 597 Variable/G gFront_R_pixelX = 0.8 // (cm) 598 Variable/G gFront_R_pixelY = 0.4 // (cm) 599 600 Variable/G gFront_T_pixelX = 0.4 // (cm) these tubes are horizontal 601 Variable/G gFront_T_pixelY = 0.8 // (cm) 602 Variable/G gFront_B_pixelX = 0.4 // (cm) 603 Variable/G gFront_B_pixelY = 0.8 // (cm) 601 604 602 605 // number of pixels in each bank (this can be modified at acquisition time, so it must be adjustable here) … … 605 608 606 609 ///// MIDDLE DETECTOR BANKS 607 Variable/G middle_LR_w = 384 //middle bank, nominal LR panel width (mm)608 Variable/G middle_LR_h = 1000609 Variable/G middle_TB_w = 500610 Variable/G middle_TB_h = 384610 Variable/G gMiddle_LR_w = 384 //middle bank, nominal LR panel width (mm) 611 Variable/G gMiddle_LR_h = 1000 612 Variable/G gMiddle_TB_w = 500 613 Variable/G gMiddle_TB_h = 384 611 614 // SDD offset of T/B (decide on units??) 612 Variable/G middle_SDDOffset = 300 // (mm)615 Variable/G gMiddle_SDDOffset = 300 // (mm) 613 616 614 617 // detector resolution (xy for each bank!) 615 Variable/G middle_L_pixelX = 0.8 // (cm) these tubes are vertical616 Variable/G middle_L_pixelY = 0.4 // (cm)617 Variable/G middle_R_pixelX = 0.8 // (cm)618 Variable/G middle_R_pixelY = 0.4 // (cm)619 620 Variable/G middle_T_pixelX = 0.4 // (cm) these tubes are horizontal621 Variable/G middle_T_pixelY = 0.8 // (cm)622 Variable/G middle_B_pixelX = 0.4 // (cm)623 Variable/G middle_B_pixelY = 0.8 // (cm)618 Variable/G gMiddle_L_pixelX = 0.8 // (cm) these tubes are vertical 619 Variable/G gMiddle_L_pixelY = 0.4 // (cm) 620 Variable/G gMiddle_R_pixelX = 0.8 // (cm) 621 Variable/G gMiddle_R_pixelY = 0.4 // (cm) 622 623 Variable/G gMiddle_T_pixelX = 0.4 // (cm) these tubes are horizontal 624 Variable/G gMiddle_T_pixelY = 0.8 // (cm) 625 Variable/G gMiddle_B_pixelX = 0.4 // (cm) 626 Variable/G gMiddle_B_pixelY = 0.8 // (cm) 624 627 625 628 626 629 627 630 //// BACK DETECTOR 628 Variable/G back_w = 320 //w and h for the back detector, (mm)629 Variable/G back_h = 320630 Variable/G back_pixelX = 0.1 // 1mm resolution (units of cm here)631 Variable/G back_pixelY = 0.1631 Variable/G gBack_w = 320 //w and h for the back detector, (mm) 632 Variable/G gBack_h = 320 633 Variable/G gBack_pixelX = 0.1 // 1mm resolution (units of cm here) 634 Variable/G gBack_pixelY = 0.1 632 635 633 636 634 637 // Generate all of the waves used for the detector and the q values 635 638 // 636 // TODO: the d imensions need to be properly defined here...639 // TODO: the detector dimensions need to be properly defined here... 637 640 // FRONT 641 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 642 638 643 Make/O/D/N=(48,256) det_FL,det_FR 639 644 Make/O/D/N=(128,48) det_FT,det_FB … … 644 649 645 650 //MIDDLE 646 // TODO: the dimensions need to be properly defined here... 651 // TODO: the detector dimensions need to be properly defined here... 652 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 653 647 654 Make/O/D/N=(48,256) det_ML,det_MR 648 655 Make/O/D/N=(128,48) det_MT,det_MB … … 653 660 654 661 // BACK 662 // TODO: the detector dimensions need to be properly defined here... 663 SetDataFolder root:Packages:NIST:VSANS:VCALC:Back 664 655 665 Make/O/D/N=(320,320) det_B 656 666 Duplicate/O det_B qTot_B,qx_B,qy_B,qz_B … … 659 669 //////////// FOR THE PANEL 660 670 671 SetDataFolder root:Packages:NIST:VSANS:VCALC 661 672 662 673 // popup strings for each tab (then use the string in the panel) … … 689 700 // limits for detector travel? or are these limits part of the panel, hard-wired there 690 701 691 692 693 694 695 696 702 697 703 // … … 775 781 776 782 // 777 Variable/G instrument = 6 // files (may) be tagged SA6 as the 6th SANS instrument778 Variable/G s12 = 54.8783 Variable/G gInstrument = 6 // files (may) be tagged SA6 as the 6th SANS instrument 784 Variable/G gS12 = 54.8 779 785 // Variable/G d_det = 0.5 780 786 // Variable/G a_pixel = 0.5 781 787 // Variable/G del_r = 0.5 782 788 // Variable/G det_width = 64.0 783 Variable/G lambda_t = 5.50784 Variable/G l2r_lower = 132.3785 Variable/G l2r_upper = 1317786 Variable/G lambda_lower = 2.5787 Variable/G lambda_upper = 20.0788 Variable/G d_upper = 25.0789 Variable/G bs_factor = 1.05790 Variable/G t1 = 0.63791 Variable/G t2 = 1.0792 Variable/G t3 = 0.75793 Variable/G l_gap = 100.0794 Variable/G g uide_width = 6.0795 Variable/G idmax = 100.0789 Variable/G gLambda_t = 5.50 790 Variable/G gL2r_lower = 132.3 791 Variable/G gL2r_upper = 1317 792 Variable/G gLambda_lower = 2.5 793 Variable/G gLambda_upper = 20.0 794 Variable/G gD_upper = 25.0 795 Variable/G gBs_factor = 1.05 796 Variable/G gT1 = 0.63 797 Variable/G gT2 = 1.0 798 Variable/G gT3 = 0.75 799 Variable/G gL_gap = 100.0 800 Variable/G gGuide_width = 6.0 801 Variable/G gIdmax = 100.0 796 802 797 803 // 798 804 // //new values, from 11/2009 --- BeamFluxReport_2009.ifn 799 Variable/G phi_0 = 2.42e13800 Variable/G b= 0.0801 Variable/G c= -0.0243805 Variable/G gPhi_0 = 2.42e13 806 Variable/G gB = 0.0 807 Variable/G gC = -0.0243 802 808 Variable/G gGuide_loss = 0.924 803 809 // 804 810 // //fwhm values (new variables) (+3, 0, -3, calibrated 2009) 805 Variable/G fwhm_narrow = 0.109806 Variable/G fwhm_mid = 0.125807 Variable/G fwhm_wide = 0.236811 Variable/G gFwhm_narrow = 0.109 812 Variable/G gFwhm_mid = 0.125 813 Variable/G gFwhm_wide = 0.236 808 814 // 809 815 // //source apertures (cm) 810 Variable/G a1_0_0 = 1.43811 Variable/G a1_0_1 = 2.54812 Variable/G a1_0_2 = 3.81813 Variable/G a1_7_0 = 2.5 // after the polarizer814 Variable/G a1_7_1 = 5.0815 Variable/G a1_7_1 = 0.95 //816 Variable/G a1_def = 5.00816 Variable/G gA1_0_0 = 1.43 817 Variable/G gA1_0_1 = 2.54 818 Variable/G gA1_0_2 = 3.81 819 Variable/G gA1_7_0 = 2.5 // after the polarizer 820 Variable/G gA1_7_1 = 5.0 821 Variable/G gA1_7_1 = 0.95 // 822 Variable/G gA1_def = 5.00 817 823 // 818 824 SetDataFolder root: -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_BinSlitAperture.ipf
r935 r947 12 12 // 13 13 // 14 Function V_fBinDetector_byRows( type)15 String type14 Function V_fBinDetector_byRows(folderStr,type) 15 String folderStr,type 16 16 17 17 SetDataFolder root:Packages:NIST:VSANS:VCALC … … 19 19 Variable pixSizeX,pixSizeY,delQx, delQy 20 20 21 WAVE inten = $(" det_"+type) // 2D detector data21 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:"+folderStr+":det_"+type) // 2D detector data 22 22 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 23 Wave qTotal = $(" qTot_"+type) // 2D q-values24 Wave qx = $(" qx_"+type)25 Wave qy = $(" qy_"+type)23 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:"+folderStr+":qTot_"+type) // 2D q-values 24 Wave qx = $("root:Packages:NIST:VSANS:VCALC:"+folderStr+":qx_"+type) 25 Wave qy = $("root:Packages:NIST:VSANS:VCALC:"+folderStr+":qy_"+type) 26 26 27 27 pixSizeX = VCALC_getPixSizeX(type) -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_DetectorBin_Middle.ipf
r936 r947 19 19 // fill the panels with fake sphere scattering data 20 20 // TODO: am I in the right data folder?? 21 SetDataFolder root:Packages:NIST:VSANS:VCALC 21 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 22 22 23 23 WAVE det_ML = det_ML … … 41 41 // view of how much of the detectors are actually collecting data 42 42 // -- I can get the separation L/R from the panel - only this "open" width is visible. 43 V_SetShadow_TopBottom(" ","MT") // TODO: -- be sure the data folder is properly set (within the function...)44 V_SetShadow_TopBottom(" ","MB")43 V_SetShadow_TopBottom("Middle","MT") // TODO: -- be sure the data folder is properly set (within the function...) 44 V_SetShadow_TopBottom("Middle","MB") 45 45 46 46 // do the q-binning for each of the panels to get I(Q) … … 87 87 M_sdd_offset = VSANS_getTopBottomSDDOffset("MT") //T/B are 30 cm farther back //TODO: make all detector parameters global, not hard-wired 88 88 89 SetDataFolder root:Packages:NIST:VSANS:VCALC 89 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 90 90 Wave det_ML,det_MR // these are (48,256) 91 91 Wave det_MT,det_MB // these are (128,48) … … 118 118 119 119 //Middle/RIGHT 120 SetDataFolder root:Packages:NIST:VSANS:VCALC 120 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 121 121 WAVE qTot_MR,qx_MR,qy_MR,qz_MR 122 122 qTot_MR = 0 … … 139 139 140 140 //Middle/TOP 141 SetDataFolder root:Packages:NIST:VSANS:VCALC 141 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 142 142 WAVE qTot_MT,qx_MT,qy_MT,qz_MT 143 143 qTot_MT = 0 … … 161 161 162 162 //Middle/BOTTOM 163 SetDataFolder root:Packages:NIST:VSANS:VCALC 163 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 164 164 WAVE qTot_MB,qx_MB,qy_MB,qz_MB 165 165 qTot_MB = 0 … … 194 194 // Display /W=(1477,44,1978,517) 195 195 196 SetDataFolder root:Packages:NIST:VSANS:VCALC 196 SetDataFolder root:Packages:NIST:VSANS:VCALC:Middle 197 197 198 198 CheckDisplayed /W=VCALC#Panels_Q det_MB … … 249 249 Proc BinAllMiddlePanels() 250 250 251 SetDeltaQ(" ","ML")252 SetDeltaQ(" ","MT")251 SetDeltaQ("Middle","ML") 252 SetDeltaQ("Middle","MT") 253 253 254 254 Variable binType … … 275 275 if(binType == 4) 276 276 /// this is for a tall, narrow slit mode 277 V_fBinDetector_byRows("M L")278 V_fBinDetector_byRows("M R")279 V_fBinDetector_byRows("M T")280 V_fBinDetector_byRows("M B")277 V_fBinDetector_byRows("Middle","ML") 278 V_fBinDetector_byRows("Middle","MR") 279 V_fBinDetector_byRows("Middle","MT") 280 V_fBinDetector_byRows("Middle","MB") 281 281 endif 282 282 End … … 415 415 // fill the panels with fake sphere scattering data 416 416 // TODO: am I in the right data folder?? 417 SetDataFolder root:Packages:NIST:VSANS:VCALC 417 SetDataFolder root:Packages:NIST:VSANS:VCALC:Back 418 418 419 419 WAVE det_B = det_B … … 463 463 464 464 // TODO (make the N along the tube length a variable, since this can be reset @ acquisition) 465 SetDataFolder root:Packages:NIST:VSANS:VCALC 465 SetDataFolder root:Packages:NIST:VSANS:VCALC:Back 466 466 Wave det_B // this is (320,320) 467 467 … … 501 501 // Display /W=(1477,44,1978,517) 502 502 503 SetDataFolder root:Packages:NIST:VSANS:VCALC 503 SetDataFolder root:Packages:NIST:VSANS:VCALC:Back 504 504 505 505 CheckDisplayed /W=VCALC#Panels_Q det_B … … 546 546 Proc BinAllBackPanels() 547 547 548 SetDeltaQ(" ","B")548 SetDeltaQ("Back","B") 549 549 550 550 Variable binType … … 557 557 if(binType == 4) 558 558 /// this is for a tall, narrow slit mode 559 V_fBinDetector_byRows("B ")559 V_fBinDetector_byRows("Back","B") 560 560 endif 561 561 -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_DetectorBinning.ipf
r936 r947 39 39 // fill the panels with fake sphere scattering data 40 40 // TODO: am I in the right data folder?? 41 SetDataFolder root:Packages:NIST:VSANS:VCALC 41 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 42 42 43 43 WAVE det_FL = det_FL … … 66 66 // -- I can get the separation L/R from the panel - only this "open" width is visible. 67 67 //TODO - make this a proper shadow - TB extent of the LR panels matters too, not just the LR separation 68 V_SetShadow_TopBottom(" ","FT") // TODO: -- be sure the data folder is properly set (within the function...)69 V_SetShadow_TopBottom(" ","FB")68 V_SetShadow_TopBottom("Front","FT") // TODO: -- be sure the data folder is properly set (within the function...) 69 V_SetShadow_TopBottom("Front","FB") 70 70 71 71 // do the q-binning for each of the panels to get I(Q) … … 83 83 String type 84 84 85 SetDataFolder root:Packages:NIST:VSANS:VCALC 85 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 86 86 87 87 // q-values and detector arrays already allocated and calculated … … 222 222 F_sdd_offset = VSANS_getTopBottomSDDOffset("FT") //T/B are 300 mm farther back //TODO: make all detector parameters global, not hard-wired 223 223 224 SetDataFolder root:Packages:NIST:VSANS:VCALC 224 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 225 225 Wave det_FL,det_FR // these are (48,256) 226 226 Wave det_FT,det_FB // these are (128,48) … … 254 254 255 255 //FRONT/RIGHT 256 SetDataFolder root:Packages:NIST:VSANS:VCALC 256 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 257 257 WAVE qTot_FR,qx_FR,qy_FR,qz_FR 258 258 qTot_FR = 0 … … 275 275 276 276 //FRONT/TOP 277 SetDataFolder root:Packages:NIST:VSANS:VCALC 277 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 278 278 WAVE qTot_FT,qx_FT,qy_FT,qz_FT 279 279 qTot_FT = 0 … … 297 297 298 298 //FRONT/BOTTOM 299 SetDataFolder root:Packages:NIST:VSANS:VCALC 299 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 300 300 WAVE qTot_FB,qx_FB,qy_FB,qz_FB 301 301 qTot_FB = 0 … … 347 347 348 348 //detector data 349 Wave det = $("root:Packages:NIST:VSANS:VCALC:"+ "det_"+type)349 Wave det = $("root:Packages:NIST:VSANS:VCALC:"+folderStr+":det_"+type) 350 350 351 351 // TODO - these are to be set from globals, not hard-wired … … 378 378 // Display /W=(1477,44,1978,517) 379 379 380 SetDataFolder root:Packages:NIST:VSANS:VCALC 380 SetDataFolder root:Packages:NIST:VSANS:VCALC:Front 381 381 382 382 CheckDisplayed /W=VCALC#Panels_Q det_FB … … 702 702 Proc BinAllFrontPanels() 703 703 704 SetDeltaQ(" ","FL")705 SetDeltaQ(" ","FT")704 SetDeltaQ("Front","FL") 705 SetDeltaQ("Front","FT") 706 706 707 707 Variable binType … … 728 728 if(binType == 4) 729 729 /// this is for a tall, narrow slit mode 730 V_fBinDetector_byRows("F L")731 V_fBinDetector_byRows("F R")732 V_fBinDetector_byRows("F T")733 V_fBinDetector_byRows("F B")730 V_fBinDetector_byRows("Front","FL") 731 V_fBinDetector_byRows("Front","FR") 732 V_fBinDetector_byRows("Front","FT") 733 V_fBinDetector_byRows("Front","FB") 734 734 endif 735 735 … … 737 737 738 738 739 //TODO -- folderStr is ignored in this function740 739 Function SetDeltaQ(folderStr,type) 741 740 String folderStr,type 742 743 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "det_"+type) // 2D detector data741 742 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + folderStr + ":det_"+type) // 2D detector data 744 743 745 744 Variable xDim,yDim,delQ … … 749 748 750 749 if(xDim<yDim) 751 WAVE qx = $("root:Packages:NIST:VSANS:VCALC:" + "qx_"+type)750 WAVE qx = $("root:Packages:NIST:VSANS:VCALC:" + folderStr + ":qx_"+type) 752 751 delQ = abs(qx[0][0] - qx[1][0])/2 753 752 else 754 WAVE qy = $("root:Packages:NIST:VSANS:VCALC:" + "qy_"+type)753 WAVE qy = $("root:Packages:NIST:VSANS:VCALC:" + folderStr + ":qy_"+type) 755 754 delQ = abs(qy[0][1] - qy[0][0])/2 756 755 endif 757 756 758 757 // set the global 759 Variable/G $("root:Packages:NIST:VSANS:VCALC:" + " delQ_"+type) = delQ758 Variable/G $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_"+type) = delQ 760 759 // Print "SET delQ = ",delQ," for ",type 761 760 … … 772 771 773 772 774 V_fDoBinning_QxQy2D( "", type)773 V_fDoBinning_QxQy2D(folderStr, type) 775 774 776 775 777 776 /// this is for a tall, narrow slit mode 778 // V_fBinDetector_byRows( type)777 // V_fBinDetector_byRows(folderStr,type) 779 778 780 779 End 780 781 781 782 782 Proc V_Graph_1D_detType(folderStr,type) … … 801 801 //Function fDoBinning_QxQy2D(inten,qx,qy,qz) 802 802 // 803 // this has been mod eified to accept different detector panels and to take arrays803 // this has been modified to accept different detector panels and to take arrays 804 804 // -- type = FL or FR or...other panel identifiers 805 805 // … … 809 809 Function V_fDoBinning_QxQy2D(folderStr,type) 810 810 String folderStr,type 811 812 // TODO: folderStr is ignored here813 folderStr = ""814 811 815 812 Variable nSets = 0 … … 818 815 Variable qVal,nq,var,avesq,aveisq 819 816 Variable binIndex,val 820 817 821 818 822 819 SetDataFolder root:Packages:NIST:VSANS:VCALC … … 829 826 case "FL": // execute if case matches expression 830 827 case "FR": 831 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_FL")832 WAVE inten = $(" det_"+type) // 2D detector data828 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_FL") 829 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+type) // 2D detector data 833 830 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 834 Wave qTotal = $(" qTot_"+type) // 2D q-values831 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+type) // 2D q-values 835 832 nSets = 1 836 833 break … … 838 835 case "FT": 839 836 case "FB": 840 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_FT")841 WAVE inten = $(" det_"+type) // 2D detector data837 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_FT") 838 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+type) // 2D detector data 842 839 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 843 Wave qTotal = $(" qTot_"+type) // 2D q-values840 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+type) // 2D q-values 844 841 nSets = 1 845 842 break … … 847 844 case "ML": 848 845 case "MR": 849 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_ML")850 WAVE inten = $(" det_"+type) // 2D detector data846 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_ML") 847 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+type) // 2D detector data 851 848 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 852 Wave qTotal = $(" qTot_"+type) // 2D q-values849 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+type) // 2D q-values 853 850 nSets = 1 854 851 break … … 856 853 case "MT": 857 854 case "MB": 858 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_MT")859 WAVE inten = $(" det_"+type) // 2D detector data855 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_MT") 856 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+type) // 2D detector data 860 857 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 861 Wave qTotal = $(" qTot_"+type) // 2D q-values858 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+type) // 2D q-values 862 859 nSets = 1 863 860 break 864 861 865 862 case "B": 866 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_B")867 WAVE inten = $(" det_"+type) // 2D detector data863 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_B") 864 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Back" + ":det_"+type) // 2D detector data 868 865 WAVE/Z iErr = $("iErr_"+type) // 2D errors -- may not exist, especially for simulation 869 Wave qTotal = $(" qTot_"+type) // 2D q-values866 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Back" +":qTot_"+type) // 2D q-values 870 867 nSets = 1 871 868 break 872 869 873 870 case "FLR": 874 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_FL")875 WAVE inten = $(" det_"+"FL") // 2D detector data871 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_FL") 872 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FL") // 2D detector data 876 873 WAVE/Z iErr = $("iErr_"+"FL") // 2D errors -- may not exist, especially for simulation 877 Wave qTotal = $(" qTot_"+"FL") // 2D q-values878 WAVE inten2 = $(" det_"+"FR") // 2D detector data874 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FL") // 2D q-values 875 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FR") // 2D detector data 879 876 WAVE/Z iErr2 = $("iErr_"+"FR") // 2D errors -- may not exist, especially for simulation 880 Wave qTotal2 = $(" qTot_"+"FR") // 2D q-values877 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FR") // 2D q-values 881 878 nSets = 2 882 879 break 883 880 884 881 case "FTB": 885 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_FT")886 WAVE inten = $(" det_"+"FT") // 2D detector data882 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_FT") 883 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FT") // 2D detector data 887 884 WAVE/Z iErr = $("iErr_"+"FT") // 2D errors -- may not exist, especially for simulation 888 Wave qTotal = $(" qTot_"+"FT") // 2D q-values889 WAVE inten2 = $(" det_"+"FB") // 2D detector data885 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FT") // 2D q-values 886 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FB") // 2D detector data 890 887 WAVE/Z iErr2 = $("iErr_"+"FB") // 2D errors -- may not exist, especially for simulation 891 Wave qTotal2 = $(" qTot_"+"FB") // 2D q-values888 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FB") // 2D q-values 892 889 nSets = 2 893 890 break 894 891 895 892 case "FLRTB": 896 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_FL")897 WAVE inten = $(" det_"+"FL") // 2D detector data893 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_FL") 894 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FL") // 2D detector data 898 895 WAVE/Z iErr = $("iErr_"+"FL") // 2D errors -- may not exist, especially for simulation 899 Wave qTotal = $(" qTot_"+"FL") // 2D q-values900 WAVE inten2 = $(" det_"+"FR") // 2D detector data896 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FL") // 2D q-values 897 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FR") // 2D detector data 901 898 WAVE/Z iErr2 = $("iErr_"+"FR") // 2D errors -- may not exist, especially for simulation 902 Wave qTotal2 = $(" qTot_"+"FR") // 2D q-values903 WAVE inten3 = $(" det_"+"FT") // 2D detector data899 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FR") // 2D q-values 900 WAVE inten3 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FT") // 2D detector data 904 901 WAVE/Z iErr3 = $("iErr_"+"FT") // 2D errors -- may not exist, especially for simulation 905 Wave qTotal3 = $(" qTot_"+"FT") // 2D q-values906 WAVE inten4 = $(" det_"+"FB") // 2D detector data902 Wave qTotal3 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FT") // 2D q-values 903 WAVE inten4 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" + ":det_"+"FB") // 2D detector data 907 904 WAVE/Z iErr4 = $("iErr_"+"FB") // 2D errors -- may not exist, especially for simulation 908 Wave qTotal4 = $(" qTot_"+"FB") // 2D q-values905 Wave qTotal4 = $("root:Packages:NIST:VSANS:VCALC:" + "Front" +":qTot_"+"FB") // 2D q-values 909 906 nSets = 4 910 907 break … … 912 909 913 910 case "MLR": 914 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_ML")915 WAVE inten = $(" det_"+"ML") // 2D detector data911 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_ML") 912 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"ML") // 2D detector data 916 913 WAVE/Z iErr = $("iErr_"+"ML") // 2D errors -- may not exist, especially for simulation 917 Wave qTotal = $(" qTot_"+"ML") // 2D q-values918 WAVE inten2 = $(" det_"+"MR") // 2D detector data914 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"ML") // 2D q-values 915 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MR") // 2D detector data 919 916 WAVE/Z iErr2 = $("iErr_"+"MR") // 2D errors -- may not exist, especially for simulation 920 Wave qTotal2 = $(" qTot_"+"MR") // 2D q-values917 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MR") // 2D q-values 921 918 nSets = 2 922 919 break 923 920 924 921 case "MTB": 925 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_MT")926 WAVE inten = $(" det_"+"MT") // 2D detector data922 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_MT") 923 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MT") // 2D detector data 927 924 WAVE/Z iErr = $("iErr_"+"MT") // 2D errors -- may not exist, especially for simulation 928 Wave qTotal = $(" qTot_"+"MT") // 2D q-values929 WAVE inten2 = $(" det_"+"MB") // 2D detector data925 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MT") // 2D q-values 926 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MB") // 2D detector data 930 927 WAVE/Z iErr2 = $("iErr_"+"MB") // 2D errors -- may not exist, especially for simulation 931 Wave qTotal2 = $(" qTot_"+"MB") // 2D q-values928 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MB") // 2D q-values 932 929 nSets = 2 933 930 break 934 931 935 932 case "MLRTB": 936 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + " delQ_ML")937 WAVE inten = $(" det_"+"ML") // 2D detector data933 NVAR delQ = $("root:Packages:NIST:VSANS:VCALC:" + "gDelQ_ML") 934 WAVE inten = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"ML") // 2D detector data 938 935 WAVE/Z iErr = $("iErr_"+"ML") // 2D errors -- may not exist, especially for simulation 939 Wave qTotal = $(" qTot_"+"ML") // 2D q-values940 WAVE inten2 = $(" det_"+"MR") // 2D detector data936 Wave qTotal = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"ML") // 2D q-values 937 WAVE inten2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MR") // 2D detector data 941 938 WAVE/Z iErr2 = $("iErr_"+"MR") // 2D errors -- may not exist, especially for simulation 942 Wave qTotal2 = $(" qTot_"+"MR") // 2D q-values943 WAVE inten3 = $(" det_"+"MT") // 2D detector data939 Wave qTotal2 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MR") // 2D q-values 940 WAVE inten3 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MT") // 2D detector data 944 941 WAVE/Z iErr3 = $("iErr_"+"MT") // 2D errors -- may not exist, especially for simulation 945 Wave qTotal3 = $(" qTot_"+"MT") // 2D q-values946 WAVE inten4 = $(" det_"+"MB") // 2D detector data942 Wave qTotal3 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MT") // 2D q-values 943 WAVE inten4 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" + ":det_"+"MB") // 2D detector data 947 944 WAVE/Z iErr4 = $("iErr_"+"MB") // 2D errors -- may not exist, especially for simulation 948 Wave qTotal4 = $(" qTot_"+"MB") // 2D q-values945 Wave qTotal4 = $("root:Packages:NIST:VSANS:VCALC:" + "Middle" +":qTot_"+"MB") // 2D q-values 949 946 nSets = 4 950 947 break … … 974 971 // note that the back panel of 320x320 (1mm res) results in 447 data points! 975 972 // - so I upped nq to 600 976 973 974 //******TODO****** -- where to put the averaged data -- right now, folderStr is forced to "" 977 975 // SetDataFolder $("root:"+folderStr) //should already be here, but make sure... 978 976 Make/O/D/N=(nq) $("root:Packages:NIST:VSANS:VCALC:"+folderStr+"iBin_qxqy"+"_"+type) -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_DetectorFrontViews.ipf
r935 r947 13 13 14 14 // Dimensions of detectors 15 NVAR F_LR_w = root:Packages:NIST:VSANS:VCALC: front_LR_w16 NVAR F_LR_h = root:Packages:NIST:VSANS:VCALC: front_LR_h17 NVAR F_TB_w = root:Packages:NIST:VSANS:VCALC: front_TB_w18 NVAR F_TB_h = root:Packages:NIST:VSANS:VCALC: front_TB_h15 NVAR F_LR_w = root:Packages:NIST:VSANS:VCALC:gFront_LR_w 16 NVAR F_LR_h = root:Packages:NIST:VSANS:VCALC:gFront_LR_h 17 NVAR F_TB_w = root:Packages:NIST:VSANS:VCALC:gFront_TB_w 18 NVAR F_TB_h = root:Packages:NIST:VSANS:VCALC:gFront_TB_h 19 19 20 NVAR M_LR_w = root:Packages:NIST:VSANS:VCALC: middle_LR_w21 NVAR M_LR_h = root:Packages:NIST:VSANS:VCALC: middle_LR_h22 NVAR M_TB_w = root:Packages:NIST:VSANS:VCALC: middle_TB_w23 NVAR M_TB_h = root:Packages:NIST:VSANS:VCALC: middle_TB_h20 NVAR M_LR_w = root:Packages:NIST:VSANS:VCALC:gMiddle_LR_w 21 NVAR M_LR_h = root:Packages:NIST:VSANS:VCALC:gMiddle_LR_h 22 NVAR M_TB_w = root:Packages:NIST:VSANS:VCALC:gMiddle_TB_w 23 NVAR M_TB_h = root:Packages:NIST:VSANS:VCALC:gMiddle_TB_h 24 24 25 NVAR B_h = root:Packages:NIST:VSANS:VCALC: back_h26 NVAR B_w = root:Packages:NIST:VSANS:VCALC: back_w25 NVAR B_h = root:Packages:NIST:VSANS:VCALC:gBack_h 26 NVAR B_w = root:Packages:NIST:VSANS:VCALC:gBack_w 27 27 28 28 … … 35 35 // these offset values are in mm !! 36 36 //in mm !! distance T/B are behind L/R - not to be confused with lateral offset 37 NVAR front_SDDOffset = root:Packages:NIST:VSANS:VCALC: front_SDDOffset38 NVAR middle_SDDOffset = root:Packages:NIST:VSANS:VCALC: middle_SDDOffset37 NVAR front_SDDOffset = root:Packages:NIST:VSANS:VCALC:gFront_SDDOffset 38 NVAR middle_SDDOffset = root:Packages:NIST:VSANS:VCALC:gMiddle_SDDOffset 39 39 40 40 -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_FrontDetectorPanelGraph.ipf
r935 r947 60 60 61 61 // Dimensions of detectors 62 NVAR F_LR_w = front_LR_w63 NVAR F_LR_h = front_LR_h64 NVAR F_TB_w = front_TB_w65 NVAR F_TB_h = front_TB_h66 67 NVAR M_LR_w = middle_LR_w68 NVAR M_LR_h = middle_LR_h69 NVAR M_TB_w = middle_TB_w70 NVAR M_TB_h = middle_TB_h71 72 NVAR B_h = back_h73 NVAR B_w = back_w62 NVAR F_LR_w = gFront_LR_w 63 NVAR F_LR_h = gFront_LR_h 64 NVAR F_TB_w = gFront_TB_w 65 NVAR F_TB_h = gFront_TB_h 66 67 NVAR M_LR_w = gMiddle_LR_w 68 NVAR M_LR_h = gMiddle_LR_h 69 NVAR M_TB_w = gMiddle_TB_w 70 NVAR M_TB_h = gMiddle_TB_h 71 72 NVAR B_h = gBack_h 73 NVAR B_w = gBack_w 74 74 75 75 76 76 // get the values from the panel 77 Variable F_LR_sep,F_TB_sep, F_SDD, F_offset , TB_SDD_offset77 Variable F_LR_sep,F_TB_sep, F_SDD, F_offset 78 78 Variable M_LR_sep,M_TB_sep, M_SDD, M_offset 79 79 Variable B_SDD, B_offset 80 80 81 // TODO -- this needs to be a global 82 TB_SDD_offset = 300 //in mm !! distance T/B are behind L/R - not to be confused with lateral offset 81 NVAR TB_SDD_offset = gFront_SDDOffset //in mm !! distance T/B are behind L/R - not to be confused with lateral offset 83 82 TB_SDD_offset /= 1000 //convert to meters for the plot 84 83 … … 287 286 288 287 // Dimensions of detectors 289 NVAR F_LR_w = front_LR_w290 NVAR F_LR_h = front_LR_h291 NVAR F_TB_w = front_TB_w292 NVAR F_TB_h = front_TB_h293 294 NVAR M_LR_w = middle_LR_w295 NVAR M_LR_h = middle_LR_h296 NVAR M_TB_w = middle_TB_w297 NVAR M_TB_h = middle_TB_h298 299 NVAR B_h = back_h300 NVAR B_w = back_w288 NVAR F_LR_w = gFront_LR_w 289 NVAR F_LR_h = gFront_LR_h 290 NVAR F_TB_w = gFront_TB_w 291 NVAR F_TB_h = gFront_TB_h 292 293 NVAR M_LR_w = gMiddle_LR_w 294 NVAR M_LR_h = gMiddle_LR_h 295 NVAR M_TB_w = gMiddle_TB_w 296 NVAR M_TB_h = gMiddle_TB_h 297 298 NVAR B_h = gBack_h 299 NVAR B_w = gBack_w 301 300 302 301 -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_UtilityProcedures.ipf
r935 r947 178 178 strswitch(type) 179 179 case "FL": 180 NVAR pixSizeX = front_L_pixelX180 NVAR pixSizeX = gFront_L_pixelX 181 181 break 182 182 case "FR": 183 NVAR pixSizeX = front_R_pixelX184 break 185 case "FT": 186 NVAR pixSizeX = front_T_pixelX183 NVAR pixSizeX = gFront_R_pixelX 184 break 185 case "FT": 186 NVAR pixSizeX = gFront_T_pixelX 187 187 break 188 188 case "FB": 189 NVAR pixSizeX = front_B_pixelX190 break 191 192 case "ML": 193 NVAR pixSizeX = middle_L_pixelX189 NVAR pixSizeX = gFront_B_pixelX 190 break 191 192 case "ML": 193 NVAR pixSizeX = gMiddle_L_pixelX 194 194 break 195 195 case "MR": 196 NVAR pixSizeX = middle_R_pixelX197 break 198 case "MT": 199 NVAR pixSizeX = middle_T_pixelX196 NVAR pixSizeX = gMiddle_R_pixelX 197 break 198 case "MT": 199 NVAR pixSizeX = gMiddle_T_pixelX 200 200 break 201 201 case "MB": 202 NVAR pixSizeX = middle_B_pixelX203 break 204 205 case "B": 206 NVAR pixSizeX = Back_pixelX202 NVAR pixSizeX = gMiddle_B_pixelX 203 break 204 205 case "B": 206 NVAR pixSizeX = gBack_pixelX 207 207 break 208 208 … … 225 225 strswitch(type) 226 226 case "FL": 227 NVAR pixSizeY = front_L_pixelY227 NVAR pixSizeY = gFront_L_pixelY 228 228 break 229 229 case "FR": 230 NVAR pixSizeY = front_R_pixelY231 break 232 case "FT": 233 NVAR pixSizeY = front_T_pixelY230 NVAR pixSizeY = gFront_R_pixelY 231 break 232 case "FT": 233 NVAR pixSizeY = gFront_T_pixelY 234 234 break 235 235 case "FB": 236 NVAR pixSizeY = front_B_pixelY237 break 238 239 case "ML": 240 NVAR pixSizeY = middle_L_pixelY236 NVAR pixSizeY = gFront_B_pixelY 237 break 238 239 case "ML": 240 NVAR pixSizeY = gMiddle_L_pixelY 241 241 break 242 242 case "MR": 243 NVAR pixSizeY = middle_R_pixelY244 break 245 case "MT": 246 NVAR pixSizeY = middle_T_pixelY243 NVAR pixSizeY = gMiddle_R_pixelY 244 break 245 case "MT": 246 NVAR pixSizeY = gMiddle_T_pixelY 247 247 break 248 248 case "MB": 249 NVAR pixSizeY = middle_B_pixelY250 break 251 252 case "B": 253 NVAR pixSizeY = Back_pixelY249 NVAR pixSizeY = gMiddle_B_pixelY 250 break 251 252 case "B": 253 NVAR pixSizeY = gBack_pixelY 254 254 break 255 255 … … 279 279 case "FT": 280 280 case "FB": 281 NVAR sdd_offset = front_SDDOffset //T/B are 300 mm farther back281 NVAR sdd_offset = gFront_SDDOffset //T/B are 300 mm farther back 282 282 break 283 283 … … 288 288 case "MT": 289 289 case "MB": 290 NVAR sdd_offset = Middle_SDDOffset //T/B are 30 cm farther back290 NVAR sdd_offset = gMiddle_SDDOffset //T/B are 30 cm farther back 291 291 break 292 292
Note: See TracChangeset
for help on using the changeset viewer.