Comments (6)
I suspect the problem is in the analysis. How are you extracting the energy spectra? Are you using an actor that filters on primary particles generated? That you have an identical number of entries in both spectra is suspicious. At the very least, there should be far fewer positrons (~32 per million decays) than brem photons (~6 % of the number of decays).
from gate.
I'm using the code you provided in plot.C:
...
// open the output file
TFile* f2 =new TFile (" fastY90Pos.root");
// and print the content
f2->ls();
// retrieve the energy spectrum histogram
TH1F * h2 = (TH1F*)f2->Get("energySpectrum");
h2->SetLineColor(kRed);
// set title and the axis legend
h2->SetTitle("fastY90 positron kinetic energy");
h2->GetXaxis()->SetTitle("Energy (MeV)");
h2->GetYaxis()->SetTitle("Counts");
h2->GetXaxis()->SetRangeUser(0. ,1);
// statistic type ( (1) write mean value, number of entries and RMS)
h2->SetStats(1);
// draw histogram using point (P) associated with error bar (E) or HIST for histogramm
c1->cd(2);
h2->Draw("HIST");
...
and the output was generated with the unmodified output.mac:
...
/gate/actor/addActor EnergySpectrumActor pos_actor
/gate/actor/pos_actor/attachTo MainBox
/gate/actor/pos_actor/addFilter particleFilter
/gate/actor/pos_actor/particleFilter/addParticle e+
/gate/actor/pos_actor/save output/fastY90Pos.root
/gate/actor/pos_actor/energySpectrum/setEmin 0 eV
/gate/actor/pos_actor/energySpectrum/setEmax 1 MeV
/gate/actor/pos_actor/energySpectrum/setNumberOfBins 100
I agree with your point on the suspicious numbers of particles being exactly the same. Even inspecting the fastY90Pos.root file directly yields the duplicate data:
from gate.
I really don't understand. I created the example plot using the code that's in the example directory, and rereading, everything appears to be in order.
One thing I do notice is that the total number of entries in your output (6504827) is exactly equal to the sum of the brem and positron entries (6501704 and 3123 respectively) that I got when I created the example. It seems like the particle filter of the actor is doing nothing and both of your outputs contain both the gammas and the positrons. The number of positrons is small, so the peak isn't visible in the plot.
from gate.
Hello,
any progress on that issue ? Should we still consider it as open ?
from gate.
I tried to reproduce this behavior on multiple machines with the same configuration (Ubuntu 16.04 LTS, Geant4.10.02.p02, Root 5.34.36) and I could not get the plots provided, i.e. the correct positron kinetic energy. Maybe someone else can try and confirm this issue?
from gate.
I re-evaluated this issue with Gate 8.0, Geant4 10.3.1 and root 6.10/00 and I cannot reproduce it any more.
With the versions mentioned it works as intended. I consider this issue closed.
from gate.
Related Issues (20)
- Issue with GateVImageVolume.cc (l.734)
- Ubuntu 20.04 Snap Version of Gate Permission Issues HOT 1
- Issue with GateSourcePhaseSpace.cc if the source is in ROOT format after installing V9.3
- Image layers on docker image for Gate9.3 doesn't match with script versions
- Attach actor to "world" volume causing summary to malfuntion. HOT 13
- Whether Gate can output the energy information of secondary particles?
- segmentation fault when using SinglesDigitizer and ROOT output HOT 1
- separating slash in branch name hinders easy plotting HOT 2
- for pPs + prompt gamma simulation with Gate
- material issue!
- setSkipEqualMaterials of ImageRegularParametrisedVolume setting to zero causing system to complaining geometry overlapping
- /gate/source/F-18/imageReader/linearTranslator/setScale 1 MBq
- Directional Keys Not Working in Gate Software
- Creating a circular source in Gate
- question about the difference of F18-source and gamma-source
- All Crystal IDs are the same with multiple crystals attached HOT 1
- Memory Leak Issue with Lu-177 Ion Source in Gate 9.0 HOT 2
- Issues with Wedge Geometry in Tungsten Slit Collimator Design in Gate 9.0
- If there is a Bug!
- questions in CT simulation
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from gate.