← Back to branch summary

~maddevelopers/mg5amcnlo/2.6.3_rwgt

« back to all changes in this revision

Viewing changes to Template/NLO/SubProcesses/add_write_info.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:
Template/NLO/SubProcesses/add_write_info.f
39
39
      external ran2
40
40
 
41
41
c Jamp amplitudes of the Born (to be filled with a call the sborn())
42
 
      double Precision amp2(maxamps), jamp2(0:maxamps)
 
42
      double Precision amp2(ngraphs), jamp2(0:ncolor)
43
43
      common/to_amps/  amp2,       jamp2
44
44
 
45
45
C iforest and other configuration info. Read once and saved.
70
70
      logical OnBW(-nexternal:0)
71
71
 
72
72
c LesHouches info
73
 
      integer maxflow
74
 
      parameter (maxflow=999)
75
73
      integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc),
76
74
     &     icolup(2,nexternal,maxflow),niprocs
77
75
c      include "leshouche.inc"
540
538
      do i=-ns,nexpart
541
539
         jpart(4,i)=icolalt(1,i)
542
540
         jpart(5,i)=icolalt(2,i)
543
 
         if(i.eq.1.or.i.eq.2) then 
 
541
         if(i.ge.1.and.i.le.nincoming) then 
544
542
            ito(i)=i            ! initial state particle
545
 
         elseif(i.ge.3) then 
 
543
         elseif(i.ge.nincoming+1) then 
546
544
            ito(i)=i+nres       ! final state particle
547
545
         elseif(i.le.-1.and.jpart(6,i).eq.2) then
548
546
            ires=ires+1
549
547
            ito(i)=2+ires       ! s-channel resonances
550
548
         else 
551
 
            ito(i)=0
 
549
            ito(i)=i
552
550
            if(i.eq.0) cycle
553
551
         endif
554
552
         if(jpart(2,i).lt.0.and.jpart(6,jpart(2,i)).ne.2) then
983
981
      implicit none
984
982
      include 'genps.inc'
985
983
      include 'nexternal.inc'
986
 
      integer    maxflow
987
 
      parameter (maxflow=999)
988
984
      integer i
989
985
      integer idup(nexternal,maxproc)
990
986
      integer mothup(2,nexternal,maxproc)