---

main.cxx

Go to the documentation of this file.

// $Id: main.cxx,v 1.1.1.1 1999/12/20 22:27:05 burnett Exp $
//
// Original author: Toby Burnett, mailto:tburnett@u.washington.edu
//
// make a histogram from a tuple
#include "analysis/Tuple.h"
#include "analysis/Histogram.h"
#include <fstream>
#include <strstream>
using namespace std;


static const char* helptext=
"\n----------------------------------------------"
"\n   makehist tuple_file histo_file item_name [min] [max] [step] [factor] [offset]"
"\n"
"\n  tuple_file  Name of standard tuple file,  or \"help\" to give advice"
"\n  histo_file  Name of histogram tuple file, or \"-\" to write to standard otput"
"\n  item_name   Name of a column in the tuple"
"\n  min         left edge of histogram [0]"
"\n  max         right edge [min+1]"
"\n  step        bin size [(max-min)/25]"
"\n  factor      multiply the item [1.0]"
"\n  offset      add to the item [0]"
"\n"
;

int main(int argc, const char** argv)
{
    int         n=0; // counts arguments
    string tuple_file("-");
    string histo_file("-");

    if( argc>++n ) tuple_file = argv[n];                // required, input file or "help"
    if (string(tuple_file).find("?") != string::npos
                || tuple_file=="help") {
        cerr << helptext << std::endl;
        return 0;
    }

    if( argc>++n ) histo_file = argv[n];                // required, output file or cout
    if (string(histo_file).find("?") != string::npos
                        && histo_file!="-") {           // "-" specifies standard output
        cerr <<  "command line problem" << std::endl;
        return 0;
    }

    // other args, all optional
    string item_name("MC_Energy");
    if( argc>++n ) item_name = argv[n];

    double xmin   = (argc>++n)?  atof(argv[n]) : 0;
    double xmax   = (argc>++n)?  atof(argv[n]) : xmax = xmin+1.0;
    double xstep  = (argc>++n)?  atof(argv[n]) : (xmax-xmin)/25;
    double factor = (argc>++n)?  atof(argv[n]) : 1.0;
    double offset = (argc>++n)?  atof(argv[n]) : 0;
    if( argc>++n ) offset= atof(argv[n]);

    // now assign input stream, either the file name or std input

    istream* infile;
    if( tuple_file== "-") {
        infile = &std::cin;
    } else {
        ifstream* realfile = new ifstream(tuple_file.c_str());
        if( !realfile->is_open()) {
            cerr << "Cannot open input file \"" << tuple_file << "\"" << endl;
            exit(1);
        }

        infile = realfile;
    }

    Tuple tup(*infile);

    Tuple::const_iterator item= tup.find(item_name);
    if( item==tup.end() ) {
        cerr << "Tuple item '" << item_name.c_str() << "' not found: header follows" << endl;
        tup.writeHeader(cerr);
        return(2);
    }

    const TupleItem& tuple_item= **item;

    // create title for histogram
    strstream hist_title;
    if (factor != 1.0) hist_title << factor << '*';
    hist_title << tuple_item.name();
    if (offset > 0.) hist_title << '+' << offset;
    else if (offset < 0.) hist_title << offset;
    hist_title << '\0';

    // create the histogram
    Histogram hist( hist_title.str(), xmin, xmax, xstep) ;

    do {
        (*infile) >> tup;
        if( ! infile->good() ) break;
        if( tuple_item() != 0) {
                hist.fill( factor*tuple_item()+offset ) ;
        }

    } while( 1 ) ;

    // now assign output stream, either the file name or std output

    ostream* outfile;
    if( histo_file== "-") {
        outfile = &std::cout;
    }
    else {
        outfile = new ofstream(histo_file.c_str());
    }
    *outfile <<  "Tuple title: \""<< tup.title() << "\"\n" ;
    hist.print(*outfile);
    return 0;
}

---