lable the tau0 output as tau_INITIAL for easier analysis later

MUSIC-fluid · Sep 30, 2024 · 8d0d729 · 8d0d729
1 parent 85eb226
commit 8d0d729
Show file tree

Hide file tree

Showing 3 changed files with 24 additions and 7 deletions.
diff --git a/src/evolve.cpp b/src/evolve.cpp
@@ -148,7 +148,8 @@ int Evolve::EvolveIt(
 
         if (it == iFreezeStart || it == iFreezeStart + 10
             || it == iFreezeStart + 30 || it == iFreezeStart + 50) {
-            grid_info.output_momentum_anisotropy_vs_etas(tau, *fpCurr);
+            grid_info.output_momentum_anisotropy_vs_etas(
+                it - iFreezeStart, tau, *fpCurr);
         }
 
         if (DATA.beastMode == 0) {

diff --git a/src/grid_info.cpp b/src/grid_info.cpp
@@ -1942,31 +1942,47 @@ void Cell_info::output_average_phase_diagram_trajectory(
 //! This function outputs system's eccentricity and momentum anisotropy
 //! as functions of eta_s
 void Cell_info::output_momentum_anisotropy_vs_etas(
-    const double tau, Fields &arena) const {
+    const int itau, const double tau, Fields &arena) const {
     ostringstream filename;
-    filename << "momentum_anisotropy_tau_" << tau << ".dat";
+    if (itau == 0) {
+        filename << "momentum_anisotropy_tau_INITIAL.dat";
+    } else {
+        filename << "momentum_anisotropy_tau_" << tau << ".dat";
+    }
     std::fstream of;
     of.open(filename.str().c_str(), std::fstream::out);
     of << "# eta_s  epsilon_p(ideal)(cos)  epsilon_p(ideal)(sin)  "
        << "epsilon_p(shear)(cos)  epsilon_p(shear)(sin)  "
        << "epsilon_p(full)(cos)  epsilon_p(full)(sin)  " << endl;
 
     ostringstream filename1;
-    filename1 << "eccentricities_evo_ed_tau_" << tau << ".dat";
+    if (itau == 0) {
+        filename1 << "eccentricities_evo_ed_tau_INITIAL.dat";
+    } else {
+        filename1 << "eccentricities_evo_ed_tau_" << tau << ".dat";
+    }
     std::fstream of1;
     of1.open(filename1.str().c_str(), std::fstream::out);
     of1 << "# eta_s  ed(GeV/fm^3)  ecc_n(cos)  ecc_n(sin) (n=1-6)  "
         << "x_o(ed)(fm)  y_o(ed)(fm)" << endl;
 
     ostringstream filename2;
-    filename2 << "eccentricities_evo_nB_tau_" << tau << ".dat";
+    if (itau == 0) {
+        filename2 << "eccentricities_evo_nB_tau_INITIAL.dat";
+    } else {
+        filename2 << "eccentricities_evo_nB_tau_" << tau << ".dat";
+    }
     std::fstream of2;
     of2.open(filename2.str().c_str(), std::fstream::out);
     of2 << "# eta_s  nB(1/fm^3)  ecc_n(cos)  ecc_n(sin) (n=1-6)  "
         << "x_o(nB)(fm)  y_o(nB)(fm)" << endl;
 
     ostringstream filename3;
-    filename3 << "meanpT_estimators_tau_" << tau << ".dat";
+    if (itau == 0) {
+        filename3 << "meanpT_estimators_tau_INITIAL.dat";
+    } else {
+        filename3 << "meanpT_estimators_tau_" << tau << ".dat";
+    }
     std::fstream of3;
     of3.open(filename3.str().c_str(), std::fstream::out);
     of3 << "# eta_s  dS/deta_s  dE/deta_s (GeV) [s] (1/fm^-3)  [r^2] (fm^2)"

diff --git a/src/grid_info.h b/src/grid_info.h
@@ -136,7 +136,7 @@ class Cell_info {
     //! This function outputs system's eccentricity and momentum anisotropy
     //! as functions of eta_s
     void output_momentum_anisotropy_vs_etas(
-        const double tau, Fields &arena) const;
+        const int itau, const double tau, Fields &arena) const;
 
     //! This function outputs hydro evolution file into memory for JETSCAPE
     void OutputEvolutionDataXYEta_memory(