Logo ROOT  
Reference Guide
mt304_fillHistos.C File Reference

Detailed Description

View in nbviewer Open in SWAN Fill histograms in parallel with automatic binning. Illustrates use of power-of-two autobin algorithm

// The number of workers
const UInt_t nWorkers = 8U;
// Reference boundaries
const Double_t xmiref = -1.;
const Double_t xmaref = 7.;
Int_t mt304_fillHistos(UInt_t nNumbers = 1001)
{
// The first, fundamental operation to be performed in order to make ROOT
// thread-aware.
// Histograms to be filled in parallel
ROOT::TThreadedObject<TH1D> h1d("h1d", "1D test histogram", 64, 0., -1.);
ROOT::TThreadedObject<TH1D> h1dr("h1dr", "1D test histogram w/ ref boundaries", 64, xmiref, xmaref);
// We define our work item
auto workItem = [&](UInt_t workerID) {
// One generator, file and ntuple per worker
TRandom3 workerRndm(workerID); // Change the seed
auto wh1d = h1d.Get();
wh1d->SetBit(TH1::kAutoBinPTwo);
auto wh1dr = h1dr.Get();
for (UInt_t i = 0; i < nNumbers; ++i) {
x = workerRndm.Gaus(3.);
wh1d->Fill(x);
wh1dr->Fill(x);
}
};
// Create the collection which will hold the threads, our "pool"
std::vector<std::thread> workers;
// Fill the "pool" with workers
for (auto workerID : ROOT::TSeqI(nWorkers)) {
workers.emplace_back(workItem, workerID);
}
// Now join them
for (auto &&worker : workers)
worker.join();
// Merge
auto fh1d = h1d.Merge();
auto fh1dr = h1dr.Merge();
// Make the canvas
auto c = new TCanvas("c", "c", 800, 800);
c->Divide(1, 2);
gStyle->SetOptStat(111110);
c->cd(1);
fh1d->DrawCopy();
c->cd(2);
fh1dr->DrawCopy();
c->Update();
return 0;
}
Date
November 2017
Author
Gerardo Ganis

Definition in file mt304_fillHistos.C.

c
#define c(i)
Definition: RSha256.hxx:119
ROOT::EnableThreadSafety
void EnableThreadSafety()
Enables the global mutex to make ROOT thread safe/aware.
Definition: TROOT.cxx:494
Int_t
int Int_t
Definition: RtypesCore.h:45
x
Double_t x[n]
Definition: legend1.C:17
gStyle
R__EXTERN TStyle * gStyle
Definition: TStyle.h:412
TRandom3
Definition: TRandom3.h:27
UInt_t
unsigned int UInt_t
Definition: RtypesCore.h:46
ROOT::TThreadedObject
A wrapper to make object instances thread private, lazily.
Definition: TThreadedObject.hxx:151
TStyle::SetOptStat
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
Definition: TStyle.cxx:1592
Double_t
double Double_t
Definition: RtypesCore.h:59
TCanvas
Definition: TCanvas.h:23
ROOT::TSeq
A pseudo container class which is a generator of indices.
Definition: TSeq.hxx:66
TH1::kAutoBinPTwo
@ kAutoBinPTwo
Use Power(2)-based algorithm for autobinning.
Definition: TH1.h:170