← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.1

« back to all changes in this revision

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

  • Committer: Marco Zaro
  • Date: 2014-01-27 16:54:10 UTC
  • mfrom: (78.124.55 MG5_aMC_2.1)
  • Revision ID: marco.zaro@gmail.com-20140127165410-5lma8c2hzbzm426j
merged with lp:~maddevelopers/madgraph5/MG5_aMC_2.1 r 267

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/LO/SubProcesses/setcuts.f
84
84
c
85
85
c
86
86
c     reading parameters
87
 
      include '../../Source/run_config.inc'
 
87
      include 'run_config.inc'
88
88
      character*20 param(maxpara),value(maxpara)
89
89
      integer npara
90
90
c
113
113
            scale=pmass(1)
114
114
            fixed_ren_scale=.true.
115
115
            fixed_fac_scale=.true.
 
116
            use_syst=.false.
116
117
         endif
117
118
c
118
119
c     set ptj and s_min if xqcut and ktscheme = 1, to improve
137
138
            write(*,*) 'Note that this might affect non-radiated jets,'
138
139
            write(*,*) 'e.g. from decays. Use cut_decays=F in run_card.'
139
140
          else if(mmjj.gt.xqcut)then
 
141
c           In principle this should never happen since the banner.py
 
142
c           is expected to correct this already.
140
143
            mmjj=0d0
141
144
            write(*,*) 'Warning! mmjj set to 0 since xqcut > 0 and'
142
145
            write(*,*) '         auto_ptj_mjj = F'
453
456
c     Compute Smin (for efficiency
454
457
c
455
458
      do i=nincoming+1,nexternal-1
456
 
        smin=0.0d0**2
 
459
      do j=nincoming+1,nexternal-1
 
460
         if(j.lt.i)then
 
461
            s_min(i,j) = s_min(j,i)
 
462
         else
 
463
            smin=0.0d0**2
457
464
 
458
465
            if(do_cuts(i).and.do_cuts(j)) then
459
466
               if(is_a_j(i).and.is_a_j(j)) s_min(j,i)=mmjj*dabs(mmjj)
472
479
     &            (idup(i,1,1)*idup(j,1,1).lt.0))
473
480
     &            s_max(j,i)=mmllmax*dabs(mmllmax)  !only on l+l- pairs (same flavour)
474
481
            endif
 
482
         endif
 
483
      enddo
475
484
      enddo
476
485
 
477
486
c