Initial data, common for all operations

Задаем рабочую директорию, устанавливаем опции для графики, задаем модель для подгонки контуров (Лорентцевским контуром).

MyPath="D:/work/HoFe3(B03)4";
SetDirectory[MyPath];
SetOptions[{Plot,ListPlot,ListLinePlot},ImageSize->600];
<<PhysicalConstants`
FileData=FileNames["*.txt"];
iName=FileData;

Pearson[i_,v_,w_,m_]:=((2 (Sqrt[-1+(2^(1/m))]) i (((1+((4 (-1+(2^(1/m))) (((-v+x))^2))/(w^2))))^(-m)) LogGamma[m])/((Sqrt[π]) w LogGamma[-(1/2)+m]));
PsdVoigt1[i_,v_,w_,m_]:=i(m(2/π)(w/(4(((x-v))^2)+(w^2)))+(1-m)(((Sqrt[4Log[2]])/((Sqrt[π])w))) (E^(-(4((Log[2])/(w^2)))(((x-v))^2))));
PsdVoigt2[i_,v_,wL_,wG_,m_]:=i(m(2/π)(wL/(4(((x-v))^2)+(wL^2)))+(1-m)((Sqrt[4Log[2]])/((Sqrt[π])wG)) (E^(-4((Log[2])/(wG^2))(((x-v))^2))));
PsdVoigt3[i_,v_,w_,m_]:=i(0.63662m(w/(4(((x-v))^2)+(w^2)))+(1-m)(((1.75415×(10^(-14)))/w)) (E^(-(4(0.6931471805599453/(w^2)))(((x-v))^2))));
Voigt[i_,v_,wL_,wG_]:=((2i Log[2] wL)/((π^(3/2))(wG^2)))(∫_(-∞)^∞)((E^(-(t^2)))/((((((Sqrt[Log[2]])(wL/wG)))^2)+((((Sqrt[4Log[2]])(((x-v))/wG)-t))^2))))\[DifferentialD]t;
Voigt1[i_,v_,wL_,wG_]:=N[((0.248i wL)/(wG^2))(∫_(-∞)^∞)((E^(-(t^2)))/(((((0.8325(wL/wG)))^2)+(((1.665(((x-v))/wG)-t))^2))))\[DifferentialD]t];
Lorentz[i_,v_,w_]:=(2/π)((i w)/(4(((x- v))^2) +(((w ))^2)));
Lorentz2[i_,v_,w_]:=(i/((((x- v))^2) +(((w/2 ))^2)));
PearsonVII[i_,v_,w_,m_]:=(i/((1+4(((((x-v)/w)))^2)((((2^(1/m))-1))^m))));
GausLorCross[i_,v_,w_,m_]:=(i/(1+m(((((x-v)/w)))^2)(E^((((1-m))/2)(((((x-v)/w)))^2)))));
nbes[ω_,T_]:=N [(1/((E^(( 100 PlanckConstant ω SpeedOfLight  )/(BoltzmannConstant T Meter Kelvin)))-1))]+1;
nbeas[ω_,T_]:=N[(1/((E^((100 PlanckConstant  ω SpeedOfLight )/(BoltzmannConstant T Meter Kelvin )))-1))];
Harmonic[i_,v_,w_]:= (2 i v^2 w x)/((v^2-x^2)^2+4 w^2 x^2);

Off[General::stop];

FIT

Sample with Harmonic

MyPath="D:/work/HoFe3(B03)4";
SetDirectory[MyPath];
FileData=FileNames["*.txt"];
iName=FileData;
(* формируем форму вывода в файл *)
ResultData={{T,i1,v1,w1,i2,v2,w2}};
(* Задаем модель для подгонки*)
If[x<0,model=nbeas[x,T] (Harmonic[i1,v1,w1]+Harmonic[i2,v2,w2])+b,
model=nbes[x,T] (Harmonic[i1,v1,w1]+Harmonic[i2,v2,w2])+c x+b];
fit={i1->1400,v1->74,w1->4,i2->220,v2->92,w2->3,c->0,b->0};
Sfrom=0;
Sto=130;
Szero=0;

Do[{Print[iName],
T=ToExpression[StringSplit[iName,{"_","K"}][[1]]];
Print[T];
FullData=ReadList[iName,Number,RecordLists->True];
Data=Select[FullData, #[[1]]>Sfrom &&#[[1]] <Sto &];

fit=FindFit[Data,{model},{
{i1,i1/.fit},{v1,v1/.fit},{w1,w1/.fit},
{i2,i2/.fit},{v2,v2/.fit},{w2,w2/.fit},
{c,c/.fit},{b,b/.fit}},{x},
{MaxIterations->5000}];

Print[Plot[{model/.fit,
If[x<0,{nbeas[x,T] Harmonic[i1,v1,w1]+b/.fit,
nbeas[x,T] Harmonic[i2,v2,w2]+b/.fit},{
nbes[x,T] Harmonic[i1,v1,w1]+b/.fit,
nbes[x,T] Harmonic[i2,v2,w2]+b/.fit}]
},{x,-Sto,Sto},
Epilog->Point[Data],
AxesOrigin->{Sfrom,Szero},
PlotRange->{0,Max[Data[[1;;,2]]]} ,
PlotStyle->{Axes->True,Thickness[0.001]},AxesLabel->{"x","y"}]];

Eva=Evaluate[model/.fit]/.{x->Data[[1;;,1]]};
Diff=Eva-Data[[1;;,2]];
Print[ListLinePlot[Diff,PlotRange->{-1000,1000},AxesOrigin->{0,0},Filling->Axis]];
tmp={T,i1,v1,w1,i2,v2,w2}/.fit;
ResultData=Append[ResultData,tmp];
Print[tmp];
},
{iName,FileData[[1;;100]]}]
TableForm[ResultData]

Sample with Lorentz

MyPath="D:/work/HoFe3(B03)4";
SetDirectory[MyPath];
FileData=FileNames["*.txt"];
iName=FileData;
(* формируем форму вывода в файл *)
ResultData={{T,i1,v1,w1,i2,v2,w2}};
(* Задаем модель для подгонки*)
model=Lorentz[i1,v1,w1]+Lorentz[i2,v2,w2]+c x+b;
fit={i1->1400,v1->74,w1->4,i2->220,v2->92,w2->3,c->0,b->0};
Sfrom=0;
Sto=130;
Szero=0;

Do[{Print[iName],
FullData=ReadList[iName,Number,RecordLists->True];
T=ToExpression[StringSplit[iName,{"_","K","."}][[1]]];
Print[T];
BoseFactor=Array[0&,{Length[FullData],2}];
Print[T];
	Eva1=(1/(Exp[((100 PlanckConstant  x SpeedOfLight )/(BoltzmannConstant T Meter Kelvin ))]-1)+1)/.{x->FullData[[1;;,1]]};
	Diff1=FullData[[1;;,2]]/Eva1;
	BoseFactor[[1;;,1]]=FullData[[1;;,1]];
	BoseFactor[[1;;,2]]=Diff1;
Data=Select[BoseFactor, #[[1]]>Sfrom &&#[[1]] <Sto &];

fit=FindFit[Data,{model},{
{i1,i1/.fit},{v1,v1/.fit},{w1,w1/.fit},
{i2,i2/.fit},{v2,v2/.fit},{w2,w2/.fit},
{c,c/.fit},{b,b/.fit}},{x},
{MaxIterations->5000}];

Print[Plot[{model/.fit,
{Lorentz[i1,v1,w1]+c/.fit+b/.fit,
Lorentz[i2,v2,w2]+c/.fit+b/.fit}
},{x,-Sto,Sto},
Epilog->Point[Data],
AxesOrigin->{Sfrom,Szero},
PlotRange->{0,Max[Data[[1;;,2]]]} ,
PlotStyle->{Axes->True,Thickness[0.001]},AxesLabel->{"x","y"}]];

Eva=Evaluate[model/.fit]/.{x->Data[[1;;,1]]};
Diff=Eva-Data[[1;;,2]];
Print[ListLinePlot[Diff,PlotRange->{-1000,1000},AxesOrigin->{0,0},Filling->Axis]];
tmp={T,i1,v1,w1,i2,v2,w2}/.fit;
ResultData=Append[ResultData,tmp];
Print[tmp];
},
{iName,FileData[[1;;100]]}]
TableForm[ResultData]