TrkCuts/TrkChi2DeflTimeCut/TrkChi2DeflTimeCut.h

/*
 * TrkChi2DeflTimeCut.h
 *
 *  Created on: 6-oct-2009
 *      Author: S. Ricciarini
 */

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

#ifndef TRKCHI2DEFLTIMECUT_H_
#define TRKCHI2DEFLTIMECUT_H_

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

#include <TH2F.h>
#include <TTimeStamp.h>

/*! @brief The tracker chi2 vs |deflection| cut.
 *
 * Events whose track has been fitted with a chi2 greater than Chi2(|eta|) are discarded.
 * Chi2(|eta|) is a function which gives the value of chi2 distribution, corresponding to a given quantile (depending on the calibration text file), as a function of |deflection| (|eta|); in current implementation it is parameterized as:
 *
 *   4 parameters p0...p3:
 *
 *   Chi2(|eta|) = p0 + ( p1 * |eta|^p2 * (1 + (p3 * |eta|)^2) )
 *   with eta expressed in GV^-1
 *
 *   5 parameters p0..p4
 *
 *   Chi2(|eta|) = p0 + ( p1 * |eta|^p2 * (1 + (p3 * |eta|)^p4) )
 *   with eta expressed in GV^-1
 *
 * The parameters are read from a calibration text file (prepared for a given quantile, e.g. 95%), whose format must be:
 * ...
 * first_day    last_day        p0      p1      p2      ...
 * with first_day and last_day expressed in the format YYMMDD
 * ...
 * where each row corresponds to a given time interval (first_day,last_day)
 *
 * For a given event the associated day is determined; the FIRST (starting from first row) time interval of the calibration file, which contains the event day, gives the parameters to be inserted in the Chi2(|eta|) for the given event.
 * If the event day does not fall in any of the intervals of the calibration file, then the event is discarded. 
 *
 */

class TrkChi2DeflTimeCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * @param cutName The cut's name.
   * @param calibFile Path for the calibration text file.
   * @param nPar Number of parameters used for the Chi2 function. Maximum allowed: 5. Default: 4.
   */
  TrkChi2DeflTimeCut(const char *cutName, const char *calibFile, int nPar=4) :
    PamCut(cutName), _calibFile(calibFile), _nPar(nPar) {
 
    _chi2CutTable.open(_calibFile);
    while (! _chi2CutTable.eof()) {
    
      _chi2CutTable >> _tstringtemp;
      _iDayFirst.push_back(_tstringtemp.Atoi());

      _chi2CutTable >> _tstringtemp;
      _iDayLast.push_back(_tstringtemp.Atoi());
      
      if (_nPar>=0 && _nPar<=5) {
        for (int iPar=0; iPar<_nPar; iPar++) {
      
          _chi2CutTable >> _tstringtemp;
          _p[iPar].push_back(_tstringtemp.Atof());
        
        }
      }
            
    }
    _chi2CutTable.close();

  }

  /*! @brief Destructor. */
  ~TrkChi2DeflTimeCut() {

  }

  /*! @brief The tracker chi2 vs |deflection| check.
   *
   *
   * @param event The event to analyze.
   * @return #CUTOK if chi2 < Chi2(|eta|).
   * @return 0 if not
   */
  int Check(PamLevel2 *event);

