TofCuts/TofNucleiZCut/TofNucleiZCut.h

/*
 * TofNucleiZCut.h
 *
 *  Created on: 21-apr-2009
 *      Author: Nicola Mori
 */

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

#ifndef NO_TOFNUCLEI

#ifndef TOFNUCLEIZCUT_H_
#define TOFNUCLEIZCUT_H_

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

#ifdef DEBUGPAMCUT
#include <TH2F.h>
#endif

/*! @enum TofNucleiZ_Return Return values for rejected events. */
enum TofNucleiZ_Return {
  TOFNUCLEIZ_OUTOFBOUNDS, ///< Discarded because charge is out of bounds in a layer.
  TOFNUCLEIZ_TOOFEWLAYERS, ///< Discarded because charge can be measured in too few layers.
  TOFNUCLEIZ_ILLEGALZ
///< Discarded because charge is not between 1 and 8
};

/*! @brief The TofNuclei charge cut.
 *
 * This cut uses the TofNuclei library to cut events depending on particle's charge value.
 *
 */

class TofNucleiZCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * The charge to look for is passed to the TofNucleiZCut object as Z. For each charge, there is an associated distribution of
   * events, with mean Zmean and standard deviation sigmaZ. These values will be used to check the events.
   *
   * @param cutName The cut's name.
   * @param Z The charge value.
   * @param lowerLimit The lower bound (in units of sigmaZ) of the desired charge range, defined as Zmean - lowerBound*sigmaZ.
   * @param upperLimit The upper bound (in units of sigmaZ) of the desired charge range, defined as Zmean + upperBound*sigmaZ.
   * @param minLayers The minimum number of layers for which Z must have been measured.
   * @param layersMask The masked layers (see #TOFLAYERS). These will be ignored throughout all the TofNucleiZCut analysis.
   * @see Check
   */
  TofNucleiZCut(const char *cutName, unsigned int Z, float lowerLimit, float upperLimit, unsigned int minLayers = 6,
      unsigned int layersMask = 0);

  /*! @brief Destructor. */
  ~TofNucleiZCut() {
    delete _tofNuclei;

#ifdef DEBUGPAMCUT

    for (UInt_t j = 0; j < 12; j++) {
      h_tof_z_beta[j]->Write();
    }

#endif

  }

  /*! @brief The TofNucleiZ check.
   *
   * The routine selects events whose value of Zmean obtained from the ToF is in the range [Zmean - sigma*lowerBound, Zmean + sigma*upperBound].
   * Zmean and sigma are defined respectively as the mean and the standard deviation of the charge distributions obtained using ToFNuclei.
   * The calibration values used in this routine (for all energies) are:
   *
   *   Z=1 |   Zmean     sigma
   *   ------------------------
   *   S11 |  0.987      0.11
   *   S12 |  0.988      0.12
   *   S21 |  0.961      0.083
   *   S22 |  0.966      0.10
   *   S31 |  0.981      0.092
   *   S32 |  0.979      0.095
   *
   *   Z=2 |   Zmean     sigma
   *   ------------------------
   *   S11 |  1.96       0.12
   *   S12 |  1.95       0.14
   *   S21 |  1.95       0.13
   *   S22 |  1.96       0.14
   *   S31 |  1.99       0.14
   *   S32 |  2.00       0.15
   *
   * In current implementation, heavier nuclei up to oxygen are partially calibrated; means are taken equal to the particle's charge (Zmean = Z)
   * and a single sigma is used for all ToF layers:
   *
   *   Z |  sigma
   *   -----------
   *   3 |  0.20
   *   4 |  0.25
   *   5 |  0.32
   *   6 |  0.40
   *   7 |  0.45
   *   8 |  0.50
   *
   * The charge of the particle is detrmined by calling ToFNuclei::Get_Charge_ToF_std_layer().
   * First of all, the number of non-masked layers giving a valid charge information is checked to be greater than minLayers; then the bound
   * on the charge value is checked for all valid planes. If at least one of the valid, non-masked planes gives a charge outside the bounds
   * the routine will return #TOFNUCLEIZ_OUTOFBOUNDS. Note that this happens even if more than minLayers give a charge inside boundaries;
   * for example, if minLayers = 5 and we have 6 good charge measurement of which only 5 are inside the boundaries, the check will fail.
   * Notice that masking a layer has the net effect of reducing the number of layers: asking minLayers = 5 and masking two planes
   * will result in return value of #TOFNUCLEIZ_TOOFEWLAYERS for all events, since the "effective ToF" has only 4 layers.
   *
   * @param event The event to analyze.
   * @return #CUTOK if the charge obtained from TofNuclei is in the range [Z - sigma*lowerBound, Z + sigma*upperBound] in at least minLayers
   *          (apart from eventual masked layers).
   * @return #TOFNUCLEIZ_OUTOFBOUNDS if charge is out of bounds in at least one non-masked layer.
   * @return #TOFNUCLEIZ_TOOFEWLAYERS if charge information is available for a number of layers lesser than minLayers.
   * @return #TRKNUCLEIZ_OUTOFBOUNDS if Z < 1 or Z > 8.
   */
  int Check(PamLevel2 *event);

private:

#ifdef DEBUGPAMCUT

  TH2F* h_tof_z_beta[12];

#endif

  unsigned int _Z;
  SimpleMatrix<float> _mean; // Index structure: [Z-1][layer]
  SimpleMatrix<float> _sigma; // Index structure: [Z-1][layer]
  static const float _sigmaZdefault[8];

  float _lowerLimit;
  float _upperLimit;
  unsigned int _minLayers;
  unsigned int _layersMask;

  ToFNuclei *_tofNuclei;
  static const unsigned int _maskArray[6];

};

