← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to madgraph/iolibs/template_files/loop/loop_matrix_standalone.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/loop/loop_matrix_standalone.inc
291
291
  ML_INIT = .FALSE.
292
292
ENDIF
293
293
 
294
 
IF(NTRY.EQ.0) THEN
295
 
 
296
294
C Setup the file paths
297
295
  CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
298
296
  CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
301
299
  CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
302
300
  CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
303
301
  CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
304
 
  
305
 
  CALL %(proc_prefix)sSET_N_EVALS(N_DP_EVAL,N_QP_EVAL)
306
 
  HELDOUBLECHECKED=.NOT.DoubleCheckHelicityFilter  
307
 
  DO J=1,NCOMB
308
 
    DO I=1,NCTAMPS
309
 
          GOODAMP(I,J)=.TRUE.
310
 
        ENDDO
311
 
  ENDDO
312
 
OPEN(1, FILE=LoopFilterFN, err=100, status='OLD',           action='READ')
313
 
  DO J=1,NCOMB
314
 
    READ(1,*,END=101) (GOODAMP(I,J),I=NCTAMPS+1,NLOOPAMPS)
315
 
  ENDDO
316
 
  GOTO 101
317
 
100  CONTINUE
318
 
  FOUNDLOOPFILTER=.FALSE.
319
 
  DO J=1,NCOMB
320
 
    DO I=NCTAMPS+1,NLOOPAMPS
321
 
          GOODAMP(I,J)=(.NOT.USELOOPFILTER)
322
 
        ENDDO
323
 
  ENDDO
324
 
101  CONTINUE
325
 
CLOSE(1)
326
 
OPEN(1, FILE=HelFilterFN, err=102, status='OLD',           action='READ')
327
 
  READ(1,*,END=103) (GOODHEL(I),I=1,NCOMB)
328
 
  GOTO 103
329
 
102  CONTINUE
330
 
  FOUNDHELFILTER=.FALSE.
331
 
  DO J=1,NCOMB
332
 
        GOODHEL(J)=.TRUE.
333
 
  ENDDO
334
 
103  CONTINUE
335
 
CLOSE(1)
 
302
 
336
303
OPEN(1, FILE=ColorNumFN, err=104, status='OLD',           action='READ')
337
304
  DO I=1,NLOOPAMPS
338
305
    READ(1,*,END=105) (CF_N(I,J),J=1,%(color_matrix_size)s)
364
331
  WRITE(*,*) 'Stopped by user request.'
365
332
  STOP
366
333
ENDIF
 
334
 
 
335
IF(NTRY.EQ.0) THEN
 
336
  CALL %(proc_prefix)sSET_N_EVALS(N_DP_EVAL,N_QP_EVAL)
 
337
  HELDOUBLECHECKED=(.NOT.DoubleCheckHelicityFilter).OR.(HelicityFilterLevel.eq.0)
 
338
  DO J=1,NCOMB
 
339
    DO I=1,NCTAMPS
 
340
          GOODAMP(I,J)=.TRUE.
 
341
        ENDDO
 
342
  ENDDO
 
343
OPEN(1, FILE=LoopFilterFN, err=100, status='OLD',           action='READ')
 
344
  DO J=1,NCOMB
 
345
    READ(1,*,END=101) (GOODAMP(I,J),I=NCTAMPS+1,NLOOPAMPS)
 
346
  ENDDO
 
347
  GOTO 101
 
348
100  CONTINUE
 
349
  FOUNDLOOPFILTER=.FALSE.
 
350
  DO J=1,NCOMB
 
351
    DO I=NCTAMPS+1,NLOOPAMPS
 
352
          GOODAMP(I,J)=(.NOT.USELOOPFILTER)
 
353
        ENDDO
 
354
  ENDDO
 
355
101  CONTINUE
 
356
CLOSE(1)
 
357
IF (HelicityFilterLevel.eq.0) then
 
358
  FOUNDHELFILTER=.TRUE.
 
359
  DO J=1,NCOMB
 
360
        GOODHEL(J)=.TRUE.
 
361
  ENDDO
 
362
  GOTO 122
 
363
ENDIF
 
364
OPEN(1, FILE=HelFilterFN, err=102, status='OLD',           action='READ')
 
365
  READ(1,*,END=103) (GOODHEL(I),I=1,NCOMB)
 
366
  GOTO 103
 
367
102  CONTINUE
 
368
  FOUNDHELFILTER=.FALSE.
 
369
  DO J=1,NCOMB
 
370
        GOODHEL(J)=.TRUE.
 
371
  ENDDO
 
372
103  CONTINUE
 
373
CLOSE(1)
 
374
122  CONTINUE
367
375
ENDIF
368
376
 
369
377
MP_DONE=.FALSE.