@@ -246,7 +246,7 @@ void HydroSourceTATB::get_hydro_energy_source(
246246 }
247247 */
248248 double y_CM = atanh ((TA - TB) / (TA + TB + Util::small_eps) * tanhYbeam_);
249- // double y_L = yL_frac_ * y_CM;
249+ // double y_L = yL_frac_ * y_CM;
250250 double y_L = compute_yL (TA, TB, y_CM, eta0_, sigma_eta_, eta_m_);
251251
252252 /*
@@ -263,7 +263,10 @@ void HydroSourceTATB::get_hydro_energy_source(
263263 double eta_envelop =
264264 eta_profile_plateau(eta_s, eta0, DATA_.eta_fall_off);
265265 */
266- double M = Util::m_N * std::sqrt (TA * TA + TB * TB + 2.0 * TA * TB * std::cosh (2.0 * ybeam_))/(Util::hbarc);
266+ double M =
267+ Util::m_N
268+ * std::sqrt (TA * TA + TB * TB + 2.0 * TA * TB * std::cosh (2.0 * ybeam_))
269+ / (Util::hbarc);
267270 /*
268271 double E_norm =
269272 tau_source
@@ -274,12 +277,14 @@ void HydroSourceTATB::get_hydro_energy_source(
274277 // double E_norm = tau_source*energy_eta_profile_normalisation_numerical(
275278 // y_CM, eta0, DATA_.eta_fall_off);
276279 // double epsilon = M_inv * eta_envelop / E_norm / dtau; // [1/fm^5]
277- double tilted_epsilon = eta_profile_plateau (eta_s, eta0_, sigma_eta_) * pow (TA, f_plus) * pow (TB, f_minus);
280+ double tilted_epsilon = eta_profile_plateau (eta_s, eta0_, sigma_eta_)
281+ * pow (TA, f_plus) * pow (TB, f_minus);
278282 double shifted_epsilon =
279283 eta_profile_plateau (eta_s - (y_CM - y_L), eta0_, sigma_eta_);
280284 double tilted_norm = energy_eta_profile_normalisation_tilted (
281- TA, TB, eta0_, eta_m_, sigma_eta_, y_CM, M, y_L)/tau_source;
282- double shifted_norm = M / C_eta_/tau_source;
285+ TA, TB, eta0_, eta_m_, sigma_eta_, y_CM, M, y_L)
286+ / tau_source;
287+ double shifted_norm = M / C_eta_ / tau_source;
283288
284289 double epsilon = (beta_ * tilted_epsilon * tilted_norm
285290 + (1 . - beta_) * shifted_epsilon * shifted_norm)
0 commit comments