- Timestamp:
- Aug 24, 2012 5:16:39 PM (11 years ago)
- Location:
- sans/Dev/trunk/NCNR_User_Procedures/Reduction/Polarization
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
sans/Dev/trunk/NCNR_User_Procedures/Reduction/Polarization/Pol_FlipperPanel.ipf
r858 r863 214 214 WAVE/T cell=$("CondCell_"+popStr) 215 215 cell[0] = "Enter cell name" 216 // SetDimLabel 0,-1,'Cell Name?',cell 216 217 217 218 Make/O/D/N=(1,8) $("Cond_"+popStr) 218 219 WAVE cond = $("Cond_"+popStr) 219 220 // set the column labels 220 SetDimLabel 1,0, UU_Trans,cond221 SetDimLabel 1,1, DU_Trans,cond222 SetDimLabel 1,2, DD_Trans,cond223 SetDimLabel 1,3, UD_Trans,cond224 SetDimLabel 1,4, Blocked,cond221 SetDimLabel 1,0,'UU_Trans?',cond 222 SetDimLabel 1,1,'DU_Trans?',cond 223 SetDimLabel 1,2,'DD_Trans?',cond 224 SetDimLabel 1,3,'UD_Trans?',cond 225 SetDimLabel 1,4,'Blocked?',cond 225 226 SetDimLabel 1,5,Pol_SM_FL,cond 226 227 SetDimLabel 1,6,Pol_SM,cond //for a mask wave, non-zero is used in the fit 227 SetDimLabel 1,7, Include,cond228 SetDimLabel 1,7,'Include?',cond 228 229 cond[0][7] = 1 //default to include the point 229 230 … … 380 381 381 382 //include this row of data? 382 if(w[selRow][% Include] == 1)383 if(w[selRow][%'Include?'] == 1) 383 384 nRowsIncluded += 1 384 385 … … 399 400 400 401 // do the calculations: 401 // cr1 = TotalCR_FromRun(w[selRow][% UU_Trans],err_cr1,1)402 // cr2 = TotalCR_FromRun(w[selRow][% DU_Trans],err_cr2,1)403 // cr3 = TotalCR_FromRun(w[selRow][% DD_Trans],err_cr3,1)404 // cr4 = TotalCR_FromRun(w[selRow][% UD_Trans],err_cr4,1)405 // cr5 = TotalCR_FromRun(w[selRow][% Blocked],err_cr5,1) //blocked beam is NOT normalized to zero attenuators402 // cr1 = TotalCR_FromRun(w[selRow][%'UU_Trans?'],err_cr1,1) 403 // cr2 = TotalCR_FromRun(w[selRow][%'DU_Trans?'],err_cr2,1) 404 // cr3 = TotalCR_FromRun(w[selRow][%'DD_Trans?'],err_cr3,1) 405 // cr4 = TotalCR_FromRun(w[selRow][%'UD_Trans?'],err_cr4,1) 406 // cr5 = TotalCR_FromRun(w[selRow][%'Blocked?'],err_cr5,1) //blocked beam is NOT normalized to zero attenuators 406 407 // Print "The Blocked CR is *NOT* rescaled to zero attenuators -- FlipperAverageButtonProc" 407 cr1 = TotalCR_FromRun(w[selRow][% UU_Trans],err_cr1,0)408 cr2 = TotalCR_FromRun(w[selRow][% DU_Trans],err_cr2,0)409 cr3 = TotalCR_FromRun(w[selRow][% DD_Trans],err_cr3,0)410 cr4 = TotalCR_FromRun(w[selRow][% UD_Trans],err_cr4,0)411 cr5 = TotalCR_FromRun(w[selRow][% Blocked],err_cr5,0) //blocked beam is normalized to zero attenuators408 cr1 = TotalCR_FromRun(w[selRow][%'UU_Trans?'],err_cr1,0) 409 cr2 = TotalCR_FromRun(w[selRow][%'DU_Trans?'],err_cr2,0) 410 cr3 = TotalCR_FromRun(w[selRow][%'DD_Trans?'],err_cr3,0) 411 cr4 = TotalCR_FromRun(w[selRow][%'UD_Trans?'],err_cr4,0) 412 cr5 = TotalCR_FromRun(w[selRow][%'Blocked?'],err_cr5,0) //blocked beam is normalized to zero attenuators 412 413 Print "The Blocked CR *IS* rescaled to zero attenuators -- FlipperAverageButtonProc" 413 414 … … 506 507 err_gam = NumberByKey("err_gamma", gCellKW, "=", ",", 0) 507 508 508 fname = FindFileFromRunNumber(w[selRow][% UU_Trans])509 fname = FindFileFromRunNumber(w[selRow][%'UU_Trans?']) 509 510 t1str = getFileCreationDate(fname) 510 511 t1 = ElapsedHours(t0Str,t1Str) 511 512 512 fname = FindFileFromRunNumber(w[selRow][% DU_Trans])513 fname = FindFileFromRunNumber(w[selRow][%'DU_Trans?']) 513 514 t2str = getFileCreationDate(fname) 514 515 t2 = ElapsedHours(t0Str,t2Str) … … 768 769 cond = 0 769 770 calc = 0 770 cellW = ""771 // cellW = "" 771 772 cond[0][7] = 1 //default to include the point 772 773 … … 867 868 868 869 // are all file numbers valid? 869 fname = FindFileFromRunNumber(w[selRow][% UU_Trans])870 fname = FindFileFromRunNumber(w[selRow][%'UU_Trans?']) 870 871 if(cmpstr(fname,"")==0) 871 DoAlert 0,"UU_Trans run "+num2str(w[selRow][% UU_Trans])+" is not a valid run number"872 DoAlert 0,"UU_Trans run "+num2str(w[selRow][%'UU_Trans?'])+" is not a valid run number" 872 873 err = 1 873 874 else … … 876 877 endif 877 878 878 fname = FindFileFromRunNumber(w[selRow][% DU_Trans])879 fname = FindFileFromRunNumber(w[selRow][%'DU_Trans?']) 879 880 if(cmpstr(fname,"")==0) 880 DoAlert 0,"DU_Trans run "+num2str(w[selRow][% DU_Trans])+" is not a valid run number"881 DoAlert 0,"DU_Trans run "+num2str(w[selRow][%'DU_Trans?'])+" is not a valid run number" 881 882 err = 1 882 883 else … … 885 886 endif 886 887 887 fname = FindFileFromRunNumber(w[selRow][% DD_Trans])888 fname = FindFileFromRunNumber(w[selRow][%'DD_Trans?']) 888 889 if(cmpstr(fname,"")==0) 889 DoAlert 0,"DD_Trans run "+num2str(w[selRow][% DD_Trans])+" is not a valid run number"890 DoAlert 0,"DD_Trans run "+num2str(w[selRow][%'DD_Trans?'])+" is not a valid run number" 890 891 err = 1 891 892 else … … 894 895 endif 895 896 896 fname = FindFileFromRunNumber(w[selRow][% UD_Trans])897 fname = FindFileFromRunNumber(w[selRow][%'UD_Trans?']) 897 898 if(cmpstr(fname,"")==0) 898 DoAlert 0,"UD_Trans run "+num2str(w[selRow][% UD_Trans])+" is not a valid run number"899 DoAlert 0,"UD_Trans run "+num2str(w[selRow][%'UD_Trans?'])+" is not a valid run number" 899 900 err = 1 900 901 else … … 903 904 endif 904 905 905 fname = FindFileFromRunNumber(w[selRow][% Blocked])906 fname = FindFileFromRunNumber(w[selRow][%'Blocked?']) 906 907 if(cmpstr(fname,"")==0) 907 DoAlert 0,"Blocked run "+num2str(w[selRow][% Blocked])+" is not a valid run number"908 DoAlert 0,"Blocked run "+num2str(w[selRow][%'Blocked?'])+" is not a valid run number" 908 909 err = 1 909 910 else -
sans/Dev/trunk/NCNR_User_Procedures/Reduction/Polarization/Pol_PolarizationPanels.ipf
r859 r863 634 634 WAVE decay = $("Decay_"+popStr) 635 635 // set the column labels 636 SetDimLabel 1,0, Trans_He_In,decay637 SetDimLabel 1,1, Trans_He_Out,decay638 SetDimLabel 1,2, Blocked,decay636 SetDimLabel 1,0,'Trans_He_In?',decay 637 SetDimLabel 1,1,'Trans_He_Out?',decay 638 SetDimLabel 1,2,'Blocked?',decay 639 639 SetDimLabel 1,3,mu_star,decay 640 640 SetDimLabel 1,4,Effective_Pol,decay 641 641 SetDimLabel 1,5,Atomic_Pol,decay 642 642 SetDimLabel 1,6,T_Major,decay 643 SetDimLabel 1,7, Include,decay //for a mask wave, non-zero is used in the fit643 SetDimLabel 1,7,'Include?',decay //for a mask wave, non-zero is used in the fit 644 644 SetDimLabel 1,8,elapsed_hr,decay 645 645 decay[0][7] = 1 //default to include the point … … 729 729 if(selRow == 0 && !overrideT0) 730 730 //find T0 731 fname = FindFileFromRunNumber(w[0][% Trans_He_In])731 fname = FindFileFromRunNumber(w[0][%'Trans_He_In?']) 732 732 t0str = getFileCreationDate(fname) 733 733 SVAR gT0 = root:Packages:NIST:Polarization:Cells:gT0 … … 758 758 // 1 for each file, return the count rate and err_CR (normalize to atten or not) 759 759 760 cr1 = TotalCR_FromRun(w[selRow][% Trans_He_In],err_cr1,0)761 cr2 = TotalCR_FromRun(w[selRow][% Trans_He_Out],err_cr2,0)762 // cr3 = TotalCR_FromRun(w[selRow][% Blocked],err_cr3,1) //blocked beam is NOT normalized to zero attenuators760 cr1 = TotalCR_FromRun(w[selRow][%'Trans_He_In?'],err_cr1,0) 761 cr2 = TotalCR_FromRun(w[selRow][%'Trans_He_Out?'],err_cr2,0) 762 // cr3 = TotalCR_FromRun(w[selRow][%'Blocked?'],err_cr3,1) //blocked beam is NOT normalized to zero attenuators 763 763 // Print "************The Blocked CR is *NOT* rescaled to zero attenuators -- CalcRowParamButton" 764 cr3 = TotalCR_FromRun(w[selRow][% Blocked],err_cr3,0) //blocked beam is normalized to zero attenuators764 cr3 = TotalCR_FromRun(w[selRow][%'Blocked?'],err_cr3,0) //blocked beam is normalized to zero attenuators 765 765 Print "************The Blocked CR *is* rescaled to zero attenuators -- CalcRowParamButton " 766 766 … … 801 801 802 802 // elapsed hours 803 fname = FindFileFromRunNumber(w[selRow][% Trans_He_In])803 fname = FindFileFromRunNumber(w[selRow][%'Trans_He_In?']) 804 804 t1str = getFileCreationDate(fname) 805 805 w[selRow][%elapsed_hr] = ElapsedHours(t0Str,t1Str) … … 1105 1105 Make/O/D/N=(num) tmp_Mask,tmp_hr,tmp_Y,tmp_err_Y,tmp_Y2 1106 1106 1107 tmp_Mask = decay[p][% Include]1107 tmp_Mask = decay[p][%'Include?'] 1108 1108 tmp_hr = decay[p][%elapsed_hr] 1109 1109 … … 1441 1441 1442 1442 // are all file numbers valid? 1443 fname = FindFileFromRunNumber(w[selRow][% Trans_He_In])1443 fname = FindFileFromRunNumber(w[selRow][%'Trans_He_In?']) 1444 1444 if(cmpstr(fname,"")==0) 1445 DoAlert 0,"Trans_He_In run "+num2str(w[selRow][% Trans_He_In])+" is not a valid run number"1445 DoAlert 0,"Trans_He_In run "+num2str(w[selRow][%'Trans_He_In?'])+" is not a valid run number" 1446 1446 err = 1 1447 1447 else … … 1450 1450 endif 1451 1451 1452 fname = FindFileFromRunNumber(w[selRow][% Trans_He_Out])1452 fname = FindFileFromRunNumber(w[selRow][%'Trans_He_Out?']) 1453 1453 if(cmpstr(fname,"")==0) 1454 DoAlert 0,"Trans_He_Out run "+num2str(w[selRow][% Trans_He_Out])+" is not a valid run number"1454 DoAlert 0,"Trans_He_Out run "+num2str(w[selRow][%'Trans_He_Out?'])+" is not a valid run number" 1455 1455 err = 1 1456 1456 else … … 1459 1459 endif 1460 1460 1461 fname = FindFileFromRunNumber(w[selRow][% Blocked])1461 fname = FindFileFromRunNumber(w[selRow][%'Blocked?']) 1462 1462 if(cmpstr(fname,"")==0) 1463 DoAlert 0,"Blocked run "+num2str(w[selRow][% Blocked])+" is not a valid run number"1463 DoAlert 0,"Blocked run "+num2str(w[selRow][%'Blocked?'])+" is not a valid run number" 1464 1464 err = 1 1465 1465 else
Note: See TracChangeset
for help on using the changeset viewer.