| 1 |
/* |
| 2 |
* FluxHistoAction.h |
| 3 |
* |
| 4 |
* Created on: 20-mag-2009 |
| 5 |
* Author: Nicola Mori |
| 6 |
*/ |
| 7 |
|
| 8 |
/*! @file FluxHistoAction.h The FluxHistoAction class declaration file */ |
| 9 |
|
| 10 |
#ifndef FLUXHISTOACTION_H_ |
| 11 |
#define FLUXHISTOACTION_H_ |
| 12 |
|
| 13 |
#include "../CollectionAction/CollectionAction.h" |
| 14 |
#include <TH1I.h> |
| 15 |
#include <stdint.h> |
| 16 |
|
| 17 |
/*! @brief An action that fills a flux histogram. |
| 18 |
* |
| 19 |
* This action reads a rigidity binning from a file and fills a flux histogram (text and ROOT format). |
| 20 |
* Currently, it fills an event-count histogram: in future implementations, it could also |
| 21 |
* read live time, efficiency and geometrical factor from external files and fill a |
| 22 |
* true flux histogram. |
| 23 |
*/ |
| 24 |
class FluxHistoAction: public CollectionAction { |
| 25 |
|
| 26 |
public: |
| 27 |
/*! @brief Constructor. |
| 28 |
* |
| 29 |
* outFileBase is the base name for output file: #Finalize will add .txt for ASCII output |
| 30 |
* and .root for ROOT output. outFileBase has to contain the path (otherwise, files will be |
| 31 |
* saved in the current directory). |
| 32 |
* The file containing the rigidity bins must be a text file. It must contain both the |
| 33 |
* lower and upper limits of the rigidity axis, so that if it contains N values it |
| 34 |
* defines a set of N-1 bins. |
| 35 |
* |
| 36 |
* @param actionName The action's name. |
| 37 |
* @param outFileBase The output file base name. |
| 38 |
* @param rigBinsFile The file containing the rigidity bins. |
| 39 |
* @param mode The mode of ROOT file creation (see documentation of TFile constructor |
| 40 |
* in ROOT's reference guide). In current implementation the text file output |
| 41 |
* does not support append. |
| 42 |
*/ |
| 43 |
FluxHistoAction(const char *actionName, TString outFileBase, TString rigBinsFile, TString mode = "UPDATE"); |
| 44 |
|
| 45 |
/*! @brief Destructor */ |
| 46 |
~FluxHistoAction() { |
| 47 |
} |
| 48 |
|
| 49 |
/*! @brief Fills histogram with the selected event. |
| 50 |
* |
| 51 |
* @param event The selected event. |
| 52 |
*/ |
| 53 |
void OnGood(PamLevel2 *event); |
| 54 |
|
| 55 |
/*! @brief Writes the histogram to the output files (ASCII and ROOT). |
| 56 |
* |
| 57 |
* The output consists of a text file and of a ROOT file where the 1-dimensional rigidity |
| 58 |
* histogram (TH1F) is saved. The format of the text output is: (lower bin limit) (upper bin limit) (\# of events). |
| 59 |
*/ |
| 60 |
void Finalize(); |
| 61 |
|
| 62 |
private: |
| 63 |
|
| 64 |
TString _outFileBase; |
| 65 |
vector<float> _bins; |
| 66 |
TH1I _rootHisto; |
| 67 |
vector<uint64_t> _textHisto; |
| 68 |
TString _mode; |
| 69 |
|
| 70 |
#ifdef DEBUGPAMCUT |
| 71 |
int _outUp, _outDown; |
| 72 |
#endif |
| 73 |
|
| 74 |
}; |
| 75 |
|
| 76 |
#endif /* FLUXHISTOACTION_H_ */ |