From 81da12aa7036052ba12bcf28467d2b5d5ef15a24 Mon Sep 17 00:00:00 2001
From: Klaus Rabbertz <klaus.rabbertz@cern.ch>
Date: Tue, 13 Apr 2021 09:07:12 +0200
Subject: [PATCH] Add kQuadraticSum option in rootout & yodaout
---
v2.3/toolkit/src/fnlo-tk-rootout.cc | 8 ++++++--
v2.3/toolkit/src/fnlo-tk-yodaout.cc | 8 ++++++--
2 files changed, 12 insertions(+), 4 deletions(-)
diff --git a/v2.3/toolkit/src/fnlo-tk-rootout.cc b/v2.3/toolkit/src/fnlo-tk-rootout.cc
index 8c7a8f7d..f74747f8 100644
--- a/v2.3/toolkit/src/fnlo-tk-rootout.cc
+++ b/v2.3/toolkit/src/fnlo-tk-rootout.cc
@@ -116,8 +116,9 @@ int main(int argc, char** argv) {
man << " Default: \"kScale1\", i.e. mur=muf=scale1," << endl;
man << " Alternatives: \"kScale2\", i.e. mur=muf=scale2," << endl;
man << " \"scale12\", i.e. mur=scale1, muf=scale2," << endl;
- man << " \"scale21\", i.e. mur=scale2, muf=scale1." << endl;
- man << " \"kProd\", i.e. mur=muf=scale1*scale2." << endl;
+ man << " \"scale21\", i.e. mur=scale2, muf=scale1," << endl;
+ man << " \"kProd\", i.e. mur=muf=scale1*scale2," << endl;
+ man << " \"kQuadraticSum\", i.e. mur=muf=sqrt(scale1^2+scale2^2)." << endl;
yell << " #" << endl;
man << "Use \"_\" to skip changing a default argument." << endl;
yell << " #" << endl;
@@ -365,6 +366,9 @@ int main(int argc, char** argv) {
} else if ( chflex == "kProd" ) {
fnlo->SetMuFFunctionalForm(kProd);
fnlo->SetMuRFunctionalForm(kProd);
+ } else if ( chflex == "kQuadraticSum" ) {
+ fnlo->SetMuFFunctionalForm(kQuadraticSum);
+ fnlo->SetMuRFunctionalForm(kQuadraticSum);
} else {
error["fnlo-tk-rootout"] << "Unknown scale choice " << chflex << ", aborted!" << endl;
}
diff --git a/v2.3/toolkit/src/fnlo-tk-yodaout.cc b/v2.3/toolkit/src/fnlo-tk-yodaout.cc
index d80c0a85..e0a6399a 100644
--- a/v2.3/toolkit/src/fnlo-tk-yodaout.cc
+++ b/v2.3/toolkit/src/fnlo-tk-yodaout.cc
@@ -119,8 +119,9 @@ int main(int argc, char** argv) {
man << " Default: \"kScale1\", i.e. mur=muf=scale1," << endl;
man << " Alternatives: \"kScale2\", i.e. mur=muf=scale2," << endl;
man << " \"scale12\", i.e. mur=scale1, muf=scale2," << endl;
- man << " \"scale21\", i.e. mur=scale2, muf=scale1." << endl;
- man << " \"kProd\", i.e. mur=muf=scale1*scale2." << endl;
+ man << " \"scale21\", i.e. mur=scale2, muf=scale1," << endl;
+ man << " \"kProd\", i.e. mur=muf=scale1*scale2," << endl;
+ man << " \"kQuadraticSum\", i.e. mur=muf=sqrt(scale1^2+scale2^2)." << endl;
man << "[np]: Apply nonperturbative corrections if available, def. = no." << endl;
man << " Alternatives: \"yes\" or \"np\"" << endl;
yell << " #" << endl;
@@ -399,6 +400,9 @@ int main(int argc, char** argv) {
} else if ( chflex == "kProd" ) {
fnlo->SetMuFFunctionalForm(kProd);
fnlo->SetMuRFunctionalForm(kProd);
+ } else if ( chflex == "kQuadraticSum" ) {
+ fnlo->SetMuFFunctionalForm(kQuadraticSum);
+ fnlo->SetMuRFunctionalForm(kQuadraticSum);
} else {
error["fnlo-tk-yodaout"] << "Unknown scale choice " << chflex << ", aborted!" << endl;
}
--
GitLab