| 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_ */ |