<div dir="ltr">Hello Maikol, <div><br></div><div>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. </div><div>As there is no doc and I don't really get what it means compared to the EnergySpectrumActor, I would suggest to remove it. </div><div><br></div><div>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.</div><div><br></div><div>thanks,</div><div>David</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 31, 2022 at 5:50 PM Maikol Salas Ramirez <<a href="mailto:mmsalas@gmail.com">mmsalas@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Gate users,</div><div><br></div><div>Recently I found the .cc and .hh files of a <b>GateDoseSpectrumActor</b>. This actor is not mentioned in <a href="http://opengate.readthedocs.io" target="_blank">opengate.readthedocs.io</a> or any webpage.</div><div>I would like to ask if someone is working on it. Is it a prototype?</div><div><br></div><div>I tested the energy calculation using the EnergySpectrumActor and the energy units differ:</div><div></div><div>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. <b>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):</b></div><div>#-------------------------------------------<br></div><div>/// Save data<br>void GateDoseSpectrumActor::SaveData()<br>{<br> GateVActor::SaveData();<br> std::ofstream DoseResponseFile;<br> OpenFileOutput(mSaveFilename, DoseResponseFile);<br> DoseResponseFile << "Incoming energy (<b>keV</b>)" << " "<br> << "Average energy deposit (<b>keV</b>)" << " "<br> << "Energy uncertainty (<b>keV</b>)" << Gateendl;<br> std::map< G4double, G4double>::iterator itermEnergy;<br> for( itermEnergy = mEnergy.begin(); itermEnergy != mEnergy.end(); itermEnergy++)<br> {<br> G4double energyIn = itermEnergy->first;<br> G4double totEnergyOut = mEnergy[energyIn];<br> G4double numPart = mNumParticPerEnergy[energyIn];<br> G4double avgEnergyOut = totEnergyOut / numPart;<br> G4double avgEnergySq = mEnergySquare[energyIn] / numPart;<br> G4double error = sqrt( (1.0 / ( numPart - 1)) * ( avgEnergySq - pow( avgEnergyOut, 2)));<br> DoseResponseFile << energyIn/<span style="color:rgb(255,0,0)"><b>keV</b> </span><< " "<br> << avgEnergyOut/<span style="color:rgb(255,0,0)"><b>keV</b> </span><< " "<br> << error/<span style="color:rgb(255,0,0)"><b>keV</b> </span> << Gateendl;<br> }</div><div>#-------------------------------------------</div><div><br></div><div>Is anyone working on this?</div><div></div><div><br></div><div><br></div><div>Best regards</div><div>Maikol<br></div><div><br></div></div>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><span style="color:rgb(0,0,0)">David Sarrut, Phd</span><br></div><div><font color="#000000">Directeur de recherche CNRS</font></div><div><font color="#000000">CREATIS, UMR CNRS 5220, Inserm U1294</font></div><div><font color="#000000">Centre de lutte contre le cancer Léon Bérard</font></div><div><font color="#000000">28 rue Laënnec, 69373 Lyon cedex 08</font></div><div><font color="#000000">Tel : 04 78 78 51 51 / 06 74 72 05 42</font></div><div><font color="#000000"><a href="http://dsarrut.github.io" target="_blank">http://dsarrut.github.io</a></font></div><div><font color="#000000">_________________________________</font></div><div><font color="#000000"> "2 + 2 = 5, for extremely large values of 2"</font></div><div><font color="#000000">_________________________________</font></div></div></div>