Manipulation of root files returned by a SecondaryProduction actor

Newbie
1

Dear Root users,

I am using GATE (an open source software dedicated to digital simulations in medical imaging and radiotherapy) to analyze the secondary particles produced after an X-rays irradiation simulation.

I use a SecondaryProductionActor and my GATE simulations give rise to .root files in which I can found the number of secondaries in function of the particle type. I know how to open the root file through the root prompt

2
2737Ă—480 61.8 KB

which gives rise to

1
11346Ă—982 31.6 KB

But I would like to display the number of particles for a specific category of secondaries, let’s say for example the e+. How can I do that ?

The

fragments->Draw(“e+”)

command seems to display all types of secondaries and the result is not really readable.
In addition, the goal is subsequently to record the number of particles produced for each type in an excel file.

If anyone could give me some pointers on how to get started, I was able to do this task for root files returned by a phaseSpace actor but failed to do so in the case of the SecondaryProductionActor…

Many thank for the help,
Regards,
Sarah

Platform : Ubuntu 20.04.6 LTS
Root version : ROOT Version: 6.30/04

2

The histogram drawing options are here:
https://root.cern/doc/master/classTHistPainter.html
I guess you tried Draw("e+") thinking this will select the bin of “e+” (if there is such a bin in the histogram), but that’s not what it does. “E” draws errors, and the “+” in “e+” is ignored (see page above).
You could just zoom in on the x-axis and display only “e+” if you find it, or find the bin number to which it corresponds and then limit the axis range with SetRangeUser, e.g., supposing it is bin 4:

h->SetRangeUser(3,4);

note that the range is from bin i-1 to bin i, if “e+” is bin i.
To find out the bin number, you can print all bin labels and numbers like this:

for (int i=1; i<=h->GetXaxis()->GetNbins(); ++i) {
  cout << i << " " << h->GetXaxis()->GetBinLabel(i) << endl;

With this, you can figure out the rest (you can use h->GetBinContent(yourbinnumber))

3

Dear dastudillo,

Thank you so much for your help!
Yes, I did try to draw the bin labelled “e+” in the histogram this way.

So know I know how to recover some bin label and its corresponding value through the root prompt :

3

In order to make a root macro that displays all bin label and corresponding values, how can I recover the variable “fragments” that corresponds to the TH1F object? This does not work :

void dataMacro()
{
    string path_HYB = "/home/sarah/Documents/ROOT/Layer1-Secondaries.root";

    TFile *f = new TFile("Layer1-Secondaries.root");
    f->ls();

    for (int i=1; i<=f.fragments->GetXaxis()->GetNbins(); ++i)
        {
            cout << i << " " << f.fragments->GetXaxis()->GetBinLabel(i) << " " << f.fragments->GetBinContent(i) << endl;
        }
}

I’m currently reading the link ROOT: THistPainter Class Reference you sent me.

Many thanks,
Sarah

4

Check out the histogram tutorials:
https://root.cern/doc/master/group__tutorial__hist.html
for instance, this one to read from a TFile:
https://root.cern/doc/master/h1ReadAndDraw_8C.html

6

It’s okey, it works fine know, thank you very much for your patience and help !

void dataMacro2()
{
   TFile *f = TFile::Open("/home/sarah/Documents/ROOT/Layer1-Secondaries.root");

   TH1 *hpx = nullptr;

   f->GetObject("fragments", hpx);

   for (int i=1; i<=hpx->GetXaxis()->GetNbins(); ++i)
        {
            if (hpx->GetBinContent(i) != 0)
            {
                cout << hpx->GetXaxis()->GetBinLabel(i) << " : " << hpx->GetBinContent(i) << endl;
            }
        }
}
1 Like