<div dir="ltr">Dear Gate users, <div><br></div><div>In Gate simulation, Ion sources are extremely slow (not really usable for large simulations) compared to the simplest back-to-back gamma source. However, back-to-back gamma can not well model the nonpure position emitters, like I-124, Ga-68, Sc-44, and many other positron-emitting isotopes. There are two simplified fast sources for Y-90 and I-124. I want to modify one of them to create my own fast source Ga-68. The builtin source fastI124 may be a good template.  I could not find any documentation about the parameters used in the source code. Could anyone help me understand the source below?  Or point me to the documentation?</div><div><br></div><div>1. What does the state number (-2,-1,0,1,2) mean? </div><div>2. what does the <font color="#ff0000">cumulative probability </font>mean? It does not seem to match the probability in the decay scheme in the attached figure.</div><div>3. What are amplitude and normalisationFactor? </div><div><br></div><div>Thanks a lot for any help.</div><div>Zhengzhi</div><div><br></div><div><br></div><div>## Source code snippet</div><div><br></div><div>## GateSimplifiedDecayTransition.hh</div><div>class GateSimplifiedDecayTransition {<br>  friend class GateSimplifiedDecay;<br> public:<br>  GateSimplifiedDecayTransition(int cs, int ns, double pr,  std::function<psd(GateSimplifiedDecayTransition*)> act, double en=0, double ampl=0, double norm=0, int Z=0):<br><span style="background-color:rgb(0,255,0)">    currentState(cs),<br>    nextState(ns),<br>    probability(pr),<br>    action(act),<br>    energy(en),<br>    amplitude(ampl),<br>    normalisationFactor(norm),<br>    atomicNumber(Z){;}</span><br></div><div><br></div><div><br></div><div>## GateFastI124.cc</div><div>void GateFastI124::InitializeFastI124()<br>{    <br>// Forces fixed energy.  Real energy will be set later<br>     m_source->SetNumberOfParticles( 1 );<br>       m_source->GetEneDist()->SetEnergyDisType( "Mono" );<br>   m_source->SetParticleTime( m_source->GetTime() );<br>       <br>      // This defines the 13 transitions forming the simplified scheme<br>  // Entries are:<br>  //             for gammas :      current state / <font color="#0000ff">next state</font> / <font color="#ff0000">cumulative probability</font> / particle emitted (gamma)    / energy<br>  //             for e+     :      current state / <font color="#0000ff">next state </font>/ <font color="#ff0000">cumulative probability</font> / particle emitted (e+)       / max energy ...<br>  //                               ... / <font color="#00ff00">amplitude</font> of majoring function / normalisation factor for energy distribution (Fermi function) / atomic number<br><br>       m_simpleDecay = new GateSimplifiedDecay();<br>    m_particleVector = new std::vector<psd>;<br>        <br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, <font color="#0000ff">1</font>, <font color="#ff0000">0.0175</font>,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     1.376   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, <font color="#0000ff">1</font>, <font color="#ff0000">0.0488</font>,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     1.509   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, <font color="#0000ff">1</font>, <font color="#ff0000">0.1189</font>,  mem_fn( &GateSimplifiedDecayTransition::issueNone )              )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, <font color="#0000ff">2</font>,<font color="#ff0000"> 0.2359</font>,  mem_fn( &GateSimplifiedDecayTransition::issuePositron),  1.534,  1.4,  0.4407471595713562, 53) );<br></div><div>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, 2, 0.3447,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     1.691   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, 2, 0.3515,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     2.283   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, 2, 0.3598,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     0.645   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0, 2, 0.6422,  mem_fn( &GateSimplifiedDecayTransition::issueNone )              )  );<br>  <br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0,-1, 0.7502,  mem_fn( &GateSimplifiedDecayTransition::issuePositron),  2.137,  1.0,  0.10072654633851122,53) );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(0,-1, 1.,      mem_fn( &GateSimplifiedDecayTransition::issueNone )              )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(1, 2, 0.8690,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     0.722   )  );<br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(1,-1, 1.0000,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     1.325   )  );<br>  <br>  m_simpleDecay->addTransition( new GateSimplifiedDecayTransition(2,-1, 1.0000,  mem_fn( &GateSimplifiedDecayTransition::issueGamma),     0.602   )  );<br></div><div><br></div><div><br></div></div>