/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus
*
* MAUS is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* MAUS is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with MAUS. If not, see .
*
*/
/* Simple code used to open a reconstructed MAUS ROOT file,
* iterate through spills and pass a spill to EVExporter
*/
#include
#include
#include
#include "DataStructure/Spill.hh"
#include "DataStructure/Data.hh"
#include "JsonCppStreamer/IRStream.hh"
#include "src/utilities/event-viewer/EVExporter.hh"
int ParseArguments(int argc, char* argv[]);
void ShowUsage(std::string arg);
std::string FILE_NAME;
int SPILL_NUM;
int SPILL_RANGE[2];
int EVENT_SELECTION;
bool VERBOSE;
int main(int argc, char* argv[]) {
SPILL_NUM = 0;
SPILL_RANGE[0] = 0;
SPILL_RANGE[1] = 0;
EVENT_SELECTION = 0;
VERBOSE = false;
if (!ParseArguments(argc, argv))
return 1;
// - start reading data:
MAUS::Data data;
irstream infile(FILE_NAME.c_str(), "Spill");
// - iterate over spills:
while (infile >> readEvent != NULL) {
infile >> branchName("data") >> data;
MAUS::Spill *spill = data.GetSpill();
if (spill != NULL && spill->GetDaqEventType() == "physics_event") {
int spillNumber = spill->GetSpillNumber();
// - instantiate EVExported and pass spill to it
if (SPILL_NUM != 0) {
if (spillNumber == SPILL_NUM) {
EventViewer::EVExporter *exp = new EventViewer::EVExporter(spill);
exp->ReadAllEvents(EVENT_SELECTION, VERBOSE);
return 0;
}
} else if (SPILL_RANGE[1] != 0) { // can not pass both --spill and --spill_range arguments
if (SPILL_RANGE[1] > SPILL_RANGE[0]) {
if (spillNumber >= SPILL_RANGE[0] && spillNumber <= SPILL_RANGE[1]) {
EventViewer::EVExporter *exp = new EventViewer::EVExporter(spill);
exp->ReadAllEvents(EVENT_SELECTION, VERBOSE);
} else if (spillNumber > SPILL_RANGE[1]) {
return 0;
}
} else {
std::cerr << std::endl << "ERROR: Invalid spill range!" << std::endl << std::endl;
return 2;
}
} else if (SPILL_NUM == 0 && SPILL_RANGE[0] == 0 && SPILL_RANGE[1] == 0) {
EventViewer::EVExporter *exp = new EventViewer::EVExporter(spill);
exp->ReadAllEvents(EVENT_SELECTION, VERBOSE);
}
}
}
infile.close();
return 0;
}
int ParseArguments(int argc, char* argv[]) {
FILE_NAME = "";
if (argc < 2) {
ShowUsage(argv[0]);
return 0;
}
for (int i = 1; i < argc; ++i) {
std::string arg = argv[i];
if (arg == "-h" || arg == "--help") {
ShowUsage(argv[0]);
return 0;
}
if (arg == "-v" || arg == "--verbose") {
VERBOSE = true;
}
if (arg == "-f" || arg == "--file")
FILE_NAME = argv[i+1];
if (arg == "-s" || arg == "--spill") {
if (i+1 < argc) {
SPILL_NUM = atoi(argv[i+1]);
} else {
std::cerr << "\n" << "WARNING: You must provide argument to --spill option" << "\n\n";
return 0;
}
}
if (arg == "-sr" || arg == "--spill_range") {
if (i+1 < argc) {
SPILL_RANGE[0] = atoi(argv[i+1]);
} else {
std::cerr << "\n" << "WARNING: You must provide argument to --spill_range option" << "\n\n";
return 0;
}
if (i+2 < argc) {
SPILL_RANGE[1] = atoi(argv[i+2]);
} else {
std::cerr << "\n" << "WARNING: You must provide second argument to --spill_range option"
<< "\n\n";
return 0;
}
}
if (arg == "-es" || arg == "--event_selection")
if (i+1 < argc) {
EVENT_SELECTION = atoi(argv[i+1]);
} else {
std::cerr << "\n" << "WARNING: You must provide argument to --event_selection option"
<< "\n\n";
return 0;
}
}
if (FILE_NAME == "") {
std::cerr << "\n" << "ERROR! You must provide INPUT_FILE argument!" << "\n\n";
ShowUsage(argv[0]);
return 0;
}
return 1;
}
void ShowUsage(std::string progName) {
std::cerr << "\n" << "USAGE: " << progName << "