Changeset 950


Ignore:
Timestamp:
May 27, 2020 1:15:29 PM (2 years ago)
Author:
ilavsky
Message:

Minor fix

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/User Procedures/Irena/IR3_SimpleFits.ipf

    r942 r950  
    10751075        W_coef[0]=Guinier_I0    //G
    10761076        W_coef[1]=Guinier_Rg            //Rg
    1077         T_Constraints[0] = {"K1 > "+num2str(Guinier_Rg/10)}
    1078         T_Constraints[1] = {"K0 > 0"}
    1079         LocalEwave[0]=(Guinier_I0/20)
    1080         LocalEwave[1]=(Guinier_Rg/20)
     1077        T_Constraints[0] = {"K0 > 0"}
     1078        T_Constraints[1] = {"K1 > "+num2str(Guinier_Rg/10)}
     1079        LocalEwave[0]=(Guinier_I0/50)
     1080        LocalEwave[1]=(Guinier_Rg/50)
    10811081        variable QminFit, QmaxFit
    10821082        QminFit = CursorAXWave[DataQstartPoint]
     
    10991099        endswitch
    11001100        if (V_FitError!=0)      //there was error in fitting
    1101                 RemoveFromGraph $("fit_"+NameOfWave(CursorAWave))
     1101                RemoveFromGraph/W=IR3J_LogLogDataDisplay /Z $("fit_"+NameOfWave(CursorAWave))
    11021102                beep
    11031103                Abort "Fitting error, check starting parameters and fitting limits"
    11041104        endif
    11051105        Wave W_sigma
     1106        W_coef =  abs(W_coef)
    11061107        string TagText, TagTextLin
    11071108        AchievedChiSquare = V_chisq/(DataQEndPoint-DataQstartPoint)
     
    11831184        w[0]=abs(w[0])
    11841185        w[1]=abs(w[1])
    1185         return w[0]/q * exp(-0.5 * q^2*w[1]^2)
     1186        return w[0]*exp(-q^2 * w[1]^2/2)/q
    11861187End
    11871188//**********************************************************************************************************
     
    16121613                        break           
    16131614                case "Guinier Sheet":                           // Guinier Sheet
    1614                         ModelLogLogInt = Guinier_I0 *exp(-ModelLogLogQ[p]^2*Guinier_Rg^2)/ModelLogLogQ^2
     1615                        ModelLogLogInt = Guinier_I0 *exp(-ModelLogLogQ[p]^2*Guinier_Rg^2)*ModelLogLogQ^(-2)
    16151616                        if(UsingLinearizedModel)
    16161617                                ModelLinLinLogInt = ln(ModelLogLogInt*ModelLogLogQ^2)   
Note: See TracChangeset for help on using the changeset viewer.