<html> <head> <style></style></head> <body class='hmmessage'><div dir='ltr'><br><div>Hi all, I want to simulate positron annihilation in water, and compare the range for diferent elements.</div><div>My problem is that I always get the same range regardless of the element I use.</div><div>Please, can you help me?</div><div>My code:</div><div><br><br> <style></style> <div dir="ltr">#=====================================================<br># VISUALISATION and VERBOSE<br>#=====================================================<br><br><br>/control/execute ./carpeta/Visualisation.mac<br>/control/execute ./carpeta/Verbose.mac<br><br>#=====================================================<br># GEOMETRY<br>#=====================================================<br><br>/gate/geometry/setMaterialDatabase ./GateMaterials.db<br><br><br>#=====================================================<br># W O R L D<br>#=====================================================<br><br>/gate/world/geometry/setXLength 30 cm<br>/gate/world/geometry/setYLength 30 cm<br>/gate/world/geometry/setZLength 30 cm<br>/gate/world/setMaterial Air<br><br>#=======================================================<br># G E O M E T R Y<br>#=======================================================<br><br>/control/execute ./carpeta/geometria.mac<br><br>#=====================================================<br># P H Y S I C S<br>#=====================================================<br><br>/gate/physics/addProcess RadioactiveDecay<br>/gate/physics/addProcess PositronAnnihilation e+<br>#/gate/physics/addProcess G4PositronAnnihilation e+<br>#/gate/physics/processes/G4PositronAnnihilation/setModel StandardModel<br><br>/gate/physics/processList Enabled<br>/gate/physics/processList Initialized<br><br><br>#=====================================================<br># INITIALISATION<br>#=====================================================<br><br>/gate/run/initialize<br><br>#=====================================================<br># D A T A O U T P U T <br>#=====================================================<br><br><br>/gate/output/root/enable<br>/gate/output/root/setFileName PuntualO<br>/gate/output/root/setRootHitFlag 1<br>/gate/output/root/setRootSinglesFlag 1<br>/gate/output/root/setRootOpticalFlag 1<br><br><br><br>#/control/execute ./archivo/Projection.mac<br><br>#=====================================================<br># MEASUREMENT SETTINGS <br>#=====================================================<br><br># A C Q U I S I T I O N<br><br>/control/execute ./carpeta/StartDAQ.mac<br><br><br>#=======================================================<br># Geometry<br>#=======================================================<br><br><br>/gate/world/daughters/name scanner<br>/gate/world/daughters/insert cylinder<br><br>#Fijamos las dimensiones del sistema.<br>#Seleccionamos el material.<br><br>/gate/scanner/geometry/setRmax 12. mm<br>/gate/scanner/geometry/setHeight 6. mm<br>/gate/scanner/setMaterial Water<br><br>/gate/scanner/attachCrystalSD<br><br>#=======================================================<br># Source<br>#=======================================================<br><br>/control/execute ./carpeta/fuenteO.mac<br><br>/gate/source/addSource Oxigen<br>/gate/source/Oxigen/setActivity 11000 becquerel<br><br>/gate/source/Oxigen/gps/particle e+<br>/gate/source/Oxigen/setForcedUnstableFlag true<br>/gate/source/Oxigen/setForcedHalfLife 122.24 s <br>/gate/source/Oxigen/gps/energytype Oxygen15<br><br><br>/gate/source/Oxigen/gps/type Plane<br>/gate/source/Oxigen/gps/shape Circle<br>/gate/source/Oxigen/gps/radius 0.5 mm<br>/gate/source/Oxigen/gps/angtype iso<br>/gate/source/Oxigen/gps/centre 0. 0. 0. cm<br><br><br><br>/gate/source/list<br> </div></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">Thank you!</div> </div></body> </html>