Skip to content

Modifying Energy and Baryon density functions to get the correct normalization #28

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
chunshen1987 merged 3 commits into from
Feb 24, 2026
Merged

Modifying Energy and Baryon density functions to get the correct normalization #28

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
255811b
Fixed the baryon density by dividing by datu
Afrid165 Feb 19, 2026
6b15123
Energy and Baryon density changed to get correct normalziation'
Afrid165 Feb 23, 2026
a6ea254
formatting
chunshen1987 Feb 24, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
19 changes: 11 additions & 8 deletions src/hydro_source_TATB.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -264,8 +264,9 @@ void HydroSourceTATB::get_hydro_energy_source(
eta_profile_plateau(eta_s, eta0, DATA_.eta_fall_off);
*/
double M =
TA * TA + TB * TB + 2 * TA * TB * Util::m_N * cosh2Ybeam_; // [1/fm^4]

Util::m_N
* std::sqrt(TA * TA + TB * TB + 2.0 * TA * TB * std::cosh(2.0 * ybeam_))
/ (Util::hbarc);
/*
double E_norm =
tau_source
Expand All @@ -276,13 +277,14 @@ void HydroSourceTATB::get_hydro_energy_source(
// double E_norm = tau_source*energy_eta_profile_normalisation_numerical(
// y_CM, eta0, DATA_.eta_fall_off);
// double epsilon = M_inv * eta_envelop / E_norm / dtau; // [1/fm^5]
double tilted_epsilon = pow(TA, f_plus) * pow(TB, f_minus);
double tilted_epsilon = eta_profile_plateau(eta_s, eta0_, sigma_eta_)
* pow(TA, f_plus) * pow(TB, f_minus);
double shifted_epsilon =
eta_profile_plateau(eta_s - (y_CM - y_L), eta0_, sigma_eta_);
double tilted_norm = tau_source
* energy_eta_profile_normalisation_tilted(
TA, TB, eta0_, eta_m_, sigma_eta_, y_CM, M, y_L);
double shifted_norm = tau_source * M / C_eta_;
double tilted_norm = energy_eta_profile_normalisation_tilted(
TA, TB, eta0_, eta_m_, sigma_eta_, y_CM, M, y_L)
/ tau_source;
double shifted_norm = M / C_eta_ / tau_source;

double epsilon = (beta_ * tilted_epsilon * tilted_norm
+ (1. - beta_) * shifted_epsilon * shifted_norm)
Expand Down Expand Up @@ -323,6 +325,7 @@ double HydroSourceTATB::get_hydro_rhob_source(
*/
res = norm_B * (1 - omega) * (TA * eta_rhob_plus + TB * eta_rhob_minus)
+ norm_B_prime * omega * TA * TB * (eta_rhob_plus + eta_rhob_minus);
res /= gridDtau_;
return (res);
}
/*
Expand All @@ -343,7 +346,7 @@ double HydroSourceTATB::eta_rhob_left_factor(const double eta) const {
*/

double HydroSourceTATB::eta_rhob_left_factor(const double eta) const {
double eta_0_nB = -std::abs(DATA_.eta_rhob_0);
double eta_0_nB = std::abs(DATA_.eta_rhob_0);
double sigma_B_plus = DATA_.eta_rhob_width_1;
double sigma_B_minus = DATA_.eta_rhob_width_2;

Expand Down