Histograms of dose profile and surface dose distribution

ROOT
1

Hello everyone!

I have information about hits in the sensitive detector: position and deposited dose corresponding to each hit which are saved in a 3D weighted histogram (dose is the weight). I have two questions regarding histograms on ROOT.

  1. I want to plot a 2D dose distribution at several depths along Z-axis (see example below):
    Surface
    Surface1274×628 173 KB

What I do is I select, for example, the first bin along the z-axis and I do a “yx” projection of my 3D histogram:

h1->GetZaxis()->SetRange(0,1)
TH2D* h2 = (TH2D*) h1->Project3D(“yx”)
h2->Draw(“colz”)

However, I am bothered by the fact that I do it in quite a cumbersome way. Is there a more nice and straightforward way to extract such 2D histogram (I am afraid I am doing it wrong selecting one bin along the z-axis etc…)

The other question is about dose profiles (see picture below). I need to get not a projection on the X axis (this is feasible) but rather a dose profile along two perpendicular lines. My idea is again select z- and y- bins in a small region around the x-axis and make a 1D projection, however, the results seem wrong to me… Is there a way to get such profiles through TH1 class commands?

Profile
Profile612×970 337 KB

Thank you very much in advance!

Kind regards,
Vasily

Please read tips for efficient and successful posting and posting code

ROOT Version: Not Provided
Platform: Not Provided
Compiler: Not Provided

2

TH2::ProjectionX
TH2::ProjectionY
TH2::ProfileX
TH2::ProfileY

3

Hello Wile,

Thank you for the prompt reply, I appreciate it much. I also tried these functions, however, my issue is that my detector is quite elongated, meaning if I use just TH2 histograms I will then sum up all the bins along the z-axis. My goal is to analyze the data only within a first few bins and plot it. That is why I thought that it would be appropriate for me to use TH3 instead of TH2 in the first place, what do you think?

Also, regarding the ProfileX/ProfileY functions: as I understood, TH2::ProfileX show the average value of Y for each bin in X, however, I need to get the profile of dose (weight in my case) along central perpendicular axis as shown before. So, I assume I would also have to work with 3D histograms in the first place…

Thanks in advance!
Kind regards,
Vasily

4

Search for methods with names “Project*” in the TH3 Class Reference.

5

Hi Wile! Thanks a lot, I’ve found my way through with the 3D projection method.

I have one more question regarding the SetProjectionX method on the ROOT canvas. I have seen a tutorial where a black window of a selected bin width pops up and draws projections on a separate screen (see image below). In my case, nothing appears and I can just scroll up down with the arrow. Do you happen to know how it can be fixed?

black
black1094×1014 71.8 KB

Cheers,
Vasily