From 3b8514708ec4bfad7303d93e1deca818ac93c1ef Mon Sep 17 00:00:00 2001 From: Tobias Neumann Date: Thu, 28 Mar 2024 12:17:01 -0700 Subject: [PATCH] Use genpt with ptbreak=5 for gen_stop (MCFM-8.3 default), this seems to work really well. --- src/Phase/gen_stop.f | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/Phase/gen_stop.f b/src/Phase/gen_stop.f index 85e2dab..8f4b1e9 100644 --- a/src/Phase/gen_stop.f +++ b/src/Phase/gen_stop.f @@ -93,8 +93,8 @@ c xjac = exp(2._dp*log(sqrts/2._dp) - 2._dp*(1-r(ijet))/r(ijet))/r(ijet)** !if ((kcase == kbq_tpq_jet) .and. (kpart == kreal) .and. usescet .and. (ipsgen == 1)) then c if (ijet == 1) then ! in 1/sqrt(pT) - pt = r(ijet)**2*sqrts/2._dp - xjac = pt*sqrts/2._dp*2._dp*r(ijet) + !pt = r(ijet)**2*sqrts/2._dp + !xjac = pt*sqrts/2._dp*2._dp*r(ijet) c else c ! 1/pt^2 c pt = 1._dp / sqrt(4._dp/sqrts**2 + (1-r(ijet))/r(ijet)) @@ -105,7 +105,7 @@ c endif !error stop "todo" !else !error stop "please explicitly setup gen_stop for this case" - !call genpt(r(ijet),ptbreak,.true.,pt,xjac) + call genpt(r(ijet),5._dp,.false.,pt,xjac) !endif wt=wt*xjac