TofCuts/TofDedxHCut/TofDedxHCut.h

/*
 * TofDedxHCut.h
 *
 *  Created on: 5 apr 2011
 *      Author: Nicola Mori
 */

/*! @file TofDedxHCut.h  The TofDedxHCut.h class definition file. */

#ifndef TOFDEDXHCUT_H_
#define TOFDEDXHCUT_H_

#include "../../PamCutBase/PamCutBase.h"

/*! @brief Cutting a band dE/dx vs. Beta for hydrogen.
 * This cut discards events where TOF dE/dx vs. Beta are outside the hydrogen "corridor". This corridor is currently
 * defined by:
 *
 * dE/dx > _p[0] / (_p[1] * beta + _p[2]) + _p[3] + _p[4] * beta
 * dE/dx < _p[5] / (_p[6] * beta + _p[7]) + _p[8] + _p[9] * beta
 *
 * where p[i] are calibration parameters. Currently, they are available for each single layer; the beta used is the beta[12] corresponding
 * to the TOF standalone track or to the TRK physical track, depending on the value of the standAlone parameter of the constructor. The
 * cut parameters have been obtained using the standalone beta. Currently, no test has been made to see if these parameters are good also
 * for beta from TRK track.
 * If the required combination of planes has no available calibration, the one for S21 will be used. Currently, the 9th reduction uses
 * a dE/dx calibration which allow to use an (almost) unique set of parameters for all layers, but the possibility to define cut parameters
 * for each layer is maintained for possible future refinements.
 * Note that there's no currently available calibration for combinations of layers, so using mean release on many layers could potentially
 * result in an undefined behavior, but given the almost constant cut parameters it should work as expected (however, THIS HAS NOT BEEN TESTED
 * AND IS NOT GUARANTEED).
 */
class TofDedxHCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * @param cutName The cut's name.
   * @param layer The layer to use. Select it using values like #S11 (for single layers) or S11+S12 (for the mean of many layers).
   * @param standAlone If true, the cut will use the standalone measure of beta; otherwise it will use the
   *                   physical track.
   * @param strict If true, a more strict cut will be applied (lower efficiency but also lower contamination).
   */
  TofDedxHCut(const char *cutName, unsigned int layer, bool standAlone = true, bool strict = false);
  /*! @brief Destructor. */
  ~TofDedxHCut() {

  }

  /*! @brief The dE/dx vs. beta helium check.
   *
   * @param event The event to analyze.
   * @return #CUTOK if event is inside helium "corridor".
   * @return 0 otherwise.
   */
  int Check(PamLevel2 *event);

private:

  unsigned int _layer;
  bool _standAlone;
  const float *_p;
  // Single layers standalone
  static const float _S11_SA_params[10], _S12_SA_params[10], _S21_SA_params[10], _S22_SA_params[10], _S31_SA_params[10],
      _S32_SA_params[10];
  // Single layers standalone, strict selection
  static const float _S11_SA_params_strict[10], _S12_SA_params_strict[10], _S21_SA_params_strict[10],
      _S22_SA_params_strict[10], _S31_SA_params_strict[10], _S32_SA_params_strict[10];
};

#endif /* TOFDEDXHCUT_H_ */
1	pam-fi	1.1	/*
2			* TofDedxHCut.h
3			*
4			* Created on: 5 apr 2011
5			* Author: Nicola Mori
6			*/
7
8			/! @file TofDedxHCut.h The TofDedxHCut.h class definition file. /
9
10	pam-fi	1.2	#ifndef TOFDEDXHCUT_H_
11			#define TOFDEDXHCUT_H_
12	pam-fi	1.1
13			#include "../../PamCutBase/PamCutBase.h"
14
15			/*! @brief Cutting a band dE/dx vs. Beta for hydrogen.
16			* This cut discards events where TOF dE/dx vs. Beta are outside the hydrogen "corridor". This corridor is currently
17			* defined by:
18			*
19	pam-fi	1.2	* dE/dx > _p[0] / (_p[1] * beta + _p[2]) + _p[3] + _p[4] * beta
20			* dE/dx < _p[5] / (_p[6] * beta + _p[7]) + _p[8] + _p[9] * beta
21	pam-fi	1.1	*
22			* 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. Currently, no test has been made to see if these parameters are good also
25			* for beta from TRK track.
26			* If the required combination of planes has no available calibration, the one for S21 will be used. Currently, the 9th reduction uses
27			* a dE/dx calibration which allow to use an (almost) unique set of parameters for all layers, but the possibility to define cut parameters
28			* for each layer is maintained for possible future refinements.
29	pam-fi	1.2	* Note that there's no currently available calibration for combinations of layers, so using mean release on many layers could potentially
30			* result in an undefined behavior, but given the almost constant cut parameters it should work as expected (however, THIS HAS NOT BEEN TESTED
31			* AND IS NOT GUARANTEED).
32	pam-fi	1.1	*/
33			class TofDedxHCut: public PamCut {
34
35			public:
36			/*! @brief Constructor.
37			*
38			* @param cutName The cut's name.
39			* @param layer The layer to use. Select it using values like #S11 (for single layers) or S11+S12 (for the mean of many layers).
40			* @param standAlone If true, the cut will use the standalone measure of beta; otherwise it will use the
41			* physical track.
42	pam-fi	1.2	* @param strict If true, a more strict cut will be applied (lower efficiency but also lower contamination).
43	pam-fi	1.1	*/
44	pam-fi	1.2	TofDedxHCut(const char *cutName, unsigned int layer, bool standAlone = true, bool strict = false);
45	pam-fi	1.1	/! @brief Destructor. /
46			~TofDedxHCut() {
47
48			}
49
50			/*! @brief The dE/dx vs. beta helium check.
51			*
52			* @param event The event to analyze.
53			* @return #CUTOK if event is inside helium "corridor".
54			* @return 0 otherwise.
55			*/
56			int Check(PamLevel2 *event);
57
58			private:
59
60			unsigned int _layer;
61			bool _standAlone;
62			const float *_p;
63			// Single layers standalone
64			static const float _S11_SA_params[10], _S12_SA_params[10], _S21_SA_params[10], _S22_SA_params[10], _S31_SA_params[10],
65			_S32_SA_params[10];
66	pam-fi	1.2	// Single layers standalone, strict selection
67			static const float _S11_SA_params_strict[10], _S12_SA_params_strict[10], _S21_SA_params_strict[10],
68			_S22_SA_params_strict[10], _S31_SA_params_strict[10], _S32_SA_params_strict[10];
69	pam-fi	1.1	};
70
71	pam-fi	1.2	#endif /* TOFDEDXHCUT_H_ */