[Gate-users] SpekCalc vs GATE x-ray energy spectrum

Emmanuel Marfo emmanuel.marfo at postgrad.otago.ac.nz
Thu Oct 12 08:29:18 CEST 2017 

Please, find the two files. The 'oldpekvsGATE' file is what I reported on earlier. The new one is a result of the implementation of your advice. I only change the number of bins from 120 to 121. It is better than before. The new file is rough because I run fewer histories. Thanks

________________________________
From: Emmanuel Marfo
Sent: Wednesday, 11 October 2017 10:04:08 p.m.
To: rachid ayad
Cc: Simon Rit; gate-users at lists.opengatecollaboration.org; Mohammad K. Ibrahim Al Matari
Subject: Re: [Gate-users] SpekCalc vs GATE x-ray energy spectrum


Thanks very much for your quick responses. I will work around the binning of the energy and see what happens. If I get them correctly, I will report back.

Best regards,

Emmanuel Marfo

________________________________
From: rachid ayad <rayad at ut.edu.sa>
Sent: Wednesday, 11 October 2017 9:44:55 p.m.
To: Emmanuel Marfo
Cc: Simon Rit; gate-users at lists.opengatecollaboration.org; Mohammad K. Ibrahim Al Matari
Subject: Re: [Gate-users] SpekCalc vs GATE x-ray energy spectrum



  Dear Emmanual, I think this is due to how we bin a histogram. Often we need to center the filled data in the histogram within the center of the bin for some specific values. For example we would like in the below example to center integer values: 0,,1, to 10 within the bin center so limits and number of bins should be given like:

      step = 1
      xmin     =  0. - step/2.
      xmax    =  10. + step/2.
      nbins   = (xmax-xmin)/step

So here you enter histograms limits xmin=-.5 and xmax=10.5 and nbins = 11 and your histograms will be centered on values 0,1,2 to 10 not shifted with one half bin if you enter xmin=0, xmax=10, and nbins = 10. Here  the first bin is centered at 0.5 and the next will be centered at 1.5 and the last at 9.5 . So I think what you did and what they did is one is shifted respect the other on how you defined your inputs to create the histogram.

Thank you, Rachid Ayad


On Wed, Oct 11, 2017 at 11:29 AM, Emmanuel Marfo <emmanuel.marfo at postgrad.otago.ac.nz<mailto:emmanuel.marfo at postgrad.otago.ac.nz>> wrote:

Please, in my simulation. I modelled the anode(rectangular geometry with thickness 1mm and cathode(source type: plane and shape: circle, source energy: monoenergetic electron(120keV)), filters(aluminium and copper), all according to the specification of our x-ray manual. I was able to generate an x-ray energy spectrum from this setup. I used both Penelope and emstandard_opt3 which apparently gives similar results.


Secondly, I used the same parameters in SpekCalc program. A program for calculating x-ray energy spectrum. I, therefore, compared the two spectra as a way of validating my spectrum from GATE. This gives the results of the image attached. If you need more info please let me know. Thanks

________________________________
From: simon.rit at gmail.com<mailto:simon.rit at gmail.com> <simon.rit at gmail.com<mailto:simon.rit at gmail.com>> on behalf of Simon Rit <simon.rit at creatis.insa-lyon.fr<mailto:simon.rit at creatis.insa-lyon.fr>>
Sent: Wednesday, 11 October 2017 7:38:45 p.m.
To: Emmanuel Marfo
Cc: gate-users at lists.opengatecollaboration.org<mailto:gate-users at lists.opengatecollaboration.org>
Subject: Re: [Gate-users] SpekCalc vs GATE x-ray energy spectrum

Hi,
I am not sure I understand what we are looking at. What was the x-ray source in your Gate simulation? I think it would help to have the full simulation.
Best regards,
Simon

On Wed, Oct 11, 2017 at 3:47 AM, Emmanuel Marfo <emmanuel.marfo at postgrad.otago.ac.nz<mailto:emmanuel.marfo at postgrad.otago.ac.nz>> wrote:

Hi GATE users,

Find attached file 'SpekCalc vs GATE x-ray energy spectrum' which compares the two spectrums. My problem is the  GATE characteristics spectrum is shifted to the right a bit when compared to the characteristics energy spectrum from SpekCalc. Please, I want to know how I can make it align with that of spekcalc. I used the energyspectrum actor definition below

### E N E R G Y S P E C T R U M ######
######################################
/gate/actor/addActor             EnergySpectrumActor espectrum
/gate/actor/espectrum/attachTo   phantom
/gate/actor/espectrum/save     Desktop/PHDthesisGATE/MARFO2017/ModellingX-raytube/XraytubeCU/output/cylinderespectrum.root
/gate/actor/espectrum/addFilter particleFilter
/gate/actor/espectrum/particleFilter/addParticle gamma
/gate/actor/espectrum/energySpectrum/setEmin           0 eV
/gate/actor/espectrum/energySpectrum/setEmax           0.12 MeV
/gate/actor/espectrum/energySpectrum/setNumberOfBins   120
/gate/actor/espectrum/energyLossHisto/setEmin          0 eV
/gate/actor/espectrum/energyLossHisto/setEmax          0.12 MeV
/gate/actor/espectrum/energyLossHisto/setNumberOfBins  120

Thanks

Best regards,

Emmanuel Marfo

University of Otago

Student

_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users


_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users



--
Dr. Rachid Ayad
Department of Physics, Faculty of Science
University Of Tabuk
P.O. Box 2072
Tabuk 71451
Kingdom Of Saudi Arabia
Tel: +966560296632



---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

هذه الرسالة و مرفقاتها (إن وجدت) تمثل وثيقة سرية قد تحتوي على معلومات تتمتع بحماية وحصانة قانونية . إذا لم تكن الشخص المعني بهذه الرسالة يجب عليك تنبيه المُرسل بخطأ وصولها إليك، و حذف الرسالة و مرفقاتها (إن وجدت) من الحاسب الآلي الخاص بك. ولا يجوز لك نسخ هذه الرسالة أو مرفقاتها (إن وجدت) أو أي جزئ منها، أو البوح بمحتوياتها لأي شخص أو استعمالها لأي غرض . علماً بأن الإفادات و الآراء التي تحويها هذه الرسالة تعبر فقط عن رأي المُرسل و ليس بالضرورة رأي جامعة تبوك  . ولا تتحمل جامعة تبوك أي مسئولية عن الأضرار الناتجة عن أي فيروسات قد يحملها هذ البريد.

This message and its attachment, if any, are confidential and may contain legally privileged information. If you are not the intended recipient, please contact the sender immediately and delete this message and its attachment , if any, from your system. You should not copy this message or disclose its contents to any other person or use it for any purpose. Statements and opinions expressed in this e-mail are those of the sender, and do not necessarily reflect those of University of Tabuk. University of Tabuk Accepts no liability for any damage caused by any virus transmitted by this email.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20171012/505e7738/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NewspekvsGATE.PNG
Type: image/png
Size: 20225 bytes
Desc: NewspekvsGATE.PNG
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20171012/505e7738/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: oldSpekvsGATE.PNG
Type: image/png
Size: 17866 bytes
Desc: oldSpekvsGATE.PNG
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20171012/505e7738/attachment-0003.png>

More information about the Gate-users mailing list