Hi, I have a question about: I have an histogram that gives the energy distribution in a detector. So, on the X axis I have the energy range and on the Y axis the total number of times each bin is filled. Now, I want to get the total energy… so I thought that the Integral() method should help me. …

The loop you suggest seems fine assuming in each bin all the counts are comming from the same value equal to the bin center. May be @moneta has a better sugestion.

Yes, It was an assumption I implicitly made since I don’t know how to consider all the possibile energies in a given bin.

if you have access to the data before the histogram filling you can compute the exact total energy.

Mh I don’t think so, because data come from a Geant4 simulation and the histos are automatically filled.

ok. maybe @moneta has more hints to provide.

I was thinking that things could be easier than expected… since my bins have all the same width, let’s say 0.1 MeV, I could simply take for each bin the total number of entries and multiply by 0.1MeV to get the total energy deposit in that bin. Then, I sum all contributes.

if they are the same width you can do what you suggest assuming the x-axis starts at 0 MeV. And as the integral is the sum of the contents of the bins may be it is enough to do 0.1*integral ?

Ok, so the first two options give the same output, which is in both cases different from the one I get using my code…

So, following your suggestion, the associated error could be computed as hist->GetMeanError()*hist->GetSumofWeights()… right?

Home | News | Documentation | Download | Privacy Notice

ROOT Forum

Get total energy from an histogram

ROOT

Wile_E_Coyote March 24, 2023, 8:15am 9

Histogram sum ROOT

Try (sensitive to user’s changes to the current x-axis range): std::cout << h->GetMean() * h->Integral() << "\n"; // mean * integral and/or (valid only for the “entire” x-axis range): std::cout << h->GetMean() * h->GetSumOfWeights() << "\n"; Another good idea would be (sensitive to user’s changes to the current x-axis range): // change the x-axis range, if desired, and then ... Double_t stats[(TH1::kNstat)]; // sums of w, w^2, w*x, w*x^2, ... h->GetStats(stats); std::cout << stats[2] << "\n…

show post in topic

Home
Categories
Guidelines
Terms of Service
Privacy Policy

Powered by Discourse, best viewed with JavaScript enabled