COSMOS v7.655  COSMOSv7655
(AirShowerMC)
gpinte.f
Go to the documentation of this file.
1  implicit none
2  real(8)::p1(4), p2(4)
3  integer irej, i, j
4  real(8):: eelec, prom, eprot, egam, sigmax, fac
5 
6  integer NMXHEP
7  parameter(nmxhep=4000)
8 
9  INTEGER NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
10  DOUBLE PRECISION PHEP,VHEP
11  COMMON /poevt1/ nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
12  & jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),
13  & vhep(4,nmxhep)
14 ! extension to standard particle data interface (PHOJET specific)
15  INTEGER IMPART,IPHIST,ICOLOR
16  COMMON /poevt2/ impart(nmxhep),iphist(2,nmxhep),icolor(2,nmxhep)
17 
18 
19 
20 ! call pho_inp(-1, irej)
21  call pho_init(-1, irej)
22 
23  prom =0.938
24 ! eprot=820.
25  eprot=prom
26  p1(1)=0
27  p1(2) = 0
28  p1(3) = sqrt(eprot**2 - prom**2)
29  p1(4) = eprot
30 
31 
32  eelec =100
33  egam = 0.45*eelec
34  p2(1) = 0
35  p2(2) = 0
36  p2(3)=-egam
37  p2(4)=egam
38 
39  call pho_setpar(1, 2212, 0, 0.d0)
40  call pho_setpar(2, 22, 0, 0.d0)
41 
42  call pho_event(-1, p1, p2,sigmax, irej)
43 
44  write(0,*) ' sigmax=',sigmax, ' ireg=',irej
45 
46 
47  do i = 1, 10
48 
49  call pho_event(3, p1, p2, fac, irej)
50  write(0,*) ' fac=',fac, ' ireg=',irej
51 
52  write(0,*) ' nhep=',nhep
53  do j = 1, nhep
54  if( isthep(j) == 1 ) then
55  write(0,'(2i5, 1p, 5g13.3)')
56  * isthep(j), idhep(j), phep(:, j)
57  endif
58  enddo
59  enddo
60 
61  call pho_event(-2, p1, p2, sigmax,irej)
62  end
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
nodes i
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5