From: Alfredo Ferrari (alfredo.ferrari@cern.ch)
Date: Mon Feb 13 2006 - 17:20:38 CET
Arggghhhh!!!
Directing particles at the center is infinitely wrong!!! Indeed you should
get infinite as a result. To get an isotropic source requires an algorithm
emitting from a spherical surface (random point) in a direction which
points inside and which forms an angle with the local normal
distributed as cos(theta) dcos(theta). The routine sflood included in
fluka does just this for a 1 cm radius (multiply by the radius you wish
the position components).
The calling sequence is:
SUBROUTINE SFLOOD ( XXX, YYY, ZZZ, UXXX, VYYY, WZZZ )
where XXX,YYY,ZZZ is the returned position (on a 1 cm radius sphere, again
multiply all three for the radius you wish) and UXXX, VYYY, WZZZ is the
direction
The fluence you generate in this way will be uniform and isotropic
everywhere inside the sphere (of course if there is no material altering
it) and numerically equal to 1/ (pi R**2) where R is the radius you used.
Ciao
Alfredo
On Mon, 13 Feb 2006, Reddell, Brandon D wrote:
> This is a multi-part message in MIME format.
>
> ------_=_NextPart_001_01C630B1.5AAE113B
> Content-Type: text/plain;
> charset="us-ascii"
> Content-Transfer-Encoding: quoted-printable
>
> Dear FLUKA users,
>
> I have a very fundamental and rather simple question that nobody seems
> to be able to answer. How do you correctly determine dose at the center
> of a solid aluminum sphere. I want to compare FLUKA 3D results with
> results from others codes, such as Shieldose2 and GEANT4.
>
> I've included a test input file as an example. The results from SCORE
> and USRBIN 1) don't agree with each other, then 2) they aren't even
> close to that from the other codes. Ultimately, I would like to get the
> dose at the center of a sphere from an isotropic source/spectrum (hence
> my source command at the bottom). For now, I am directing all energy to
> the center -- the answer should be relatively close to that from an
> isotropic source.
>
> Can any body help me determine dose in Rad or Gy at the center of a
> solid sphere --- and provide an answer that consistent with what is
> published in literature?
>
> Sincerely,
>
> Brandon
>
> <<test1.inp>>=20
>
> ------_=_NextPart_001_01C630B1.5AAE113B
> Content-Type: application/octet-stream;
> name="test1.inp"
> Content-Transfer-Encoding: base64
> Content-Description: test1.inp
> Content-Disposition: attachment;
> filename="test1.inp"
>
> VElUTEUNCnNwaGVyZQ0KREVGQVVMVFMgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICAgICAgICAgIFBSRUNJU0lPDQoqMjM0NTY3ODkgMTIzNDU2Nzg5
> IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2
> Nzg5DQpCRUFNICAgICAgLTAuNSAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgUFJPVE9ODQpCRUFNUE9TICAgMC4wMDAwICAgIDAuMDAwMCAgICAt
> MC4wNTAwDQpHRU9CRUdJTiAgICAgICAgICAgIDAuMDEgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICAgQ09NQklOQVQNCioyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2
> Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkNCiAg
> ICAgICAgICAgICAgICAgICAgICAgICAgIGEgc3BoZXJlIA0KICBTUEggICAgMSAgICAgICAwLjAg
> ICAgICAgMC4wICAgICAgIDAuMCAgKzIwMDAwLjANCiAgU1BIICAgIDIgICAgICAgMC4wICAgICAg
> IDAuMCAgICAgICAwLjAgICArOTAwMC4wIA0KICBTUEggICAgMyAgICAgICAwLjAgICAgICAgMC4w
> ICAgICAgIDAuMCAgIDAuNzYwNTcgIA0KICBSUFAgICAgNCAgICAtMC4wMDUgICAgKzAuMDA1ICAg
> IC0wLjAwNSAgICArMC4wMDUgICAgLTAuMDA1ICAgICswLjAwNSAgDQogIEVORA0KKiBibGFjayBo
> b2xlDQogICAgMSAgICA2ICAgICAxICAgIC0yDQoqIHZhY3V1bQ0KICAgIDIgICAgNiAgICAgMiAg
> ICAtMw0KKiBBbHVtaW51bSBMYXllciANCiAgICAzICAgIDYgICAgIDMgICAgLTQgICAgIA0KKiBE
> ZXRlY3RvciANCiAgICA0ICAgIDYgICAgIDQNCiAgRU5EDQpHRU9FTkQgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgREVCVUcgDQoq
> MjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4
> OSAxMjM0NTY3ODkgMTIzNDU2Nzg5DQoqIEVsZW1lbnRzDQpNQVRFUklBTCAgMTMuMCAgICAgIDI2
> Ljk4MCAgICAyLjYyOTU5ICAgMTAuMCAgICAgICAgICAgICAgICAgICAgICAgICAgQUxVTUlOVU0N
> CioNCioyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIz
> NDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkNCiogIEV4dGVybmFsIEJsYWNrIEhvbGUNCkFTU0lH
> Tk1BVCAxLjAgICAgICAgIDEuMA0KKiAgVmFjdXVtDQpBU1NJR05NQVQgMi4wICAgICAgICAyLjAN
> CiogIEFsdW1pbnVtDQpBU1NJR05NQVQgMTAuMCAgICAgICAzLjANCkFTU0lHTk1BVCAxMC4wICAg
> ICAgIDQuMA0KKg0KTE9XLU5FVVQgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgMC4w
> DQpFTUYgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgDQpFTUZDVVQgICAgICAgICAgIDAuMCAgICAgICAwLjAgICAgICAgMC4w
> ICAgICAgIDMuMCAgICAgICA0LjAgICAgICAgICAgUFJPRC1DVVQNCioyMzQ1Njc4OSAxMjM0NTY3
> ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0
> NTY3ODkNClNDT1JFICAgICAyMDguMCAgICAgMjExLjAgICAgIDIwMS4wDQoqICoqKiogU2Nvcmlu
> ZyAqKioqDQoqMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5
> IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5DQpVU1JCSU4gICAgICAgICAxMC4wICAgICAy
> MDguMCAgICAgLTUwLjAgICAgIDAuMDA1ICAgICAwLjAwNSAgICAgMC4wMDUgZGV0ZWN0b3INClVT
> UkJJTiAgICAgICAtMC4wMDUgICAgLTAuMDA1ICAgIC0wLjAwNSAgICAgIDIwLjAgICAgICAyMC4w
> ICAgICAgMjAuMCAmDQoqDQpSQU5ET01JWkUgIDEuICAgICAgIDAuMTA2MkUrMDQNCioyMzQ1Njc4
> OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1Njc4OSAxMjM0NTY3ODkgMTIzNDU2Nzg5IDEyMzQ1
> Njc4OSAxMjM0NTY3ODkNClNPVVJDRSAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICANClNUQVJUICAgICAwLjEwMDBFKzA3DQpTVE9QDQo=
>
> ------_=_NextPart_001_01C630B1.5AAE113B--
>
-- +----------------------------------------------------------------------------+ | Alfredo Ferrari || Tel.: +41.22.767.6119 | | CERN-AB || Fax.: +41.22.767.7555 | | 1211 Geneva 23 || e-mail: Alfredo.Ferrari@cern.ch | | Switzerland || Alfredo.Ferrari@mi.infn.it | +----------------------------------------------------------------------------+
This archive was generated by hypermail 2.1.6 : Mon Feb 13 2006 - 23:29:48 CET