Ignore:
Timestamp:
Sep 10, 2009 3:45:30 PM (13 years ago)
Author:
srkline
Message:

in SASCALC - changed some of the logic to allow 2D simulation

in MultiScatter_MC:

  • changed results wave to report coherent multiple scattering counters correctly (also changed XOP)
  • added save of 1D sim data

in U_CALC:

  • adjusted display for optimal view on windows
  • 7 angle ranges now (added defaults, propagated where all controls are modified)
  • added the angle axis correctly now as a free axis, with a hook linked to the bottom axis
  • removed XS and trans calculations, as integration of ran_dev isn't well

in USANS_EmptyWaves:

  • double precision waves
File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r550 r551  
    1515// 
    1616// - more presets 
    17 // - fast ways to increase/decrease number of points 
    18 // - fast ways to increase/decrease counting time 
    19 // 
    20 // - plot as countrate, not absolute scale 
    21 // - 3e-5 cutoff 
    22 // - ? don't plot lowest angle range (but needs to be in the count time) 
     17 
     18// 
     19// X plot as countrate, not absolute scale 
     20// X 3e-5 cutoff 
     21// X ? don't plot lowest angle range (but needs to be in the count time) 
    2322// - need empty beam and empty cell count rate vs. aperture (Cd vs. Gd?) 
    2423// 
    … …  
    2827// or at least compare to what the empty cell would give in the same count time 
    2928// 
    30 // model the direct beam?? currently the "red" region from -1 to 0.6 is almost entirely 
     29// -? model the direct beam?? currently the "red" region from -1 to 0.6 is almost entirely 
    3130// the primary beam, so it's a bit artificial (qmin is really ~ 3e-5) 
    3231// 
    3332 
    3433// 
    35 // Need T_wide, T_rock, I peak for proper absolute scaling 
     34// Need T_wide, T_rock, I peak for proper absolute scaling, but I don't know if it's really important 
     35// to be able to simlulate to this extent 
    3636// 
    3737 
    … …  
    103103        Variable/G gCtTime6 = 2000 
    104104        Variable/G gIncr6 = 2    
     105         
     106        Variable/G gAngLow7 = 50 
     107        Variable/G gAngHigh7 = 95 
     108        Variable/G gNumPts7 = 10 
     109        Variable/G gCtTime7 = 0 
     110        Variable/G gIncr7 = 5    
    105111         
    106112        // results, setup values 
    … …  
    287293        Variable pos=strsearch(platform,"WINDOWS",0) 
    288294        if(pos >= 0)            //windows 
    289                 Display /W=(55,44,955,544) /K=1 
     295                Display /W=(55,44,855,450) /K=1 
    290296        else            //mac 
    291297                Display /W=(55,44,1055,544) /K=1 
    … …  
    297303        ControlBar/L 500 
    298304         
    299         GroupBox group0,pos={5,1},size={493,159},title="Instrument Setup" 
    300         GroupBox group1,pos={5,165},size={240,147},title="Sample Setup" 
    301         GroupBox group2,pos={5,325},size={259,147},title="Results" 
     305        GroupBox group0,pos={5,1},size={493,177},title="Instrument Setup" 
     306        GroupBox group1,pos={5,183},size={240,147},title="Sample Setup" 
     307        GroupBox group2,pos={5,343},size={259,147},title="Results" 
    302308         
    303309        PopupMenu popup0,pos={17,19},size={165,20},title="Sample Aperture Diam (in)" 
    … …  
    316322        SetVariable setvar1d,pos={left+299,top},size={100,15},title="# points",value= gNumPts1 
    317323        SetVariable setvar1e,pos={left+399,top},size={93,15},title="count (s)",value= gCtTime1 
    318         SetVariable setvar1a,labelBack=(65535,32768,32768) 
     324//      SetVariable setvar1a,labelBack=(65535,32768,32768)              //old rainbow 
     325        SetVariable setvar1a,labelBack=(49858,65535,65535) 
    319326         
    320327        pt += inc 
    … …  
    324331        SetVariable setvar2d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts2 
    325332        SetVariable setvar2e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime2 
    326         SetVariable setvar2a labelBack=(65535,65533,32768) 
     333//      SetVariable setvar2a labelBack=(65535,65533,32768)              //old rainbow 
     334        SetVariable setvar2a labelBack=(21074,8995,21074) 
    327335         
    328336        pt += inc 
    … …  
    332340        SetVariable setvar3d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts3 
    333341        SetVariable setvar3e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime3 
    334         SetVariable setvar3a labelBack=(32769,65535,32768) 
     342//      SetVariable setvar3a labelBack=(32769,65535,32768)              //old rainbow 
     343        SetVariable setvar3a labelBack=(0,60652,60652) 
    335344         
    336345        pt += inc 
    … …  
    340349        SetVariable setvar4d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts4 
    341350        SetVariable setvar4e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime4 
    342         SetVariable setvar4a labelBack=(32768,65535,65535) 
     351//      SetVariable setvar4a labelBack=(32768,65535,65535)              //old rainbow 
     352        SetVariable setvar4a labelBack=(0,51400,0) 
    343353         
    344354        pt += inc 
    … …  
    348358        SetVariable setvar5d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts5 
    349359        SetVariable setvar5e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime5 
    350         SetVariable setvar5a labelBack=(32768,54615,65535) 
     360//      SetVariable setvar5a labelBack=(32768,54615,65535)              //old rainbow 
     361        SetVariable setvar5a labelBack=(59367,49344,0) 
    351362         
    352363        pt += inc 
    … …  
    356367        SetVariable setvar6d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts6 
    357368        SetVariable setvar6e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime6 
    358         SetVariable setvar6a labelBack=(44253,29492,58982) 
     369//      SetVariable setvar6a labelBack=(44253,29492,58982)              //old rainbow 
     370        SetVariable setvar6a labelBack=(54998,0,0) 
     371 
     372        pt += inc 
     373        SetVariable setvar7a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow7 
     374        SetVariable setvar7b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh7 
     375        SetVariable setvar7c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr7 
     376        SetVariable setvar7d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts7 
     377        SetVariable setvar7e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime7 
     378//      SetVariable setvar7a labelBack=(44253,29492,58982)              //old rainbow 
     379        SetVariable setvar7a labelBack=(39321,21845,51657) 
    359380 
    360381// the action procedures and limits/increments 
    … …  
    365386        SetVariable setvar5a proc=ThetaMinSetVarProc 
    366387        SetVariable setvar6a proc=ThetaMinSetVarProc 
     388        SetVariable setvar7a proc=ThetaMinSetVarProc 
    367389 
    368390// 
    … …  
    373395        SetVariable setvar5b proc=ThetaMaxSetVarProc 
    374396        SetVariable setvar6b proc=ThetaMaxSetVarProc 
     397        SetVariable setvar7b proc=ThetaMaxSetVarProc 
    375398// 
    376399        SetVariable setvar1c proc=IncrSetVarProc,limits={0.01,0.1,0.01} 
    … …  
    380403        SetVariable setvar5c proc=IncrSetVarProc,limits={0.5,5,1} 
    381404        SetVariable setvar6c proc=IncrSetVarProc,limits={1,10,2} 
     405        SetVariable setvar7c proc=IncrSetVarProc,limits={1,20,2} 
    382406// 
    383407        SetVariable setvar1d proc=NumPtsSetVarProc,limits={2,50,1} 
    … …  
    387411        SetVariable setvar5d proc=NumPtsSetVarProc,limits={2,50,1} 
    388412        SetVariable setvar6d proc=NumPtsSetVarProc,limits={2,50,1} 
     413        SetVariable setvar7d proc=NumPtsSetVarProc,limits={2,50,1} 
    389414// 
    390415        SetVariable setvar1e proc=CtTimeSetVarProc,limits={-1,50000,1} 
    … …  
    394419        SetVariable setvar5e proc=CtTimeSetVarProc,limits={-1,50000,100} 
    395420        SetVariable setvar6e proc=CtTimeSetVarProc,limits={-1,50000,100} 
    396          
    397         Button button0,pos={left+280,180},size={60,20},fColor=(65535,65535,0),proc=U_SimPlotButtonProc,title="Plot" 
    398         CheckBox check0_2,pos={left+280,230},size={60,14},title="CountRate?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR 
    399         CheckBox check0_3,pos={left+280,250},size={60,14},title="Noise?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_AddNoise 
    400         CheckBox check0_4 title="Show EMP?",pos={left+280,270},proc=ShowEMPCheckProc,value=0 
     421        SetVariable setvar7e proc=CtTimeSetVarProc,limits={-1,50000,100} 
     422         
     423        Button button0,pos={left+280,200},size={130,20},fColor=(65535,65535,0),proc=U_SimPlotButtonProc,title="Simulate USANS" 
     424        CheckBox check0_2,pos={left+280,250},size={60,14},title="CountRate?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR 
     425        CheckBox check0_3,pos={left+280,270},size={60,14},title="Noise?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_AddNoise 
     426        CheckBox check0_4 title="Show EMP?",pos={left+280,290},proc=ShowEMPCheckProc,value=0 
    401427 
    402428 
    … …  
    405431        ThetaEditMode(2)                //checked on startup 
    406432         
    407 //      instrument setup 
    408         SetVariable U_setvar0_1,pos={left+20,211},size={160,15},title="Thickness (cm)" 
     433//      sample setup 
     434        SetVariable U_setvar0_1,pos={left+20,231},size={160,15},title="Thickness (cm)" 
    409435        SetVariable U_setvar0_1,limits={0,inf,0.1},value= root:Packages:NIST:USANS:Globals:U_Sim:gThick  
    410         SetVariable U_setvar0_3,pos={left+20,235},size={160,15},title="Sample Transmission" 
     436        SetVariable U_setvar0_3,pos={left+20,255},size={160,15},title="Sample Transmission" 
    411437        SetVariable U_setvar0_3,limits={0,1,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:gSamTrans 
    412         PopupMenu U_popup0,pos={left+20,185},size={165,20},proc=Sim_USANS_ModelPopMenuProc,title="Model" 
     438        PopupMenu U_popup0,pos={left+20,205},size={165,20},proc=Sim_USANS_ModelPopMenuProc,title="Model" 
    413439        PopupMenu U_popup0,mode=1,value= #"U_FunctionPopupList()" 
    414         SetVariable setvar0,pos={left+20,259},size={120,15},title="Empty Level",disable=2 
     440        SetVariable setvar0,pos={left+20,279},size={120,15},title="Empty Level",disable=2 
    415441        SetVariable setvar0,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_EmptyLevel 
    416         SetVariable setvar0_1,pos={left+20,284},size={120,15},title="Bkg Level",disable=2 
     442        SetVariable setvar0_1,pos={left+20,304},size={120,15},title="Bkg Level",disable=2 
    417443        SetVariable setvar0_1,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_BkgLevel 
    418444         
    … …  
    420446         
    421447// a box for the results 
    422         SetVariable totalTime,pos={left+20,350},size={150,15},title="Count time (h:m)",value= gTotTimeStr 
    423 //      ValDisplay valdisp0,pos={338,210},size={220,13},title="Total detector counts" 
    424 //      ValDisplay valdisp0,limits={0,0,0},barmisc={0,1000},value= root:Packages:NIST:USANS:Globals:U_Sim:g_1DTotCts 
    425         ValDisplay valdisp0_2,pos={left+20,380},size={220,13},title="Fraction of beam scattered" 
     448        SetVariable totalTime,pos={left+20,370},size={150,15},title="Count time (h:m)",value= gTotTimeStr 
     449        ValDisplay valdisp0_2,pos={left+20,395},size={220,13},title="Fraction of beam scattered" 
    426450        ValDisplay valdisp0_2,limits={0,0,0},barmisc={0,1000},value= root:Packages:NIST:USANS:Globals:U_Sim:g_1DFracScatt 
    427         ValDisplay valdisp0_3,pos={left+20,410},size={220,13},title="Estimated transmission" 
     451        ValDisplay valdisp0_2,disable=1 
     452        ValDisplay valdisp0_3,pos={left+20,420},size={220,13},title="Estimated transmission" 
    428453        ValDisplay valdisp0_3,limits={0,0,0},barmisc={0,1000},value=root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans 
    429         Button button1,pos={left+20,440},size={50,20},proc=U_SaveButtonProc,title="Save" 
     454        ValDisplay valdisp0_3,disable=1 
     455        Button button1,pos={left+20,440},size={150,20},proc=U_SaveButtonProc,title="Save Simulated Data" 
    430456 
    431457         
    … …  
    550576End 
    551577 
    552 // hard-wired for 6 controls 
     578// hard-wired for 7 controls 
    553579// return value is in seconds 
    554580// global display string is set with hrs:min 
    … …  
    557583        Variable ii,num,totTime=0 
    558584        String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:" 
    559         num=6 
     585        num=7 
    560586         
    561587        for(ii=1;ii<=num;ii+=1) 
    … …  
    585611End 
    586612 
    587 // changes edit mode of the theta min/max boxes for simplified setup 
     613// changes edit mode of the theta min/max boxes and increment for simplified setup 
    588614// val = 2 = disable 
    589615// val = 0 = edit enabled 
    … …  
    593619        SetVariable setvar1a,win=UCALC,disable=val 
    594620        SetVariable setvar1b,win=UCALC,disable=val 
     621        SetVariable setvar1c,win=UCALC,disable=val 
    595622         
    596623        SetVariable setvar2a,win=UCALC,disable=val 
    597624        SetVariable setvar2b,win=UCALC,disable=val 
     625        SetVariable setvar2c,win=UCALC,disable=val 
    598626         
    599627        SetVariable setvar3a,win=UCALC,disable=val 
    600628        SetVariable setvar3b,win=UCALC,disable=val 
     629        SetVariable setvar3c,win=UCALC,disable=val 
    601630         
    602631        SetVariable setvar4a,win=UCALC,disable=val 
    603632        SetVariable setvar4b,win=UCALC,disable=val 
     633        SetVariable setvar4c,win=UCALC,disable=val 
    604634         
    605635        SetVariable setvar5a,win=UCALC,disable=val 
    606636        SetVariable setvar5b,win=UCALC,disable=val 
     637        SetVariable setvar5c,win=UCALC,disable=val 
    607638         
    608639        SetVariable setvar6a,win=UCALC,disable=val 
    609640        SetVariable setvar6b,win=UCALC,disable=val 
    610          
     641        SetVariable setvar6c,win=UCALC,disable=val 
     642         
     643        SetVariable setvar7a,win=UCALC,disable=val 
     644        SetVariable setvar7b,win=UCALC,disable=val 
     645        SetVariable setvar7c,win=UCALC,disable=val 
    611646        return(0) 
    612647End 
    … …  
    735770        SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder 
    736771         
    737         num = 6                 //# of angular ranges 
     772        num = 7                 //# of angular ranges 
    738773         
    739774        // only try to plot ranges with non-zero count times 
    … …  
    832867                Imon = GetUSANSBeamIntensity()                          //based on the aperture size, select the beam intensity 
    833868                Print "imon=",imon 
     869                 
    834870                // calculate the scattering cross section simply to be able to estimate the transmission 
    835                  
     871                // unfortunately, this calculation is useless in the USANS range. So although it's calculated, the 
     872                // results are never reported. 
     873                // -- the main issue is that the integration range is optimized for SANS, and is not useful for USANS 
     874                // -there are only a few points in the USANS range, and extending it lower caused issues before. 
    836875                CalculateRandomDeviate(funcUnsmeared,$coefStr,wavelength,"root:Packages:NIST:SAS:ran_dev",sig_sas) 
    837876                 
    … …  
    10691108// add SIM data to the graph if it exists and is not already on the graph 
    10701109// 
     1110// 
     1111// ** currently, I have changed the graph to not display the angle at the top. I have not yet 
     1112// found a way to keep the scaling of the two axes in-sync when the empty cell data is added to the 
     1113// graph (versus Q) 
     1114// 
    10711115Function GraphSIM() 
    10721116 
    … …  
    10901134 
    10911135// lines for the no-noise result vs angle 
    1092                 AppendToGraph/T Smeared_inten vs angle 
     1136//              AppendToGraph/T Smeared_inten vs angle 
     1137                AppendToGraph Smeared_inten vs Sim_USANS_q 
    10931138                ModifyGraph rgb(Smeared_inten)=(1,12815,52428) 
    10941139                ModifyGraph mode(Smeared_inten)=4,marker(Smeared_inten)=19,msize(Smeared_inten)=2 
    … …  
    10981143                AppendToGraph Sim_USANS_i vs Sim_USANS_q 
    10991144                ModifyGraph mode(Sim_USANS_i)=3,marker(Sim_USANS_i)=19,msize(Sim_USANS_i)=4 
    1100                 ModifyGraph zColor(Sim_USANS_i)={setNumber,1,6,Rainbow,0}                       //force the colors from 1->6 
     1145                //don't reverse the dbZ21 or the highest angle will have "invisible" light blue error bars 
     1146                ModifyGraph zColor(Sim_USANS_i)={setNumber,1,7,dBZ21,0}                         //better for 7 colors 
     1147//              ModifyGraph zColor(Sim_USANS_i)={setNumber,1,7,Rainbow,0}                       //force the colors from 1->7 
    11011148                ModifyGraph useMrkStrokeRGB(Sim_USANS_i)=1 
    11021149                ErrorBars/T=0 Sim_USANS_i Y,wave=(Sim_USANS_s,Sim_USANS_s) 
    … …  
    11041151                ModifyGraph log=1 
    11051152                ModifyGraph mirror(left)=1 
     1153                // if no top axis, then 
     1154//              ModifyGraph mirror(bottom)=1 
    11061155                ModifyGraph grid=2 
    11071156                ModifyGraph standoff=0 
    11081157                 
    11091158                // to make sure that the scales are the same (but fails on zoom) 
    1110 //              SetAxis bottom 2e-06,0.003 
    1111 //              SetAxis top 2e-06/5.55e-5,0.003/5.55e-5 
     1159                NewFreeAxis/O/T top_angle 
     1160                ModifyFreeAxis/Z top_angle,master= bottom,hook= TransformAngleAxisHook 
     1161                ModifyGraph lblPos(top_angle)=50,freePos(top_angle)=0 
     1162                ModifyGraph log(top_angle)=1 
    11121163                 
    11131164                SetDrawEnv linefgc= (39321,1,1),dash= 3,linethick= 3 
    11141165                SetDrawEnv xcoord= bottom 
    11151166                DrawLine 3e-05,0.01,3e-05,0.99 
    1116  
    1117                 Label top "Angle" 
     1167                 
     1168                Label top_angle "Angle" 
    11181169                Label bottom "Q (1/A)" 
    11191170                Label left "Counts or Count Rate" 
    … …  
    11251176End 
    11261177 
     1178Function TransformAngleAxisHook(s) 
     1179        STRUCT WMAxisHookStruct &s 
     1180 
     1181        s.max= s.max/5.55e-5 
     1182        s.min= s.min/5.55e-5 
     1183         
     1184        return 0 
     1185End 
    11271186 
    11281187//fakes a folder with loaded 1-d usans data, no calculation of the matrix 
    … …  
    12831342                        NVAR gIncr6 = gIncr6 
    12841343                         
    1285                          
     1344                        NVAR gAngLow7 = gAngLow7 
     1345                        NVAR gAngHigh7 = gAngHigh7 
     1346                        NVAR gNumPts7 = gNumPts7 
     1347                        NVAR gCtTime7 = gCtTime7 
     1348                        NVAR gIncr7 = gIncr7 
    12861349                         
    12871350                        strswitch(popStr)       // string switch 
    … …  
    13241387                                        gIncr6 = 2 
    13251388                                         
     1389                                        gAngLow7 = 50 
     1390                                        gAngHigh7 = 95 
     1391                                        gNumPts7 = 10 
     1392                                        gCtTime7 = 0 
     1393                                        gIncr7 = 5 
     1394                                         
    13261395                                        break                                            
    13271396                                case "Medium Count":     
    … …  
    13621431                                        gCtTime6 = 1200 
    13631432                                        gIncr6 = 2       
    1364                                                                          
     1433                                         
     1434                                        gAngLow7 = 50 
     1435                                        gAngHigh7 = 95 
     1436                                        gNumPts7 = 10 
     1437                                        gCtTime7 = 0 
     1438                                        gIncr7 = 5 
     1439                                                                 
    13651440                                        break 
    13661441                                case "Long Count":       
    … …  
    14021477                                        gIncr6 = 2 
    14031478                                         
     1479                                        gAngLow7 = 50 
     1480                                        gAngHigh7 = 95 
     1481                                        gNumPts7 = 10 
     1482                                        gCtTime7 = 3000 
     1483                                        gIncr7 = 5 
     1484                                         
    14041485                                        break 
    14051486                                default:                         
    … …  
    14531534                        String baseName="SIMUL" 
    14541535                         
    1455                         num = 6                 //# of angular ranges 
     1536                        num = 7                 //# of angular ranges 
    14561537                         
    14571538                        // only try to plot ranges with non-zero count times 
Note: See TracChangeset for help on using the changeset viewer.