[Gate-users] Fwd: histogram-source

[Uwe Pietrzyk]

Dear Andrew,

this is an eMail, which I have sent to Johannes Bender in reply
to his similar question concerning ion sources. The name of the
file is in fact GateTrajectoryNavigator.cc. I have attached our version
to this eMail.

But as you said, I was wrong and took the wrong file, as there is
a file GateTrajectoryNavigator.cc also in version 6. So we should just
need a similar correction as for version 5.

Hope this is helpful

Kind regards,

Uwe
-----
Begin forwarded message:

> From: "Pietrzyk, Uwe" <u.pietrzyk at fz-juelich.de>
> Date: 30. Oktober 2009 09:07:47 MEZ
> To: Johannes Bender <JohannesBender at gmx.de>
> Cc: "gate-users at lists.healthgrid.org" <gate-users at lists.healthgrid.org>
> Subject: Re: [Gate-users] histogram-source
> 
> Hi Johannes,
> 
> for the ION-Source to work properly experts informed me, that a part in
> .../source/management/src/GateTrajectory.cc needs to be changed.
> It actually is only the line of code from Jan De Neenhouwer. See below.
> Also note, that depending on the type of ION-source, some of the emitted
> particles simply do not travel very far. Concerning electron capture,
> I cannot tell, whether it is implemented and if not, how to do this.
> The general particle source (gps), however, I would guess should
> have this mode included, but please check in the file that are
> anyhow required by GEANT: RadioactiveDecay3.2 ...
> 
> I hope this is helpful.
> 
> Kind regards,
> 
> Uwe
> 
> 
> void GateTrajectoryNavigator::SetIonID()
> {
>   G4int n_trajectories = m_trajectoryContainer->entries();
>   for (G4int iTrj=0; iTrj<n_trajectories; iTrj++) {
>     G4Trajectory* trj = (G4Trajectory*)((*m_trajectoryContainer)[iTrj]);
>     if (trj->GetTrackID() == 1) {
> //       G4cout << "GateTrajectoryNavigator::SetIonID: PDG " << trj->GetPDGEncoding() 
> // 	     << " Name <" << trj->GetParticleName() << ">" << G4endl;
> //      if (trj->GetPDGEncoding() == 0) { // 0 == Ion -----> to be changed !! 
>       if (trj->GetPDGEncoding() > 1000000000) { // Correction from Jan De Beenhouwer 
> 
> 	m_ionID = 1;
>       } else {
> 	m_ionID = 0;
>       }
>     }
>   }
> }
> 
> 
> 
> 
> On 29.10.2009, at 17:57, Johannes Bender wrote:
> 
>> Hallo,
>> 
>> I tried to check, wether the implemented data for the ion-source is correct.
>> 
>> I tried to used the histogram soucre with data from TOI, but whenever I start my simulation, gate is starting up, but nothing is visible in the visualisation-window.
>> 
>> The last line, displayed in the terminal is: 
>> 
>> /gate/decay/setPrimaryDecayTimeGeneration 0
>> 
>> btw: does somebody know how to implement electron capture?
>> 
>> Sincerely yours,
>> 
>> Johannes Bender 
>> TU Dresden, Germany
>> 
> 

/----------------------------------------------------------------------
Prof. Dr. Uwe Pietrzyk
Diplom-Physiker
Institut für Neurowissenschaften und Medizin / INM-1
Forschungszentrum Jülich GmbH
52425 Jülich
Tel: +49-2461-61 4763
Fax: +49-2461-61 2820 + 8229
Mobil: +49-162-103-4871
eMail: U.Pietrzyk at fz-juelich.de
http://www.fz-juelich.de/inm
&
Bergische Universität Wuppertal
Fachbereich C - Physik / D.08.01
Gaußstr. 20
42097 Wuppertal
Tel: +49-(0)202-439-3523
Fax: +49-(0)202-439-2811
eMail: Uwe.Pietrzyk at physik.uni-wuppertal.de
http://www.medizinphysik.uni-wuppertal.de
----------------------------------------------------------------------/