29 int main(
int argc,
char **argv)
41 JParser<> zap(
"Example program to histogram trigger efficiency.");
49 catch(
const exception &error) {
50 FATAL(error.what() << endl);
58 TH1D h0(
"h0", NULL, 100, 0.5, 100.5);
59 TH1D h1(
"h1", NULL, 100, 0.5, 100.5);
60 TH1D h2(
"h2", NULL, 100, 0.5, 100.5);
71 const Evt*
event = ps;
73 const int n =
event->mc_hits.size() - count_if(event->mc_hits.begin(), event->mc_hits.cend(), &
is_noise);
75 h1.Fill((Double_t)
n, 1.0);
82 STATUS(
"event: " << setw(10) << in.getCounter() <<
'\r');
DEBUG(endl);
84 const Evt*
event = in.next();
86 const int n =
event->mc_hits.size() - count_if(event->mc_hits.begin(), event->mc_hits.cend(), &
is_noise);
88 h0.Fill((Double_t)
n, 1.0);
93 for (Int_t i = 1; i <= h1.GetNbinsX(); ++i) {
95 const Double_t y1 = h1.GetBinContent(i);
96 const Double_t y0 = h0.GetBinContent(i);
100 const Double_t y3 = y1 / y0;
101 const Double_t w3 = sqrt(y1 * (y0 - y1) / (y0*y0*y0));
103 h2.SetBinContent(i, y3);
104 h2.SetBinError (i, w3);
107 ERROR(
"Bin " << h0.GetName() <<
"[" << i <<
"] empty." << endl);
General purpose messaging.
#define DEBUG(A)
Message macros.
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
ROOT TTree parameter settings of various packages.
int main(int argc, char **argv)
Synchronously read DAQ events and Monte Carlo events (and optionally other events).
Utility class to parse command line options.
virtual bool hasNext() override
Check availability of next element.
counter_type getCounter() const
Get counter.
bool is_noise(const Hit &hit)
Verify hit origin.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
KM3NeT DAQ data structures and auxiliaries.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
General purpose class for multiple pointers.
Auxiliary class to synchronously read DAQ events and Monte Carlo events (and optionally other events)...
virtual bool hasNext() override
Check availability of next element.
virtual const multi_pointer_type & next() override
Get next element.