← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.4_array

« back to all changes in this revision

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

  • Committer: olivier-mattelaer
  • Date: 2019-10-16 11:38:38 UTC
  • mfrom: (282.1.28 2.6.7)
  • Revision ID: olivier-mattelaer-20191016113838-6u4j5u3hn06rq2wc
pass to 2.6.7

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/LO/SubProcesses/reweight.f
403
403
c     Color singlet
404
404
         ipart(1,imo)=ipart(1,ida1)
405
405
         ipart(2,imo)=ipart(1,ida2)
 
406
      else if (abs(get_color(idmo)).eq.3.and.abs(get_color(idda1)).eq.3.and.abs(get_color(idda2)).eq.3)then
 
407
c     EPSILON^IJK color structure -- handle it as photon emission --
 
408
         ipart(1,imo)=ipart(1,ida1)
 
409
         ipart(2,imo)=0
 
410
      else if (abs(get_color(idmo)).eq.6.and.abs(get_color(idda1)).eq.3.and.abs(get_color(idda2)).eq.3)then
 
411
c     sextet -> (anti-)quark (anti-)quark': use both, but take hardest as 1
 
412
        if(p(1,ipart(1,ida1))**2+p(2,ipart(1,ida1))**2.gt.
 
413
     $     p(1,ipart(1,ida2))**2+p(2,ipart(1,ida2))**2) then
 
414
          ipart(1,imo)=ipart(1,ida1)
 
415
          ipart(2,imo)=ipart(1,ida2)
 
416
        else
 
417
          ipart(1,imo)=ipart(1,ida2)
 
418
          ipart(2,imo)=ipart(1,ida1)
 
419
        endif
406
420
      else
407
421
         write(*,*) idmo,'>', idda1, idda2, 'color', get_color(idmo),'>', get_color(idda1), get_color(idda2)
408
422
         write(*,*) "failed for ipartupdate. Please retry without MLM/default dynamical scale"