Changeset 800 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/NCNR_Utils.ipf

Timestamp:

May 4, 2011 9:47:33 AM (11 years ago)

Author:

srkline

Message:

Changed the 2D resolution calculation to include gravity. This seems to work, but will still require some testing. It requires no modifications to the smearing calculation, still using parallel and perpendicular directions.

Added a Fake Update feature to the RealTime? update. There are specific, separate instructions for how to set this up. Usef (possibly) for summer schools or demos.

Adjusted the 2D MonteCarlo? simulation to a corrected beam center. The Gauss Peak was not symmetric around the old beam center.

Corrected the AAO resolution smearing functions to work with cursors.

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/NCNR_Utils.ipf (modified) (2 diffs)

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

@@ -307 +307 @@
         
         // this is really sigma_Q_parallel
-        SigmaQX = kap*kap/12* (3*(S1/L1)^2 + 3*(S2/LP)^2 + (proj_DDet/L2)^2 + (sin(phi))^2*8*(a_val_L2)^2*lambda^4*lambdaWidth^2)
+        SigmaQX = kap*kap/12 * (3*(S1/L1)^2 + 3*(S2/LP)^2 + (proj_DDet/L2)^2 + (sin(phi))^2*8*(a_val_L2)^2*lambda^4*lambdaWidth^2)
         SigmaQX += inQ*inQ*v_lambda
 
@@ -313 +313 @@
         proj_DDet = DDet*sin(phi) + DDet*cos(phi)               //not necessary, since DDet is the same in both X and Y directions
 
-        SigmaQY = 3*(S1/L1)^2 + 3*(S2/LP)^2 + (proj_DDet/L2)^2 + (cos(phi))^2*8*(a_val_L2)^2*lambda^4*lambdaWidth^2
-        SigmaQY *= kap*kap/12
+        SigmaQY = kap*kap/12 * (3*(S1/L1)^2 + 3*(S2/LP)^2 + (proj_DDet/L2)^2 + (cos(phi))^2*8*(a_val_L2)^2*lambda^4*lambdaWidth^2)
         
         SigmaQX = sqrt(SigmaQX)
         SigmaQy = sqrt(SigmaQY)
         
+
+/////////////////////////////////////////////////
+///     
+        ////// this is all experimental, inclusion of gravity effect into the parallel component
+        //
+        Variable yg_d,acc,sdd,ssd,lambda0,DL_L,sig_l
+        Variable var_qlx,var_qly,var_ql,qx,qy,sig_perp,sig_para
+//      G = 981.  //!   ACCELERATION OF GRAVITY, CM/SEC^2
+        acc = vz_1              //      3.956E5 //!     CONVERT WAVELENGTH TO VELOCITY CM/SEC
+        SDD = L2                //1317
+        SSD = L1                //1627          //cm
+        lambda0 = lambda                //              15
+        DL_L = lambdaWidth              //0.236
+        SIG_L = DL_L/sqrt(6)
+        YG_d = -0.5*G*SDD*(SSD+SDD)*(LAMBDA0/acc)^2
+//      Print "DISTANCE BEAM FALLS DUE TO GRAVITY (CM) = ",YG
+        
+        sig_perp = kap*kap/12 * (3*(S1/L1)^2 + 3*(S2/LP)^2 + (proj_DDet/L2)^2)
+        sig_perp = sqrt(sig_perp)
+        
+        
+        FindQxQy(inQ,phi,qx,qy)
+        
+        VAR_QLY = SIG_L^2 * (QY+4*PI*YG_d/(SDD*LAMBDA0))^2
+        VAR_QLX = (SIG_L*QX)^2
+        VAR_QL = VAR_QLY + VAR_QLX  //! WAVELENGTH CONTRIBUTION TO VARIANCE
+        sig_para = (sig_perp^2 + VAR_QL)^0.5
+        
+//return values PBR     
+        SigmaQX = sig_para
+        SigmaQy = sig_perp
+        
+//      
         
         results = "success"