23 int main(
int argc,
char **argv)
33 JParser<> zap(
"Program to print PDF of Cherenkov light from elongated EM-shower.");
40 catch(
const exception &error) {
41 FATAL(error.what() << endl);
45 typedef JPolint0Function1D_t JFunction1D_t;
46 typedef JMAPLIST<JPolint0FunctionalMap,
47 JPolint0FunctionalMap,
48 JPolint0FunctionalMap,
49 JPolint0FunctionalGridMap,
50 JPolint0FunctionalGridMap>::maplist JMapList_t;
51 typedef JPDFTable<JFunction1D_t, JMapList_t> JPDF_t;
57 NOTICE(
"loading input from file " << inputFile <<
"... " << flush);
59 pdf.load(inputFile.c_str());
63 catch(
const JException& error) {
64 FATAL(error.what() << endl);
67 int number_of_counts = 0;
68 int number_of_errors = 0;
70 for (JPDF_t::super_const_iterator
i = pdf.super_begin();
i != pdf.super_end(); ++
i) {
72 const JFunction1D_t&
f1 =
i.getValue();
74 for (JFunction1D_t::const_iterator p = f1.begin(); p != f1.end(); ++p) {
76 number_of_counts += 1;
78 if (p->getY() != p->getY()) {
81 cout <<
FIXED(15,5) <<
i->first <<
' ';
82 cout <<
FIXED(15,5) <<
i->second->first <<
' ';
83 cout <<
FIXED(15,5) <<
i->second->second->first <<
' ';
84 cout <<
FIXED(15,5) <<
i->second->second->second->first <<
' ';
85 cout <<
FIXED(15,5) <<
i->second->second->second->second->first <<
' ';
86 cout <<
FIXED(15,5) << p->getX() <<
' ';
87 cout <<
FIXED(15,5) << p->getY() << endl;
90 number_of_errors += 1;
95 cout <<
"Number of errors / counts: " << number_of_errors <<
" / " << number_of_counts << endl;
97 if (number_of_errors != 0) {
98 FATAL(
"Number of errors " << number_of_errors << endl);
Utility class to parse command line options.
int main(int argc, char *argv[])
Auxiliary data structure for floating point format specification.
Various implementations of functional maps.
Numbering scheme for PDF types.
const JPolynome f1(1.0, 2.0, 3.0)
Function.
I/O formatting auxiliaries.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
General purpose messaging.
Utility class to parse command line options.