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 deflFactor 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	pam-fi	1.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	pam-fi	1.2	* 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	pam-fi	1.1	*
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	pam-fi	1.3	* @param sigmaDeflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
33	pam-fi	1.1	*/
34	pam-fi	1.2	TrkSigmaDeflCut(const char *cutName, std::vector<float> binning, float sigmaDeflFactor = 1.) :
35			PamCut(cutName), _binning(binning), _sigmaDeflFactor(sigmaDeflFactor) {
36	pam-fi	1.1	}
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	pam-fi	1.3	* @param deflFactor The factor to multiply sigma_deflection by before comparing it to the deflection.
46	pam-fi	1.1	*/
47	pam-fi	1.2	TrkSigmaDeflCut(const char cutName, const char rigBinListFileName, float deflFactor = 1.);
48	pam-fi	1.1
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	pam-fi	1.2	float _sigmaDeflFactor;
68	pam-fi	1.1
69			};
70
71			#endif /* TRKSIGMADEFLCUT_H_ */