Changeset 148 for sans/SASCalc

Timestamp:

Aug 15, 2007 2:33:02 PM (16 years ago)

Author:

srkline

Message:

Tried, in at least an approximate way to account for lenses using SASCALC

• added the correct resolution calculation as in the reduction macros

--

• flux is not really correct as the lens transmission has not been included
• the lenses can be 'used' in unphysical cases, may need to be more stringent

File:

• sans/SASCalc/trunk/SASCALC.ipf (modified) (13 diffs)

sans/SASCalc/trunk/SASCALC.ipf

@@ -10 +10 @@
 // 19 MAR 2007 SRK - corrections added for projected BS diameter at anode plane
 // 11 APR 2007 SRK - default aperture offset of 5 cm added to match VAX implementation
+// nn AUG 2007 SRK - added defulat sample aperture size
+//                                               added option of lenses, approximated beamDiam=sourceDiam, BSdiam=1"
+//                                               Lens flux, trans is not corrected for lens/prism transmission
+//                                               Lenses can still be inserted in incorrect cases, and are not automatically taken out
 //
 // calculate what q-values you get based on the instruments settings
@@ -89 +93 @@
         String/G root:SAS:gApPopStr = "1/16\";1/8\";3/16\";1/4\";5/16\";3/8\";7/16\";1/2\";9/16\";5/8\";11/16\";3/4\";other;"
         Variable/G root:SAS:gSamApOther = 10            //non-standard aperture diameter, in mm
+        Variable/G root:SAS:gUsingLenses = 0            //0=no lenses, 1=lenses(or prisms)
         
         //tick labels for SDD slider
@@ -278 +283 @@
         Slider SC_Slider_2,pos={394,21},size={47,65},proc=OffsetSliderProc,live=0,ticks=4
         Slider SC_Slider_2,limits={0,25,1},variable= root:SAS:gOffset,thumbColor= (1,16019,65535)
-        CheckBox checkNG3,pos={20,16},size={36,14},proc=SelectInstrumentCheckProc,title="NG3"
+        CheckBox checkNG3,pos={20,18},size={36,14},proc=SelectInstrumentCheckProc,title="NG3"
         CheckBox checkNG3,value=1,mode=1
-        CheckBox checkNG7,pos={66,16},size={36,14},proc=SelectInstrumentCheckProc,title="NG7"
+        CheckBox checkNG7,pos={66,18},size={36,14},proc=SelectInstrumentCheckProc,title="NG7"
         CheckBox checkNG7,value=0,mode=1
         CheckBox checkChamber,pos={172,48},size={57,14},proc=TableCheckProc,title="Chamber"
@@ -307 +312 @@
         SetVariable setvar0_3,pos={140,94},size={110,15},title="Diam (mm)",disable=1
         SetVariable setvar0_3,limits={0,100,0.1},value= root:SAS:gSamApOther,proc=SourceApOtherSetVarProc
+        
+        CheckBox checkLens,pos={6,155},size={44,14},proc=LensCheckProc,title="Lenses?"
+        CheckBox checkLens,value=root:SAS:gUsingLenses
+        
         // set up a fake dependency to trigger recalculation
         //root:SAS:gCalculate := ReCalculateInten(root:SAS:gTouched)
@@ -392 +401 @@
         
         if(event %& 0x1)        // bit 0, value set
+                if(sliderValue != 0)
+                        LensCheckProc("",0)             //make sure lenses are deselected
+                endif
                 sourceToSampleDist()            //updates the SSD global and wave
                 //change the sourceAp popup, SDD range, etc
@@ -408 +420 @@
 
         if(event %& 0x1)        // bit 0, value set
+                if(sliderValue < 1300)
+                        LensCheckProc("",0)             //make sure lenses are deselected
+                endif
                 sampleToDetectorDist()          //changes the SDD and wave (DetDist is the global)
                 ReCalculateInten(1)
@@ -471 +486 @@
         ReCalculateInten(1)
 End
