16 |
* This cut discards events where TOF dE/dx vs. Beta are outside the helium "corridor". This corridor is currently |
* This cut discards events where TOF dE/dx vs. Beta are outside the helium "corridor". This corridor is currently |
17 |
* defined by: |
* defined by: |
18 |
* |
* |
19 |
* dE/dx > 1 / (0.2 * (beta - 0.05)) - 2. |
* dE/dx > 1 / (p[0] * pow(beta, p[4]) + p[1]) + p[2] + p[3]*beta |
20 |
* dE/dx < 1 / (0.25 * (beta - 0.1)) |
* dE/dx < 1 / (p[5] * pow((*beta), p[9]) + p[6]) + p[7] + p[8]*(*beta); |
|
* dE/dx > 11. * beta - 7.7 |
|
21 |
* |
* |
22 |
* dE/dx is the mean of the energy releases on S11, S12, S21, S22. |
* where p[i] are calibration parameters. Currently, they are available for each single layer; the beta used is the beta[12] corresponding |
23 |
|
* to the TOF standalone track or to the TRK physical track, depending on the value of the standAlone parameter of the constructor. The |
24 |
|
* cut parameters have been obtained using the standalone beta, but they have proven to be good also for beta from tracker, so a single |
25 |
|
* set of parameter will be used for both cases. |
26 |
|
* If the required combination of planes has no available calibration, the one for S21 will be used (but it probably won't work well...). |
27 |
|
* In particular, there's no currently available calibration for combinations of layers. |
28 |
*/ |
*/ |
29 |
class TofDedxHeCut: public PamCut { |
class TofDedxHeCut: public PamCut { |
30 |
|
|
32 |
/*! @brief Constructor. |
/*! @brief Constructor. |
33 |
* |
* |
34 |
* @param cutName The cut's name. |
* @param cutName The cut's name. |
35 |
|
* @param layer The layer to use. Select it using values like #S11 (for single layers) or S11+S12 (for the mean of many layers). |
36 |
* @param standAlone If true, the cut will use the standalone measure of beta; otherwise it will use the |
* @param standAlone If true, the cut will use the standalone measure of beta; otherwise it will use the |
37 |
* physical track. |
* physical track. |
38 |
*/ |
*/ |
39 |
TofDedxHeCut(const char *cutName, bool standAlone = true) : |
TofDedxHeCut(const char *cutName, unsigned int layer, bool standAlone = true); |
|
PamCut(cutName), _standAlone(standAlone) { |
|
|
} |
|
40 |
/*! @brief Destructor. */ |
/*! @brief Destructor. */ |
41 |
~TofDedxHeCut() { |
~TofDedxHeCut() { |
42 |
|
|
43 |
} |
} |
44 |
|
|
45 |
/*! @brief The dE/dx vs. beta Helium check. |
/*! @brief The dE/dx vs. beta helium check. |
46 |
* |
* |
47 |
* @param event The event to analyze. |
* @param event The event to analyze. |
48 |
* @return #CUTOK if event is inside helium "corridor". |
* @return #CUTOK if event is inside helium "corridor". |
52 |
|
|
53 |
private: |
private: |
54 |
|
|
55 |
|
unsigned int _layer; |
56 |
bool _standAlone; |
bool _standAlone; |
57 |
|
const float *_p; |
58 |
|
// Single layers standalone |
59 |
|
static const float _S11_SA_params[10], _S12_SA_params[10], _S21_SA_params[10], _S22_SA_params[10], _S31_SA_params[10], |
60 |
|
_S32_SA_params[10]; |
61 |
}; |
}; |
62 |
|
|
63 |
#endif /* TOFDEDXHECUT_H_ */ |
#endif /* TOFDEDXHECUT_H_ */ |