1 / 8

Radionuclides in TGeo

Radionuclides in TGeo. M.Gheata 28/09/06. Radionuclides table. Added support for radionuclides as elements that may compose a TGeo material/mixture: TGeoElementRN class Database of ~3500 radionuclides connected via decay channels ( NUBASE2003 ) Table of radionuclides loaded on demand:

fbarnum
Download Presentation

Radionuclides in TGeo

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Radionuclides in TGeo M.Gheata 28/09/06

  2. Radionuclides table • Added support for radionuclides as elements that may compose a TGeo material/mixture: TGeoElementRN class • Database of ~3500 radionuclides connected via decay channels (NUBASE2003) • Table of radionuclides loaded on demand: • gGeoManager->GetElementTable()->GetElementRN(Int_t A, Int_t Z, Int_t iso=0) ; iso = isomeric level

  3. Radionuclides and their decays • Nuclide properties: A, Z, Half-life, mass excess, isomeric number and level, natural abundency, toxicity • Decay channels: TGeoDecayChannel objects in a list stored by each radionuclide • alpha, beta decays, neutron/proton emission, electron capture, isomeric transitions, spontaneous fission • Decay properties: decay type, branching ratio, Q value, parent and daughter elements, delta iso.

  4. Decay chains • Having the decay chain: P->Q->R->S->... the time evolution of the numbers of nuclei of each type is given by the Bateman equations: • The general solution: • The coefficients aij depend on li and branching ratios • Computed only for the case P0≠0, Q0=R0=…=0 • General case – superposition of solutions as above • A functionality of TGeoElementRN class is to find the evolution of the population resulting from an initial radionuclide: • TGeoElementRN::FillPopulation(TObjArray *pop, Double_t precision, Double_t normalization=1.) • The method will fill the input array with all elements in the decay tree and attach to each element a TGeoBatemanSol object

  5. Radioactive materials and mixtures • Can be made in the same way as for normal ones, but using radionuclides as components: TGeoMaterial *mat = new TGeoMaterial("C14", pElemC14, density); TGeoMixture *mix = new TGeoMixture("mix", nelem, density); mix->AddElement(pElem, weigth); ... • Time evolution of the components of a material/mixture: FillMaterialEvolution(TObjArray *pop, Double_t precision) • Material/mixture evolved at a given moment: DecayMaterial(Double_t time, Double_t precision)

  6. Example of usage • See $ROOTSYS/tutorials/RadioNuclides.C Material evolution

  7. Example (cont) Mixture evolution

  8. Conclusion • Support for radionuclides added in TGeo • Possibility to create materials and mixtures based on these • Evolution in time according Bateman equations • Possibility to produce a material derived from a radioactive one at a moment in time • Possibility to compute activity for a given decay (API will be provided soon)

More Related