Histo1DActions/TrkSigmaDeflHistoAction/TrkSigmaDeflHistoAction.cpp

/*
 * TrkSigmaDeflHistoAction.cpp
 *
 *  Created on: 2009-06-17
 *      Author: S. Ricciarini
 */

/*! @file TrkSigmaDeflHistoAction.cpp The TrkSigmaDeflHistoAction class implementation file. */

#include "TrkSigmaDeflHistoAction.h"

TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, TString mode,
    float sigmaDeflFactor, bool outRoot, bool outText, TString title) :
  Histo1DAction<Float_t> (actionName, title, outFileBase, mode, outRoot, outText), _sigmaDeflFactor(sigmaDeflFactor) {
}

void TrkSigmaDeflHistoAction::OnGood(PamLevel2 *event) {

  float sigmaDefl = pow(event->GetTrack(0)->GetTrkTrack()->coval[4][4], 0.5);

  // Bin filling: sigmaDefl must be less than (<) 1/(bin maximum)

  UInt_t i = 0;
  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++;
  }

}
1	pam-fi	1.1	/*
2			* TrkSigmaDeflHistoAction.cpp
3			*
4			* Created on: 2009-06-17
5			* Author: S. Ricciarini
6			*/
7
8			/! @file TrkSigmaDeflHistoAction.cpp The TrkSigmaDeflHistoAction class implementation file. /
9
10			#include "TrkSigmaDeflHistoAction.h"
11
12	pam-fi	1.2	TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, TString mode,
13			float sigmaDeflFactor, bool outRoot, bool outText, TString title) :
14	pam-fi	1.3	Histo1DAction<Float_t> (actionName, title, outFileBase, mode, outRoot, outText), _sigmaDeflFactor(sigmaDeflFactor) {
15	pam-fi	1.1	}
16
17			void TrkSigmaDeflHistoAction::OnGood(PamLevel2 *event) {
18
19			float sigmaDefl = pow(event->GetTrack(0)->GetTrkTrack()->coval[4][4], 0.5);
20
21			// Bin filling: sigmaDefl must be less than (<) 1/(bin maximum)
22
23			UInt_t i = 0;
24	pam-fi	1.2	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	pam-fi	1.1	_histo[i]++;
26			i++;
27			}
28
29			}