Re: VOXEL GEOMETRY (with attachment)

tp1.mi.infn.it [192.84.138.69]); Fri, 16 Apr 2010 18:46:24 +0200 (CEST)
Sender: owner-fluka-discuss_at_mi.infn.it

--8323328-1499803886-1271433471=3D:6764
Content-Type: text/plain; charset=3D"ISO-8859-7"; format=3Dflowed
Content-Transfer-Encoding: 8BIT
Content-ID: <alpine.LFD.2.00.1004161800021.6764_at_localhost.localdomain>

Dear Antonello,

may be I misunderstand your actual needs, but I wonder if indeed what you=
=20
want is just to get an energy deposition map over your ionization chamber.=
=20
If this is the case, you have simply to input an USRBIN detector, defining=
=20
a (cartesian or cylindrical) grid which will give you the energy=20
deposition density [GeV/cm3] averaged over each cell (if this=20
size-adjustable scoring cell is what you call voxel).
On the other hand, the voxel geometry capability of FLUKA is intended to=20
implement a density/material map, typically coming from a CT scan and=20
imported after suitable pre-processing by the fortran tool you attached.=20
If this is really your aim, then I suggest you take a careful look at the=
=20
beginner course dedicated lecture, available at=20
https://www.fluka.org/free_download/course/mumbai2009/Lectures/Voxels_and_m=
ed_App1009.pdf
Of course in that case you realize that you must have an input "Hounsfield=
=20
Unit" matrix.

Hope this helps

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Fri, 16 Apr 2010, ANTONELLO CARLONI wrote:

>=20
>
> Dear Fluka team,
> sorry to bother again you with stupid questions but as I told you
> I am really a beginner in programming.
> I need to subdvid into voxels a ionization chamber and then
> detect the release of energy on it by a proton beam.
> On the fluka manual i read that I must write down compile and
> then link the Writegolem program that i attached herewithin.
> My questions are:
> #1: What do the constants IREG and KREG stand for and, therefore,
> what kind of value should i give them.....
> #2: How do I have to change the program since i do not have a
> matrix of values in input such as the one named
> 'ascii_segm_golem' which is being called to be read in the
> command:
> OPEN(UNIT=3D30,FILE=3D=A2ascii_segm_golem=A2,STATUS=3D=A2OLD=A2)
> !?
> Thank you so much in advance for your kind assistance,
> Antonello Carloni
>=20
>=20
>=20
>
--8323328-1499803886-1271433471=3D:6764--

-------
Received on Sat Apr 17 2010 - 00:20:24 CEST