← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/RWGTNLO/setrun.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/RWGTNLO/setrun.f
1
 
C-------------------------------------------------
2
 
C   GET_PDFUP
3
 
C   Convert MadEvent pdf name to LHAPDF number
4
 
C-------------------------------------------------
5
 
 
6
 
      subroutine get_pdfup(pdfin,pdfgup,pdfsup,lhaid)
7
 
      implicit none
8
 
 
9
 
      character*(*) pdfin
10
 
      integer mpdf
11
 
      integer npdfs,i,pdfgup(2),pdfsup(2),lhaid
12
 
 
13
 
      parameter (npdfs=16)
14
 
      character*7 pdflabs(npdfs)
15
 
      data pdflabs/
16
 
     $   'none',
17
 
     $   'mrs02nl',
18
 
     $   'mrs02nn',
19
 
     $   'cteq4_m',
20
 
     $   'cteq4_l',
21
 
     $   'cteq4_d',
22
 
     $   'cteq5_m',
23
 
     $   'cteq5_d',
24
 
     $   'cteq5_l',
25
 
     $   'cteq5m1',
26
 
     $   'cteq6_m',
27
 
     $   'cteq6_l',
28
 
     $   'cteq6l1',     
29
 
     $   'nn23lo',
30
 
     $   'nn23lo1',
31
 
     $   'nn23nlo'/
32
 
      integer numspdf(npdfs)
33
 
      data numspdf/
34
 
     $   00000,
35
 
     $   20250,
36
 
     $   20270,
37
 
     $   19150,
38
 
     $   19170,
39
 
     $   19160,
40
 
     $   19050,
41
 
     $   19060,
42
 
     $   19070,
43
 
     $   19051,
44
 
     $   10000,
45
 
     $   10041,
46
 
     $   10042,
47
 
     $   246800,
48
 
     $   247000,
49
 
     $   244800/
50
 
 
51
 
 
52
 
      if(pdfin.eq."lhapdf") then
53
 
        write(*,*)'using LHAPDF'
54
 
        do i=1,2
55
 
           pdfgup(i)=-1
56
 
           pdfsup(i)=lhaid
57
 
        enddo
58
 
        return
59
 
      endif
60
 
 
61
 
      
62
 
      mpdf=-1
63
 
      do i=1,npdfs
64
 
        if(pdfin(1:len_trim(pdfin)) .eq. pdflabs(i))then
65
 
          mpdf=numspdf(i)
66
 
        endif
67
 
      enddo
68
 
 
69
 
      if(mpdf.eq.-1) then
70
 
        write(*,*)'ERROR: pdf ',pdfin,' not implemented in get_pdfup.'
71
 
        write(*,*)'known pdfs are'
72
 
        write(*,*) pdflabs
73
 
        open(unit=26,file='../../error',status='unknown')
74
 
        write(26,*)'ERROR: pdf ',pdfin,' not implemented in get_pdfup.'
75
 
        write(26,*)'known pdfs are'
76
 
        write(26,*) pdflabs
77
 
        stop 1
78
 
c$$$        write(*,*)'using ',pdflabs(12)
79
 
c$$$        mpdf=numspdf(12)
80
 
      endif
81
 
 
82
 
      do i=1,2
83
 
        pdfgup(i)=-1
84
 
        pdfsup(i)=mpdf
85
 
      enddo
86
 
 
87
 
      return
88
 
      end