|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
subroutine vegaslogo(vegasopen)
implicit double precision(a-h, o-z)
implicit integer(i-n)
#include "invegas.h"
#include "bcvegpy_set_par.inc"
common/vegasinf/number,nitmx
logical vegasopen
if(vegasopen) then
write(*,'(a,i5)') 'the used bin number in vegas:',NVEGBIN
write(3,'(a,i5)') 'the used bin number in vegas:',NVEGBIN
write(*,'(a,i9,a,i3)') 'using vegas: number in each iter.='
& ,number,' iter.=',nitmx
write(3,'(a,i9,a,i3)') 'using vegas: number in each iter.='
& ,number,' iter.=',nitmx
write(*,*)
write(3,*)
write(*,'(a)')
& '......................................................'
write(*,'(a)')
& '.............. end of initialization .................'
write(*,'(a)')
& '......................................................'
write(3,'(a)')
& '......................................................'
write(3,'(a)')
& '.............. end of initialization .................'
write(3,'(a)')
& '......................................................'
write(*,*)
write(3,*)
write(*,'(a)') '.....waiting...... vegas running......'
write(3,'(a)') '.....waiting...... vegas running......'
write(*,*)
write(3,*)
end if
end
c*******************************************************************
subroutine vegasend(vegasopen,ievntdis,usegrade)
implicit double precision(a-h, o-z)
implicit integer(i-n)
logical vegasopen,usegrade
common/mixevnt/xbcsec(8),imix,imixtype
if(vegasopen) then
write(*,'(a)') '...... end of vegas, grade generated ......'
write(3,'(a)') '...... end of vegas, grade generated ......'
if(ievntdis.eq.1) then
write(*,'(a)')'ok, to get the evnt number distributions'
write(3,'(a)')'ok, to get the evnt number distributions'
else
write(*,'(a)')'ok, to get the differential distributions'
write(3,'(a)')'ok, to get the differential distributions'
end if
write(*,*)
write(3,*)
else
write(*,'(a)') '...... vegas has not been used ......'
write(3,'(a)') '...... vegas has not been used ......'
if(usegrade) then
write(*,'(a)') 'using existed grade to generated events.....'
write(3,'(a)') 'using existed grade to generated events.....'
write(*,'(a)')'(existed) grade has not been further improved'
write(3,'(a)')'(existed) grade has not been further improved'
else
write(*,'(a)') 'no grade is adopted, using trivial MC.'
write(3,'(a)') 'no grade is adopted, using trivial MC.'
if(imix.eq.1) then
write(*,'(a)') 'not an appreciable way for mixed events'
write(3,'(a)') 'not an appreciable way for mixed events'
write(*,'(a)') 'the precision is lower than other methods'
write(3,'(a)') 'the precision is lower than other methods'
end if
end if
write(*,*)
write(3,*)
write(*,'(a)')
& '......................................................'
write(*,'(a)')
& '.............. end of initialization .................'
write(*,'(a)')
& '......................................................'
write(3,'(a)')
& '......................................................'
write(3,'(a)')
& '.............. end of initialization .................'
write(3,'(a)')
& '......................................................'
write(*,*)
write(3,*)
if(ievntdis.eq.1) then
write(*,'(a)')'ok, to get the evnt number distributions.....'
write(3,'(a)')'ok, to get the evnt number distributions.....'
else
write(*,'(a)')'ok, to get the differential distributions....'
write(3,'(a)')'ok, to get the differential distributions....'
end if
end if
end
|