Changeset 583 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/U_CALC.ipf

Timestamp:

Oct 27, 2009 2:07:14 PM (13 years ago)

Author:

srkline

Message:

Added a calculation of the total fraction of multiple coherent scattering. This value is now displayed on the 1D sim panel and on the UCALC panel. If the total fraction of multiple coherent scattering (relative to the fraction coherently scattered) is greater than 0.10, then the field turns red as a warning.

Changed the numbering and suffix scheme for saving simulated VAX raw data files.

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/U_CALC.ipf (modified) (4 diffs)

sans/Dev/trunk/NCNR_User_Procedures/Reduction/USANS/U_CALC.ipf

@@ -134 +134 @@
         Variable/G g_1DFracScatt= 0     // ??
         Variable/G g_1DEstTrans = 0     // ? can I calculate this?
+        Variable/G g_MultScattFraction = 0
+
 
 // not on panel yet
@@ -311 +313 @@
         
         GroupBox group0,pos={5,1},size={493,177},title="Instrument Setup"
-        GroupBox group1,pos={5,183},size={240,147},title="Sample Setup"
-        GroupBox group2,pos={5,343},size={259,147},title="Results"
+        GroupBox group1,pos={5,183},size={240,110},title="Sample Setup"
+        GroupBox group2,pos={5,310},size={240,170},title="Results"
         
         PopupMenu popup0,pos={17,19},size={165,20},title="Sample Aperture Diam (in)"
@@ -455 +457 @@
         
 // a box for the results
-        SetVariable totalTime,pos={left+20,370},size={150,15},title="Count time (h:m)",value= gTotTimeStr
-        ValDisplay valdisp0_2,pos={left+20,395},size={220,13},title="Fraction of beam scattered"
-        ValDisplay valdisp0_2,limits={0,0,0},barmisc={0,1000},value= root:Packages:NIST:USANS:Globals:U_Sim:g_1DFracScatt
-        ValDisplay valdisp0_2,disable=1
-        ValDisplay valdisp0_3,pos={left+20,420},size={220,13},title="Estimated transmission"
-        ValDisplay valdisp0_3,limits={0,0,0},barmisc={0,1000},value=root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans
-        ValDisplay valdisp0_3,disable=1
-        Button button1,pos={left+20,400},size={150,20},proc=U_SavePanelProc,title="Save PNG"
-        Button button2,pos={left+20,430},size={150,20},proc=U_ConfigTextProc,title="Config Text"
-        Button button3,pos={left+20,460},size={150,20},proc=U_SaveButtonProc,title="Save Simulated Data"
+        SetVariable totalTime,pos={left+20,335},size={150,15},title="Count time (h:m)",value= gTotTimeStr
+        SetVariable totalTime,limits={-inf,inf,0},noedit=1
+        SetVariable valdisp0_2,pos={left+20,365},size={190,15},title="Multiple Coherent Scattering"
+        SetVariable valdisp0_2,limits={-inf,inf,0},noedit=1,value=g_MultScattFraction
+//      ValDisplay valdisp0_3,pos={left+20,420},size={220,13},title="Estimated transmission"
+//      ValDisplay valdisp0_3,limits={0,0,0},barmisc={0,1000},value=root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans
+//      ValDisplay valdisp0_3,disable=1
+        Button button1,pos={left+20,390},size={150,20},proc=U_SavePanelProc,title="Save PNG"
+        Button button2,pos={left+20,420},size={150,20},proc=U_ConfigTextProc,title="Config Text"
+        Button button3,pos={left+20,450},size={150,20},proc=U_SaveButtonProc,title="Save Simulated Data"
 
 // help, done buttons
@@ -906 +908 @@
 //                      sprintf abortStr,"sig_sas = %g. Please check that the model coefficients have a zero background, or the low q is well-behaved.",sig_sas
 //              endif
+                // calculate the multiple scattering fraction for display (10/2009)
+                NVAR mScat = root:Packages:NIST:USANS:Globals:U_Sim:g_MultScattFraction
+                Variable nMax=10,tau
+                mScat=0
+                tau = thick*sig_sas
+                // this sums the normalized scattering P', so the result is the fraction of multiply coherently scattered
+                // neutrons out of those that were scattered
+                for(ii=2;ii<nMax;ii+=1)
+                        mScat += tau^(ii)/factorial(ii)
+        //              print tau^(ii)/factorial(ii)
+                endfor
                 estTrans = exp(-1*thick*sig_sas)                //thickness and sigma both in units of cm
+                mscat *= (estTrans)/(1-estTrans)
+                
+                if(mScat > 0.1)
+                        SetVariable valdisp0_2 win=UCALC,labelBack=(65535,32768,32768)
+                else
+                        SetVariable valdisp0_2 win=UCALC,labelBack=0
+                endif
+
+
                 Print "Sig_sas = ",sig_sas