Skip to content
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

Add mass over temperature ratio check output #175

Merged
merged 4 commits into from
Dec 13, 2024
Merged

Add mass over temperature ratio check output #175

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b1381d3
add mass over temperature ratio check output
lisabiermann Nov 4, 2024
3bd4f9d
Rename tmp in a meaningful way
lisabiermann Dec 11, 2024
6aa71c9
Merge branch 'develop' into highT_validity_check
lisabiermann Dec 11, 2024
3ac2b2b
Add max ratio output and unit test
lisabiermann Dec 13, 2024
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
13 changes: 11 additions & 2 deletions include/BSMPT/transition_tracer/transition_tracer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -181,14 +181,23 @@ class TransitionTracer

/**
* @brief Store the list of bounce solutions
*
*/
std::vector<BounceSolution> ListBounceSolution;

/**
* @brief output data storage
*/
output output_store;

/**
* @brief CheckMassRatio
* @param input
* @param vev
* @param temp
*/
void CheckMassRatio(const user_input &input,
const std::vector<double> &vec,
const double &temp) const;
};

} // namespace BSMPT
} // namespace BSMPT
67 changes: 66 additions & 1 deletion src/transition_tracer/transition_tracer.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,9 @@
new_transition_data.crit_false_vev =
pair.false_phase.Get(pair.crit_temp).point;

CheckMassRatio(
input, new_transition_data.crit_false_vev, pair.crit_temp);

BounceSolution bounce(input.modelPointer,
mintracer,
pair,
Expand Down Expand Up @@ -172,6 +175,10 @@
.Get(new_transition_data.nucl_approx_temp.value_or(
EmptyValue))
.point;

CheckMassRatio(input,
new_transition_data.nucl_approx_false_vev,
bounce.GetNucleationTempApprox());
}
else
{
Expand All @@ -195,6 +202,9 @@
pair.false_phase
.Get(new_transition_data.nucl_temp.value_or(EmptyValue))
.point;
CheckMassRatio(input,
new_transition_data.nucl_false_vev,
bounce.GetNucleationTemp());
}
else
{
Expand All @@ -218,6 +228,9 @@
pair.false_phase
.Get(new_transition_data.perc_temp.value_or(EmptyValue))
.point;
CheckMassRatio(input,
new_transition_data.perc_false_vev,
bounce.GetPercolationTemp());
}
else
{
Expand All @@ -243,6 +256,9 @@
.Get(
new_transition_data.compl_temp.value_or(EmptyValue))
.point;
CheckMassRatio(input,
new_transition_data.compl_false_vev,
bounce.GetCompletionTemp());
}
else
{
Expand Down Expand Up @@ -489,4 +505,53 @@
{
}

} // namespace BSMPT
void TransitionTracer::CheckMassRatio(const user_input &input,
const std::vector<double> &vec,
const double &temp) const
{
std::stringstream ss;
std::vector<double> massOverTempSq, tmp;
vollous marked this conversation as resolved.
Show resolved Hide resolved
massOverTempSq = input.modelPointer->HiggsMassesSquared(
input.modelPointer->MinimizeOrderVEV(vec), temp) /
std::pow(temp, 2);
tmp = input.modelPointer->GaugeMassesSquared(
input.modelPointer->MinimizeOrderVEV(vec), temp) /
std::pow(temp, 2);

massOverTempSq.insert(massOverTempSq.end(), tmp.begin(), tmp.end());

int color = 0;
for (auto el : massOverTempSq)
{
if (el > 0.25) // m/T > 0.5
{
color = 1;
if (el > 1) // m/T > 1.0
{
color = 2;
break;
}
}
}

if (color == 0)
{
ss << "\n\033[1;92mm^2(vev_false, T = " << std::to_string(temp)
<< ") / T^2 = " << vec_to_string(massOverTempSq) << "\033[0m\n";

Check warning on line 540 in src/transition_tracer/transition_tracer.cpp

View check run for this annotation

Codecov / codecov/patch

src/transition_tracer/transition_tracer.cpp#L539-L540 

Added lines #L539 - L540 were not covered by tests
}
else if (color == 1)
{
ss << "\n\033[1;93mm^2(vev_false, T = " << std::to_string(temp)
<< ") / T^2 = " << vec_to_string(massOverTempSq) << "\033[0m\n";
}
else
{
ss << "\n\033[1;91mm^2(vev_false, T = " << std::to_string(temp)
<< ") / T^2 = " << vec_to_string(massOverTempSq) << "\033[0m\n";

Check warning on line 550 in src/transition_tracer/transition_tracer.cpp

View check run for this annotation

Codecov / codecov/patch

src/transition_tracer/transition_tracer.cpp#L549-L550 

Added lines #L549 - L550 were not covered by tests
}

Logger::Write(LoggingLevel::TransitionDetailed, ss.str());
return;
}

} // namespace BSMPT