--- PamCut/CollectionActions/Histo1DActions/TrkSigmaDeflHistoAction/TrkSigmaDeflHistoAction.cpp 2009/09/25 15:36:43 1.1 +++ PamCut/CollectionActions/Histo1DActions/TrkSigmaDeflHistoAction/TrkSigmaDeflHistoAction.cpp 2010/05/06 17:09:55 1.2.2.1 @@ -9,9 +9,9 @@ #include "TrkSigmaDeflHistoAction.h" -TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, - TString mode, bool outRoot, bool outText, TString title) : - Histo1DAction (actionName, title, outFileBase, mode, outRoot, outText) { +TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, TString mode, + float sigmaDeflFactor, bool outRoot, bool outText, TString title) : + Histo1DAction (actionName, title, outFileBase, mode, outRoot, outText), _sigmaDeflFactor(sigmaDeflFactor) { } void TrkSigmaDeflHistoAction::OnGood(PamLevel2 *event) { @@ -21,7 +21,7 @@ // Bin filling: sigmaDefl must be less than (<) 1/(bin maximum) UInt_t i = 0; - 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 _histo[i]++; i++; }