~maddevelopers/mg5amcnlo/2.6.6_bug_1813292

« back to all changes in this revision

Viewing changes to Template/NLO/SubProcesses/reweight.f

  • Committer: olivier-mattelaer
  • Date: 2017-05-26 07:48:55 UTC
  • mfrom: (271.1.33 2.5.5)
  • Revision ID: olivier-mattelaer-20170526074855-r463wfxlom110fiu
passĀ theĀ 2.5.5

Show diffs side-by-side

added added

removed removed

Lines of Context:
73
73
         if (iimode.eq.2) then
74
74
            gamma=gamma+CF*alphasq0**2*log(Q1/q0)*kappa/(2d0*pi**2) ! A2
75
75
         endif
76
 
         if (iipdg.gt.NF) then ! include mass effects
 
76
         if (abs(iipdg).gt.NF) then ! include mass effects
77
77
            qom=q0/mass(iipdg)
78
78
            gamma=gamma+CF*alphasq0/pi/2d0*( 0.5d0 - qom*atan(1d0/qom) -
79
79
     $           (1d0-0.5d0*qom**2)*log(1d0+1d0/qom**2) )
490
490
      if((ipart(1,ida1).ge.1.and.ipart(1,ida1).le.2).or.
491
491
     $   (ipart(1,ida2).ge.1.and.ipart(1,ida2).le.2))then
492
492
c     Check if ida1 is outgoing parton or ida2 is outgoing parton
493
 
         if(.not.islast.and.ipart(1,ida2).ge.1
494
 
     $        .and.ipart(1,ida2).le.2.and.isparton(idda1).or.ipart(1
495
 
     $        ,ida1).ge.1.and.ipart(1
496
 
     $        ,ida1).le.2.and.isparton(idda2))then
 
493
         if((.not.islast) .and. ((ipart(1,ida2).ge.1 .and. ipart(1
 
494
     $        ,ida2).le.2 .and. isparton(idda1)) .or. (ipart(1
 
495
     $        ,ida1).ge.1 .and. ipart(1,ida1).le.2 .and.
 
496
     $        isparton(idda2))))then
497
497
           ispartonvx=.true.
498
498
        else
499
499
           ispartonvx=.false.
1222
1222
     $     ,fks_configs,lmaxconfigs),jlast(2)
1223
1223
      real*8 q2bck(2)
1224
1224
      common /to_rw/jlast,njetstore,iqjetstore,q2bck
1225
 
      integer maxflow
1226
 
      parameter (maxflow=999)
1227
1225
      integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc),
1228
1226
     &     icolup(2,nexternal,maxflow),niprocs
1229
1227
      common /c_leshouche_inc/idup,mothup,icolup,niprocs