← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/LO/SubProcesses/unwgt.f
478
478
      double precision beta, get_betaz
479
479
      double precision ebi(0:3), ebo(0:3)
480
480
      double precision ptcltmp(nexternal), pdum(0:3)
481
 
 
 
481
      integer beam_number
 
482
      
482
483
      integer idup(nexternal,maxproc,maxsproc)
483
484
      integer mothup(2,nexternal)
484
485
      integer icolup(2,nexternal,maxflow,maxsproc)
705
706
         else
706
707
            write(s_buff(3), '(a)') '<asrwt>0</asrwt>'
707
708
         endif
708
 
         if(n_pdfrw(1).gt.0)then
 
709
         beam_number = 1
 
710
         if (flip) then
 
711
            beam_number =2
 
712
         endif
 
713
         
 
714
         if(n_pdfrw(1).gt.0.and.abs(lpp(1)).eq.1)then
709
715
            if(2*n_pdfrw(1).lt.10) then
710
 
               write(cfmt,'(a,I1,a,I1,a)') '(a,I3,',
 
716
               write(cfmt,'(a,I1,a,I1,a)') '(a,I1,a,I3,',
711
717
     $              n_pdfrw(1),'I9,',2*n_pdfrw(1),'E15.8,a)'
712
718
            else
713
 
               write(cfmt,'(a,I1,a,I2,a)') '(a,I3,',
 
719
               write(cfmt,'(a,I1,a,I2,a)') '(a,I1,a,I3,',
714
720
     $              n_pdfrw(1),'I9,',2*n_pdfrw(1),'E15.8,a)'
715
721
            endif
716
 
            write(s_buff(4), cfmt) '<pdfrwt beam="1">',
 
722
            
 
723
            write(s_buff(4), cfmt) '<pdfrwt beam="', beam_number, '">',
717
724
     $           n_pdfrw(1),(i_pdgpdf(i,1),i=1,n_pdfrw(1)),
718
725
     $           (s_xpdf(i,1),i=1,n_pdfrw(1)),
719
726
     $           (s_qpdf(i,1),i=1,n_pdfrw(1)),
720
727
     $           '</pdfrwt>'
721
728
         else
722
 
            write(s_buff(4), '(a)') '<pdfrwt beam="1">0</pdfrwt>'
723
 
         endif
724
 
         if(n_pdfrw(2).gt.0)then
 
729
            write(s_buff(4), '(a,I1,a)') '<pdfrwt beam="',
 
730
     $       beam_number,'">0</pdfrwt>'
 
731
         endif
 
732
         beam_number = 2
 
733
         if (flip) then
 
734
            beam_number = 1
 
735
         endif
 
736
         if(n_pdfrw(2).gt.0.and.abs(lpp(2)).eq.1)then
725
737
            if(2*n_pdfrw(2).lt.10) then
726
 
               write(cfmt,'(a,I1,a,I1,a)') '(a,I3,',
 
738
               write(cfmt,'(a,I1,a,I1,a)') '(a,I1,a,I3,',
727
739
     $              n_pdfrw(2),'I9,',2*n_pdfrw(2),'E15.8,a)'
728
740
            else
729
 
               write(cfmt,'(a,I1,a,I2,a)') '(a,I3,',
 
741
               write(cfmt,'(a,I1,a,I2,a)') '(a,I1,a,I3,',
730
742
     $              n_pdfrw(2),'I9,',2*n_pdfrw(2),'E15.8,a)'
731
743
            endif
732
 
            write(s_buff(5), cfmt) '<pdfrwt beam="2">',
 
744
            write(s_buff(5), cfmt) '<pdfrwt beam="',beam_number,'">',
733
745
     $           n_pdfrw(2),(i_pdgpdf(i,2),i=1,n_pdfrw(2)),
734
746
     $           (s_xpdf(i,2),i=1,n_pdfrw(2)),
735
747
     $           (s_qpdf(i,2),i=1,n_pdfrw(2)),