bug in FL NC with IC when computed in the FONLL scheme fixed

scarrazza · Dec 14, 2015 · 47f6454 · 47f6454
1 parent 1530ca8
commit 47f6454
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/src/DIS/ComputeDISOperators.f b/src/DIS/ComputeDISOperators.f
@@ -668,7 +668,7 @@ subroutine ComputeDISOperators(Q)
                         if(Nf_FF.lt.4)then
 *     FFNS
                            do pt=0,ipt_max_IC
-                              CLnsp(4) = CLnsp(4) + as(pt)
+                              CLnsp(4) = CLnsp(4) + cFLIC * as(pt)
      1                             * ( c0(4)
      2                             * SCLmNC(jgrid,ixi(4),
      3                             3,pt,alpha,beta)
@@ -683,10 +683,9 @@ subroutine ComputeDISOperators(Q)
      5                             * SC2mNC(jgrid,ixi(4)+1,
      6                             3,pt,alpha,beta) )
                            enddo
-                           CLnsp(4) = cFLIC * CLnsp(4)
 *     FFN0
                            do pt=0,ipt_max_IC
-c                              CLnsp(4) = CLnsp(4) - as(pt)
+c                              CLnsp(4) = CLnsp(4) - cFLIC * as(pt)
 c     1                             *  damp(4) * ( c0(4)
 c     2                             * SCLm0NC(jgrid,ixi(4),
 c     3                             3,pt,alpha,beta)

diff --git a/src/DIS/RSLintegralsDIS.f b/src/DIS/RSLintegralsDIS.f
@@ -848,7 +848,7 @@ subroutine RSLintegralsDIS(beta,alpha)
                      SC2mNC(igrid,ixi,k,2,beta,alpha) =
      1                    sig1 * SC2mNC(igrid,ixi,k,2,beta,alpha)
      2                    + sig2 * SC2m0NC(igrid,ixi,k,2,beta,alpha)
-                      SCLmNC(igrid,ixi,k,2,beta,alpha) =
+                     SCLmNC(igrid,ixi,k,2,beta,alpha) =
      1                    sig1 * SCLmNC(igrid,ixi,k,2,beta,alpha)
      2                    + sig2 * SCLm0NC(igrid,ixi,k,2,beta,alpha)
                   enddo