CaloCuts/CaloGeomCut/CaloGeomCut.h

/*
 * CaloGeomCut.h
 *
 *  Created on: 19-mar-2009
 *      Author: Sergio Ricciarini
 */

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

#ifndef NO_CALOAXIS

#ifndef CALOGEOMCUT_H_
#define CALOGEOMCUT_H_

#include "../../PamCutBase/PamCutBase.h"
#include <CaloAxis.h>

/*! @brief The geometric cut using the calorimeter track.
 *
 * This cut checks if the track obtained from the calorimeter (with different methods) is inside the fiducial acceptance.
 * The current implementation uses the CaloAxis objects; to save computing time, the class
 * assumes that the track is externally computed for each event and stored in
 * CaloAxis objects; pointers to these objects are passed as arguments to the constructor.
 * The Check method will then ignore the PamLevel2 *event and assume that the
 * current content of the CaloAxis objects are relative to the current event.
 * It is an user's task to ensure that these assumptions are fulfilled every time
 * Check or ApplyCut are called.
 *
 * CUT DEPENDENCIES: CaloTrackCut for the existence of the track (it can also provide the CaloAxis objects).
 */

class CaloGeomCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * The CaloAxis arguments are pointers to objects which contain the calorimeter
   * track information for current event.
   *
   * @param cutName The cut's name.
   * @param xCaloAxis The pointer to the CaloAxis object for X axis.
   * @param yCaloAxis The pointer to the CaloAxis object for Y axis.
   * @param iMethod index of the method used to track: 0 means straight line;
   *                1 means curved line (taking into account magnetic field and beta
   *                from TOF stand-alone with 3 parameters specified
   *                below) starting from the top of the calorimeter and going backward
   *                (straight line inside the calorimeter).
   *                Various combinations of TOF beta quality parameters are possible:
   *                default, low-quality beta (used for beta[12]): 10.,10.,20.;
   *                medium-quality beta: 5.,15.,4; high-quality beta: 3.,20.,3.
   * @param xTolCaloTrack The tolerance for X view which defines the fiducial acceptance.
   * @param yTolCaloTrack The tolerance for Y view which defines the fiducial acceptance.
   * @param mass The particle's mass (used only when iMethod = 1). Default is #H_MASS.
   * @param resMax Default: 10. (equivalent to beta[12])
   * @param qualCut Default: 10. (equivalent to beta[12])
   * @param chi2Cut Default: 20. (equivalent to beta[12])
   */
  CaloGeomCut(const char *cutName, CaloAxis *xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
      Float_t xTolCaloTrack = 0.7, Float_t yTolCaloTrack = 0.7, Float_t mass = H_MASS, Float_t resMax = 10.,
      Float_t qualCut = 10., Float_t chi2Cut = 20.);

  /*! @brief Destructor. */
  ~CaloGeomCut();

  /*! @brief The geometry check using the calorimeter's track.
   *
   * @param event The event to analyze.
   * @return #CUTOK if charge released in plane 22 is greater than 0 (from CaloAxis) for both X and Y
   * @return 0 otherwise
   */
  int Check(PamLevel2 *event);

private:
  CaloAxis *_xCaloAxis, *_yCaloAxis;
  UInt_t _iMethod;
  Float_t _xTolCaloTrack, _yTolCaloTrack;
  Float_t _mass;
  Float_t _resMax, _qualCut, _chi2Cut;

  Trajectory* _traj;
  static const Int_t _nPoint;
  vector<Float_t> _zIn; // MUST be Float_t
  static const Float_t _zCaloTop; // cm MUST be Float_t
};

