bin/BCVEGPY/updatafile.F

0001       subroutine updatafile
0002       implicit double precision(a-h, o-z)
0003         implicit integer(i-n)
0004
0005 c...data files.

0006 c...record usual data: pt--bc p_t; rap--bc rapidity; shat--

0007 c...the subprocess c.m.energy; y--rapidity but equal to

0008 c...ln((x1-x2)/(x1+x2)); pseta--bc pseudo-rapidity;

0009 c...z--(2(k1+k2).p_bc)/(square of the subprocess energy);

0010
0011         open(unit=4,file='data/pt.dat',status='unknown')
0012         open(unit=5,file='data/rap.dat',status='unknown')
0013         open(unit=16,file='data/shat.dat',status='unknown')
0014         open(unit=7,file='data/y.dat',status='unknown')
0015         open(unit=8,file='data/pseta.dat',status='unknown')
0016       open(unit=9,file='data/z.dat',status='unknown')
0017
0018 c...record y-result with different pt cut: 5,...,100 gev.

0019       open(unit=21,file='data/pt005y.dat',status='unknown')
0020       open(unit=22,file='data/pt020y.dat',status='unknown')
0021       open(unit=23,file='data/pt035y.dat',status='unknown')
0022       open(unit=24,file='data/pt050y.dat',status='unknown')
0023       open(unit=25,file='data/pt100y.dat',status='unknown')
0024
0025 c...record p_t result with different y cut: |y|<0.5,...,2.5

0026         open(unit=31,file='data/y05pt.dat',status='unknown')
0027         open(unit=32,file='data/y10pt.dat',status='unknown')
0028         open(unit=33,file='data/y15pt.dat',status='unknown')
0029         open(unit=34,file='data/y20pt.dat',status='unknown')
0030         open(unit=35,file='data/y25pt.dat',status='unknown')
0031
0032 c...record shat-result with different y cut: |y|<0.5,...,2.5

0033         open(unit=41,file='data/y05shat.dat',status='unknown')
0034         open(unit=42,file='data/y10shat.dat',status='unknown')
0035         open(unit=43,file='data/y15shat.dat',status='unknown')
0036         open(unit=44,file='data/y20shat.dat',status='unknown')
0037         open(unit=45,file='data/y25shat.dat',status='unknown')
0038
0039 c...record shat-result with different p_t cut: p_t<5,...,100 gev

0040         open(unit=51,file='data/pt005shat.dat',status='unknown')
0041         open(unit=52,file='data/pt020shat.dat',status='unknown')
0042         open(unit=53,file='data/pt035shat.dat',status='unknown')
0043         open(unit=54,file='data/pt050shat.dat',status='unknown')
0044         open(unit=55,file='data/pt100shat.dat',status='unknown')
0045
0046 c...ya means |y|.

0047 c...record p_t result with different y cut: 0.5<|y|<1.0,...

0048         open(unit=62,file='data/ya05-10pt.dat',status='unknown')
0049         open(unit=63,file='data/ya10-15pt.dat',status='unknown')
0050         open(unit=64,file='data/ya15-20pt.dat',status='unknown')
0051         open(unit=65,file='data/ya20-25pt.dat',status='unknown')
0052
0053         end