1 |
/* |
2 |
* TrkDedxHCut.h |
3 |
* |
4 |
* Created on: 15-mar-2010 |
5 |
* Author: S. Ricciarini |
6 |
*/ |
7 |
|
8 |
/*! @file TrkDedxHCut.h The TrkDedxHCut class definition file */ |
9 |
|
10 |
#ifndef TRKDEDXHCUT_H_ |
11 |
#define TRKDEDXHCUT_H_ |
12 |
|
13 |
#include "../../PamCutBase/PamCutBase.h" |
14 |
|
15 |
/*! @brief The dE/dX cut to select H and anti-H. |
16 |
* This cut discards all the events whose rigidity modulus and dE/dX (both given by TRK) are outside a pre-defined band (lower band: by Nico De Simone; upper band: by Nicola Mori). |
17 |
* Here the rigidity is defined as pc/Ze (GV) where Z is the particle charge (WITH SIGN) and p the momentum modulus: therefore the rigidity can be positive or negative. |
18 |
* CUT DEPENDENCIES: TrkPhysSin for object trkTrack |
19 |
*/ |
20 |
class TrkDedxHCut: public PamCut { |
21 |
|
22 |
public: |
23 |
|
24 |
/*! @brief Constructor. |
25 |
* |
26 |
* @param cutName The cut's name. |
27 |
*/ |
28 |
TrkDedxHCut(const char *cutName, const char *trkAlg) : |
29 |
PamCut(cutName), _trkAlg(trkAlg) { |
30 |
|
31 |
} |
32 |
|
33 |
/*! @brief Destructor. */ |
34 |
~TrkDedxHCut() { |
35 |
|
36 |
} |
37 |
|
38 |
/*! @brief The rigidity check. |
39 |
* |
40 |
* @param event The event to analyze. |
41 |
* @return #CUTOK if the (rigidity modulus , dE/dX) point is inside the pre-defined band. |
42 |
* @return 0 otherwise. |
43 |
*/ |
44 |
int Check(PamLevel2 *event); |
45 |
|
46 |
private: |
47 |
|
48 |
const char *_trkAlg; |
49 |
|
50 |
// lower cut line |
51 |
static const int _nLow = 23; |
52 |
static const Double_t _xLow[]; |
53 |
static const Double_t _yLow[]; |
54 |
|
55 |
}; |
56 |
#endif /* TRKDEDXHCUT_H_ */ |