Changeset 472 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/SASCALC.ipf

Timestamp:

Feb 3, 2009 3:57:08 PM (14 years ago)

Author:

srkline

Message:

A large number of changes and fixes:

--168/169/170: panels and windows are now at least on-screen for all packages. Tested
with 1024x768 resolution.
-- closed ticket 176 which was a question about resampling data to generate error estimates
on fitted parameters. Useful for reflectometry, not needed for SANS.
--157: bug of low Q power law extrapolation in Invariant fixed by avoiding q==0
--178/180: Tr/Tw? notification in USANS. log/lin checkbox for display.
--167: saveData checkbox for USANS not behaving well. turns off/on better now.
--197: changed all (?) 1D writing routines to enforce 26 characters as the maximum length
to make sure that file loading will never exceed 31 characters

-- lots of changes to MonteCarlo? and SASCALC

• SASCALC now enforces *exact* lens conditions, rather than a close approximation
• improved MonteCarlo? interface panel
• added writer for simlulated VAX binary data file
• can save 2D as ABS or raw counts
• can freeze w/no offset

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/SASCALC.ipf (modified) (24 diffs)

sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/SASCALC.ipf

@@ -17 +17 @@
 //                                               Lens flux, trans is not corrected for lens/prism transmission
 //                                               Lenses can still be inserted in incorrect cases, and are not automatically taken out
+// 27 JAN 2009 SRK - Changed behavior of Lens checkbox. Now, it SETS parameters as needed for proper
+//                                               configuration. 17.2 can be typed in for lens/prism on NG3. Invalid conditions
+//                                               will automatically uncheck the box
 //
 // calculate what q-values you get based on the instruments settings
@@ -25 +28 @@
 //
 // To Do:
-// - add in instrument conditions for lens/(lens+prism) configurations
-// - proper resolution calculation for lens/prism
-//
 //
 // Optional:
@@ -61 +61 @@
                 ReCalculateInten(1)             //will use defaults
         Endif
-        
-        DoWindow/F MC_SASCALC
-        if(V_flag==0)
-                MC_SASCALC()
-        endif
+
+// now a checkbox as needed
+//      DoWindow/F MC_SASCALC
+//      if(V_flag==0)
+//              MC_SASCALC()
+//      endif
 End
 
@@ -108 +109 @@
         
         // for the MC simulation
+        Variable/G root:Packages:NIST:SAS:gRanDateTime=datetime
         Variable/G root:Packages:NIST:SAS:gImon = 10000
         Variable/G root:Packages:NIST:SAS:gThick = 0.1
@@ -115 +117 @@
         Variable/G root:Packages:NIST:SAS:gCntTime = 1
         Variable/G root:Packages:NIST:SAS:gDoMonteCarlo = 0
+        Variable/G root:Packages:NIST:SAS:gRawCounts = 0
+        Variable/G root:Packages:NIST:SAS:gSaveIndex = 100
+        String/G root:Packages:NIST:SAS:gSavePrefix = "SIMUL"
         Make/O/D/N=10 root:Packages:NIST:SAS:results = 0
-        Make/O/T/N=10 root:Packages:NIST:SAS:results_desc = {"total X-section (1/cm)","SAS X-section (1/cm)","number that scatter","number that reach detector","avg # times scattered","fraction single coherent","fraction double coherent","fraction multiple scattered","fraction transmitted","detector counts w/beamstop"}
+        Make/O/T/N=10 root:Packages:NIST:SAS:results_desc = {"total X-section (1/cm)","SAS X-section (1/cm)","number that scatter","number that reach detector","avg # times scattered","fraction single coherent","fraction double coherent","fraction multiple scattered","fraction transmitted","detector counts w/o beamstop"}
         
         //tick labels for SDD slider
@@ -132 +137 @@
         //for plotting
         Variable/G root:Packages:NIST:SAS:gFreezeCount=1                //start the count at 1 to keep Jeff happy
+        Variable/G root:Packages:NIST:SAS:gDoTraceOffset=1              // (1==Yes, offset 2^n), 0==turn off the offset
+        
 End
 
@@ -283 +290 @@
 
         PauseUpdate; Silent 1           // building window...
