| 9 |
|
|
| 10 |
#include "TrkSigmaDeflHistoAction.h" |
#include "TrkSigmaDeflHistoAction.h" |
| 11 |
|
|
| 12 |
TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, |
TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, TString mode, |
| 13 |
TString mode, bool outRoot, bool outText, TString title) : |
float sigmaDeflFactor, bool outRoot, bool outText, TString title) : |
| 14 |
Histo1DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText) { |
Histo1DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _sigmaDeflFactor(sigmaDeflFactor) { |
| 15 |
} |
} |
| 16 |
|
|
| 17 |
void TrkSigmaDeflHistoAction::OnGood(PamLevel2 *event) { |
void TrkSigmaDeflHistoAction::OnGood(PamLevel2 *event) { |
| 21 |
// Bin filling: sigmaDefl must be less than (<) 1/(bin maximum) |
// Bin filling: sigmaDefl must be less than (<) 1/(bin maximum) |
| 22 |
|
|
| 23 |
UInt_t i = 0; |
UInt_t i = 0; |
| 24 |
while ((sigmaDefl < 1. / _bins[i + 1]) && (i <= _bins.size() - 1)) { // rigidity bin i maximum = _binning[i+1]. Exit while loop when i exceeds the bin number = _binning.size()-1 |
while ((sigmaDefl * _sigmaDeflFactor < 1. / _bins[i + 1]) && (i <= _bins.size() - 1)) { // rigidity bin i maximum = _binning[i+1]. Exit while loop when i exceeds the bin number = _binning.size()-1 |
| 25 |
_histo[i]++; |
_histo[i]++; |
| 26 |
i++; |
i++; |
| 27 |
} |
} |