30int main(
int argc,
char **argv)
49 JParser<> zap(
"Auxiliary program to project 2D histograms.");
51 zap[
'f'] =
make_field(inputFile,
"<input file>:<object name>");
52 zap[
'P'] =
make_field(project,
"projection") =
' ',
'x',
'X',
'y',
'Y';
57 zap[
'F'] =
make_field(format,
"format, e.g. \"%s %i\" or \"%s %f %f\"") =
"";
58 zap[
'O'] =
make_field(option,
"option, see TH2::Projection(X|Y)") =
"";
63 catch(
const exception &error) {
64 FATAL(error.what() << endl);
67 const bool px = (project ==
'x' || project ==
'X' || !Y.empty());
68 const bool py = (project ==
'y' || project ==
'Y' || !X.empty());
71 FATAL(
"Invalid operation: "
72 << (px ?
"" :
"no") <<
" X projection " <<
" and "
73 << (py ?
"" :
"no") <<
" Y projection " << endl);
80 if (px) { format +=
"px"; }
81 if (py) { format +=
"py"; }
83 if (X.empty() && Y.empty())
91 for (vector<JRootObjectID>::const_iterator input = inputFile.begin(); input != inputFile.end(); ++input) {
93 DEBUG(
"Input: " << *input << endl);
95 TDirectory* dir = getDirectory(*input);
98 ERROR(
"File: " << input->getFullFilename() <<
" not opened." << endl);
102 const TRegexp regexp(input->getObjectName());
104 TIter iter(dir->GetListOfKeys());
106 for (TKey* key; (key = (TKey*) iter.Next()) != NULL; ) {
108 const TString tag(key->GetName());
110 DEBUG(
"Key: " << tag <<
" match = " << tag.Contains(regexp) << endl);
114 if (tag.Contains(regexp) && isTObject(key)) {
116 TH2* h2 =
dynamic_cast<TH2*
>(key->ReadObj());
124 for (Int_t i = (overflow ? 0 : 1); i <= h2->GetYaxis()->GetNbins() + (overflow ? 1 : 0); ++i) {
125 listOfObjects.push_back(h2->ProjectionX(TString::Format(format.c_str(), h2->GetName(), i),
132 for (Int_t i = 0; i != (Int_t) Y.size(); ++i) {
133 listOfObjects.push_back(h2->ProjectionX(TString::Format(format.c_str(), h2->GetName(), Y[i].getLowerLimit(), Y[i].getUpperLimit()),
134 h2->GetYaxis()->FindBin(Y[i].getLowerLimit()),
135 h2->GetYaxis()->FindBin(Y[i].getUpperLimit()) - 1, option.c_str()));
143 for (Int_t i = (overflow ? 0 : 1); i <= h2->GetXaxis()->GetNbins() + (overflow ? 1 : 0); ++i) {
144 listOfObjects.push_back(h2->ProjectionY(TString::Format(format.c_str(), h2->GetName(), i), i, i, option.c_str()));
149 for (Int_t i = 0; i != (Int_t) X.size(); ++i) {
150 listOfObjects.push_back(h2->ProjectionY(TString::Format(format.c_str(), h2->GetName(), X[i].getLowerLimit(), X[i].getUpperLimit()),
151 h2->GetXaxis()->FindBin(X[i].getLowerLimit()),
152 h2->GetXaxis()->FindBin(X[i].getUpperLimit()) - 1, option.c_str()));
161 if (!listOfObjects.empty()) {
165 for (vector<TObject*>::const_iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {