Changeset 824 for sans/Dev/trunk/NCNR_User_Procedures/Analysis/Models/NewModels_2010/Two_Yukawa_v40.ipf
 Timestamp:
 Jul 20, 2011 9:40:32 AM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

sans/Dev/trunk/NCNR_User_Procedures/Analysis/Models/NewModels_2010/Two_Yukawa_v40.ipf
r773 r824 61 61 xwave_1yuk = alog(log(qmin) + x*((log(qmax)log(qmin))/num)) 62 62 Make/O/D coef_1yuk = {0.1,50,1,10} 63 make/o/t parameters_1yuk = {"volume fraction","Radius (A)","scale, K"," charge, Z"}63 make/o/t parameters_1yuk = {"volume fraction","Radius (A)","scale, K","Decay constant, Z"} 64 64 Edit parameters_1yuk,coef_1yuk 65 65 Variable/G root:g_1yuk … … 113 113 xwave_2yuk = alog(log(qmin) + x*((log(qmax)log(qmin))/num)) 114 114 Make/O/D coef_2yuk = {0.2,50,6,10,1,2} 115 make/o/t parameters_2yuk = {"volume fraction","Radius (A)","scale, K1"," charge, Z1","scale, K2","charge, Z2"}115 make/o/t parameters_2yuk = {"volume fraction","Radius (A)","scale, K1","Decay constant, Z1","scale, K2","Decay constant, Z2"} 116 116 Edit parameters_2yuk,coef_2yuk 117 117 Variable/G root:g_2yuk … … 162 162 End 163 163 164 // with the regular 2yukawa plotted, this uses the coefficients to plot g(r) 165 // 166 //  no dependency is created, it would just slow things down. So you'll 167 // need to rerun this every time. 168 // 169 // gr is scaled to dimensionless distance, r/diameter 170 // 171 Macro Plot_2Yukawa_Gr() 172 //if the regular 2yukawa procedure is already plotted 173 //  then append it to thte graph yourself 174 Duplicate/O ywave_2yuk ywave_2yuk_Igor 175 176 fTwoYukawa(coef_2yuk,xwave_2yuk,ywave_2yuk_Igor) 177 178 DoWindow/F Gr_plot 179 if(V_flag==0) 180 Display gr 181 DoWindow/C Gr_plot 182 Modifygraph log=0 183 SetAxis bottom 0,10 184 Modifygraph lsize=2 185 Label left "g(r)";DelayUpdate 186 Label bottom "dimensionless distance (r/diameter)" 187 legend 188 endif 189 190 191 192 End 193 194 195 164 196 // 165 197 Function fTwoYukawa(cw,xw,yw) : FitFunc … … 190 222 yw = SqTwoYukawa(xw*radius*2, Z1, Z2, K1, K2, phi, a, b, c1, c2, d1, d2) 191 223 // printf("%g %g\n",q,sq) 224 else 225 yw = 1000 //return a really bogus answer, as Yun suggests 192 226 endif 193 227 endif … … 1602 1636 gr[0] = 0 1603 1637 1638 SetScale/P x,0,dr, gr 1639 1604 1640 // Killwaves/Z temp 1605 1641
Note: See TracChangeset
for help on using the changeset viewer.