← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

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

  • Committer: Rikkert Frederix
  • Date: 2019-10-17 08:36:05 UTC
  • mfrom: (275.7.31 2.6.7)
  • Revision ID: frederix@physik.uzh.ch-20191017083605-bg78bda5jiq8u614
merge with public 2.6.7

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/LO/SubProcesses/setcuts.f
664
664
         endif
665
665
      enddo
666
666
 
667
 
c     ensure symmetry of s_min(i,j)
 
667
c     ensure symmetry of s_min(i,j) + mass effect
668
668
      do i=nincoming+1,nexternal-1
669
669
        do j=nincoming+1,nexternal-1
670
670
           if(j.lt.i)then
671
 
               s_min(i,j) = max(s_min(j,i),s_min(i,j))
 
671
               s_min(i,j) = max(s_min(j,i), s_min(i,j), (pmass(i)+pmass(j))**2)
672
672
            endif
673
673
        enddo
674
674
      enddo
 
675
 
 
676
c     ensure that smin is bigger than mass threshold
 
677
      smin_p = 0d0
 
678
      do i=nincoming+1, nexternal
 
679
         smin_p = smin_p + pmass(i)
 
680
      enddo
 
681
      smin = max(smin, smin_p**2)
675
682
      write(*,*) "Define smin to", smin
676
683
 
677
684
c Check that results are consistent among all the grouped subprocesses