5 #include "../FleshHist/Zprivate1.h" 6 #include "../FleshHist/Zprivate3.h" 8 real*8 Et, wx, wy, wz, firstdep
14 write(0,*)
' file name is ' 16 open(fnodat, file=file, form=
"unformatted")
32 read(fnodat,end=100) bufc, (buf(
i),
i=1, bufc)
36 *
'(6i3, 1pE11.3, 0p,f6.1,1p2E11.3,0p, 2f8.4,f10.6)')
38 * buf(
i).
charge, buf(
i).ridx, buf(
i).faiidx,
39 * buf(
i).rinmu, buf(
i).fai, buf(
i).ek,
40 * buf(
i).
t, buf(
i).wx, buf(
i).wy, buf(
i).wz
43 *
'(6i3, 1pE11.3, 0p,f6.1,1p2E11.3,0p, 2f8.4,f10.6,1pE11.3)')
45 * buf(
i).
charge, buf(
i).ridx, buf(
i).faiidx,
46 * buf(
i).rinmu, buf(
i).fai, buf(
i).ek,
47 * buf(
i).
t, buf(
i).wx, buf(
i).wy, buf(
i).wz,
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
*Zfirst p fm *Zfirst p Zfirst p code
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
*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
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode