ideal hydro - nwE0.0001T0.34Tf0.15

oheader = "hydro3_2/final/nwE0.0001T0.34Tf0.15/"

<< /int/users/paulrom/hydro3/runs/getHBT_new.m;

\(\(massarr = 
        Flatten[Import["/int/users/paulrom/hydro3_2/decays/masses.dat"]];\)\)

\(\(namesarr = 
        Flatten[Import["/int/users/paulrom/hydro3_2/decays/names.dat"]];\)\)

Tf = 0.15

Rside

ra = {}; 
    Do[cc = Table[{qq, idCside[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 
            0.05, 0.005}]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRside.dat", %]

Rout

ra = {}; 
    Do[cc = Table[{qq, idCout[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 
            0.05, 0.005}]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRout.dat", %]

Rlong

ra = {}; 
    Do[cc = Table[{qq, idClong[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 
            0.05, 0.005}]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRlong.dat", %]

Transpose[{%65[[All, 
            1]], %65[[All, 2]]/%63[[All, 2]]}]

Export[header <> "results/HBTRoORs.dat", %]

viscous hydro - C3.0E0.08T0.34Tf0.18

oheader = "hydro3_2/final/C3.0E0.08T0.34Tf0.18/"

<< /int/users/paulrom/hydro3/runs/getHBT_new.m;

\(\(massarr = 
        Flatten[Import["/int/users/paulrom/hydro3_2/decays/masses.dat"]];\)\)

\(\(namesarr = 
        Flatten[Import["/int/users/paulrom/hydro3_2/decays/names.dat"]];\)\)

Tf = 0.18

Rside

ra = {}; 
    Do[cc = Table[{qq, Cside[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 0.1, 
            0.01}]; 
ListPlot[cc, PlotJoined -> True]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRside.dat", %]

Rout

ra = {}; 
    Do[cc = Table[{qq, Cout[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 0.05, 
            0.005}]; ListPlot[cc, PlotJoined -> True]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRout.dat", %]

Rlong

ra = {}; 
    Do[cc = Table[{qq, Clong[i/10. , qq, Tf, massarr[[1]]]}, {qq, 0, 0.1, 
            0.005}]; ListPlot[cc, PlotJoined -> True]; 
      r = x /. FindRoot[
            Normal[Fit[cc, {1, x, x^2, x^3, x^4}, x]] == 0.5, {x, 0, 
              0.05}]; 
Print[i/10. , "  :  ", 
        r,  " radius: ", 
        Sqrt[Log[2]]/r 0.1973]; 
ra = 
        Append[ra, {i/10. , Sqrt[Log[2]]/r 0.1973}], {i, 1, 12}];  
    Print[ra]

Export[header <> "results/HBTRlong.dat", %]

Transpose[{%75[[All, 
            1]], %75[[All, 2]]/%73[[All, 2]]}]

Export[header <> "results/HBTRoORs.dat", %]