260 |
/*! @brief The ROOT histogram. */ |
/*! @brief The ROOT histogram. */ |
261 |
TH2 *_rootHisto; |
TH2 *_rootHisto; |
262 |
|
|
263 |
|
TString _outFileBase; |
264 |
|
TString _mode; |
265 |
|
TString _title, _xLabel, _yLabel; |
266 |
|
|
267 |
private: |
private: |
268 |
|
|
269 |
vector<HistoType> _xUnderflow, _xOverflow, _yUnderflow, _yOverflow; |
vector<HistoType> _xUnderflow, _xOverflow, _yUnderflow, _yOverflow; |
270 |
HistoType _xUnderYUnderflow, _xOverYOverflow, _xUnderYOverflow, _xOverYUnderflow; |
HistoType _xUnderYUnderflow, _xOverYOverflow, _xUnderYOverflow, _xOverYUnderflow; |
|
TString _outFileBase; |
|
|
TString _mode; |
|
|
TString _title, _xLabel, _yLabel; |
|
271 |
bool _outRoot; |
bool _outRoot; |
272 |
bool _outText; |
bool _outText; |
|
|
|
273 |
void _CreateHisto(); |
void _CreateHisto(); |
274 |
void _InitHistos(); |
void _InitHistos(); |
275 |
}; |
}; |
326 |
/* The row index (first) corresponds to the position on the vertical (Y) axis. */ |
/* The row index (first) corresponds to the position on the vertical (Y) axis. */ |
327 |
_histo.Resize(_yBins.size() - 1, _xBins.size() - 1); |
_histo.Resize(_yBins.size() - 1, _xBins.size() - 1); |
328 |
|
|
329 |
_xUnderflow.resize(_yBins.size()); |
_xUnderflow.resize(_yBins.size() - 1); |
330 |
_xOverflow.resize(_yBins.size()); |
_xOverflow.resize(_yBins.size() - 1); |
331 |
_yUnderflow.resize(_xBins.size()); |
_yUnderflow.resize(_xBins.size() - 1); |
332 |
_yOverflow.resize(_xBins.size()); |
_yOverflow.resize(_xBins.size() - 1); |
333 |
|
|
334 |
} |
} |
335 |
|
|
343 |
template<class HistoType> |
template<class HistoType> |
344 |
Histo2DAction<HistoType>::Histo2DAction(const char *actionName, TString title, TString outFileBase, TString mode, |
Histo2DAction<HistoType>::Histo2DAction(const char *actionName, TString title, TString outFileBase, TString mode, |
345 |
bool outRoot, bool outText) : |
bool outRoot, bool outText) : |
346 |
CollectionAction(actionName), _xBins(0), _yBins(0), _histo(0, 0), _rootHisto(NULL), _xUnderflow(0), _xOverflow(0), |
CollectionAction(actionName), _xBins(0), _yBins(0), _histo(0, 0), _rootHisto(NULL), _outFileBase(outFileBase), _mode( |
347 |
_yUnderflow(0), _yOverflow(0), _outFileBase(outFileBase), _mode(mode), _title(title), _xLabel(""), _yLabel(""), |
mode), _title(title), _xLabel(""), _yLabel(""), _xUnderflow(0), _xOverflow(0), _yUnderflow(0), _yOverflow(0), |
348 |
_outRoot(outRoot), _outText(outText) { |
_xUnderYUnderflow((HistoType) 0), _xOverYOverflow((HistoType) 0), _xUnderYOverflow((HistoType) 0), |
349 |
|
_xOverYUnderflow((HistoType) 0), _outRoot(outRoot), _outText(outText) { |
350 |
|
|
351 |
} |
} |
352 |
|
|
528 |
else { |
else { |
529 |
if (yBin == (int) _yBins.size()) { |
if (yBin == (int) _yBins.size()) { |
530 |
_xUnderYOverflow += (HistoType) weight; |
_xUnderYOverflow += (HistoType) weight; |
531 |
|
cout << "_xUnderYOverflow: " << GetXUnderYOverflow() << ", weight: " << weight << ", (HistoType) weight: " |
532 |
|
<< (HistoType) weight << endl; |
533 |
return; |
return; |
534 |
} |
} |
535 |
else { |
else { |
586 |
ofstream outTextFile((_outFileBase + ".txt").Data(), ios_base::out); |
ofstream outTextFile((_outFileBase + ".txt").Data(), ios_base::out); |
587 |
for (unsigned int i = 0; i < _histo.GetNRows(); i++) { |
for (unsigned int i = 0; i < _histo.GetNRows(); i++) { |
588 |
for (unsigned int j = 0; j < _histo.GetNCols(); j++) { |
for (unsigned int j = 0; j < _histo.GetNCols(); j++) { |
589 |
outTextFile << _histo[i][j] << " "; |
outTextFile << setw(7) << _histo[i][j] << " "; |
590 |
} |
} |
591 |
outTextFile << "\n"; |
outTextFile << "\n"; |
592 |
} |
} |
|
outTextFile << endl; |
|
593 |
outTextFile.close(); |
outTextFile.close(); |
594 |
} |
} |
595 |
} |
} |