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)
 * multiplied by a constant parameter 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).
   * @param sigmaDeflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
   */
  TrkSigmaDeflCut(const char *cutName, std::vector<float> binning, float sigmaDeflFactor = 1.) :
    PamCut(cutName), _binning(binning), _sigmaDeflFactor(sigmaDeflFactor) {
  }

  /*! @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).
   * @param sigmaDeflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
   */
  TrkSigmaDeflCut(const char *cutName, const char* rigBinListFileName, float deflFactor = 1.);

  /*! @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;
  float _sigmaDeflFactor;

};

#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)
18	* multiplied by a constant parameter is not less than the minimum deflection for the bin which
19	* contains the event deflection modulus (deflection of the first TRK track).
20	*
21	*/
22	class TrkSigmaDeflCut: public PamCut {
23
24	public:
25	/*! @brief Constructor.
26	*
27	* @param cutName The cut's name.
28	* @param binning A vector containing the histogram binning in rigidity.
29	* Note that it must contain both the upper and lower limits, and
30	* that the elements must be ordered (ie., lowest rigidity value in the first
31	* element and so on).
32	* @param sigmaDeflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
33	*/
34	TrkSigmaDeflCut(const char *cutName, std::vector<float> binning, float sigmaDeflFactor = 1.) :
35	PamCut(cutName), _binning(binning), _sigmaDeflFactor(sigmaDeflFactor) {
36	}
37
38	/*! @brief Constructor.
39	*
40	* @param cutName The cut's name.
41	* @param rigBinListFileName A text file containing the histogram binning.
42	* Note that it must contain both the upper and lower limits, and
43	* that the elements must be ordered (ie., lowest rigidity value in the first
44	* element and so on).
45	* @param sigmaDeflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
46	*/
47	TrkSigmaDeflCut(const char cutName, const char rigBinListFileName, float deflFactor = 1.);
48
49	/! @brief Destructor. /
50	~TrkSigmaDeflCut() {
51	}
52
53	/*! @brief The sigma_deflection check.
54	*
55	* The event is discarded if its deflection is less than the lower
56	* 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.
57	*
58	* @param event The event to analyze.
59	* @return #CUTOK if the sigma_deflection condition is satisfied.
60	* @return 0 otherwise.
61	*/
62	int Check(PamLevel2 *event);
63
64	private:
65
66	std::vector<float> _binning;
67	float _sigmaDeflFactor;
68
69	};
70
71	#endif /* TRKSIGMADEFLCUT_H_ */