Changeset 598 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/MultScatter_MonteCarlo_2D.ipf

Timestamp:

Nov 6, 2009 5:50:45 PM (13 years ago)

Author:

srkline

Message:

Procedures for mostly my use that enable me to (crudely) script the MC simulations so that they can be done in a big batch overnight.

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/MultScatter_MonteCarlo_2D.ipf (modified) (7 diffs)

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

@@ -1098 +1098 @@
         Endif
         
+        NVAR actSimTime = root:Packages:NIST:SAS:g_actSimTime
+        String str = ""
+        sprintf str,"%30s\t\t%g seconds\r","MonteCarlo Simulation time = ",actSimTime
+                
         Open refNum as fullpath+".txt"
                 wfprintf refNum, "%30s\t\t%g\r",results_desc,results
+                FBinWrite refNum,str
                 FStatus refNum
                 FSetPos refNum,V_logEOF
@@ -1239 +1244 @@
         t0 *= imon/1000/ThreadProcessorCount                    //projected time, in seconds (using threads for the calculation)
 
+        Print "Estimated Simulation time (s) = ",t0
+        
 // to correct for detector efficiency, send only the fraction of neutrons that are actually counted     
         NVAR detectorEff = root:Packages:NIST:SAS:g_detectorEff
+        NVAR actSimTime = root:Packages:NIST:SAS:g_actSimTime
+        NVAR SimTimeWarn = root:Packages:NIST:SAS:g_SimTimeWarn
 
         inputWave[0] = imon     * detectorEff                   //reset number of input neutrons before full simulation
         
-        if(t0>10)
+        if(t0>SimTimeWarn)
                 sprintf str,"The simulation will take approximately %d seconds.\r- Proceed?",t0
                 DoAlert 1,str
@@ -1268 +1277 @@
         t0 = (stopMSTimer(-2) - t0)*1e-6
         Printf  "MC sim time = %g seconds\r",t0
+        actSimTime = t0
         
         trans = results[8]                      //(n1-n2)/n1
@@ -1307 +1317 @@
         if(!rawCts)                     //go ahead and do the linear scaling
                 linear_data = linear_data / kappa
+                linear_data /= detectorEff
         endif           
         data = linear_data
@@ -1781 +1792 @@
         sigave[10,50;10] = 10*sigave[p]
 
-        // convert to absolute scale
+        // convert to absolute scale, remembering to un-correct for the detector efficiency
         if(doABS)
                 Variable kappa = thick*(pixSize/sdd)^2*trans*iMon
                 aveint /= kappa
                 sigave /= kappa
+                aveint /= detectorEff
+                sigave /= detectorEff
         endif
                                 
@@ -1837 +1850 @@
 //Function Script_2DMC()
 //
+//
+//      NVAR SimTimeWarn = root:Packages:NIST:SAS:g_SimTimeWarn
+//      SimTimeWarn = 36000                     //sets the threshold for the warning dialog to 10 hours
 //      STRUCT WMButtonAction ba
 //      ba.eventCode = 2                        //fake mouse click on button
 //      
 //      NVAR detDist = root:Packages:NIST:SAS:gDetDist
-//      
+//
 //      detDist = 200           //set directly in cm
 //      MC_DoItButtonProc(ba)
@@ -1854 +1870 @@
 //      SaveAsVAXButtonProc("",runIndex=107,simLabel="this is run 107, SDD = 400")
 //      
+//
+// SimTimeWarn = 10             //back to 10 seconds for manual operation
 //      return(0)
 //end