| 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 |
|
|
} |