[Gate-users] GateDoseSpectrumActor

[David Sarrut]

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/76029154/attachment.html>