private:

    virtual double _GetChi2(PamLevel2 *event);

    const char* _calibFile;
    int _nPar;

    TString _tstringtemp;
    
    std::vector<Int_t> _iDayFirst;
    std::vector<Int_t> _iDayLast;
    std::vector<Double_t> _p[5];

    ifstream _chi2CutTable;
    
    Double_t _pSel[5];
    
    TTimeStamp _time;

};
#endif /* TRKCHI2DEFLTIMECUT_H_ */
1	/*
2	* TrkChi2DeflTimeCut.h
3	*
4	* Created on: 6-oct-2009
5	* Author: S. Ricciarini
6	*/
7
8	/! @file TrkChi2DeflTimeCut.h The TrkChi2DeflTimeCut class definition file /
9
10	#ifndef TRKCHI2DEFLTIMECUT_H_
11	#define TRKCHI2DEFLTIMECUT_H_
12
13	#include "../../PamCutBase/PamCutBase.h"
14
15	#include <TH2F.h>
16	#include <TTimeStamp.h>
17
18	/*! @brief The tracker chi2 vs \|deflection\| cut.
19	*
20	* Events whose track has been fitted with a chi2 greater than Chi2(\|eta\|) are discarded.
21	* Chi2(\|eta\|) is a function which gives the value of chi2 distribution, corresponding to a given quantile (depending on the calibration text file), as a function of \|deflection\| (\|eta\|); in current implementation it is parameterized as:
22	*
23	* 4 parameters p0...p3:
24	*
25	* Chi2(\|eta\|) = p0 + ( p1 * \|eta\|^p2 * (1 + (p3 * \|eta\|)^2) )
26	* with eta expressed in GV^-1
27	*
28	* 5 parameters p0..p4
29	*
30	* Chi2(\|eta\|) = p0 + ( p1 * \|eta\|^p2 * (1 + (p3 * \|eta\|)^p4) )
31	* with eta expressed in GV^-1
32	*
33	* The parameters are read from a calibration text file (prepared for a given quantile, e.g. 95%), whose format must be:
34	* ...
35	* first_day last_day p0 p1 p2 ...
36	* with first_day and last_day expressed in the format YYMMDD
37	* ...
38	* where each row corresponds to a given time interval (first_day,last_day)
39	*
40	* For a given event the associated day is determined; the FIRST (starting from first row) time interval of the calibration file, which contains the event day, gives the parameters to be inserted in the Chi2(\|eta\|) for the given event.
41	* If the event day does not fall in any of the intervals of the calibration file, then the event is discarded.
42	*
43	*/
44
45	class TrkChi2DeflTimeCut: public PamCut {
46
47	public:
48	/*! @brief Constructor.
49	*
50	* @param cutName The cut's name.
51	* @param calibFile Path for the calibration text file.
52	* @param nPar Number of parameters used for the Chi2 function. Maximum allowed: 5. Default: 4.
53	*/
54	TrkChi2DeflTimeCut(const char cutName, const char calibFile, int nPar=4) :
55	PamCut(cutName), _calibFile(calibFile), _nPar(nPar) {
56
57	_chi2CutTable.open(_calibFile);
58	while (! _chi2CutTable.eof()) {
59
60	_chi2CutTable >> _tstringtemp;
61	_iDayFirst.push_back(_tstringtemp.Atoi());
62
63	_chi2CutTable >> _tstringtemp;
64	_iDayLast.push_back(_tstringtemp.Atoi());
65
66	if (_nPar>=0 && _nPar<=5) {
67	for (int iPar=0; iPar<_nPar; iPar++) {
68
69	_chi2CutTable >> _tstringtemp;
70	_p[iPar].push_back(_tstringtemp.Atof());
71
72	}
73	}
74
75	}
76	_chi2CutTable.close();
77
78	}
79
80	/! @brief Destructor. /
81	~TrkChi2DeflTimeCut() {
82
83	}
84
85	/*! @brief The tracker chi2 vs \|deflection\| check.
86	*
87	*
88	* @param event The event to analyze.
89	* @return #CUTOK if chi2 < Chi2(\|eta\|).
90	* @return 0 if not
91	*/
92	int Check(PamLevel2 *event);
93
94	private:
95
96	virtual double _GetChi2(PamLevel2 *event);
97
98	const char* _calibFile;
99	int _nPar;
100
101	TString _tstringtemp;
102
103	std::vector<Int_t> _iDayFirst;
104	std::vector<Int_t> _iDayLast;
105	std::vector<Double_t> _p[5];
106
107	ifstream _chi2CutTable;
108
109	Double_t _pSel[5];
110
111	TTimeStamp _time;
112
113	};
114	#endif /* TRKCHI2DEFLTIMECUT_H_ */