+
+// if I make the graph a subwindow in a panel, it breaks the "append 1d" from the wrapper       
+//      NewPanel/W=(5,44,463,570)/K=1 as "SASCALC"
+//      DoWindow/C SASCALC
+//      ModifyPanel cbRGB=(49151,53155,65535)
+//
+//      
+//      String fldrSav0= GetDataFolder(1)
+//      SetDataFolder root:Packages:NIST:SAS:
+//      
+//      Display/HOST=#/W=(5,200,463,570) aveint vs qval
+//      ModifyGraph mode=3
+//      ModifyGraph marker=19
+//      ModifyGraph rgb=(0,0,0)
+//      Modifygraph log=1
+//      Modifygraph grid=1
+//      Modifygraph mirror=2
+//      ModifyGraph msize(aveint)=2
+//      ErrorBars/T=0 aveint Y,wave=(sigave,sigave)
+//      Label bottom, "Q (1/A)"
+//      Label left, "Relative Intensity"
+//      legend
+//      
+//      RenameWindow #,G_aveint
+//      SetActiveSubwindow ##
+
+//// end panel commands
+
+/// draw as a graph
         String fldrSav0= GetDataFolder(1)
         SetDataFolder root:Packages:NIST:SAS:
-        Display /W=(5,44,463,570)/K=1 aveint vs qval as "SASCALC"
+        
+        Display/W=(5,44,463,570)/K=1  aveint vs qval as "SASCALC"
         DoWindow/C SASCALC
         ModifyGraph cbRGB=(49151,53155,65535)
@@ -299 +336 @@
         Label left, "Relative Intensity"
         legend
+        
+        ControlBar/T 200
+
         SetDataFolder fldrSav0
-
-
-        ControlBar 200
         
         Slider SC_Slider,pos={11,46},size={150,45},proc=GuideSliderProc,live=0
@@ -344 +381 @@
         CheckBox checkLens,value=root:Packages:NIST:SAS:gUsingLenses
         
-//      CheckBox checkSim,pos={20,165},size={44,14},proc=SimCheckProc,title="Simulation?"
-//      CheckBox checkSim,value=0
+        CheckBox checkSim,pos={6,175},size={44,14},proc=SimCheckProc,title="MC Simulation?"
+        CheckBox checkSim,value=0
         
         // set up a fake dependency to trigger recalculation
@@ -431 +468 @@
         Variable event  // bit field: bit 0: value set, 1: mouse down, 2: mouse up, 3: mouse moved
         
+        Variable recalc=0
+        
         if(event %& 0x1)        // bit 0, value set
-                if(sliderValue != 0)
-                        LensCheckProc("",0)             //make sure lenses are deselected
+                if(cmpstr(ctrlName,"") != 0)            //here by direct action, so do LensCheck and recalculate
+                        recalc=1
+                        LensCheckProc("",2)             //make sure lenses are deselected
                 endif
                 sourceToSampleDist()            //updates the SSD global and wave
                 //change the sourceAp popup, SDD range, etc
                 UpdateControls()
-                ReCalculateInten(1)
+                ReCalculateInten(recalc)
         endif
         return 0
@@ -450 +490 @@
         Variable event  // bit field: bit 0: value set, 1: mouse down, 2: mouse up, 3: mouse moved
 
+        Variable recalc=0
         if(event %& 0x1)        // bit 0, value set
-                if(sliderValue < 1300)
-                        LensCheckProc("",0)             //make sure lenses are deselected
+                if(cmpstr(ctrlName,"") != 0)
+                        recalc=1
+                        LensCheckProc("",2)             //make sure lenses are only selected for valid configurations
                 endif
                 sampleToDetectorDist()          //changes the SDD and wave (DetDist is the global)
-                ReCalculateInten(1)
+                ReCalculateInten(recalc)
         endif
 
@@ -493 +535 @@
                 initNG7()
         endif
+        LensCheckProc("",2)             //check if lenses are still valid (they won't be)
         UpdateControls()
         ReCalculateInten(1)
@@ -520 +563 @@
 
 //lenses (or prisms) in/out changes resolution
+//
+// passing in a checked == 2 will do a check of the configuration without
+// affecting the box state
+//
+// if an invalid configuration is detected, the box is unchecked
+// and the lens flag is set to zero (=out).
+//
+// When necessary controls are "popped", a ctrlName="" is passed to signal the control
+// to NOT recalculate the intensity, and to NOT (recursively) call LensCheckProc again
+//
+// currently, the 17.2 A for lens/prism @ ng3 must be typed in
+//
 Function LensCheckProc(ctrlName,checked) : CheckBoxControl
         String ctrlName
@@ -531 +586 @@
         NVAR lam = root:Packages:NIST:SAS:gLambda
         NVAR dist = root:Packages:NIST:SAS:gDetDist
-        
-        if( (Ng != 0) || (lam < 8) || (dist < 1300) )
+        NVAR instrument = root:Packages:NIST:SAS:instrument
+        Wave rw=root:Packages:NIST:SAS:realsRead
+        
+        // directly uncheck the box, just set the flag and get out
+        if(checked == 0)
                 lens = 0
                 CheckBox checkLens,value=0
