Get number of entries within certain range from TOY generated data sets

[Debashis_sahoo]

Hello,

I generate some TOY data sets from a PDF within a certain range. Now I want to get the number of data points with different sub-ranges from the generated data sets.

The generated toy data distribution is shown here.

.

Now I want to get how many events are there in between (1,1.3) and (1.3,1.8) etc. can anyone let me know how to calculate the number of entries within the subranges.

Thanks

[Wile_E_Coyote]

TH1::Integral

[Debashis_sahoo]

Hi,
Thank You for the prompt answer. But in that way, I need to convert from Roodataset to TH1. Can you let me know how to convert?

Thanks

[couet]

@Debashis_sahoo

I need to convert from Roodataset to TH1.

RooFit: saving RooDataHist as ROOT TH1 histogram?

Hello - I would like to save my generated data from RooDataSet and RooDataHist in an output text or ROOT file. For RooDataSet I see a simple way: RooDataSet* tmodel_data=tmodel.generate(t,nevtot) ; tmodel_data->write(Form(“datagen/tmodeldata_%d.dat”,niter)) (a question: could above two lines be accomplished in a single call?) But I could not figure out how to (in a few lines) save RooDataHist* tmodel_hist=tmodel.generateBinned(t,nevtot) ; as a ROOT TH1 histogram …:sunglasses: Is there a f…

[StephanH]

RooDataSet::sumEntries(const char* cutSpec, const char* cutRange = 0)

You can either set ranges on the observable, and name them one by one, or you just list cuts such as myDataSet.sumEntries("0 < x && x < 5").
Replace x with the name of your observable.

[Debashis_sahoo]

Hi Stephan,

Thank You! This works directly without converting to TH1. But I try to convert to TH1. After converting,the Integral() is not working properly when I specify limits. But Integral() without mentioning the range gives the total number of entries.

cout<<“Sum of entries…”<< data->sumEntries(“1.3 < x && x < 1.5”)<<endl;

TH1 h1 = (TH1)data->createHistogram(“h1”,x);

Double_t calc;
calc = h1->Integral(1.3,1.5);
cout<<"Integral …= "<<calc<<endl;

Output:

: warning: implicit conversion from ‘double’ to ‘Int_t’ (aka ‘int’) changes value from 1.5 to 1 [-Wliteral-conversion]
calc = h1->Integral(1.3,1.5);
~~~~~~~~ ^~~

Output:
Sum of entries…4
Integral …= 1

[StephanH]

I don’t recommend converting to TH1, as you will lose information when things are sorted into bins. Nevertheless,

@Debashis_sahoo

: warning: implicit conversion from ‘double’ to ‘Int_t’ (aka ‘int’) changes value from 1.5 to 1 [-Wliteral-conversion]
calc = h1->Integral(1.3,1.5);

The integral function takes bin numbers, not x values.

[Wile_E_Coyote]

TH1::FindFixBin