← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to madgraph/iolibs/template_files/madevent_symmetry.f

  • Committer: Rikkert Frederix
  • Date: 2021-09-09 15:51:40 UTC
  • mfrom: (78.75.502 3.2.1)
  • Revision ID: frederix@physik.uzh.ch-20210909155140-rg6umfq68h6h47cf
merge with 3.2.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
madgraph/iolibs/template_files/madevent_symmetry.f
47
47
c     DATA
48
48
c
49
49
      integer iforest(2,-max_branch:-1,lmaxconfigs)
 
50
      integer tstrategy(lmaxconfigs)
50
51
      integer sprop(maxsproc,-max_branch:-1,lmaxconfigs)
51
52
      integer tprid(-max_branch:-1,lmaxconfigs)
52
53
      include 'configs.inc'
103
104
         if (abs(lpp(1)) .eq. 3) m1 = 0.000511d0
104
105
         if (abs(lpp(2)) .eq. 3) m2 = 0.000511d0
105
106
         if (mass_ion(1).ge.0d0) m1 = mass_ion(1)
106
 
         if (mass_ion(2).ge.0d0) m1 = mass_ion(2)
 
107
         if (mass_ion(2).ge.0d0) m2 = mass_ion(2)
107
108
         if(ebeam(1).lt.m1.and.lpp(1).ne.9) ebeam(1)=m1
108
109
         if(ebeam(2).lt.m2.and.lpp(2).ne.9) ebeam(2)=m2
109
110
         pi1(0)=ebeam(1)
155
156
c     Arguments
156
157
c
157
158
      integer mapconfig(0:lmaxconfigs),use_config(0:lmaxconfigs)
158
 
      double precision prwidth(-max_branch:-1,lmaxconfigs)  !Propagotor width
 
159
      double precision prwidth(-max_branch:-1,lmaxconfigs) !Propagotor width
159
160
      integer iforest(2,-max_branch:-1,lmaxconfigs)
160
161
      integer sprop(maxsproc,-max_branch:-1,lmaxconfigs)
161
162
      integer jcomp
200
201
               enddo
201
202
               do j=1,nexternal-3
202
203
c                  write(*,*) 'Width',prwidth(-j,i),j,ic                 
203
 
                  if (prwidth(-j,i) .gt. 0d0 .and. iforest(1,-j,i).ne.1) then
 
204
            if (.not.(iforest(1,-j,i) .eq. 1 .or. prwidth(-j,i).le.0.or.
 
205
     &                       (nincoming.eq.2.and.iforest(1,-j,i) .eq. 2))) then
204
206
                     nbw=nbw+1
205
207
c                     write(*,*) 'Got bw',-nbw,j
206
208
c                    JA 4/8/11 don't treat forced BW differently
354
356
c     chain and could potentially conflict
355
357
c
356
358
      i=1
357
 
      do while (i .lt. nexternal-2 .and. itree(1,-i) .ne. 1)
 
359
      do while (i .lt. nexternal-2 .and. .not. (
 
360
     &    itree(1,-i) .eq. 1.or.(nincoming.eq.2.and.itree(1,-i).eq.2)))
358
361
         xmass(-i) = xmass(itree(1,-i))+xmass(itree(2,-i))
359
362
         mtot=mtot-xmass(-i)
360
363
         if (prwidth(-i,iconfig) .gt. 0d0) then