106 |
*/ |
*/ |
107 |
void SetYAxis(TString label, unsigned int nBins, float min, float max, bool logBinning = false); |
void SetYAxis(TString label, unsigned int nBins, float min, float max, bool logBinning = false); |
108 |
|
|
109 |
|
/*! @brief Sets the ROOT histogram's title. |
110 |
|
* |
111 |
|
* @param title The histogram title as it will appear on the histogram itself. |
112 |
|
*/ |
113 |
|
void SetTitle(TString &title) { |
114 |
|
_title = title; |
115 |
|
} |
116 |
/*! @brief Sets up the histogram |
/*! @brief Sets up the histogram |
117 |
* |
* |
118 |
* This routine effectively prepares the histogram, after the desired parameters has been set by #SetXAxis() and #SetYAxis(). |
* This routine effectively prepares the histogram, after the desired parameters has been set by #SetXAxis() and #SetYAxis(). |
119 |
* |
* |
120 |
* @param events Pointer to PamLevel2 events (unused). |
* @param events Pointer to PamLevel2 events (unused). |
121 |
*/ |
*/ |
|
|
|
|
/*! @brief Sets the ROOT histogram's title. */ |
|
|
void SetTitle(TString &title){ |
|
|
_title = title; |
|
|
} |
|
122 |
void Setup(PamLevel2 *events) { |
void Setup(PamLevel2 *events) { |
123 |
CollectionAction::Setup(events); |
CollectionAction::Setup(events); |
124 |
_InitHistos(); |
_InitHistos(); |
143 |
return _histo; |
return _histo; |
144 |
} |
} |
145 |
|
|
146 |
/*! Fills the ROOT and the vector histogram. */ |
/*! Fills the ROOT and the vector histogram. |
147 |
|
* |
148 |
|
* @param xValue The value of the X coordinate associated to the event. |
149 |
|
* @param yValue The value of the Y coordinate associated to the event. |
150 |
|
* @param weight The weight which will be applied to the event. |
151 |
|
*/ |
152 |
void Fill(double xValue, double yValue, double weight = 1.); |
void Fill(double xValue, double yValue, double weight = 1.); |
153 |
|
|
154 |
/*! @brief Gets the X overflow histogram. |
/*! @brief Gets the X overflow histogram. |
243 |
|
|
244 |
protected: |
protected: |
245 |
|
|
246 |
|
/*! @brief The vector containing the limits of the X bins(from lower to higher). */ |
247 |
std::vector<float> _xBins; |
std::vector<float> _xBins; |
248 |
|
/*! @brief The vector containing the limits of the Y bins(from lower to higher). */ |
249 |
std::vector<float> _yBins; |
std::vector<float> _yBins; |
250 |
|
/*! @brief A matrix containing the value of the histogram for each X-Y bin. */ |
251 |
SimpleMatrix<HistoType> _histo; |
SimpleMatrix<HistoType> _histo; |
252 |
|
/*! @brief The ROOT histogram. */ |
253 |
TH2 *_rootHisto; |
TH2 *_rootHisto; |
254 |
|
|
255 |
private: |
private: |
273 |
_rootHisto = NULL; |
_rootHisto = NULL; |
274 |
} |
} |
275 |
|
|
276 |
|
// Specializations for _CreateHistos(). See Histo2DAction.cpp |
277 |
|
template<> |
278 |
|
void Histo2DAction<Int_t>::_CreateHisto(); |
279 |
|
|
280 |
|
template<> |
281 |
|
void Histo2DAction<Float_t>::_CreateHisto(); |
282 |
|
|
283 |
|
template<> |
284 |
|
void Histo2DAction<Double_t>::_CreateHisto(); |
285 |
|
|
286 |
template<class HistoType> |
template<class HistoType> |
287 |
void Histo2DAction<HistoType>::_InitHistos() { |
void Histo2DAction<HistoType>::_InitHistos() { |
288 |
|
|
289 |
_CreateHisto(); |
_CreateHisto(); |
290 |
|
|
291 |
if (_xBins.size() < 2){ |
if (_xBins.size() < 2) { |
292 |
_xBins.resize(2); |
_xBins.resize(2); |
293 |
_xBins[0] = 0.; |
_xBins[0] = 0.; |
294 |
_xBins[1] = 1.; |
_xBins[1] = 1.; |
295 |
} |
} |
296 |
|
|
297 |
if (_yBins.size() < 2){ |
if (_yBins.size() < 2) { |
298 |
_yBins.resize(2); |
_yBins.resize(2); |
299 |
_yBins[0] = 0.; |
_yBins[0] = 0.; |
300 |
_yBins[1] = 1.; |
_yBins[1] = 1.; |
335 |
|
|
336 |
template<class HistoType> |
template<class HistoType> |
337 |
Histo2DAction<HistoType>::~Histo2DAction() { |
Histo2DAction<HistoType>::~Histo2DAction() { |
338 |
|
|
339 |
|
delete _rootHisto; |
340 |
|
_rootHisto = NULL; |
341 |
} |
} |
342 |
|
|
343 |
template<class HistoType> |
template<class HistoType> |
525 |
return; |
return; |
526 |
} |
} |
527 |
else { |
else { |
528 |
if (yBin == (int)_yBins.size()) { |
if (yBin == (int) _yBins.size()) { |
529 |
_xUnderYOverflow += (HistoType) weight; |
_xUnderYOverflow += (HistoType) weight; |
530 |
return; |
return; |
531 |
} |
} |
536 |
} |
} |
537 |
} |
} |
538 |
|
|
539 |
if (xBin == (int)_xBins.size()) { |
if (xBin == (int) _xBins.size()) { |
540 |
if (yBin == -1) { |
if (yBin == -1) { |
541 |
_xOverYUnderflow += (HistoType) weight; |
_xOverYUnderflow += (HistoType) weight; |
542 |
return; |
return; |
543 |
} |
} |
544 |
else { |
else { |
545 |
if (yBin ==(int) _yBins.size()) { |
if (yBin == (int) _yBins.size()) { |
546 |
_xOverYOverflow += (HistoType) weight; |
_xOverYOverflow += (HistoType) weight; |
547 |
return; |
return; |
548 |
} |
} |
558 |
return; |
return; |
559 |
} |
} |
560 |
|
|
561 |
if (yBin == (int)_yBins.size()) { |
if (yBin == (int) _yBins.size()) { |
562 |
_yOverflow[xBin] += (HistoType) weight; |
_yOverflow[xBin] += (HistoType) weight; |
563 |
return; |
return; |
564 |
} |
} |