← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

  • Committer: Rikkert Frederix
  • Date: 2015-08-19 12:50:52 UTC
  • mto: This revision was merged to the branch mainline in revision 260.
  • Revision ID: frederix@physik.uzh.ch-20150819125052-18lm4roak6iaevml
Fix by Marius Wiesemann in his stuff for the treatment of the bottom
Yukawa. The user need to uncomment this code explicitly to be able to
use it.

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/NLO/SubProcesses/fks_singular.f
1190
1190
         iwgt=1
1191
1191
         wgt_wo_pdf=(wgt(1,i) + wgt(2,i)*log(mu2_r/mu2_q) + wgt(3,i)
1192
1192
     &        *log(mu2_f/mu2_q))*g_strong(i)**QCDpower(i)
1193
 
     &        *rwgt_muR_dep_fac(sqrt(mu2_r))
 
1193
     &        *rwgt_muR_dep_fac(sqrt(mu2_r),sqrt(mu2_r))
1194
1194
         wgts(iwgt,i)=xlum * wgt_wo_pdf
1195
1195
         do j=1,iproc
1196
1196
            parton_iproc(j,i)=parton_iproc(j,i) * wgt_wo_pdf
1200
1200
c *_wgt_mint variable should be directly passed to the mint-integrator
1201
1201
c and not be part of the plots nor computation of the cross section.
1202
1202
            virt_wgt_mint=virt_wgt_mint*xlum*g_strong(i)**QCDpower(i)
1203
 
     &           *rwgt_muR_dep_fac(sqrt(mu2_r))
 
1203
     &           *rwgt_muR_dep_fac(sqrt(mu2_r),sqrt(mu2_r))
1204
1204
            born_wgt_mint=born_wgt_mint*xlum*g_strong(i)**QCDpower(i)
1205
 
     &           /(8d0*Pi**2)*rwgt_muR_dep_fac(sqrt(mu2_r))
 
1205
     &           /(8d0*Pi**2)*rwgt_muR_dep_fac(sqrt(mu2_r),sqrt(mu2_r))
1206
1206
         endif
1207
1207
      enddo
1208
1208
      call cpu_time(tAfter)
1314
1314
     &              /mu2_q)+wgt(3,i)*log(mu2_f(kf)/mu2_q))*g(kr)
1315
1315
     &              **QCDpower(i)
1316
1316
               wgts(iwgt,i)=wgts(iwgt,i)
1317
 
     &              *rwgt_muR_dep_fac(sqrt(mu2_r(kr)))
 
1317
     &              *rwgt_muR_dep_fac(sqrt(mu2_r(kr)),sqrt(mu2_r(1)))
1318
1318
            enddo
1319
1319
         enddo
1320
1320
      enddo
1399
1399
     &                 *veto_compensating_factor_new
1400
1400
               endif
1401
1401
               wgts(iwgt,i)=wgts(iwgt,i)
1402
 
     &              *rwgt_muR_dep_fac(sqrt(mu2_r(ks)))
 
1402
     &              *rwgt_muR_dep_fac(sqrt(mu2_r(ks)),sqrt(mu2_r(1)))
1403
1403
               wgts(iwgt,i)=wgts(iwgt,i)*veto_multiplier_new(ks,kh)
1404
1404
            enddo
1405
1405
         enddo
1453
1453
c add the weights to the array
1454
1454
            wgts(iwgt,i)=xlum * (wgt(1,i) + wgt(2,i)*log(mu2_r/mu2_q) +
1455
1455
     &           wgt(3,i)*log(mu2_f/mu2_q))*g_strong(i)**QCDpower(i)
1456
 
            wgts(iwgt,i)=wgts(iwgt,i)*rwgt_muR_dep_fac(sqrt(mu2_r))
 
1456
            wgts(iwgt,i)=wgts(iwgt,i)*
 
1457
     &                       rwgt_muR_dep_fac(sqrt(mu2_r),sqrt(mu2_r))
1457
1458
         enddo
1458
1459
      enddo
1459
1460
      call InitPDF(izero)