← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to madgraph/iolibs/template_files/auto_dsig_v4.inc

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

Show diffs side-by-side

added added

removed removed

Lines of Context:
madgraph/iolibs/template_files/auto_dsig_v4.inc
10
10
C             pp    4 momentum of external particles
11
11
C             wgt   weight from Monte Carlo
12
12
C             imode 0 run, 1 init, 2 reweight, 
13
 
C                   3 finalize, 4 only PDFs
 
13
C                   3 finalize, 4 only PDFs,
 
14
C                   5 squared amplitude only (never
 
15
C                     generate events)
14
16
C     Output:
15
17
C             Amplitude squared and summed
16
18
C ****************************************************
89
91
 20     WRITE(*,*)'Error opening symfact.dat. No symmetry factor used.'
90
92
        RETURN
91
93
      ENDIF
92
 
C   Only run if IMODE is 0
93
 
      IF(IMODE.NE.0.AND.IMODE.NE.4) RETURN
 
94
C     Continue only if IMODE is 0, 4 or 5
 
95
      IF(IMODE.NE.0.AND.IMODE.NE.4.and.IMODE.NE.5) RETURN
94
96
 
95
97
%(passcuts_begin)s
96
98
%(pdf_lines)s
98
100
            DSIG%(proc_id)s = PD(0)
99
101
            RETURN
100
102
         ENDIF
101
 
         CALL SMATRIX%(proc_id)s(PP,DSIGUU)                                              
 
103
         CALL SMATRIX%(proc_id)s(PP,DSIGUU)
 
104
                 IF (IMODE.EQ.5) THEN
 
105
            IF (DSIGUU.LT.1D199) THEN            
 
106
                      DSIG%(proc_id)s = DSIGUU*CONV
 
107
                        else
 
108
                      DSIG%(proc_id)s = 0.0d0
 
109
                        endif
 
110
                        RETURN
 
111
                 ENDIF
102
112
C     Select a flavor combination (need to do here for right sign)
103
113
         CALL RANMAR(R)
104
114
         IPSEL=0
115
125
             DSIGUU=0D0
116
126
             DSIG%(proc_id)s=0D0
117
127
         ENDIF
 
128
C       Generate events only if IMODE is 0.
118
129
        IF(IMODE.EQ.0.AND.DABS(DSIG%(proc_id)s).GT.0D0)THEN
119
130
C       Call UNWGT to unweight and store events
120
131
           CALL UNWGT(PP,DSIG%(proc_id)s*WGT,%(numproc)d)
121
132
        ENDIF
122
133
%(passcuts_end)s
123
134
      END
 
135
## if(read_write_good_hel) {
 
136
C
 
137
C     Functionality to handling grid
 
138
C
 
139
%(read_write_good_hel)s
 
140
 
 
141
## }
 
 
b'\\ No newline at end of file'