runPDSPThinSliceFit.c++

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <sstream>
#include <string>

#include "cetlib/filepath_maker.h"
#include "fhiclcpp/intermediate_table.h"
#include "fhiclcpp/make_ParameterSet.h"
#include "fhiclcpp/ParameterSet.h"

#include "PDSPThinSliceFitter.h"

#include "TFile.h"

int main(int argc, char ** argv){

  std::string fcl_file;
  std::string output_file;
  // Options to run
  for (int iArg = 1; iArg < argc; iArg++) {
    if (!strcasecmp(argv[iArg],"-c")) {
     fcl_file = argv[++iArg];
    }
    if (!strcasecmp(argv[iArg],"-o")) {
      output_file = argv[++iArg];
    }
    if (!strcasecmp(argv[iArg],"-h")) {
      std::cout << "Usage: runPDSPThinSliceFit -c fclfile.fcl " << 
                    "-o outputfile.root " << std::endl;
      return 1;
    }
  }

  protoana::PDSPThinSliceFitter * fit
      = new protoana::PDSPThinSliceFitter(fcl_file, output_file);
  fit->InitializeMCSamples();
  fit->BuildMCSamples();
  fit->RunFitAndSave();

  return 0;
}