32 int main(
int argc,
char **argv)
43 JParser<> zap(
"Program to add multi-dimensional histograms of shower light");
51 catch(
const exception &error) {
52 FATAL(error.what() << endl);
59 JMAPLIST<JHistogramMap_t,
62 JHistogramGridMap_t>::maplist> JMultiHistogram_t;
64 typedef JPDFTransformer<4, abscissa_type> JFunction4DTransformer_t;
69 h1.transformer.reset(
new JFunction4DTransformer_t());
75 NOTICE(
"loading input from file " << *file_name <<
"... " << flush);
79 JFileStreamReader
in(file_name->c_str());
86 p1.transformer.reset(
new JFunction4DTransformer_t());
88 for(JMultiHistogram_t* p : { &p0, &p1 }) {
98 double integral1 = 0.0;
99 double integral2 = 0.0;
101 for(JMultiHistogram_t::super_iterator
102 i0 = h0.super_begin(), i1 =
h1.super_begin(),
103 j0 = p0.super_begin(), j1 = p1.super_begin(); i1 !=
h1.super_end(); ++i0, ++i1, ++j0, ++j1) {
105 integral1 += i0.getValue().getIntegral();
107 i0.getValue().add(j0.getValue());
108 i1.getValue().add(j1.getValue());
110 integral2 += i0.getValue().getIntegral();
115 for(JMultiHistogram_t* p : { &h0, &
h1 }) {
126 catch(
const JException& error) {
127 FATAL(error.what() << endl);
133 NOTICE(
"Storing, " << flush);
135 for (
const JMultiHistogram_t* p : { &h0, &
h1 }) {
Utility class to parse command line options.
then for HISTOGRAM in h0 h1
Various implementations of functional maps.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
General purpose messaging.
Utility class to parse command line options.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
int main(int argc, char *argv[])