← Back to branch summary

~maddevelopers/mg5amcnlo/3.1.1_lo_vectorization

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/super_auto_dsig.f

  • Committer: olivier-mattelaer
  • Date: 2022-05-20 18:53:15 UTC
  • Revision ID: olivier-mattelaer-20220520185315-5zvm5gdhs4qdcnyp
fix the creation of model for NLO mode

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/super_auto_dsig.f
50
50
      COMMON /TO_SYMCONF/ SYMCONF
51
51
 
52
52
 
53
 
      DOUBLE PRECISION XDUM, XSDUM
 
53
      DOUBLE PRECISION XDUM, XSDUM, DUM
54
54
 
55
55
      INTEGER LMAPPED
56
56
 
87
87
              IF(IMIRROR.EQ.1.OR.MIRRORPROCS(IPROC))THEN
88
88
C               Calculate PDF weight for all subprocesses
89
89
                XSDUM =  DSIGPROC(PP,J,IPROC,IMIRROR,SYMCONF,CONFSUB
90
 
     $           ,1D0,4)
 
90
     $           ,DUM,4)
91
91
                SELPROC(IMIRROR,IPROC,J)= SELPROC(IMIRROR,IPROC,J) +
92
92
     $            XSDUM
93
93
                IF(MC_GROUPED_SUBPROC) THEN
299
299
      IPROC_GLOBAL = IPROC
300
300
      ICONFIG = ICONF
301
301
 
302
 
C     set the running scale (MLM not working)
 
302
C     set the running scale 
303
303
C     and update the couplings accordingly
304
 
      CALL UPDATE_SCALE_COUPLING(ALL_P, ALL_WGT, ALL_Q2FACT, NB_PAGE)
305
 
 
 
304
      CALL UPDATE_SCALE_COUPLING_VEC(ALL_P, ALL_WGT, ALL_Q2FACT,
 
305
     $  NB_PAGE)
306
306
 
307
307
      IF(GROUPED_MC_GRID_STATUS.EQ.0) THEN
308
308
C       If we were in the initialization phase of the grid for MC over
459
459
      DATA LAST_ICONF/-1/
460
460
      COMMON/TO_LAST_ICONF/LAST_ICONF
461
461
 
 
462
      DOUBLE PRECISION DUM
462
463
      LOGICAL INIT_MODE
463
464
      COMMON /TO_DETERMINE_ZERO_HEL/INIT_MODE
464
465
C     ----------
588
589
              IF(IMIRROR.EQ.1.OR.MIRRORPROCS(IPROC))THEN
589
590
C               Calculate PDF weight for all subprocesses
590
591
                SELPROC(IMIRROR,IPROC,J)=DSIGPROC(PP,J,IPROC,IMIRROR
591
 
     $           ,SYMCONF,CONFSUB,1D0,4)
 
592
     $           ,SYMCONF,CONFSUB,DUM,4)
592
593
                IF(MC_GROUPED_SUBPROC) THEN
593
594
                  CALL MAP_3_TO_1(J,IPROC,IMIRROR,MAXSPROC,2,LMAPPED)
594
595
                  CALL DS_ADD_ENTRY('PDF_convolution',LMAPPED
615
616
        RETURN
616
617
      ENDIF
617
618
 
 
619
 
618
620
C     Perform the selection
619
621
      CALL RANMAR(R)
620
622
 
701
703
 
702
704
        IF(IPROC.EQ.0) RETURN
703
705
 
 
706
 
704
707
C       Update weigth w.r.t SELPROC normalized to selection probability
705
708
 
706
709
        WGT=WGT*(SUMPROB/SELPROC(IMIRROR,IPROC,ICONF))
723
726
C        the call DSIGPROC just below.
724
727
        ALLOW_HELICITY_GRID_ENTRIES = .FALSE.
725
728
      ENDIF
 
729
 
726
730
C     Call DSIGPROC to calculate sigma for process
727
731
      DSIG=DSIGPROC(PP,ICONF,IPROC,IMIRROR,SYMCONF,CONFSUB,WGT,IMODE)
728
732
C     Reset ALLOW_HELICITY_GRID_ENTRIES
771
775
      INCLUDE 'maxamps.inc'
772
776
      INCLUDE 'coupl.inc'
773
777
      INCLUDE 'run.inc'
 
778
      INCLUDE 'vector.inc'
774
779
C     
775
780
C     ARGUMENTS 
776
781
C     
851
856
 
852
857
      DSIGPROC=0D0
853
858
 
 
859
C     not needed anymore ... can be removed ... set for debugging only
 
860
C        
854
861
      IF (.NOT.PASSCUTS(P1)) THEN
855
862
        STOP 1
856
863
      ENDIF
 
864
C     set the running scale 
 
865
C     and update the couplings accordingly
 
866
      IF (NB_PAGE.LE.1) THEN
 
867
        CALL UPDATE_SCALE_COUPLING(PP, WGT)
 
868
      ENDIF
 
869
 
 
870
 
 
871
 
 
872
 
857
873
      IF (IMODE.EQ.0D0.AND.NB_PASS_CUTS.LT.2**12)THEN
858
874
        NB_PASS_CUTS = NB_PASS_CUTS + 1
859
875
      ENDIF