Hello everyone, I’m trying to follow the example of Interpolation function that you can find here: https://root.cern.ch/function-interpolation . My code is: #include <stdint.h> #include <stdio.h> #include "TFile.h" #include "TH1F.h" #include "TF1.h" #include <vector> #include <iostream> #include "…

You at least need some error checking: ... TFile *f=new TFile("run_4/result.root" ,"READ"); //The file might not have opened correctly. if (!f->IsOpen()) { cerr<<"ERROR: The file didn't open!\n"; return -1; } const Int_t ni = xmax-xmin/6.; Double_t xi[ni], yi[ni]; const In…

Thanks very much for your help, It breaks anyway, but I guess it’s more complete now, always good to know how to better code. I will look into the error hint,thanks!

I found the line where I have a problem, and I corrected some more errors, just in case someone is interested in helping me, the line is yi[i] = f1(xi[i]); And the corrected code is: #include <stdio.h> #include "TFile.h" #include "TH1F.h" #include "TF1.h" #include <iostream> #include "Math/Polyno…

[image] Deb: for(int i_bin=1; i_bin<200; i_bin++){ if(h->GetBinContent(i_bin)!=0 && n==0){ binmin=i_bin; n++; } if(h->GetBinContent(i_bin)==0 && h->GetBinContent(i_bin+5)==0 && n>0){ binmax=i_bin; i_bin=200; } } What about this piece of code? You don’…

[image] ksmith: You don’t check that the histogram actually has that number of bins, and the i_bin+5 is bound to be problematic. Maybe this for loop works better for you? for(int i_bin=1; i_bin<200 && i_bin<=h->GetNbinsX()-5; i_bin++){ After the above correction and removing the following bl…

Thank you very much, I knew the number of bins in the histogram but I didn’t think about the -5, thanks. Except it doesn’t work yet, at this point I think it’s something in my rootfile, something that has to do with my histos I guess. Anyway thanks again, have a nice day. Debora

Home | News | Documentation | Download

ROOT Forum

ROOT::Interpolation: Busy flag cleared

ROOT

root

ksmith August 29, 2017, 1:22pm 2

Just FYI:

Getting the same histograms after cuts ROOT

Remark on the way to post code snippets: When you post pieces of code on the forum enclose it with ``` . That makes your post much more readable (I systematically edit your posts that way) . So do: ``` // this is a C++ piece of code. ``` And it will display as: // this is a C++ piece of code.

show post in topic

Home
Categories
Guidelines
Terms of Service
Privacy Policy

Powered by Discourse, best viewed with JavaScript enabled