← Back to branch summary

~maddevelopers/mg5amcnlo/WWW5_caching

« back to all changes in this revision

Viewing changes to users/mardelcourt/PROC_129738/PROC_129738/SubProcesses/initcluster.f

  • Committer: John Doe
  • Date: 2013-03-25 20:27:02 UTC
  • Revision ID: john.doe@gmail.com-20130325202702-5sk3t1r8h33ca4p4
first clean version

Show diffs side-by-side

added added

removed removed

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