#endif /* TOFNUCLEIZCUT_H_ */

#endif /* NO_TOFNUCLEI */
1	/*
2	* TofNucleiZCut.h
3	*
4	* Created on: 21-apr-2009
5	* Author: Nicola Mori
6	*/
7
8	/! @file TofNucleiZCut.h The TofNucleiZCut class definition file /
9
10	#ifndef NO_TOFNUCLEI
11
12	#ifndef TOFNUCLEIZCUT_H_
13	#define TOFNUCLEIZCUT_H_
14
15	#include "../../PamCutBase/PamCutBase.h"
16	#include <ToFNuclei.h>
17
18	#ifdef DEBUGPAMCUT
19	#include <TH2F.h>
20	#endif
21
22	/! @enum TofNucleiZ_Return Return values for rejected events. /
23	enum TofNucleiZ_Return {
24	TOFNUCLEIZ_OUTOFBOUNDS, ///< Discarded because charge is out of bounds in a layer.
25	TOFNUCLEIZ_TOOFEWLAYERS, ///< Discarded because charge can be measured in too few layers.
26	TOFNUCLEIZ_ILLEGALZ
27	///< Discarded because charge is not between 1 and 8
28	};
29
30	/*! @brief The TofNuclei charge cut.
31	*
32	* This cut uses the TofNuclei library to cut events depending on particle's charge value.
33	*
34	*/
35
36	class TofNucleiZCut: public PamCut {
37
38	public:
39	/*! @brief Constructor.
40	*
41	* The charge to look for is passed to the TofNucleiZCut object as Z. For each charge, there is an associated distribution of
42	* events, with mean Zmean and standard deviation sigmaZ. These values will be used to check the events.
43	*
44	* @param cutName The cut's name.
45	* @param Z The charge value.
46	* @param lowerLimit The lower bound (in units of sigmaZ) of the desired charge range, defined as Zmean - lowerBound*sigmaZ.
47	* @param upperLimit The upper bound (in units of sigmaZ) of the desired charge range, defined as Zmean + upperBound*sigmaZ.
48	* @param minLayers The minimum number of layers for which Z must have been measured.
49	* @param layersMask The masked layers (see #TOFLAYERS). These will be ignored throughout all the TofNucleiZCut analysis.
50	* @see Check
51	*/
52	TofNucleiZCut(const char *cutName, unsigned int Z, float lowerLimit, float upperLimit, unsigned int minLayers = 6,
53	unsigned int layersMask = 0);
54
55	/! @brief Destructor. /
56	~TofNucleiZCut() {
57	delete _tofNuclei;
58
59	#ifdef DEBUGPAMCUT
60
61	for (UInt_t j = 0; j < 12; j++) {
62	h_tof_z_beta[j]->Write();
63	}
64
65	#endif
66
67	}
68
69	/*! @brief The TofNucleiZ check.
70	*
71	* The routine selects events whose value of Zmean obtained from the ToF is in the range [Zmean - sigmalowerBound, Zmean + sigmaupperBound].
72	* Zmean and sigma are defined respectively as the mean and the standard deviation of the charge distributions obtained using ToFNuclei.
73	* The calibration values used in this routine (for all energies) are:
74	*
75	* Z=1 \| Zmean sigma
76	* ------------------------
77	* S11 \| 0.987 0.11
78	* S12 \| 0.988 0.12
79	* S21 \| 0.961 0.083
80	* S22 \| 0.966 0.10
81	* S31 \| 0.981 0.092
82	* S32 \| 0.979 0.095
83	*
84	* Z=2 \| Zmean sigma
85	* ------------------------
86	* S11 \| 1.96 0.12
87	* S12 \| 1.95 0.14
88	* S21 \| 1.95 0.13
89	* S22 \| 1.96 0.14
90	* S31 \| 1.99 0.14
91	* S32 \| 2.00 0.15
92	*
93	* In current implementation, heavier nuclei up to oxygen are partially calibrated; means are taken equal to the particle's charge (Zmean = Z)
94	* and a single sigma is used for all ToF layers:
95	*
96	* Z \| sigma
97	* -----------
98	* 3 \| 0.20
99	* 4 \| 0.25
100	* 5 \| 0.32
101	* 6 \| 0.40
102	* 7 \| 0.45
103	* 8 \| 0.50
104	*
105	* The charge of the particle is detrmined by calling ToFNuclei::Get_Charge_ToF_std_layer().
106	* First of all, the number of non-masked layers giving a valid charge information is checked to be greater than minLayers; then the bound
107	* on the charge value is checked for all valid planes. If at least one of the valid, non-masked planes gives a charge outside the bounds
108	* the routine will return #TOFNUCLEIZ_OUTOFBOUNDS. Note that this happens even if more than minLayers give a charge inside boundaries;
109	* for example, if minLayers = 5 and we have 6 good charge measurement of which only 5 are inside the boundaries, the check will fail.
110	* Notice that masking a layer has the net effect of reducing the number of layers: asking minLayers = 5 and masking two planes
111	* will result in return value of #TOFNUCLEIZ_TOOFEWLAYERS for all events, since the "effective ToF" has only 4 layers.
112	*
113	* @param event The event to analyze.
114	* @return #CUTOK if the charge obtained from TofNuclei is in the range [Z - sigmalowerBound, Z + sigmaupperBound] in at least minLayers
115	* (apart from eventual masked layers).
116	* @return #TOFNUCLEIZ_OUTOFBOUNDS if charge is out of bounds in at least one non-masked layer.
117	* @return #TOFNUCLEIZ_TOOFEWLAYERS if charge information is available for a number of layers lesser than minLayers.
118	* @return #TRKNUCLEIZ_OUTOFBOUNDS if Z < 1 or Z > 8.
119	*/
120	int Check(PamLevel2 *event);
121
122	private:
123
124	#ifdef DEBUGPAMCUT
125
126	TH2F* h_tof_z_beta[12];
127
128	#endif
129
130	unsigned int _Z;
131	SimpleMatrix<float> _mean; // Index structure: [Z-1][layer]
132	SimpleMatrix<float> _sigma; // Index structure: [Z-1][layer]
133	static const float _sigmaZdefault[8];
134
135	float _lowerLimit;
136	float _upperLimit;
137	unsigned int _minLayers;
138	unsigned int _layersMask;
139
140	ToFNuclei *_tofNuclei;
141	static const unsigned int _maskArray[6];
142
143	};
144
145	#endif /* TOFNUCLEIZCUT_H_ */
146
147	#endif /* NO_TOFNUCLEI */