#endif /* CALOGEOMCUT_H_ */
#endif /* NO_CALOAXIS */
1	pam-fi	1.1	/*
2			* CaloGeomCut.h
3			*
4			* Created on: 19-mar-2009
5	pam-fi	1.2	* Author: Sergio Ricciarini
6	pam-fi	1.1	*/
7
8			/! @file CaloGeomCut.h The CaloGeomCut class definition file /
9
10	pam-fi	1.4	#ifndef NO_CALOAXIS
11
12	pam-fi	1.1	#ifndef CALOGEOMCUT_H_
13			#define CALOGEOMCUT_H_
14
15			#include "../../PamCutBase/PamCutBase.h"
16	pam-fi	1.4	#include <CaloAxis.h>
17	pam-fi	1.1
18	pam-fi	1.2	/*! @brief The geometric cut using the calorimeter track.
19	pam-fi	1.1	*
20	pam-fi	1.2	* This cut checks if the track obtained from the calorimeter (with different methods) is inside the fiducial acceptance.
21	pam-fi	1.1	* The current implementation uses the CaloAxis objects; to save computing time, the class
22			* assumes that the track is externally computed for each event and stored in
23			* CaloAxis objects; pointers to these objects are passed as arguments to the constructor.
24			* The Check method will then ignore the PamLevel2 *event and assume that the
25			* current content of the CaloAxis objects are relative to the current event.
26			* It is an user's task to ensure that these assumptions are fulfilled every time
27			* Check or ApplyCut are called.
28			*
29			* CUT DEPENDENCIES: CaloTrackCut for the existence of the track (it can also provide the CaloAxis objects).
30			*/
31
32			class CaloGeomCut: public PamCut {
33
34			public:
35			/*! @brief Constructor.
36			*
37			* The CaloAxis arguments are pointers to objects which contain the calorimeter
38			* track information for current event.
39			*
40			* @param cutName The cut's name.
41			* @param xCaloAxis The pointer to the CaloAxis object for X axis.
42			* @param yCaloAxis The pointer to the CaloAxis object for Y axis.
43	pam-fi	1.2	* @param iMethod index of the method used to track: 0 means straight line;
44			* 1 means curved line (taking into account magnetic field and beta
45	pam-fi	1.3	* from TOF stand-alone with 3 parameters specified
46	pam-fi	1.2	* below) starting from the top of the calorimeter and going backward
47			* (straight line inside the calorimeter).
48	pam-fi	1.3	* Various combinations of TOF beta quality parameters are possible:
49			* default, low-quality beta (used for beta[12]): 10.,10.,20.;
50			* medium-quality beta: 5.,15.,4; high-quality beta: 3.,20.,3.
51	pam-fi	1.1	* @param xTolCaloTrack The tolerance for X view which defines the fiducial acceptance.
52			* @param yTolCaloTrack The tolerance for Y view which defines the fiducial acceptance.
53	pam-fi	1.2	* @param mass The particle's mass (used only when iMethod = 1). Default is #H_MASS.
54			* @param resMax Default: 10. (equivalent to beta[12])
55			* @param qualCut Default: 10. (equivalent to beta[12])
56			* @param chi2Cut Default: 20. (equivalent to beta[12])
57	pam-fi	1.1	*/
58	pam-fi	1.5	CaloGeomCut(const char cutName, CaloAxis xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
59			Float_t xTolCaloTrack = 0.7, Float_t yTolCaloTrack = 0.7, Float_t mass = H_MASS, Float_t resMax = 10.,
60			Float_t qualCut = 10., Float_t chi2Cut = 20.);
61
62	pam-fi	1.1	/! @brief Destructor. /
63	pam-fi	1.5	~CaloGeomCut();
64	pam-fi	1.1
65			/*! @brief The geometry check using the calorimeter's track.
66			*
67			* @param event The event to analyze.
68			* @return #CUTOK if charge released in plane 22 is greater than 0 (from CaloAxis) for both X and Y
69			* @return 0 otherwise
70			*/
71			int Check(PamLevel2 *event);
72
73			private:
74			CaloAxis _xCaloAxis, _yCaloAxis;
75	pam-fi	1.2	UInt_t _iMethod;
76			Float_t _xTolCaloTrack, _yTolCaloTrack;
77			Float_t _mass;
78			Float_t _resMax, _qualCut, _chi2Cut;
79	pam-fi	1.5
80			Trajectory* _traj;
81			static const Int_t _nPoint;
82			vector<Float_t> _zIn; // MUST be Float_t
83			static const Float_t _zCaloTop; // cm MUST be Float_t
84	pam-fi	1.1	};
85
86			#endif /* CALOGEOMCUT_H_ */
87	pam-fi	1.4	#endif /* NO_CALOAXIS */