<html> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> </head> <body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> <div>Hi Johannes,</div> <div> </div> <div>for the ION-Source to work properly experts informed me, that a part in</div> <div>.../source/management/src/GateTrajectory.cc needs to be changed.</div> <div>It actually is only the line of code from Jan De Neenhouwer. See below.</div> <div>Also note, that depending on the type of ION-source, some of the emitted</div> <div>particles simply do not travel very far. Concerning electron capture,</div> <div>I cannot tell, whether it is implemented and if not, how to do this.</div> <div>The general particle source (gps), however, I would guess should</div> <div>have this mode included, but please check in the file that are</div> <div>anyhow required by GEANT: RadioactiveDecay3.2 ...</div> <div> </div> <div>I hope this is helpful.</div> <div> </div> <div>Kind regards,</div> <div> </div> <div>Uwe</div> <div> </div> <div> </div> <div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; "> void GateTrajectoryNavigator::SetIonID()</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; "> {</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">   G4int n_trajectories = m_trajectoryContainer->entries();</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">   for (G4int iTrj=0; iTrj<n_trajectories; iTrj++) {</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">     G4Trajectory* trj = (G4Trajectory*)((*m_trajectoryContainer)[iTrj]);</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">     if (trj->GetTrackID() == 1) {</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; color: rgb(0, 140, 0); "> //       G4cout << "GateTrajectoryNavigator::SetIonID: PDG " << trj->GetPDGEncoding() </div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; color: rgb(0, 140, 0); "> //      << " Name <" << trj->GetParticleName() << ">" << G4endl;</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; color: rgb(0, 140, 0); "> //      if (trj->GetPDGEncoding() == 0) { // 0 == Ion -----> to be changed !! </div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; color: rgb(0, 140, 0); ">       if (trj->GetPDGEncoding() > 1000000000) { // Correction from Jan De Beenhouwer </div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; min-height: 13px; "> </div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; "> m_ionID = 1;</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">       } else {</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; "> m_ionID = 0;</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">       }</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">     }</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; ">   }</div> <div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Menlo; "> }</div> </div> <div> </div> <div> </div> <div> </div> <div> <div>On 29.10.2009, at 17:57, Johannes Bender wrote:</div> <blockquote type="cite"> <div>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 (I'm sorry for the long listing) # S O U R C E /gate/source/addSource 86YGamma /gate/source/86YGamma/setActivity 100000 Bq /gate/source/86YGamma/gps/centre 0. 0. 0. cm /gate/source/86YGamma/gps/particle gamma /gate/source/86YGamma/gps/energytype Epn /gate/source/86YGamma/gps/emin 132.34 keV /gate/source/86YGamma/gps/emax 3877 keV /gate/source/86YGamma/gps/type Volume /gate/source/86YGamma/gps/shape Sphere /gate/source/86YGamma/gps/radius 1 mm /gate/source/86YGamma/gps/angtype iso /gate/source/86YGamma/gps/mintheta    0. deg /gate/source/86YGamma/gps/maxtheta    180. deg /gate/source/86YGamma/gps/minphi      0. deg /gate/source/86YGamma/gps/maxphi      360. deg /gate/source/list /gate/source/86YGamma/gps/histpoint 132.34 0.17 /gate/source/86YGamma/gps/histpoint 144.5 0.03 /gate/source/86YGamma/gps/histpoint 182.34 0.11 /gate/source/86YGamma/gps/histpoint 187.87 1.26 /gate/source/86YGamma/gps/histpoint 190.79 1.01 /gate/source/86YGamma/gps/histpoint 209.8 0.4 /gate/source/86YGamma/gps/histpoint 235.4 0.4 /gate/source/86YGamma/gps/histpoint 237.9 0.13 /gate/source/86YGamma/gps/histpoint 252.05 0.37 /gate/source/86YGamma/gps/histpoint 256.4 0.07 /gate/source/86YGamma/gps/histpoint 264.53 0.54 /gate/source/86YGamma/gps/histpoint 307 3.46 /gate/source/86YGamma/gps/histpoint 331.08 0.83 /gate/source/86YGamma/gps/histpoint 355.07 0.1 /gate/source/86YGamma/gps/histpoint 370.28 0.82 /gate/source/86YGamma/gps/histpoint 380.4 0.45 /gate/source/86YGamma/gps/histpoint 383.02 3.63 /gate/source/86YGamma/gps/histpoint 425.97 0.31 /gate/source/86YGamma/gps/histpoint 439.5 0.2 /gate/source/86YGamma/gps/histpoint 443.14 16.9 /gate/source/86YGamma/gps/histpoint 444.18 0.64 /gate/source/86YGamma/gps/histpoint 448.1 0.07 /gate/source/86YGamma/gps/histpoint 469.24 0.3 /gate/source/86YGamma/gps/histpoint 503 0.09 /gate/source/86YGamma/gps/histpoint 512.42 /gate/source/86YGamma/gps/histpoint 515.18 4.89 /gate/source/86YGamma/gps/histpoint 580.57 4.79 /gate/source/86YGamma/gps/histpoint 608.29 2.01 /gate/source/86YGamma/gps/histpoint 619.1 0.21 /gate/source/86YGamma/gps/histpoint 627.72 32.6 /gate/source/86YGamma/gps/histpoint 634.78 0.09 /gate/source/86YGamma/gps/histpoint 644.8 2.2 /gate/source/86YGamma/gps/histpoint 645.9 9.2 /gate/source/86YGamma/gps/histpoint 648.6 /gate/source/86YGamma/gps/histpoint 689.29 0.17 /gate/source/86YGamma/gps/histpoint 702.2 0.25 /gate/source/86YGamma/gps/histpoint 703.34 15.4 /gate/source/86YGamma/gps/histpoint 709.9 2.62 /gate/source/86YGamma/gps/histpoint 719.17 0.22 /gate/source/86YGamma/gps/histpoint 740.81 1.36 /gate/source/86YGamma/gps/histpoint 767.63 2.4 /gate/source/86YGamma/gps/histpoint 768.3 0.32 /gate/source/86YGamma/gps/histpoint 777.35 22.4 /gate/source/86YGamma/gps/histpoint 783.56 0.26 /gate/source/86YGamma/gps/histpoint 826.02 3.3 /gate/source/86YGamma/gps/histpoint 833.7 1.5 /gate/source/86YGamma/gps/histpoint 835.7 4.4 /gate/source/86YGamma/gps/histpoint 882.96 0.25 /gate/source/86YGamma/gps/histpoint 887.4 0.44 /gate/source/86YGamma/gps/histpoint 955.35 1.04 /gate/source/86YGamma/gps/histpoint 971.43 0.27 /gate/source/86YGamma/gps/histpoint 1017.93 0.18 /gate/source/86YGamma/gps/histpoint 1024.04 3.79 /gate/source/86YGamma/gps/histpoint 1076.64 83 /gate/source/86YGamma/gps/histpoint 1087.6 0.04 /gate/source/86YGamma/gps/histpoint 1092.68 0.69 /gate/source/86YGamma/gps/histpoint 1102.02 0.2 /gate/source/86YGamma/gps/histpoint 1133.3 0.3 /gate/source/86YGamma/gps/histpoint 1142.3 0.1 /gate/source/86YGamma/gps/histpoint 1150.3 /gate/source/86YGamma/gps/histpoint 1153.01 30.5 /gate/source/86YGamma/gps/histpoint 1154 /gate/source/86YGamma/gps/histpoint 1163.03 1.18 /gate/source/86YGamma/gps/histpoint 1253.11 1.53 /gate/source/86YGamma/gps/histpoint 1270.16 0.65 /gate/source/86YGamma/gps/histpoint 1283.96 0.29 /gate/source/86YGamma/gps/histpoint 1294.9 0.29 /gate/source/86YGamma/gps/histpoint 1296.03 0.54 /gate/source/86YGamma/gps/histpoint 1327.5 0.09 /gate/source/86YGamma/gps/histpoint 1349.15 2.95 /gate/source/86YGamma/gps/histpoint 1404.8 0.18 /gate/source/86YGamma/gps/histpoint 1415.2 0.33 /gate/source/86YGamma/gps/histpoint 1507.86 0.35 /gate/source/86YGamma/gps/histpoint 1533.19 0.22 /gate/source/86YGamma/gps/histpoint 1535.67 0.12 /gate/source/86YGamma/gps/histpoint 1564.4 0.18 /gate/source/86YGamma/gps/histpoint 1696.25 0.64 /gate/source/86YGamma/gps/histpoint 1711.6 0.17 /gate/source/86YGamma/gps/histpoint 1724.15 0.55 /gate/source/86YGamma/gps/histpoint 1790.9 1 /gate/source/86YGamma/gps/histpoint 1801.7 1.65 /gate/source/86YGamma/gps/histpoint 1854.38 17.2 /gate/source/86YGamma/gps/histpoint 1920.72 20.8 /gate/source/86YGamma/gps/histpoint 1969.1 0.05 /gate/source/86YGamma/gps/histpoint 2017.1 0.13 /gate/source/86YGamma/gps/histpoint 2088.09 0.25 /gate/source/86YGamma/gps/histpoint 2108.9 0.05 /gate/source/86YGamma/gps/histpoint 2180.8 0.03 /gate/source/86YGamma/gps/histpoint 2291.8 0.12 /gate/source/86YGamma/gps/histpoint 2482.08 0.12 /gate/source/86YGamma/gps/histpoint 2555.3 0.03 /gate/source/86YGamma/gps/histpoint 2567.97 2.25 /gate/source/86YGamma/gps/histpoint 2610.11 1.24 /gate/source/86YGamma/gps/histpoint 2641.9 0.17 /gate/source/86YGamma/gps/histpoint 2790 0.01 /gate/source/86YGamma/gps/histpoint 2794.9 0.21 /gate/source/86YGamma/gps/histpoint 2827.7 0.06 /gate/source/86YGamma/gps/histpoint 2862 0.01 /gate/source/86YGamma/gps/histpoint 2865.9 0.38 /gate/source/86YGamma/gps/histpoint 2997.6 0.01 /gate/source/86YGamma/gps/histpoint 3069.7 0.12 /gate/source/86YGamma/gps/histpoint 3334 0.12 /gate/source/86YGamma/gps/histpoint 3642 0.04 /gate/source/86YGamma/gps/histpoint 3877 0.05 /gate/source/addSource 86YPositron /gate/source/86YPositron/setActivity 100000 Bq /gate/source/86YPositron/gps/centre 0. 0. 0. cm /gate/source/86YPositron/gps/particle e+ /gate/source/86YPositron/gps/energytype Epn /gate/source/86YPositron/gps/emin 249.11 keV /gate/source/86YPositron/gps/emax 3141.32 keV /gate/source/86YPositron/gps/type Volume /gate/source/86YPositron/gps/shape Sphere /gate/source/86YPositron/gps/radius 1 mm /gate/source/86YPositron/gps/angtype iso /gate/source/86YPositron/gps/mintheta    0. deg /gate/source/86YPositron/gps/maxtheta    180. deg /gate/source/86YPositron/gps/minphi      0. deg /gate/source/86YPositron/gps/maxphi      360. deg /gate/source/list /gate/source/86YPositron/gps/histpoint 249.11 0 /gate/source/86YPositron/gps/histpoint 275.57 0 /gate/source/86YPositron/gps/histpoint 292.03 0.04 /gate/source/86YPositron/gps/histpoint 346.5 0 /gate/source/86YPositron/gps/histpoint 386.94 0.22 /gate/source/86YPositron/gps/histpoint 443.09 0.01 /gate/source/86YPositron/gps/histpoint 452.31 0.26 /gate/source/86YPositron/gps/histpoint 531.17 0.06 /gate/source/86YPositron/gps/histpoint 573.02 0.31 /gate/source/86YPositron/gps/histpoint 855.93 0.18 /gate/source/86YPositron/gps/histpoint 900.38 1.1 /gate/source/86YPositron/gps/histpoint 926.7 0.04 /gate/source/86YPositron/gps/histpoint 1032.77 1.9 /gate/source/86YPositron/gps/histpoint 1162.22 1.33 /gate/source/86YPositron/gps/histpoint 1220.64 11.9 /gate/source/86YPositron/gps/histpoint 1339.72 0.69 /gate/source/86YPositron/gps/histpoint 1429.5 0.05 /gate/source/86YPositron/gps/histpoint 1545.16 5.6 /gate/source/86YPositron/gps/histpoint 1575.81 0.05 /gate/source/86YPositron/gps/histpoint 1736.09 1.7 /gate/source/86YPositron/gps/histpoint 1988.26 3.6 /gate/source/86YPositron/gps/histpoint 2363.83 0.9 /gate/source/86YPositron/gps/histpoint 3141.32 2 _______________________________________________ Gate-users mailing list <a href="mailto:Gate-users@lists.healthgrid.org">Gate-users@lists.healthgrid.org</a> http://lists.healthgrid.org/mailman/listinfo/gate-users </div> </blockquote> </div> <div> <div>/---------------------------------------------------------------------- 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: <a href="mailto:U.Pietrzyk@fz-juelich.de">U.Pietrzyk@fz-juelich.de</a> <a href="http://www.fz-juelich.de/inm">http://www.fz-juelich.de/inm</a> & 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: <a href="mailto:Uwe.Pietrzyk@physik.uni-wuppertal.de">Uwe.Pietrzyk@physik.uni-wuppertal.de</a> <a href="http://www.medizinphysik.uni-wuppertal.de/">http://www.medizinphysik.uni-wuppertal.de</a> ----------------------------------------------------------------------/ </div> </div> ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ </body> </html>