[Gate-users] GateDoseSpectrumActor

Maikol Salas Ramirez mmsalas at gmail.com
Tue Feb 1 10:17:07 CET 2022


Hello David,

Thank you for the answer. I re-check my calculation, and I found out
that the units are well (I mean the EnergySpectrumActor and
GateDoseSpectrumActor agree in the deposited energy calculations).

To check it, I attached both actors to the same logical volume and
calculated the deposited energy for comparison. In the case of
EnergySpectrumActor, I summed the energy spectrum, and for
GateDoseSpectrumActor, it provides the total energy (sum) plus the
uncertainty directly.

The GateDoseSpectrumActor has the advantage that it provides the total
deposited energy with its uncertainty. This information is impossible to
obtain with the EnergySpectrumActor, because it does not give the square of
the sum of the energies.

I think the GateDoseSpectrumActor could be useful for the community. But,
it is not documented in the opengate.readthedocs.io, and the C++ classes
are probably incomplete (probably the goal was the absorbed dose
calculation, but its development stops with the energy calculation). At
this moment, the GateDoseSpectrumActor help me resolve my
task straightforwardly. I think it is better to keep it. Hopefully, it will
be a chance to improve it in the future.

I am finalizing some details of the simulation (blood irradiation with
radionuclides); if it would interest you, I can share it for the GATE
examples with some documentation.

Best regards,
Maikol

El mar, 1 feb 2022 a las 8:05, David Sarrut (<
david.sarrut at creatis.insa-lyon.fr>) escribió:

> Hello Maikol,
>
> I don't know this actor and, looking at the history, it was started before
> we started the git version (2015), so the original author is lost.
> As there is no doc and I don't really get what it means compared to the
> EnergySpectrumActor, I would suggest to remove it.
>
> However, if you think it could be useful for you, please let us know. The
> idea would be to propose some docs and use case.
>
> thanks,
> David
>
> On Mon, Jan 31, 2022 at 5:50 PM Maikol Salas Ramirez <mmsalas at gmail.com>
> wrote:
>
>> Dear Gate users,
>>
>> Recently I found the .cc and .hh files of a *GateDoseSpectrumActor*.
>> This actor is not mentioned in opengate.readthedocs.io or any webpage.
>> I would like to ask if someone is working on it. Is it a prototype?
>>
>> I tested the energy calculation using the EnergySpectrumActor and the
>> energy units differ:
>> I simulated a Ac-225 source using the Ion source definition. I scored the
>> deposited energy with both actors (EnergySpectrumActor and
>> GateDoseSpectrumActor). The deposited energy differs by a factor of 1000. *The
>> only justification that I can find is that, in the GateDoseSpectrumActor,
>> the deposited energy is expressed in keV/(1000 events) OR (eV/event) and
>> there is a typo in the code (the units shoud be eV instead of keV):*
>> #-------------------------------------------
>> /// Save data
>> void GateDoseSpectrumActor::SaveData()
>> {
>>   GateVActor::SaveData();
>>   std::ofstream DoseResponseFile;
>>   OpenFileOutput(mSaveFilename, DoseResponseFile);
>>   DoseResponseFile << "Incoming energy (*keV*)" << " "
>>                    << "Average energy deposit (*keV*)" << " "
>>                    << "Energy uncertainty (*keV*)"  << Gateendl;
>>   std::map< G4double, G4double>::iterator itermEnergy;
>>   for( itermEnergy = mEnergy.begin(); itermEnergy != mEnergy.end();
>> itermEnergy++)
>>   {
>>     G4double energyIn = itermEnergy->first;
>>     G4double totEnergyOut = mEnergy[energyIn];
>>     G4double numPart = mNumParticPerEnergy[energyIn];
>>     G4double avgEnergyOut = totEnergyOut / numPart;
>>     G4double avgEnergySq  = mEnergySquare[energyIn] / numPart;
>>     G4double error = sqrt( (1.0 / ( numPart - 1)) * ( avgEnergySq - pow(
>> avgEnergyOut, 2)));
>>     DoseResponseFile << energyIn/*keV* << " "
>>                      << avgEnergyOut/*keV* << " "
>>                      << error/*keV*  << Gateendl;
>>   }
>> #-------------------------------------------
>>
>> Is anyone working on this?
>>
>>
>> Best regards
>> Maikol
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
>
>
> --
> David Sarrut, Phd
> Directeur de recherche CNRS
> CREATIS, UMR CNRS 5220, Inserm U1294
> Centre de lutte contre le cancer Léon Bérard
> 28 rue Laënnec, 69373 Lyon cedex 08
> Tel : 04 78 78 51 51 / 06 74 72 05 42
> http://dsarrut.github.io
> _________________________________
>  "2 + 2 = 5,  for extremely large values of 2"
> _________________________________
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220201/d3f0f826/attachment-0001.html>


More information about the Gate-users mailing list