Docker-in-Docker (DinD) capabilities of public runners deactivated. More info

Commit 740b2faf authored by Grégoire Uhlrich's avatar Grégoire Uhlrich
Browse files

Fixed orderExternalFermions options

parent c244dd03
...@@ -200,35 +200,99 @@ public: ...@@ -200,35 +200,99 @@ public:
); );
void rotateFields( void rotateFields(
std::initializer_list<std::string> fields, std::vector<std::string> const &fields,
std::initializer_list<std::string> newFields, std::vector<std::string> const &newFields,
std::vector<std::vector<csl::Expr>> const &rotation, std::vector<std::vector<csl::Expr>> const &rotation,
bool diagonalizeMasses = false, bool diagonalizeMasses = false,
int nMassLessFields = 0 int nMassLessFields = 0
); );
void rotateFields( void rotateFields(
std::initializer_list<std::string> fields1, std::vector<std::string> const &fields1,
std::initializer_list<std::string> newFields1, std::vector<std::string> const &newFields1,
std::vector<std::vector<csl::Expr>> const &rotation1, std::vector<std::vector<csl::Expr>> const &rotation1,
std::initializer_list<std::string> fields2, std::vector<std::string> const &fields2,
std::initializer_list<std::string> newFields2, std::vector<std::string> const &newFields2,
std::vector<std::vector<csl::Expr>> const &rotation2, std::vector<std::vector<csl::Expr>> const &rotation2,
int nMassLessFields = 0 int nMassLessFields = 0
); );
void rotateFields( void rotateFields(
std::initializer_list<std::string> fields, std::initializer_list<std::string> const &fields,
bool diagonalizeMasses = false, std::initializer_list<std::string> const &newFields,
int nMassLessFields = 0 std::vector<std::vector<csl::Expr>> const &rotation,
bool diagonalizeMasses = false,
int nMassLessFields = 0
)
{
return rotateFields(
std::vector<std::string>(fields),
std::vector<std::string>(newFields),
rotation,
diagonalizeMasses,
nMassLessFields
);
}
void rotateFields(
std::initializer_list<std::string> const &fields1,
std::initializer_list<std::string> const &newFields1,
std::vector<std::vector<csl::Expr>> const &rotation1,
std::initializer_list<std::string> const &fields2,
std::initializer_list<std::string> const &newFields2,
std::vector<std::vector<csl::Expr>> const &rotation2,
int nMassLessFields = 0
)
{
return rotateFields(
std::vector<std::string>(fields1),
std::vector<std::string>(newFields1),
rotation1,
std::vector<std::string>(fields2),
std::vector<std::string>(newFields2),
rotation2,
nMassLessFields
);
}
void rotateFields(
std::vector<std::string> const &fields,
bool diagonalizeMasses = false,
int nMassLessFields = 0
); );
void birotateFields( void birotateFields(
std::initializer_list<std::string> fields1, std::vector<std::string> const &fields1,
std::initializer_list<std::string> fields2, std::vector<std::string> const &fields2,
int nMassLessFields = 0 int nMassLessFields = 0
); );
void rotateFields(
std::initializer_list<std::string> const &fields,
bool diagonalizeMasses = false,
int nMassLessFields = 0
)
{
return rotateFields(
std::vector<std::string>(fields),
diagonalizeMasses,
nMassLessFields
);
}
void birotateFields(
std::initializer_list<std::string> const &fields1,
std::initializer_list<std::string> const &fields2,
int nMassLessFields = 0
)
{
return birotateFields(
std::vector<std::string>(fields1),
std::vector<std::string>(fields2),
nMassLessFields
);
}
void applyUnitaryCondition( void applyUnitaryCondition(
std::vector<std::vector<csl::Expr>> const& unitary std::vector<std::vector<csl::Expr>> const& unitary
); );
......
...@@ -39,7 +39,7 @@ namespace mty::option { ...@@ -39,7 +39,7 @@ namespace mty::option {
/////////////////////////////////////////////////// ///////////////////////////////////////////////////
inline bool simplifyAmplitudes = true; inline bool simplifyAmplitudes = true;
inline bool orderExternalFermions = true; inline bool orderExternalFermions = false;
inline bool discardLowerOrders = true; inline bool discardLowerOrders = true;
inline bool evaluateFermionTraces = true; inline bool evaluateFermionTraces = true;
inline bool excludeTadpoles = true; inline bool excludeTadpoles = true;
......
...@@ -537,6 +537,7 @@ WilsonSet Model::computeWilsonCoefficients( ...@@ -537,6 +537,7 @@ WilsonSet Model::computeWilsonCoefficients(
FeynOptions feynOptions FeynOptions feynOptions
) )
{ {
feynOptions.orderExternalFermions = true;
if (order == TreeLevel) { if (order == TreeLevel) {
return computeWilsonCoefficients_default( return computeWilsonCoefficients_default(
order, insertions, feynOptions order, insertions, feynOptions
......
...@@ -735,8 +735,8 @@ void ModelBuilder::birotateFields( ...@@ -735,8 +735,8 @@ void ModelBuilder::birotateFields(
} }
void ModelBuilder::rotateFields( void ModelBuilder::rotateFields(
std::initializer_list<std::string> fields, std::vector<std::string> const &fields,
std::initializer_list<std::string> newFields, std::vector<std::string> const &newFields,
std::vector<std::vector<csl::Expr>> const &rotation, std::vector<std::vector<csl::Expr>> const &rotation,
bool diagonalizeMasses, bool diagonalizeMasses,
int nMassLessFields int nMassLessFields
...@@ -763,11 +763,11 @@ void ModelBuilder::rotateFields( ...@@ -763,11 +763,11 @@ void ModelBuilder::rotateFields(
} }
void ModelBuilder::rotateFields( void ModelBuilder::rotateFields(
std::initializer_list<std::string> fields1, std::vector<std::string> const &fields1,
std::initializer_list<std::string> newFields1, std::vector<std::string> const &newFields1,
std::vector<std::vector<csl::Expr>> const &rotation1, std::vector<std::vector<csl::Expr>> const &rotation1,
std::initializer_list<std::string> fields2, std::vector<std::string> const &fields2,
std::initializer_list<std::string> newFields2, std::vector<std::string> const &newFields2,
std::vector<std::vector<csl::Expr>> const &rotation2, std::vector<std::vector<csl::Expr>> const &rotation2,
int nMassLessFields int nMassLessFields
) )
...@@ -798,9 +798,9 @@ void ModelBuilder::rotateFields( ...@@ -798,9 +798,9 @@ void ModelBuilder::rotateFields(
} }
void ModelBuilder::rotateFields( void ModelBuilder::rotateFields(
std::initializer_list<std::string> fields, std::vector<std::string> const &fields,
bool diagonalizeMasses, bool diagonalizeMasses,
int nMassLessFields int nMassLessFields
) )
{ {
std::vector<mty::Particle> fields_p; std::vector<mty::Particle> fields_p;
...@@ -812,9 +812,9 @@ void ModelBuilder::rotateFields( ...@@ -812,9 +812,9 @@ void ModelBuilder::rotateFields(
} }
void ModelBuilder::birotateFields( void ModelBuilder::birotateFields(
std::initializer_list<std::string> fields1, std::vector<std::string> const &fields1,
std::initializer_list<std::string> fields2, std::vector<std::string> const &fields2,
int nMassLessFields int nMassLessFields
) )
{ {
std::vector<mty::Particle> fields1_p; std::vector<mty::Particle> fields1_p;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment