11 if (idhist.gt.maxhist)
return
13 if (histo(idhist)) call hdelet(idhist)
16 +
title =
'Edep in absorber '//char(ichar(
'0')+(idhist))
17 if ((idhist.gt.maxabs).and.(idhist.lt.2*maxabs))
18 +
title =
'Edep long. profile in absorber '
19 + //char(ichar(
'0')+(idhist-maxabs))//
' (MeV/event)'
20 if (idhist.eq.(2*maxabs+1))
21 +
title =
'energy flow (MeV/event)'
22 if (idhist.eq.(2*maxabs+2))
23 +
title =
'lateral energy leak (MeV/event)'
25 call hbook1(idhist,
title,nbbins,valmin,valmax,0.)
27 histo(idhist) = .
true.
28 binwidth(idhist) = (valmax-valmin)/nbbins
29 if (valunit.le.0.) valunit = 1.
30 histunit(idhist) = valunit