+                rw[28]=0                //flag for lenses out
+                ReCalculateInten(1)
                 return(0)
         endif
-        lens = checked
-        
-        ReCalculateInten(1)
-End
-
-////simulation control panel
-//Function SimCheckProc(ctrlName,checked) : CheckBoxControl
-//      String ctrlName
-//      Variable checked
-//
-//      if(checked)
-//              DoWindow/F MC_SASCALC
-//              if(V_flag==0)
-//                      Execute "MC_SASCALC()"
-//              endif
-//      endif
-//      return(0)
-//End
+        
+        // check the box, enforce the proper conditions
+        if(checked == 1)
+                lens = checked
+                if(instrument == 3)
+                        dist = 1317
+                        DetDistSliderProc("",1317,1)
+                        
+                        lam = 8.4
+                        LambdaSetVarProc("",8.4,"8.4","") 
+
+                        ng=0
+                        GuideSliderProc("",0,1)         //this updates the controls to the new # of guides
+                        
+                        PopupMenu popup0,mode=1,popvalue="1.43 cm"              //first item in source aperture menu
+                        
+                        PopupMenu popup0_2,mode=2               //deltaLambda
+                        ControlInfo popup0_2
+                        DeltaLambdaPopMenuProc("",0,S_value)                    //zero as 2nd param skips recalculation
+                else
+                        dist = 1531
+                        DetDistSliderProc("",1531,1)
+                        
+                        lam = 8.09
+                        LambdaSetVarProc("",8.09,"8.09","") 
+                        
+                        ng=0
+                        GuideSliderProc("",0,1)
+                        PopupMenu popup0,mode=1,popvalue="1.43 cm"              //first item
+                        
+                        PopupMenu popup0_2,mode=2               //deltaLambda
+                        ControlInfo popup0_2
+                        DeltaLambdaPopMenuProc("",0,S_value)                    //zero as 2nd param skips recalculation
+                endif
+                rw[28]=1                //flag for lenses in (not the true number, but OK)
+                ReCalculateInten(1)
+        endif
+        
+        // this is my internal check to see if conditions are still valid
+        // I'll uncheck as needed
+        if(checked == 2)
+
+                // source aperture must be 1.43 cm diameter
+                // number of guides must be zero
+                Variable a1 = sourceApertureDiam()
+                if(a1 != 1.43  || Ng !=0)
+                        lens = 0
+                        CheckBox checkLens,value=0
+                        rw[28]=0                //flag for lenses out
+                        return(0)
+                endif
+        
+                // instrument specific distance requirements
+                if(instrument == 3 && dist != 1317)
+                        lens = 0
+                        CheckBox checkLens,value=0
+                        rw[28]=0                //flag for lenses out
+                        return(0)
+                endif
+        
+                if(instrument == 7 && dist != 1531)
+                        lens = 0
+                        CheckBox checkLens,value=0
+                        rw[28]=0                //flag for lenses out
+                        return(0)
+                endif
+                
+                // instrument specific wavelength requirements
+                if(instrument == 3 && !(lam == 8.4 || lam == 17.2) )
+                        lens = 0
+                        CheckBox checkLens,value=0
+                        rw[28]=0                //flag for lenses out
+                        return(0)
+                endif
+                
+                if(instrument == 7 && lam != 8.09 )
+                        lens = 0
+                        CheckBox checkLens,value=0
+                        rw[28]=0                //flag for lenses out
+                        return(0)
+                endif
+                
+        endif
+
+        return(1)               //return value not used
+End
+
+//simulation ccontrols as a control bar that toggles on/off to the right
+Function SimCheckProc(ctrlName,checked) : CheckBoxControl
+        String ctrlName
+        Variable checked
+
+        if(checked)
+                //MoveWindow/W=SASCALC 5,44,763,570             //to resize
+                // draw the controls
+                DoWindow/F MC_SASCALC
+                if(V_flag==0)
+                        Execute "MC_SASCALC()"
+                        AutoPositionWindow/M=1/R=SASCALC MC_SASCALC
+                endif
+        else
+                //get rid of the controls
+                DoWindow MC_SASCALC
+                if(V_flag!=0)
+                        KillWindow MC_SASCALC
+                endif
+                //MoveWindow/W=SASCALC 5,44,463,570             //to resize
+                //KillWindow SASCALC#T_results
+        endif
+
+        return(0)
+End
 
 // change the source aperture
@@ -565 +719 @@
         Variable a1 = sourceApertureDiam()              //sets the new value in the wave
         
