SetDirectory[NotebookDirectory[]];

Import FINESSE calculation

run1=Import["lcgt2drun1.dat","Data"];
run1fit=Fit[run1,Table[f^n,{n,0,10}],f];

P1 = ListLogLogPlot[run1, Joined -> True]

P2 = LogLogPlot[run1fit, {f, 1, 10^4}, PlotStyle -> Hue[.4]]

Show[P1, P2]

run2=Import["lcgt2drun2.dat","Data"];
run2fit=Fit[run2,Table[f^n,{n,0,10}],f];

run3=Import["lcgt2drun3.dat","Data"];
run3fit=Fit[run3,Table[f^n,{n,0,10}],f];

run4=Import["lcgt2drun4.dat","Data"];
run4fit=Fit[run4,Table[f^n,{n,0,10}],f];

run5=Import["lcgt2drun5.dat","Data"];
run5fit=Fit[run5,Table[f^n,{n,0,10}],f];

Other noise

hb=1.05457266912510183 (10^(-34));(* reduced Planck constant (J*s) *)
ω0=1.8 (10^15);
c=2.99792458000000117`*^8;(* speed of light (m/s) *)
Ω=2 π f;
kB=1.38*10^(-23);(* Boltzmann's constant (J/K) *)
g=9.8;(* acceleration of the gravity (m/s^2) *)
Gravity=6.673*10^(-11);(* Gravitational constant (m^3/Kg/s^2) *)
solarmass=1.989*10^(30);(* mass of Sun (Kg) *)
Mpc=365.2422*24*60*60*c*3.26*10^6;(* 1 Mpc (m) *)

m=30;
L=3000;
I0=825;
λ=1064 (10^(-9));
T=0.004;
ρ=0.92;

F=((2π)/T);
γ=((T c)/(4L));

Irse=((I0 (1-(ρ^2)))/(1+(ρ^2)+2ρ Cos[2ArcTan[Ω/γ]]));

hRP=((16 F)/(m L (Ω^2)))(Sqrt[((hb Irse)/(π c λ))((γ^2)/((γ^2)+(Ω^2)))]);

hshot=(Sqrt[((hb λ)/(4π c Irse))((γ^2)+(Ω^2))]);

hseis=Sqrt[(2*(10^(-17))*(f^(-6.5)))^2+(2*(10^(-16))*(f^(-8)))^2];

hmirrorcoaETM=Import["hmirrorcoaETMsap.m"];
hmirrorcoaITM=Import["hmirrorcoaITMsap.m"];
hmirrorstrETM=Import["hmirrorstrETMsap.m"];
hmirrorstrITM=Import["hmirrorstrITMsap.m"];
hmirrorthermoETM=Import["hmirrorthermoETMsap.m"];
hmirrorthermoITM=Import["hmirrorthermoITMsap.m"];
hsuspETM=Import["hsuspETMsap2.m"];
hsuspITM=Import["hsuspITMsap.m"];

hmirror=Sqrt[hmirrorstrITM^2+hmirrorthermoITM^2+hmirrorstrETM^2+hmirrorthermoETM^2+hmirrorcoaITM^2+hmirrorcoaETM^2];
hsusp=Sqrt[hsuspITM^2+hsuspETM^2];

LogLogPlot[{hseis,hsusp,hmirror,Sqrt[(hRP^2)+(((run1fit/L*Sqrt[2]))^2)],Sqrt[(hRP^2)+(hshot^2)]},{f,10,1000},PlotRange->{5*(10^(-25)),(10^(-21))}]

FINESSE shot noise is wrong by sqrt[2] so we should mupltiply rt2 to runfit.

Inspiral range

m1=1.4*solarmass;(* mass of a compact object (Kg) *)
m2=1.4*solarmass;(* mass of other compace object (Kg) *)
mtotal=m1+m2;(* total moss of binary system (Kg) *)
mreduce=m1*m2/mtotal;(* reduced mass of binary system (Kg) *)
fmin=20;(* low cut off frequency of integral (Hz) *)
fmax=c^3/(6^(3/2)*Pi*Gravity*mtotal);(* high cut off frequency of integral (Hz) *)
SNR=8;(* threshold of Signal to Noise Ratio *)

Nini=100;

integ1=(∑_(ff=1)^Nini)(1/(f^(7/3)*((run1fit/L*Sqrt[2])^2+hRP^2+hseis^2+hmirror^2+hsusp^2))*(fmax-fmin)/Nini/.f->fmin+(fmax-fmin)*(ff-0.5)/Nini);

integ2=(∑_(ff=1)^Nini)(1/(f^(7/3)*((run2fit/L*Sqrt[2])^2+hRP^2+hseis^2+hmirror^2+hsusp^2))*(fmax-fmin)/Nini/.f->fmin+(fmax-fmin)*(ff-0.5)/Nini);

integ3=(∑_(ff=1)^Nini)(1/(f^(7/3)*((run3fit/L*Sqrt[2])^2+hRP^2+hseis^2+hmirror^2+hsusp^2))*(fmax-fmin)/Nini/.f->fmin+(fmax-fmin)*(ff-0.5)/Nini);

integ4=(∑_(ff=1)^Nini)(1/(f^(7/3)*((run4fit/L*Sqrt[2])^2+hRP^2+hseis^2+hmirror^2+hsusp^2))*(fmax-fmin)/Nini/.f->fmin+(fmax-fmin)*(ff-0.5)/Nini);

integ5=(∑_(ff=1)^Nini)(1/(f^(7/3)*((run5fit/L*Sqrt[2])^2+hRP^2+hseis^2+hmirror^2+hsusp^2))*(fmax-fmin)/Nini/.f->fmin+(fmax-fmin)*(ff-0.5)/Nini);

Range1=(1/SNR)/Sqrt[2]*2*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ1]/Mpc;

Range2=(1/SNR)/Sqrt[2]*2*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ2]/Mpc;

Range3=(1/SNR)/Sqrt[2]*2*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ3]/Mpc;

Range4=(1/SNR)/Sqrt[2]*2*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ4]/Mpc;

Range5=(1/SNR)/Sqrt[2]*2*(Gravity/c^3)^(5/6)*c*(5/96*mreduce)^(1/2)*(mtotal/Pi^2)^(1/3)*Sqrt[4*integ5]/Mpc;

Result

{Range1,Range2,Range3,Range4,Range5}

Quit[]