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 |
|
|
TrkSigmaDeflHistoAction::TrkSigmaDeflHistoAction(const char *actionName, TString outFileBase, |
13 |
|
|
TString mode, bool outRoot, bool outText, TString title) : |
14 |
|
|
Histo1DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText) { |
15 |
|
|
} |
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 |
|
|
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 |
25 |
|
|
_histo[i]++; |
26 |
|
|
i++; |
27 |
|
|
} |
28 |
|
|
|
29 |
|
|
} |