+        // if LensCheckProc is calling this, don't call it again.
+        // if this is called by a direct pop, then check
+        if(cmpstr(ctrlName,"") != 0)
+                LensCheckProc("",2)             //make sure lenses are only selected for valid configurations
+        endif
+        
         ReCalculateInten(popnum)                //skip the recalculation if I pass in a zero
 End
@@ -588 +748 @@
         
         sampleToDetectorDist()
-        
+        if(cmpstr(ctrlName,"") != 0)
+                LensCheckProc("",2)             //make sure lenses are only selected for valid configurations
+        endif
         ReCalculateInten(1)
 End
@@ -615 +777 @@
 
         WAVE rw=root:Packages:NIST:SAS:realsRead
+        Variable recalc=0
+        
         rw[26] = str2num(varStr)
-        ReCalculateInten(1)
+        if(cmpstr(ctrlName,"") != 0)
+                recalc=1
+                LensCheckProc("",2)             //make sure lenses are only selected for valid configurations
+        endif
+        ReCalculateInten(recalc)
         return(0)
 End
@@ -663 +831 @@
                 return(0)
         endif
+//      Print "recalculate"
         
         // update the wave with the beamstop diameter here, since I don't know what
@@ -784 +953 @@
 
                 t0 = (stopMSTimer(-2) - t0)*1e-6
-                Printf  "MC sim time = %g seconds\r\r",t0
+                Printf  "MC sim time = %g seconds\r",t0
                 
                 trans = results[8]                      //(n1-n2)/n1
@@ -794 +963 @@
                 
                 linear_data[xCtr][yCtr] = 0                     //snip out the transmitted spike
-                Print "counts on detector not transmitted = ",sum(linear_data,-inf,inf)
+//              Print "counts on detector not transmitted = ",sum(linear_data,-inf,inf)
 
                 // or simulate a beamstop
@@ -805 +974 @@
                 
                 linear_data *= tmp_mask
-                Print "counts on detector not behind beamstop = ",sum(linear_data,-inf,inf)
                 results[9] = sum(linear_data,-inf,inf)
+                //              Print "counts on detector not behind beamstop = ",results[9]
                 
                 // convert to absolute scale
@@ -812 +981 @@
 //              kappa = beamInten*pi*r1*r1*thick*(pixSize/sdd)^2*trans*(iMon/beaminten)
                 kappa = thick*(pixSize/sdd)^2*trans*iMon
-
-                linear_data = linear_data / kappa
-                
+                
+                //use kappa to get back to counts => linear_data = round(linear_data*kappa)
+                Note/K linear_data ,"KAPPA="+num2str(kappa)+";"
+                
+                NVAR rawCts = root:Packages:NIST:SAS:gRawCounts
+                if(!rawCts)                     //go ahead and do the linear scaling
+                        linear_data = linear_data / kappa
+                endif           
                 data = linear_data
                 
                 // re-average the 2D data
                 S_CircularAverageTo1D("SAS")
-                // multiply either estimate by beamstop shadowing
-                
-//              aveint *= fSubS
-
+                
                 // put the new result into the simulation folder
                 Fake1DDataFolder(qval,aveint,sigave,sigmaQ,qbar,fSubs,"Simulation")     
@@ -873 +1044 @@
         Duplicate/O sigave,$("sigave_"+num2str(ct))
         Appendtograph $("aveint_"+num2str(ct)) vs $("qval_"+num2str(ct))
-        ModifyGraph mode=3
-        ModifyGraph marker=19
+        ModifyGraph mode($("aveint_"+num2str(ct)))=3
+        ModifyGraph marker($("aveint_"+num2str(ct)))=19
         ModifyGraph msize($("aveint_"+num2str(ct)))=2
         ErrorBars/T=0 $("aveint_"+num2str(ct)) Y,wave=($("sigave_"+num2str(ct)),$("sigave_"+num2str(ct)))
@@ -911 +1082 @@
         endswitch
         
+        NVAR doTraceOffset = root:Packages:NIST:SAS:gDoTraceOffset
         NVAR offset = root:Packages:NIST:SAS:gModelOffsetFactor
-        offset = 2^ct
-        //multiply by current offset (>=1)
-        Wave inten = $("aveint_"+num2str(ct))
-        inten *= offset
-//      Print "new offset = ",offset
+        if(doTraceOffset)
+                offset = 2^ct
+                //multiply by current offset (>=1)
+                Wave inten = $("aveint_"+num2str(ct))
+                inten *= offset
+                //      Print "new offset = ",offset
+        endif
         
         ct +=1