+
+
+//lenses (or prisms) in/out changes resolution
+Function LensCheckProc(ctrlName,checked) : CheckBoxControl
+        String ctrlName
+        Variable checked
+
+        // don't let the box get checked if the conditions are wrong
+        // lambda != 8.09,8.4,17.2
+        // Ng != 0
+        NVAR lens = root:SAS:gUsingLenses
+        NVAR ng = root:SAS:gNg
+        NVAR lam = root:SAS:gLambda
+        NVAR dist = root:SAS:gDetDist
+        
+        if( (Ng != 0) || (lam < 8) || (dist < 1300) )
+                lens = 0
+                CheckBox checkLens,value=0
+                return(0)
+        endif
+        lens = checked
+        
+        ReCalculateInten(1)
+End
+
 
 // change the source aperture
@@ -1075 +1115 @@
         String results
         results ="Failure"
+        
+        NVAR usingLenses = root:SAS:gUsingLenses
 
         //rename for working variables,  these must be gotten from global
@@ -1103 +1145 @@
 
         v_lambda = lambdaWidth^2/6.0
-        v_b = 0.25*(S1*L2/L1)^2 +0.25*(S2*L2/lp)^2
+        
+        if(usingLenses==1)                      //SRK 2007
+                v_b = 0.25*(S1*L2/L1)^2 +0.25*(2/3)*(lambdaWidth/lambda)^2*(S2*L2/lp)^2         //correction to 2nd term
+        else
+                v_b = 0.25*(S1*L2/L1)^2 +0.25*(S2*L2/lp)^2              //original form
+        endif
+        
         v_d = (DDet/2.3548)^2 + del_r^2/12.0
         vz = vz_1 / lambda
         yg = 0.5*g*L2*(L1+L2)/vz^2
-        v_g = 2.0*yg^2*v_lambda
+        v_g = 2.0*(2.0*yg^2*v_lambda)                                   //factor of 2 correction, B. Hammouda, 2007
 
         r0 = L2*tan(2.0*asin(lambda*inQ/(4.0*Pi) ))
@@ -1212 +1260 @@
         NVAR instrument = instrument
         NVAR L2diff = L2diff
+        NVAR lens = root:SAS:gUsingLenses
         
         
@@ -1263 +1312 @@
         sprintf temp,"Sample Aperture to Sample Position =\t%.2f cm\r", L2Diff
         str += temp
+        if(lens==1)
+                sprintf temp,"Lenses are IN\r"
+        else
+                sprintf temp,"Lenses are OUT\r"
+        endif
+        str += temp 
         
    setDataFolder root:
@@ -1395 +1450 @@
 
 //direction = one of "vertical;horizontal;maximum;"
+// all of this is bypassed if the lenses are in
+//
 Function beamDiameter(direction)
         String direction
 
+        NVAR lens = root:SAS:gUsingLenses
+        if(lens)
+                return sourceApertureDiam()
+        endif
+        
     Variable l1 = sourceToSampleDist()
     Variable l2 //= sampleAperToDetDist()
@@ -1436 +1498 @@
 //on NG3 and NG7, allowable sizes are 1,2,3,4" diameter
 //will return values larger than 4.0*2.54 if a larger beam is needed
+//
+// - in an approximate way, account for lenses
 Function beamstopDiam()
 
-    Variable bm = beamDiameter("maximum")
-    Variable bs=0.0
-    
-    do
-        bs += 1
-    while ( (bs*2.54 < bm) || (bs > 30.0))                      //30 = ridiculous limit to avoid inf loop
- 
+        NVAR yesLens = root:SAS:gUsingLenses
+        Variable bm=0
+        Variable bs=0.0
+   
+        if(yesLens)
+                //bm = sourceApertureDiam()             //ideal result, not needed
+                bs = 1                                                          //force the diameter to 1"
+        else
+                bm = beamDiameter("maximum")
+                do
+                bs += 1
+           while ( (bs*2.54 < bm) || (bs > 30.0))                       //30 = ridiculous limit to avoid inf loop
+        endif
+
         //update the wave
         WAVE rw=root:SAS:realsRead