<html> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> </head> <body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> Thanks, Devin. That’s much more clear. <div class=""><br class=""> </div> <div class="">I’ve never used anything but the default policy. However at first sight, it looks like you’re doing everything right. Maybe someone with more experience can chime in, but it’s also possible this is a bug. If you’re comfortable with C++, I’d suggest you take a quick look at the area of the code where the error message is produced and see if anything jumps at you.</div> <div class=""><br class=""> </div> <div class="">Marc</div> <div class=""> <div apple-content-edited="true" class=""><br class=""> <br class=""> </div> <div> <blockquote type="cite" class=""> <div class="">On Jun 3, 2015, at 10:04, Devin VanElburg <<a href="mailto:vanelbud@tbh.net" class="">vanelbud@tbh.net</a>> wrote:</div> <br class="Apple-interchange-newline"> <div class=""> <meta name="Author" content="Novell GroupWise WebAccess" class=""> <div style="font-family: Tahoma, sans-serif; font-size: 13px; " class="">Hello Marc, <br class=""> <br class=""> I'll try and explain a little better what I'm trying to do. I have defined my system as a PETscanner and attached the daughter volumes as follows:<br class=""> <br class=""> /gate/systems/PETscanner/level1/attach headblock<br class=""> /gate/systems/PETscanner/level2/attach crystalboard<br class=""> /gate/systems/PETscanner/level3/attach crystal<br class=""> /gate/systems/PETscanner/layer0/attach lyso<br class=""> /gate/lyso/attachCrystalSD<br class=""> <br class=""> where 'crystal' is repeated as a cubicArray (24x24) and 'headblock' is repeated twice with a ring repeater. I now have two opposing lyso crystals with the attached CrystalSD. When using the default readout policy:<br class=""> <br class=""> /gate/digitizer/Singles/readout/setPolicy TakeEnergyWinner<br class=""> <br class=""> My simulation runs fine, and the root output I receive is expected (511keV peak as I am using a Na-22 source for positron annihilation) for both Singles and Coincidences. However, when I try:<br class=""> <br class=""> /gate/digitizer/Singles/readout/setPolicy TakeEnergyCentroid<br class=""> <br class=""> I get a break with the following message:<br class=""> <br class=""> [systems/PETscanner::FindComponent]:<br class=""> Could not find the request system component 'crystal'!<br class=""> <br class=""> -------- EEEE ------- G4Exception-START -------- EEEE -------<br class=""> *** G4Exception : ProcessPulseList<br class=""> issued by : GateReadout::ProcessPulseList<br class=""> Failed to get the array component corresponding to the crystal. Abort.<br class=""> <br class=""> *** Fatal Exception *** core dump ***<br class=""> -------- EEEE -------- G4Exception-END --------- EEEE -------<br class=""> <br class=""> <br class=""> *** G4Exception: Aborting execution ***<br class=""> Aborted (core dumped)<br class=""> <br class=""> <br class=""> What I take from this is that I need to use a system that I can attach a volume to a so-called 'crystal' level. I've been playing with cylindricalPET and redefined my system as follows:<br class=""> <br class=""> /gate/systems/cylindricalPET/rsector/attach headblock<br class=""> /gate/systems/cylindricalPET/module/attach crystalboard<br class=""> /gate/systems/cylindricalPET/crystal/attach crystal<br class=""> /gate/systems/cylindricalPET/layer0/attach lyso<br class=""> /gate/lyso/attachCrystalSD<br class=""> <br class=""> that way, I have attached my 'crystal' volume to a 'crystal' level in the system hierarchy as I assume the error message suggested. Even still, I get the break segmentation violation. I also get the same error when switching to TakeEnergyCentroid for the cylindricalPET example that GATE supplies. The wiki says to ensure the Depth is taken just above the crystal, which I did. In fact I tried running it with the depth at all levels and they all failed. I'm at a loss, I can't seem to get the TakeEnergyCentroid policy to work despite my best efforts!<br class=""> <br class=""> By the way, ignore my comment with my previous email concerning the number of plates with a ring repeater. That was a readout depth issue and it's working now.<br class=""> <br class=""> I'd appreciate any help I can get here!<br class=""> <br class=""> Regards,<br class=""> Devin Van Elburg<br class=""> <br class=""> ------------------------<br class=""> Physics w/ c/ Biomedical Physics<br class=""> Lakehead University<br class=""> <br class=""> <br class=""> <br class=""> <br class=""> >>> Marc Chamberland <marcchamberland@cmail.carleton.ca class="">06/02/15 5:00 PM >>><br class=""> Hi Devin, <div class=""><br class=""> </div> <div class="">I’m not sure I understand what you’re trying to do and what failed to work. Can you provide more details?</div> <div class=""><br class=""> </div> <div class="">Also, did you attach a crystalSD to your plates, as outlined here: <a href="http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.1:Attaching_the_sensitive_detectors#The_crystalSD" class="">http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.1:Attaching_the_sensitive_detectors#The_crystalSD</a></div> <div class=""><br class=""> </div> <div class="">Cheers,<br class=""> <div class="">Marc<br class=""> <br class=""> <br class=""> <br class=""> __________________________<br class=""> <br class=""> Marc Chamberland, PhD<br class=""> Department of Physics<br class=""> Carleton University<br class=""> Ottawa (ON) </div> <br class=""> <div class=""> <blockquote type="cite" class=""> <div class="">On Jun 2, 2015, at 16:49, Devin VanElburg <<a href="mailto:vanelbud@tbh.net" class="">vanelbud@tbh.net</a>> wrote:</div> <br class="Apple-interchange-newline"> <div class=""> <meta name="Author" content="Novell GroupWise WebAccess" class=""> <div style="font-family: Tahoma, sans-serif; font-size: 13px; " class="">Hello gate users,<br class=""> <br class=""> I am simulating a PET system with two plates. I would like to use the takeEnergyCentroid option for the digitizer's readout. However, I need to attach a volume to a 'crystal' of a system (according to the error message I keep getting). Thus a PETscanner system cannot work since it uses level1-4 and layers for system attachment. The I tried a cylindricalPET system with a ring repeater of 2, and that also didn't work. Interestingly, cylindricalPET works for a 4+ plate ring system, but no singles or coincidences in the output for 2-3 plates (with only the ring repeat number changed). Can anyone offer any insight as to what I should do?<br class=""> <br class=""> Regards,<br class=""> Devin Van Elburg<br class=""> <br class=""> ------------------------<br class=""> Physics w/ c/ Biomedical Physics<br class=""> Lakehead University<br class=""> </div> _______________________________________________<br class=""> Gate-users mailing list<br class=""> <a href="mailto:Gate-users@lists.opengatecollaboration.org" class="">Gate-users@lists.opengatecollaboration.org</a><br class=""> <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" class="">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></div> </blockquote> </div> <br class=""> </div> </marcchamberland@cmail.carleton.ca></div> </div> </blockquote> </div> <br class=""> </div> </body> </html>