← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

  • Committer: olivier-mattelaer
  • Date: 2018-06-08 15:09:50 UTC
  • mto: This revision was merged to the branch mainline in revision 278.
  • Revision ID: olivier-mattelaer-20180608150950-i6ok8tzmjd7fd7la
fixing the DIS case

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/LO/SubProcesses/genps.f
279
279
         cm_rap=.5d0*dlog((p0+p3)/(p0-p3))
280
280
         set_cm_rap=.true.
281
281
c        Set shat
282
 
         s(-nbranch) = m2**2+2*xbk(1)*ebeam(1) *
283
 
     $                 (ebeam(2)+sqrt(ebeam(2)**2-m2**2))
 
282
         s(-nbranch) = x(ndim)*stot         
284
283
      elseif (abs(lpp(2)) .ge. 1) then
285
284
         call sample_get_x(sjac,x(ndim),ndim,mincfig,0d0,1d0)
286
285
         xbk(2) = x(ndim)
 
286
 
287
287
c        Set CM rapidity for use in the rap() function
288
288
         p0=ebeam(1)+xbk(2)*ebeam(2)
289
289
         p3=sqrt(ebeam(1)**2-m1**2)-xbk(2)*ebeam(2)
290
290
         cm_rap=.5d0*dlog((p0+p3)/(p0-p3))
291
291
         set_cm_rap=.true.
292
292
c        Set shat
293
 
         s(-nbranch) = m1**2+2*(ebeam(1)+sqrt(ebeam(1)**2-m1**2))
294
 
     $                 * xbk(2)*ebeam(2)
 
293
         s(-nbranch) =  x(ndim)*stot         
295
294
      else
296
295
c        Set CM rapidity for use in the rap() function
297
296
         p0=ebeam(1) + ebeam(2)