← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/SubProcesses/reweight.f

  • Committer: Johan Alwall
  • Date: 2012-11-18 05:41:03 UTC
  • mfrom: (234.1.7 1.5.5)
  • Revision ID: johan@alwall.net-20121118054103-26xyfyez2n1d0732
Merged up to v. 1.5.5: Fixed bug in pdf reweighting giving asymmetric gamma+X distribution, fixed single top scale with maxjetflavor=4, fixed Aloha bug, check gauge bug, and improvement for LSF cluster.

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/SubProcesses/reweight.f
413
413
      data first/.true./
414
414
      data nwarning/0/
415
415
 
416
 
      logical isqcd,isjet,isparton,cluster
 
416
      logical isqcd,isjet,isparton,cluster,isjetvx
417
417
      double precision alphas
418
 
      external isqcd, isjet, isparton, cluster, alphas
 
418
      external isqcd, isjet, isparton, cluster, isjetvx, alphas
419
419
 
420
420
      setclscales=.true.
421
421
 
479
479
            q2fact(2)=q2bck(2)
480
480
         endif
481
481
      endif
482
 
      jfirst(1)=0
483
 
      jfirst(2)=0
484
482
 
485
483
      do i=1,2
486
484
         ibeam(i)=ishft(1,i-1)
 
485
         jfirst(i)=0
487
486
         jlast(i)=0
488
487
         jcentral(i)=0
489
 
         partonline(i)=isjet(ipdgcl(ibeam(i),igraphs(1),iproc))
 
488
         partonline(i)=isparton(ipdgcl(ibeam(i),igraphs(1),iproc))
490
489
         qcdline(i)=isqcd(ipdgcl(ibeam(i),igraphs(1),iproc))
491
490
      enddo
492
491
 
502
501
c             f1(x1,pt2E) is given by DSIG, just need to set scale.
503
502
                 ibeam(j)=imocl(n)
504
503
                 if(jfirst(j).eq.0)then
505
 
                    if(isjet(ipdgcl(imocl(n),igraphs(1),iproc))) then
 
504
                    if(isjetvx(imocl(n),idacl(n,1),idacl(n,2),
 
505
     $                 ipdgcl(1,igraphs(1),iproc),ipart,n.eq.nexternal-2)) then
506
506
                       jfirst(j)=n
507
507
                    else
508
508
                       jfirst(j)=-1