← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to Template/NLO/Source/setrun.f

  • Committer: Rikkert Frederix
  • Date: 2021-09-09 15:51:40 UTC
  • mfrom: (78.75.502 3.2.1)
  • Revision ID: frederix@physik.uzh.ch-20210909155140-rg6umfq68h6h47cf
merge with 3.2.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/NLO/Source/setrun.f
29
29
      integer nevents
30
30
      character*7 event_norm
31
31
      common /event_normalisation/event_norm
32
 
      integer iappl
33
 
      common /for_applgrid/ iappl
 
32
      logical pineappl
 
33
      common /for_pineappl/ pineappl
34
34
      integer idum
35
35
      logical              fixed_order,nlo_ps
36
36
      common /c_fnlo_nlops/fixed_order,nlo_ps
38
38
      double precision D
39
39
      common/to_dj/D
40
40
c Include all the parameters set in the run_card.dat
 
41
 
41
42
      include 'run_card.inc'
42
43
c Change shower_MC string to upper case
43
44
      call to_upper(shower_MC)
65
66
      q2fact(2) = muF2_ref_fixed**2      ! fact scale**2 for pdf2     
66
67
      scalefact=muR_over_ref
67
68
      ellissextonfact=QES_over_ref
 
69
 
 
70
      call fill_needed_splittings()! this sub is inside fks_singular.f
 
71
 
68
72
c check that the event normalization input is reasoble
69
73
      buff = event_norm 
70
74
      call case_trap2(buff) ! requires a string of length 20 at least
110
114
          write(*,*) 'The default order of alpha_s running is fixed to '
111
115
     &         ,nloop
112
116
      endif
113
 
      if (nlo_ps) then
 
117
      if (nlo_ps.or.pineappl) then
114
118
C Fill common block for Les Houches init info
115
119
         do i=1,2
116
120
            if(lpp(i).eq.1.or.lpp(i).eq.2) then
121
125
               idbmup(i)=11
122
126
            elseif(lpp(i).eq.-3) then
123
127
               idbmup(i)=-11
 
128
            elseif(lpp(i).eq.4) then
 
129
               idbmup(i)=13
 
130
            elseif(lpp(i).eq.-4) then
 
131
               idbmup(i)=-13               
124
132
            elseif(lpp(i).eq.0) then
125
133
               open (unit=71,status='old',file='initial_states_map.dat')
126
134
               read (71,*,err=100)idum,idum,idbmup(1),idbmup(2)
130
138
            endif
131
139
            ebmup(i)=ebeam(i)
132
140
         enddo
133
 
         call get_pdfup(pdlabel,pdfgup,pdfsup,lhaid)
 
141
         if (abs(lpp(1)).eq.1 .or. abs(lpp(2)).eq.1)
 
142
     $       call get_pdfup(pdlabel,pdfgup,pdfsup,lhaid)
134
143
      endif
135
144
c Fill the nmemPDF(i) array with the number of PDF error set. This we
136
145
c get from LHAPDF.
155
164
      character*(*) pdfin
156
165
      integer mpdf
157
166
      integer npdfs,i,pdfgup(2),pdfsup(2),lhaid
158
 
      parameter (npdfs=16)
 
167
      parameter (npdfs=20)
159
168
      character*7 pdflabs(npdfs)
160
169
      data pdflabs/ 'none', 'mrs02nl', 'mrs02nn', 'cteq4_m', 'cteq4_l',
161
170
     $     'cteq4_d', 'cteq5_m', 'cteq5_d', 'cteq5_l', 'cteq5m1',
162
171
     $     'cteq6_m', 'cteq6_l', 'cteq6l1', 'nn23lo', 'nn23lo1',
163
 
     $     'nn23nlo'/
 
172
     $     'nn23nlo', 'ct14q00', 'ct14q07', 'ct14q14', 'ct14q21'/
164
173
      integer numspdf(npdfs)
165
174
      data numspdf/ 00000, 20250, 20270, 19150, 19170, 19160, 19050,
166
175
     $     19060, 19070, 19051, 10000, 10041, 10042, 246800, 247000,
167
 
     $     244800/
 
176
     $     244800, 666666, 666666, 666666, 666666/
168
177
      if(pdfin.eq."lhapdf") then
169
178
         write(*,*)'using LHAPDF'
170
179
         do i=1,2