COSMOS v7.655  COSMOSv7655
(AirShowerMC)
readBin.f
Go to the documentation of this file.
1  include '../FleshHist/Zprivate1.h'
2  character*120 file
3  real*8 Et, wx, wy, wz
4  integer EventNo
5  integer*2 code
6 
7  read(*,'(a)') file
8  write(0,*) ' file name is '
9  write(0,*) file
10  open(fnodat, file=file, form="unformatted")
11  read(fnodat)
12  * eventno, code, et, wx, wy, wz
13  write(*,'("i ", i3, i4, g13.4,3f11.7)')
14  * eventno, code, et, wx, wy, wz
15 
16  do while(.true.)
17  read(fnodat,end=100) bufc, buf
18  do i = 1, bufc
19 #if KeepWeight != yes
20  write(*,
21  * '(6i3, 1pE11.3, 0p,f6.1,1p2E11.3,0p, 2f8.4,f10.6)')
22  * buf(i).ldep, buf(i).code, buf(i).subcode,
23  * buf(i).charge, buf(i).ridx, buf(i).faiidx,
24  * buf(i).rinmu, buf(i).fai, buf(i).ek,
25  * buf(i).t, buf(i).wx, buf(i).wy, buf(i).wz
26 #else
27  write(*,
28  * '(6i3, 1pE11.3, 0p,f6.1,1p2E11.3,0p, 2f8.4,f10.6,1pE11.3')
29  * buf(i).ldep, buf(i).code, buf(i).subcode,
30  * buf(i).charge, buf(i).ridx, buf(i).faiidx,
31  * buf(i).rinmu, buf(i).fai, buf(i).ek,
32  * buf(i).t, buf(i).wx, buf(i).wy, buf(i).wz,
33  * buf(i).wgt
34 #endif
35 
36  enddo
37  enddo
38  100 continue
39  end
i
nodes i
Definition: cstdatmos0.fNew2.h:93
true
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7
code
*Zfirst p fm *Zfirst p Zfirst p code
Definition: ZavoidUnionMap.h:1
charge
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
Definition: ZavoidUnionMap.h:1
t
nodes t
Definition: cstdatmos0.fNew2.h:84
wgt
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst wgt
Definition: ZavoidUnionMap.h:1
subcode
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode
Definition: ZavoidUnionMap.h:1