TrkCuts/TrkSigmaDeflCut/TrkSigmaDeflCut.h

/*
 * TrkSigmaDeflCut.h
 *
 *  Created on: 23-apr-2009
 *      Author: S. Ricciarini
 */

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

#ifndef TRKSIGMADEFLCUT_H_
#define TRKSIGMADEFLCUT_H_

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

/*! @brief The cut on TRK sigma_deflection.
 *
 * This cut discards an event if its sigma_deflection (defined as sqrt(coval[4][4]), in GV^-1) is not less than the minimum deflection for the bin which contains the event deflection modulus (deflection of the first TRK track).
 *
 */
class TrkSigmaDeflCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * @param cutName The cut's name.
   * @param binning A vector containing the histogram binning in rigidity.
   *                Note that it must contain both the upper and lower limits, and
   *                that the elements must be ordered (ie., lowest rigidity value in the first
   *                element and so on).
   */
  TrkSigmaDeflCut(const char *cutName, std::vector<float> binning) :
    PamCut(cutName), _binning(binning) {
  }

  /*! @brief Constructor.
   *
   * @param cutName The cut's name.
   * @param rigBinListFileName A text file containing the histogram binning.
   *                Note that it must contain both the upper and lower limits, and
   *                that the elements must be ordered (ie., lowest rigidity value in the first
   *                element and so on).
   */
  TrkSigmaDeflCut(const char *cutName, const char* rigBinListFileName);

  /*! @brief Destructor. */
  ~TrkSigmaDeflCut() {
  }

  /*! @brief The sigma_deflection check.
   *
   * The event is discarded if its deflection is less than the lower
   * histogram limit or higher than the upper limit. The event is also discarded if its sigma_deflection is not less than the minimum deflection for the bin which contains the event deflection modulus.
   *
   * @param event The event to analyze.
   * @return #CUTOK if the sigma_deflection condition is satisfied.
   * @return 0 otherwise.
   */
  int Check(PamLevel2 *event);


private:

  std::vector<float> _binning;

};

#endif /* TRKSIGMADEFLCUT_H_ */
1	/*
2	* TrkSigmaDeflCut.h
3	*
4	* Created on: 23-apr-2009
5	* Author: S. Ricciarini
6	*/
7
8	/! @file TrkSigmaDeflCut.h The TrkSigmaDeflCut class definition file /
9
10	#ifndef TRKSIGMADEFLCUT_H_
11	#define TRKSIGMADEFLCUT_H_
12
13	#include "../../PamCutBase/PamCutBase.h"
14
15	/*! @brief The cut on TRK sigma_deflection.
16	*
17	* This cut discards an event if its sigma_deflection (defined as sqrt(coval[4][4]), in GV^-1) is not less than the minimum deflection for the bin which contains the event deflection modulus (deflection of the first TRK track).
18	*
19	*/
20	class TrkSigmaDeflCut: public PamCut {
21
22	public:
23	/*! @brief Constructor.
24	*
25	* @param cutName The cut's name.
26	* @param binning A vector containing the histogram binning in rigidity.
27	* Note that it must contain both the upper and lower limits, and
28	* that the elements must be ordered (ie., lowest rigidity value in the first
29	* element and so on).
30	*/
31	TrkSigmaDeflCut(const char *cutName, std::vector<float> binning) :
32	PamCut(cutName), _binning(binning) {
33	}
34
35	/*! @brief Constructor.
36	*
37	* @param cutName The cut's name.
38	* @param rigBinListFileName A text file containing the histogram binning.
39	* Note that it must contain both the upper and lower limits, and
40	* that the elements must be ordered (ie., lowest rigidity value in the first
41	* element and so on).
42	*/
43	TrkSigmaDeflCut(const char cutName, const char rigBinListFileName);
44
45	/! @brief Destructor. /
46	~TrkSigmaDeflCut() {
47	}
48
49	/*! @brief The sigma_deflection check.
50	*
51	* The event is discarded if its deflection is less than the lower
52	* histogram limit or higher than the upper limit. The event is also discarded if its sigma_deflection is not less than the minimum deflection for the bin which contains the event deflection modulus.
53	*
54	* @param event The event to analyze.
55	* @return #CUTOK if the sigma_deflection condition is satisfied.
56	* @return 0 otherwise.
57	*/
58	int Check(PamLevel2 *event);
59
60
61	private:
62
63	std::vector<float> _binning;
64
65	};
66
67	#endif /* TRKSIGMADEFLCUT_H_ */