From bc93697f5bc55588a656322485b3deb59911e00b Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Fri, 25 Jun 2021 14:10:46 +0200 Subject: [PATCH 1/5] save --- higgsfodptN.cpp | 8 ++++---- runcards/Higgs-FO-as-N.yaml | 4 ++-- src/integration.cpp | 27 ++++++++++++++------------- src/partonic.cpp | 2 +- 4 files changed, 21 insertions(+), 20 deletions(-) diff --git a/higgsfodptN.cpp b/higgsfodptN.cpp index c3dc132..fea74da 100644 --- a/higgsfodptN.cpp +++ b/higgsfodptN.cpp @@ -145,18 +145,18 @@ int main(int argc, char* argv[]) { if (sectype == "hadronic") { results = higgshard.higgsdpt(pt,y1,y2); } else if (sectype == "partonic") { - results = higgspart.partonichiggsdpt(pt,nn); + results = higgspart.partonichiggsdpt(pt,std::pow(10,nn)); } else { std::cout << "Error in entry sectype!" << std::endl; exit(EXIT_FAILURE); } // Generate some output logs & write to output file - printf("N=%e: dHdpt = %e +- %e. \n", nn, results[0], results[1]); + printf("N=%e: dHdpt = %e +- %e. %e\n", std::pow(10,nn), results[0], results[1], results[2]); output_file.setf(std::ios_base::scientific); - output_file << nn << std::setw(space) + output_file << std::pow(10,nn) << std::setw(space) << results[0] << std::setw(space) - << results[1] << "\n"; + << "\n"; output_file.flush(); nn += Nbin; diff --git a/runcards/Higgs-FO-as-N.yaml b/runcards/Higgs-FO-as-N.yaml index 635397b..6a15e03 100644 --- a/runcards/Higgs-FO-as-N.yaml +++ b/runcards/Higgs-FO-as-N.yaml @@ -11,8 +11,8 @@ order : 1 # N bins, now x-bins Nmin : -12 -Nmax : -3 -Nbin : 0.1 +Nmax : -1 +Nbin : 0.5 # relevant scales in GeV mh : 125 # mass of colorless final state diff --git a/src/integration.cpp b/src/integration.cpp index a6b631d..14da3bb 100644 --- a/src/integration.cpp +++ b/src/integration.cpp @@ -50,18 +50,19 @@ int integration::CubaIntegrator( static const int nvec = 1; // static const int spin = -1; static const int maxpass = 5; - static const int mineval = 0; - static const int nnew = 1000; - static const int nstart = 1000; - static const int nbatch = 1000; - static const int nincrease = 500; + static const int mineval = 1000; + static const int nnew = 1e5; + static const int nmin = 1; + static const int nstart = 1e6; + static const int nbatch = 1e4; + static const int nincrease = 5e5; static const int ldxgiven = ndim; - static const int maxeval = 100000; + static const int maxeval = 1e9; static const int verbose = verbos; static const double border = 0.1; static const double maxchisq = 10.; - static const double epsabs = 1e-15; - static const double flatness = 25.; + static const double epsabs = 1e-100; + static const double flatness = 1; static const double mindeviation = 0.25; const char *statefile = NULL; @@ -82,7 +83,7 @@ int integration::CubaIntegrator( break; case 1: Suave(ndim, ncomp, func, par, nvec, prec, epsabs, verbose, seed, mineval, - maxeval, nnew, 5, flatness, statefile, NULL, &nregions, &neval, + maxeval, nnew, nmin, flatness, statefile, NULL, &nregions, &neval, fail, ris, err, proba); break; case 2: @@ -119,7 +120,7 @@ double integration::IntegrateDeltaPartonic(int method, double(Func)(double, void *), void *pp, double *error) { int fail; - double prec = 1e-8; + double prec = 1e-3; double *res = NULL; double *err = NULL; double *prb = NULL; @@ -156,7 +157,7 @@ double integration::IntegrateDistrPartonic(int method, double(Func)(double, double, void *), void *pp, double *error) { int fail; - double prec = 1e-8; + double prec = 1e-3; double *res = NULL; double *err = NULL; double *prb = NULL; @@ -198,7 +199,7 @@ double integration::IntegrateDeltaHadronic(int method, double(Func)(double, double, void *), void *pp, double *error) { int fail; - double prec = 1e-8; + double prec = 1e-3; double *res = NULL; double *err = NULL; double *prb = NULL; @@ -241,7 +242,7 @@ double integration::IntegrateDistrHadronic(int method, void *), void *pp, double *error) { int fail; - double prec = 1e-8; + double prec = 1e-3; double *res = NULL; double *err = NULL; double *prb = NULL; diff --git a/src/partonic.cpp b/src/partonic.cpp index d031704..f6d0043 100644 --- a/src/partonic.cpp +++ b/src/partonic.cpp @@ -88,7 +88,7 @@ std::vector CrossHiggs::partonichiggsdpt(double pt, double nn) { // y2: upper boundary of the rapidity yh // ////////////////////////////////////////////////// - int method = 0; + int method = 1; double deltres, delterr; double distres = 0, disterr = 0; std::vector results; From 1df5c5a27fb399e2142d4b288b5285dca94a14fe Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Sat, 26 Jun 2021 23:58:34 +0200 Subject: [PATCH 2/5] save --- src/higgsptpartonic.cpp | 46 ++++++++++++++++++----------------------- 1 file changed, 20 insertions(+), 26 deletions(-) diff --git a/src/higgsptpartonic.cpp b/src/higgsptpartonic.cpp index 0a95cd1..46fa4d2 100644 --- a/src/higgsptpartonic.cpp +++ b/src/higgsptpartonic.cpp @@ -988,9 +988,16 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, za = -th / (QQ2 - th); double jac10 = QQ2max / sqrt(std::pow(sh + MH2 - QQ2, 2) - 4. * sh * mt2); - double a10factor = + + double apolefactor = (log(qq) / qq + log(QQ2max * za / -th) / qq) / QQ2max * (-th / za); - double b10factor = 1. / qq / QQ2max * (-th / za); + double adeltafactor = + 0.5 * std::pow(log(-QQ2max / th), 2) / QQ2max * (-th / za); + double a10factor = apolefactor - adeltafactor; + double bdeltafactor = log(-QQ2max / th) / QQ2max * (-th / za); + double bpolefactor = 1. / qq / QQ2max * (-th / za); + double b10factor = bpolefactor - bdeltafactor; + shnew = za * sh; thnew = th; @@ -1039,20 +1046,11 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, } break; } - double adeltafactor = - 0.5 * std::pow(log(-QQ2max / th), 2) / QQ2max * (-th / za); - double bdeltafactor = log(-QQ2max / th) / QQ2max * (-th / za); - - double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor + - b10 * jac10 * bdeltafactor; - double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor + - a10 * jac10 * adeltafactor; + double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor; + double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor; double cfinal = c1 * jac1; double nonsingularfinal = nonsingular * jac1; - - - // nonsingular = 0; double result = afinal + bfinal + cfinal + nonsingularfinal; // dsigma/pt² to dsigma/dpt @@ -1206,9 +1204,14 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, za = -th / (QQ2 - th); double jac10 = QQ2max / sqrt(std::pow(sh + MH2 - QQ2, 2) - 4. * sh * mt2); - double a10factor = + double apolefactor = (log(qq) / qq + log(QQ2max * za / -th) / qq) / QQ2max * (-th / za); - double b10factor = 1. / qq / QQ2max * (-th / za); + double adeltafactor = + 0.5 * std::pow(log(-QQ2max / th), 2) / QQ2max * (-th / za); + double a10factor = apolefactor - adeltafactor; + double bdeltafactor = log(-QQ2max / th) / QQ2max * (-th / za); + double bpolefactor = 1. / qq / QQ2max * (-th / za); + double b10factor = bpolefactor - bdeltafactor; shnew = za * sh; thnew = th; @@ -1257,20 +1260,11 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, } break; } - double adeltafactor = - 0.5 * std::pow(log(-QQ2max / th), 2) / QQ2max * (-th / za); - double bdeltafactor = log(-QQ2max / th) / QQ2max * (-th / za); - - double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor + - b10 * jac10 * bdeltafactor; - double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor + - a10 * jac10 * adeltafactor; + double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor; + double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor; double cfinal = c1 * jac1; double nonsingularfinal = nonsingular * jac1; - - - // nonsingular = 0; double result = afinal + bfinal + cfinal + nonsingularfinal; // dsigma/pt² to dsigma/dpt From fa94d6a289b1d220036fd1a60956e83eb3c42cad Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Sun, 27 Jun 2021 09:11:08 +0200 Subject: [PATCH 3/5] save --- src/higgsptpartonic.cpp | 3 +++ 1 file changed, 3 insertions(+) diff --git a/src/higgsptpartonic.cpp b/src/higgsptpartonic.cpp index 46fa4d2..bd193d9 100644 --- a/src/higgsptpartonic.cpp +++ b/src/higgsptpartonic.cpp @@ -1194,6 +1194,9 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, } break; } + // double a10factor = log(1. - za) / (1. - za); + // double b10factor = 1. / (1. - za); + // Here we deal with the za=1 part of the plus distribution. za=1 corresponds // to setting QQ2=0. QQ2 = 0; From 89c74af06e24121b2d52b8cd0cd62ef1ca612b12 Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Mon, 28 Jun 2021 16:44:03 +0200 Subject: [PATCH 4/5] fix x-space agreement with HE --- src/higgsptpartonic.cpp | 21 +++++++++++---------- 1 file changed, 11 insertions(+), 10 deletions(-) diff --git a/src/higgsptpartonic.cpp b/src/higgsptpartonic.cpp index bd193d9..6f175e4 100644 --- a/src/higgsptpartonic.cpp +++ b/src/higgsptpartonic.cpp @@ -860,7 +860,7 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, // multiplied by either delta(Q) or delta(Q^2) (both singular // // and non-Singular) // //////////////////////////////////////////////////////////////// - double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0; + double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0, d10 = 0; zz2=1; //unused integration variable double qq = zz1; // qq = QQ2/QQ2max is an integration variable used to @@ -988,7 +988,6 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, za = -th / (QQ2 - th); double jac10 = QQ2max / sqrt(std::pow(sh + MH2 - QQ2, 2) - 4. * sh * mt2); - double apolefactor = (log(qq) / qq + log(QQ2max * za / -th) / qq) / QQ2max * (-th / za); double adeltafactor = @@ -998,12 +997,12 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, double bpolefactor = 1. / qq / QQ2max * (-th / za); double b10factor = bpolefactor - bdeltafactor; + double d10factor = -th/za/QQ2max; shnew = za * sh; thnew = th; uhnew = za*(uh-MH2)+MH2; - coeff(pt, uh, th, sh, MH2); switch (CHANNEL) { @@ -1015,6 +1014,7 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, (1. / th * pgg(za) * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2) + za / th * big1 * log((QQ2 + pt * pt) * za / (-th)) + za / th * big2); + d10 += 1. / th * beta0 * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2); } break; case (1): // gq-channel { @@ -1049,9 +1049,10 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor; double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor; double cfinal = c1 * jac1; + double dfinal = d10 * jac10 * d10factor; double nonsingularfinal = nonsingular * jac1; - double result = afinal + bfinal + cfinal + nonsingularfinal; + double result = afinal + bfinal + cfinal + dfinal + nonsingularfinal; // dsigma/pt² to dsigma/dpt result *= 2. * pt; @@ -1076,7 +1077,7 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, // multiplied by either delta(Q) or delta(Q^2) (both singular // // and non-Singular) // //////////////////////////////////////////////////////////////// - double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0; + double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0, d10 = 0; zz2=1; //unused integration variable double qq = zz1; // qq = QQ2/QQ2max is an integration variable used to @@ -1194,9 +1195,6 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, } break; } - // double a10factor = log(1. - za) / (1. - za); - // double b10factor = 1. / (1. - za); - // Here we deal with the za=1 part of the plus distribution. za=1 corresponds // to setting QQ2=0. QQ2 = 0; @@ -1216,11 +1214,12 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, double bpolefactor = 1. / qq / QQ2max * (-th / za); double b10factor = bpolefactor - bdeltafactor; + double d10factor = -th/za/QQ2max; + shnew = za * sh; thnew = th; uhnew = za*(uh-MH2)+MH2; - coeff(pt, uh, th, sh, MH2); switch (CHANNEL) { @@ -1232,6 +1231,7 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, (1. / th * pgg(za) * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2) + za / th * big1 * log((QQ2 + pt * pt) * za / (-th)) + za / th * big2); + d10 += 1. / th * beta0 * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2); } break; case (1): // gq-channel { @@ -1266,9 +1266,10 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor; double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor; double cfinal = c1 * jac1; + double dfinal = d10 * jac10 * d10factor; double nonsingularfinal = nonsingular * jac1; - double result = afinal + bfinal + cfinal + nonsingularfinal; + double result = afinal + bfinal + cfinal + dfinal + nonsingularfinal; // dsigma/pt² to dsigma/dpt result *= 2. * pt; From 42e1e963015bda1c14715b183b3dfe8bb3065575 Mon Sep 17 00:00:00 2001 From: RoyStegeman Date: Thu, 1 Jul 2021 16:24:04 +0000 Subject: [PATCH 5/5] Committing clang-format changes --- src/higgsptpartonic.cpp | 57 +++++++++++++++++++++-------------------- 1 file changed, 29 insertions(+), 28 deletions(-) diff --git a/src/higgsptpartonic.cpp b/src/higgsptpartonic.cpp index 6f175e4..598d2c1 100644 --- a/src/higgsptpartonic.cpp +++ b/src/higgsptpartonic.cpp @@ -692,7 +692,7 @@ double HiggsDpTpartonic::deltapartonic(double pt, double nn, double zz) { /////////////////////////////////////////////////////////////// double xx = nn; double result = 0; - zz=1; // unused integration variable + zz = 1; // unused integration variable // This function calculates the terms of the cross section proportional to // delta(Q²) @@ -740,8 +740,8 @@ double HiggsDpTpartonic::deltapartonic(double pt, double nn, double zz) { // part of the xsec that is also calculated in this function. case (0): // gg-channel { - // dominant contribution in small-x are the terms multiplied by uu. - // specifically the logs (not Li2) dominate uu + // dominant contribution in small-x are the terms multiplied by uu. + // specifically the logs (not Li2) dominate uu double uu = 0.5 * std::pow(log(uh / th), 2) + std::pow(M_PI, 2) / 3. - log(sh / MH2) * log(-th / MH2) - log(sh / MH2) * log(-uh / MH2) - @@ -749,13 +749,12 @@ double HiggsDpTpartonic::deltapartonic(double pt, double nn, double zz) { std::pow(log(MH2 / (MH2 - th)), 2) + std::pow(log(MH2 / (MH2 - uh)), 2) + 2. * Li2(1. - MH2 / sh) + 2. * Li2(MH2 / (MH2 - th)) + 2. * Li2(MH2 / (MH2 - uh)); - uu = - 0.5 * std::pow(log(uh / th), 2) + std::pow(M_PI, 2) / 3. - - log(sh / MH2) * log(-th / MH2) - log(sh / MH2) * log(-uh / MH2) - - log(-uh / MH2) * log(-th / MH2) + std::pow(log(MH2 / sh), 2) + - std::pow(log(MH2 / (MH2 - th)), 2) + - std::pow(log(MH2 / (MH2 - uh)), 2) + 2. * Li2(1. - MH2 / sh) + - 2. * Li2(MH2 / (MH2 - th)) + 2. * Li2(MH2 / (MH2 - uh)); + uu = 0.5 * std::pow(log(uh / th), 2) + std::pow(M_PI, 2) / 3. - + log(sh / MH2) * log(-th / MH2) - log(sh / MH2) * log(-uh / MH2) - + log(-uh / MH2) * log(-th / MH2) + std::pow(log(MH2 / sh), 2) + + std::pow(log(MH2 / (MH2 - th)), 2) + + std::pow(log(MH2 / (MH2 - uh)), 2) + 2. * Li2(1. - MH2 / sh) + + 2. * Li2(MH2 / (MH2 - th)) + 2. * Li2(MH2 / (MH2 - uh)); double de = 1.5 * beta0 * (log(-MUR2 / th) + log(-MUR2 / uh)) + 67. / 6. - 5. / 9. * NF; result += aass / (2. * M_PI) * @@ -862,7 +861,7 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, //////////////////////////////////////////////////////////////// double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0, d10 = 0; - zz2=1; //unused integration variable + zz2 = 1; // unused integration variable double qq = zz1; // qq = QQ2/QQ2max is an integration variable used to // integrate out rapidity double xx = nn; // xx = Q²/sh @@ -891,15 +890,16 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, // double a1factor = log(1. - za) / (1. - za); // double b1factor = 1. / (1. - za); // Alternatively one can write `a1factor' and `b1factor' in terms of qq: - double a1factor = (log(qq)/qq+log(QQ2max*za/-th)/qq)/QQ2max*(-th/za); - double b1factor = 1./qq/QQ2max*(-th/za); + double a1factor = + (log(qq) / qq + log(QQ2max * za / -th) / qq) / QQ2max * (-th / za); + double b1factor = 1. / qq / QQ2max * (-th / za); coeff(pt, uh, th, sh, MH2); REG(pt, uh, th, sh, MH2); double shnew = za * sh; double thnew = th; - double uhnew = za*(uh-MH2)+MH2; + double uhnew = za * (uh - MH2) + MH2; switch (CHANNEL) { // a1:: (log(1-za)/(1-za))+ terms @@ -997,11 +997,11 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, double bpolefactor = 1. / qq / QQ2max * (-th / za); double b10factor = bpolefactor - bdeltafactor; - double d10factor = -th/za/QQ2max; + double d10factor = -th / za / QQ2max; shnew = za * sh; thnew = th; - uhnew = za*(uh-MH2)+MH2; + uhnew = za * (uh - MH2) + MH2; coeff(pt, uh, th, sh, MH2); @@ -1014,7 +1014,8 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, (1. / th * pgg(za) * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2) + za / th * big1 * log((QQ2 + pt * pt) * za / (-th)) + za / th * big2); - d10 += 1. / th * beta0 * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2); + d10 += 1. / th * beta0 * log(-MUF2 * za / th) * + gg0(shnew, thnew, uhnew, MH2); } break; case (1): // gq-channel { @@ -1069,9 +1070,8 @@ double HiggsDpTpartonic::distrpartonic(double pt, double nn, double zz1, //==============================================================================================// - double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, - double zz2) { + double zz2) { //////////////////////////////////////////////////////////////// // This function computes the remaining terms that are not // // multiplied by either delta(Q) or delta(Q^2) (both singular // @@ -1079,7 +1079,7 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, //////////////////////////////////////////////////////////////// double nonsingular = 0, a1 = 0, b1 = 0, c1 = 0, a10 = 0, b10 = 0, d10 = 0; - zz2=1; //unused integration variable + zz2 = 1; // unused integration variable double qq = zz1; // qq = QQ2/QQ2max is an integration variable used to // integrate out rapidity double xx = nn; // xx = Q²/sh @@ -1108,15 +1108,16 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, // double a1factor = log(1. - za) / (1. - za); // double b1factor = 1. / (1. - za); // Alternatively one can write `a1factor' and `b1factor' in terms of qq: - double a1factor = (log(qq)/qq+log(QQ2max*za/-th)/qq)/QQ2max*(-th/za); - double b1factor = 1./qq/QQ2max*(-th/za); + double a1factor = + (log(qq) / qq + log(QQ2max * za / -th) / qq) / QQ2max * (-th / za); + double b1factor = 1. / qq / QQ2max * (-th / za); coeff(pt, uh, th, sh, MH2); REG(pt, uh, th, sh, MH2); double shnew = za * sh; double thnew = th; - double uhnew = za*(uh-MH2)+MH2; + double uhnew = za * (uh - MH2) + MH2; switch (CHANNEL) { // a1:: (log(1-za)/(1-za))+ terms @@ -1214,11 +1215,11 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, double bpolefactor = 1. / qq / QQ2max * (-th / za); double b10factor = bpolefactor - bdeltafactor; - double d10factor = -th/za/QQ2max; + double d10factor = -th / za / QQ2max; shnew = za * sh; thnew = th; - uhnew = za*(uh-MH2)+MH2; + uhnew = za * (uh - MH2) + MH2; coeff(pt, uh, th, sh, MH2); @@ -1231,7 +1232,8 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, (1. / th * pgg(za) * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2) + za / th * big1 * log((QQ2 + pt * pt) * za / (-th)) + za / th * big2); - d10 += 1. / th * beta0 * log(-MUF2 * za / th) * gg0(shnew, thnew, uhnew, MH2); + d10 += 1. / th * beta0 * log(-MUF2 * za / th) * + gg0(shnew, thnew, uhnew, MH2); } break; case (1): // gq-channel { @@ -1266,7 +1268,7 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, double bfinal = b1 * jac1 * b1factor - b10 * jac10 * b10factor; double afinal = a1 * jac1 * a1factor - a10 * jac10 * a10factor; double cfinal = c1 * jac1; - double dfinal = d10 * jac10 * d10factor; + double dfinal = d10 * jac10 * d10factor; double nonsingularfinal = nonsingular * jac1; double result = afinal + bfinal + cfinal + dfinal + nonsingularfinal; @@ -1285,4 +1287,3 @@ double HiggsDpTpartonic::distrpartoniccross(double pt, double nn, double zz1, } //==============================================================================================// -