← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to Template/NLO/FixedOrderAnalysis/dbook.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/FixedOrderAnalysis/dbook.f
124
124
c     initialization
125
125
      IMPLICIT NONE
126
126
      include 'dbook.inc'
127
 
c     APPLgrid commons
128
 
      include "reweight_appl.inc"
129
 
      include "appl_common.inc"
130
 
      integer iappl
131
 
      common /for_applgrid/ iappl
 
127
c     PineAPPL commons
 
128
      include "reweight_pineappl.inc"
 
129
      include "pineappl_common.inc"
 
130
      logical pineappl
 
131
      common /for_pineappl/ pineappl
132
132
c
133
133
c     LOCAL
134
134
c
139
139
      NHIST=0
140
140
      DO 1, I=1,NPLOTS             
141
141
   1  BOOK(I)=' NO'
142
 
C     Initialize the number of bins of the aMCfast grids
143
 
      if(iappl.ne.0) appl_obs_nbins = 0
 
142
C     Initialize the number of bins of the amcblast grids
 
143
      if(pineappl) appl_obs_nbins = 0
144
144
      END  
145
145
 
146
146
 
1001
1001
      implicit none
1002
1002
      integer n,n_by4
1003
1003
      double precision var,www
1004
 
c     APPLgrid commons
1005
 
      include "reweight_appl.inc"
1006
 
      include "appl_common.inc"
1007
 
      integer iappl
1008
 
      common /for_applgrid/ iappl
 
1004
c     PineAPPL commons
 
1005
      include "reweight_pineappl.inc"
 
1006
      include "pineappl_common.inc"
 
1007
      logical pineappl
 
1008
      common /for_pineappl/ pineappl
1009
1009
      integer j
1010
 
      if(iappl.ne.0)then
 
1010
      if(pineappl)then
1011
1011
         do j=1,nh_obs
1012
1012
            if(n.eq.ih_obs(j))then
1013
1013
               appl_obs_num   = j
1014
1014
               appl_obs_histo = var 
1015
 
c     Fill the reference APPLgrid histograms
1016
 
               call APPL_fill_ref
1017
 
c     Fill the APPLgrid files
 
1015
c     Fill the reference PineAPPL histograms
 
1016
c     Fill the PineAPPL files
1018
1017
               call APPL_fill
1019
1018
            endif
1020
1019
         enddo
1029
1028
      integer i,n,n_by4
1030
1029
      character*(*) string
1031
1030
      double precision del,xl,xu
1032
 
c     APPLgrid commons
1033
 
      include "reweight_appl.inc"
1034
 
      include "appl_common.inc"
1035
 
      integer iappl
1036
 
      common /for_applgrid/ iappl
1037
 
c     Initialize the grids only if the switch "iappl" is different from zero
 
1031
c     PineAPPL commons
 
1032
      include "reweight_pineappl.inc"
 
1033
      include "pineappl_common.inc"
 
1034
      logical pineappl
 
1035
      common /for_pineappl/ pineappl
 
1036
c     Initialize the grids only if the switch "pineappl" is different from zero
1038
1037
c     and if the title string containes the word "central" and does not contain
1039
1038
c     the word "Born". 
1040
 
      if(iappl.ne.0.and.index(string,"central").ne.0.and.
 
1039
      if(pineappl.and.index(string,"central").ne.0.and.
1041
1040
     1                  index(string,"Born").eq.0)then
1042
1041
c     Observable parameters
1043
1042
c     Compute number of bins and edges only if they have not been given by the user.
1051
1050
         appl_obs_min = appl_obs_bins(0)
1052
1051
         appl_obs_max = appl_obs_bins(appl_obs_nbins)
1053
1052
         if(abs(appl_obs_max-xu).gt.0.00000001d0)then
1054
 
            write(*,*) 'APPLgrid Histogram: ', 
 
1053
            write(*,*) 'PineAPPL Histogram: ', 
1055
1054
     1                 'Change of the upper limit:',xu,'-->',
1056
1055
     2                  appl_obs_max
1057
1056
         endif
1058
 
c     Initialize APPLgrid routines
 
1057
c     Initialize PineAPPL routines
1059
1058
         call APPL_init
1060
1059
c     Keep track of the position of this histogram
1061
1060
         nh_obs = nh_obs + 1
1072
1071
      implicit none
1073
1072
      integer n,n_by4,m_by4
1074
1073
      character*(*) string1,string2,string3
1075
 
c     APPLgrid commons
1076
 
      include "reweight_appl.inc"
1077
 
      include "appl_common.inc"
1078
 
      integer iappl
1079
 
      common /for_applgrid/ iappl
 
1074
c     PineAPPL commons
 
1075
      include "reweight_pineappl.inc"
 
1076
      include "pineappl_common.inc"
 
1077
      logical pineappl
 
1078
      common /for_pineappl/ pineappl
1080
1079
      integer j
1081
 
      if(iappl.ne.0)then
 
1080
      if(pineappl)then
1082
1081
         do j=1,nh_obs
1083
1082
            if(n.eq.ih_obs(j))then
1084
1083
               appl_obs_num = j
1085
 
               call APPL_fill_ref_out
1086
1084
               call APPL_term
1087
1085
            endif
1088
1086
         enddo
1098
1096
      implicit none
1099
1097
      integer n,n_by4,lr,lh,m_by4
1100
1098
      character*(*) string1,string2,string3
1101
 
c     APPLgrid commons
1102
 
      include "reweight_appl.inc"
1103
 
      include "appl_common.inc"
1104
 
      integer iappl
1105
 
      common /for_applgrid/ iappl
 
1099
c     PineAPPL commons
 
1100
      include "reweight_pineappl.inc"
 
1101
      include "pineappl_common.inc"
 
1102
      logical pineappl
 
1103
      common /for_pineappl/ pineappl
1106
1104
      integer j
1107
 
      if(iappl.ne.0)then
 
1105
      if(pineappl)then
1108
1106
         do j=1,nh_obs
1109
1107
            if(n.eq.ih_obs(j))then
1110
1108
               appl_obs_num = j
1111
 
               call APPL_fill_ref_out
1112
1109
               call APPL_term
1113
1110
            endif
1114
1111
         enddo