3 |
* |
* |
4 |
* Created on: 13-mar-2009 |
* Created on: 13-mar-2009 |
5 |
* Author: Nicola Mori, S. Ricciarini |
* Author: Nicola Mori, S. Ricciarini |
6 |
|
* Reworked on: 21-jul-2011 |
7 |
|
* Author: Nicola Mori, Valeria Di Felice |
8 |
*/ |
*/ |
9 |
|
|
10 |
/*! @file TofDedxS1Cut.h The TofDedxS1Cut class definition file */ |
/*! @file TofDedxS1Cut.h The TofDedxS1Cut class definition file */ |
15 |
#include "../../PamCutBase/PamCutBase.h" |
#include "../../PamCutBase/PamCutBase.h" |
16 |
|
|
17 |
|
|
18 |
/*! @brief The TofQual cut. |
/*! @brief A cut on energy release on S1. |
|
* This cut is done for each S1 layer: if there is one and only one hit paddle for the layer, then the dE/dX of the paddle (referred to a specified track) must be less than a fixed threshold. NOTE: 'hit paddle' means that both TDC signals are present AND they are both 'good' (according to TOF definition of 'good' TDC signal) |
|
19 |
* |
* |
20 |
|
* This cut discards events whose mean dE/dx on S11+S12 is greater than a specified functional form. |
21 |
|
* Currently, this functional form is: |
22 |
|
* |
23 |
|
* <dE/dx> < (3.102 - 2.6*TMath::Log(R))*(R<1.05) + 3.*(R=1.05) |
24 |
|
* |
25 |
|
* Since this cut uses the rigidity, it is necessary to enforce the existence of a physical track |
26 |
|
* using #TrkPhsinCut before using this cut. |
27 |
*/ |
*/ |
28 |
|
|
29 |
class TofDedxS1Cut: public PamCut { |
class TofDedxS1Cut: public PamCut { |
31 |
public: |
public: |
32 |
/*! @brief Constructor. |
/*! @brief Constructor. |
33 |
* |
* |
|
* The parameters are: maximum dE/dX for S1 layers and the TOF-index for the track which is necessary to input to derive the dE/dX for S1 layers. |
|
|
* |
|
34 |
* @param cutName The cut's name. |
* @param cutName The cut's name. |
|
* @param maxDedxS1 The maximum dE/dX for each S1 layer, for which the event is accepted |
|
|
* @param trackTof Pointer to the TOF-index indicating the track to be used when evaluating the dE/dX for S1 layers (default value for TOF_index is 0, corresponding to the TOF stand-alone track [no TRK information to build the TOF stand-alone track]). The pointer to the TOF_index for TRK-track as given by TrkPhSinCut can be used. |
|
35 |
* |
* |
36 |
*/ |
*/ |
37 |
TofDedxS1Cut(const char *cutName, float maxDedxS1, int *trackTof=NULL): |
TofDedxS1Cut(const char *cutName): |
38 |
PamCut(cutName), _maxDedxS1(maxDedxS1), _trackTof(trackTof) { |
PamCut(cutName){ |
39 |
} |
} |
40 |
/*! @brief Destructor. */ |
/*! @brief Destructor. */ |
41 |
~TofDedxS1Cut() { |
~TofDedxS1Cut() { |
44 |
/*! @brief The TofDedxS1 check. |
/*! @brief The TofDedxS1 check. |
45 |
* |
* |
46 |
* @param event The event to analyze. |
* @param event The event to analyze. |
47 |
* @return #CUTOK if for each layer of S1 with one and only one hit paddle, the paddle dE/dX is less than a fixed threshold |
* @return #CUTOK if <dE/dx> is lesser that the threshold function. |
48 |
* @return 0 otherwise |
* @return 0 otherwise |
49 |
*/ |
*/ |
50 |
int Check(PamLevel2 *event); |
int Check(PamLevel2 *event); |
51 |
|
|
52 |
private: |
private: |
53 |
|
|
|
float _maxDedxS1; |
|
|
int *_trackTof; |
|
|
int _nHitPadPair[3]; // Number of hit paddles on each pair of layers |
|
|
|
|
54 |
}; |
}; |
55 |
|
|
56 |
#endif /* TOFDEDXS1CUT_H_ */ |
#endif /* TOFDEDXS1CUT_H_ */ |