[Gate-users] Fluence Actor

[Simon Rit]

1/ You miss line 249:
GateScatterOrderTrackInformation * info =
dynamic_cast<GateScatterOrderTrackInformation
*>(step->GetTrack()->GetUserInformation());
So it needs to have scatter order information to compute the images. To do
that, use the ScatterOrderTrackInformationActor.

The ScatterOrderActor was outdated, I have made a small change:
https://github.com/OpenGATE/Gate/commit/700a0ef1b5b5cf58776302d198f15ea9cec22017

I have then modified the
examples/example_CT/fixedForcedDetectionCT/mac/monteCarloCT.mac to
illustrate how to get the scatter per order:
https://github.com/OpenGATE/Gate/commit/914b106477ee40d32491fd9b99b62dbcabd49d41
You still have to uncomment the
/gate/actor/DetectorFluence/scatterOrderFilename
to make it happen.

2/ No, I don't know a better solution.

On Thu, Dec 1, 2016 at 3:25 PM, Triltsch, Nicolas <nicolas.triltsch at tum.de>
wrote:

> Hey Simon,
>
> ----1----
>
> Now I have time to pick up this topic. Honestly, I didn't really
> understand your reply. I see that GateScatterOrderTrackInformation.hh is
> a package which is included in the GateFluenceActor.cc file. I thought by
> setting mSeparateProcessFilename='scatter' in that file, I might get the
> different *rayleigh* and *compton scatter images*, because it states in
> lines 182-200:
>
> -   /* Printing compton or rayleigh or fluorescence scatter images */
> -   if (*mSeparateProcessFilename* != "")
> -     {
> -     /* Saving separately process images (e.g. Compton, Rayleigh...) */
> -     std::map<G4String, GateImage*>::iterator it = mProcesses.end();
> -     for (unsigned int i = 0; i < mProcessName.size(); i++)
> -       {
> -       it = mProcesses.find(mProcessName[i]);
> -       if (it != mProcesses.end())
> -         {
> -         std::stringstream filenamestream;
> -         filenamestream << mSeparateProcessFilename << "_" << mProcessName[i] << ".mhd";
> -         sprintf(filename, filenamestream.str().c_str(), rID);
> -         mProcesses[mProcessName[i]]->Write((G4String) filename);
> -         }
> -       it = mProcesses.end();
> -       }
> -     }
> -   }
>
> But that didn't seem to work out. What can I do to get the separate
> scatter images? Or am I supposed to use the GateScatterOrderTrackInformationActor
> to obtain that information. I am stuck here, as the .cc scripting is quite
> difficult to understand for me.
>
> ---2---
>
> Is there a way to obtain the* flatfield* with the Fluence Actor while
> simulating with a phantom? I need the flatfield to normalize the images
> later. What I am doing right now, I remove the phantom and start the
> simulation again. It works, but takes two times the simulation time and I
> have to run the script twice.
>
>
> Thanks again and have a good day!
>
> Nico
>
>
>
> On 11/21/2016 07:48 AM, Simon Rit wrote:
>
> Hi,
> I think that this is meant to be used in combination with
> GateScatterOrderTrackInformation. Every time an interaction occurs in the
> volume to which GateScatterOrderTrackInformation is attached, this
> information will be added to a user structure of the particle track. When
> recording the fluence with the GateFluenceActor, it will check the same
> structure of the particle and record the particle in separate images. This
> has been developped to measure the amount of Compton vs Rayleigh scatter
> and I think it currently works for photons only since only Compton,
> Rayleigh and Fluorescence are tracked.
> I hope this helps,
> Simon
>
> On Sat, Nov 19, 2016 at 7:17 PM, Triltsch, Nicolas <
> <nicolas.triltsch at tum.de>nicolas.triltsch at tum.de> wrote:
>
>> Hi everyone,
>>
>> I really enjoy this support, as it was a great help so far. A new
>> question poped up my mind when dealing with the so called FluenceActor.
>>
>> As it is written in the Wiki, the option enableScatter prints out the
>> .mhd image for scattering. I assume that Compton and Rayleigh scatter are
>> combined here. But it should be possible to get the two scatter images
>> separate using the command
>> */gate/actor/DetectorFluence/separateProcessFilename*. In the Wiki it
>> says to do the following:
>>
>> "Set the file name for the different scatter x-rays processes that hit
>> the detector (printf format with runId as a single parameter)"
>>
>> But I don't understand that! Any help would be great?
>>
>> Cheers,
>>
>> Nico
>>
>>
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20161205/15534342/attachment-0001.html>