← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to Template/SubProcesses/initcluster.f

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/SubProcesses/initcluster.f
 
1
      subroutine initcluster()
 
2
 
 
3
      implicit none
 
4
 
 
5
      include 'message.inc'
 
6
      include 'run.inc'
 
7
      include 'nexternal.inc'
 
8
      include 'cluster.inc'
 
9
C
 
10
C     SPECIAL CUTS
 
11
C
 
12
      real*8 xptj,xptb,xpta,xptl,xmtc
 
13
      real*8 xetamin,xqcut,deltaeta
 
14
      common /to_specxpt/xptj,xptb,xpta,xptl,xmtc,xetamin,xqcut,deltaeta
 
15
 
 
16
      integer i,j
 
17
      logical filmap, cluster
 
18
      external filmap, cluster
 
19
 
 
20
c     
 
21
c     check whether y_cut is used -> set scale to y_cut*S
 
22
c
 
23
 
 
24
c      if (ickkw.le.0) return
 
25
      if (ickkw.le.0.and.xqcut.le.0d0.and.fixed_ren_scale.and.fixed_fac_scale) return
 
26
 
 
27
c      if(ickkw.eq.2.and.xqcut.le.0d0)then
 
28
c        write(*,*)'Must set qcut > 0 for ickkw = 2'
 
29
c        write(*,*)'Exiting...'
 
30
c        stop
 
31
c      endif
 
32
 
 
33
c      if(xqcut.gt.0d0)then
 
34
c      if(ickkw.eq.2)then
 
35
c        scale = xqcut
 
36
c        q2fact(1) = scale**2    ! fact scale**2 for pdf1
 
37
c        q2fact(2) = scale**2    ! fact scale**2 for pdf2
 
38
c        fixed_ren_scale=.true.
 
39
c        fixed_fac_scale=.true.
 
40
c      endif
 
41
c   
 
42
c     initialize clustering map
 
43
c         
 
44
      if (.not.filmap()) then
 
45
        write(*,*) 'cuts.f: cluster map initialization failed'
 
46
        stop
 
47
      endif
 
48
      if (btest(mlevel,3)) then
 
49
        do i=1,ishft(1,nexternal+1)
 
50
          write(*,*) 'prop ',i,' in'
 
51
          do j=1,id_cl(i,0)
 
52
            write(*,*) '  graph ',id_cl(i,j)
 
53
          enddo
 
54
        enddo
 
55
        write(*,*)'ok'
 
56
      endif
 
57
      igscl(0)=0
 
58
       
 
59
      RETURN
 
60
      END
 
61