← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/NLO/SubProcesses/montecarlocounter.f

merge with 2.0.0beta4 revision 216

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/NLO/SubProcesses/montecarlocounter.f
497
497
c
498
498
 
499
499
      if (softtest.or.colltest) then
500
 
         tiny=1d-8
 
500
         tiny=1d-6
501
501
      else
502
 
         tiny=1d-6
 
502
         tiny=1d-4
503
503
      endif
504
504
 
505
505
      if(pp(0,1).le.0.d0)then
1209
1209
c
1210
1210
 
1211
1211
      if (softtest.or.colltest) then
1212
 
         tiny=1d-8
 
1212
         tiny=1d-6
1213
1213
      else
1214
 
         tiny=1d-6
 
1214
         tiny=1d-4
1215
1215
      endif
1216
1216
 
1217
1217
      if(pp(0,1).le.0.d0)then
1915
1915
c
1916
1916
 
1917
1917
      if (softtest.or.colltest) then
1918
 
         tiny=1d-8
 
1918
         tiny=1d-6
1919
1919
      else
1920
 
         tiny=1d-6
 
1920
         tiny=1d-4
1921
1921
      endif
1922
1922
 
1923
1923
      if(pp(0,1).le.0.d0)then
2710
2710
c
2711
2711
 
2712
2712
      if (softtest.or.colltest) then
2713
 
         tiny=1d-8
 
2713
         tiny=1d-6
2714
2714
      else
2715
 
         tiny=1d-6
 
2715
         tiny=1d-4
2716
2716
      endif
2717
2717
 
2718
2718
      if(pp(0,1).le.0.d0)then
3438
3438
c
3439
3439
 
3440
3440
      if (softtest.or.colltest) then
3441
 
         tiny=1d-8
 
3441
         tiny=1d-6
3442
3442
      else
3443
 
         tiny=1d-6
 
3443
         tiny=1d-4
3444
3444
      endif
3445
3445
 
3446
3446
      if(pp(0,1).le.0.d0)then
3575
3575
            call assign_scalemax(shat,xi_i_fks,upper_scale)
3576
3576
            if(sqrt(xi(npartner)).gt.upper_scale)lzone(npartner)=.false.
3577
3577
         endif
3578
 
c z limits, following  pages 14-16 of hep-ph/0408302
 
3578
c z limits, following pages 11 and 14-16 of hep-ph/0408302
3579
3579
         if(ileg.le.2)then
3580
3580
            zplus=1-sqrt(xi(npartner)/z(npartner)/shat)*
3581
3581
     &           (sqrt(1+xi(npartner)/(4*z(npartner)*shat))-sqrt(xi(npartner)/(4*z(npartner)*shat)))
3582
3582
            if(z(npartner).gt.zplus)lzone(npartner)=.false.
3583
3583
         endif
3584
 
c z limits for the 'constrained' definition, following
3585
 
c strictly page 354 of hep-ph/0603175
 
3584
c
3586
3585
         if(ileg.gt.2)then
3587
3586
            if(ileg.eq.3)then
3588
3587
               xmm2=xm12
3608
3607
               ma2=0d0
3609
3608
            endif
3610
3609
            ma=sqrt(ma2)
3611
 
            mbeff=sqrt(xmm2)
 
3610
            mbeff=0d0
3612
3611
            mceff=0d0
3613
3612
            en_fks=sqrt(s)*(1-x)/2.d0
3614
3613
            en_mother=en_fks+sqrt(xmm2+veckn_ev**2)