-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathforpaw.f
61 lines (40 loc) · 1.3 KB
/
forpaw.f
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
*CMZ : 27/05/2004 16.27.44 by H1 Saclay
*-- Author : Unknown 17/12/2003
*--------------------------------
subroutine ntinit
*--------------------------------
include 'forpaw.common'
* Initialization of ntuple and histos
REAL hmemor
PARAMETER( NWPAWC=500000 )
COMMON/pawc/hmemor(NWPAWC)
* ----------------------------------------------------------
*---Create the Ntuple:
IDNT = 1
CALL HBNT(IDNT,'essai',' ')
*---Define the ntuple :
NTID = 1
call hbname(NTID,'dvcsvar1',xntp,'xntp,'//
+ 'qntp,yntp,phintp,tntp,'//
+ 'plilab(4),plolab(4),ppilab(4),ppolab(4),'//
+ 'pvglab(4),prglab(4),pvisr(4),'//
+ 'pliprf(4),ploprf(4),ppiprf(4),ppoprf(4),'//
+ 'pvgprf(4),prgprf(4),'//
+ 'plibel(4),plobel(4),ppibel(4),ppobel(4),'//
+ 'pvgbel(4),prgbel(4),'//
+ 'wbhntp,wdvcsntp,wintntp,wtotntp,'//
+ 'ym,phibel,resol,phibelrec,phibelgen')
return
end
*---------------------------------------------
subroutine ntend
*---------------------------------------------
* close the file and write out ntuple and histo :
* call hrout(1,icycle,' ')
* call hrout(500,icycle,' ')
call hrout(0,icycle,' ')
call hrend('toto')
LUNNT = 31
close(LUNNT)
return
end