TrkCuts/TrkChi2DeflTimeCut/TrkChi2DeflTimeCut.cpp

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

#include "TrkChi2DeflTimeCut.h"

double TrkChi2DeflTimeCut::_GetChi2(PamLevel2 *event) {
  return event->GetTrack(0)->GetTrkTrack()->chi2;
}

int TrkChi2DeflTimeCut::Check(PamLevel2 *event) {

  double chi2 = _GetChi2(event);
  if (chi2 < 0.)
    return 0;
  Double_t etaMod = fabs(event->GetTrack(0)->GetTrkTrack()->GetDeflection());

  _time.Set(event->GetOrbitalInfo()->absTime, kFALSE, 0, kFALSE);
  // TTimestamp::GetDate() returns date in format YYYYMMDD so to compare it
  // with 6-digits dates (YYMMDD) we must subtract 20000000

  Int_t iDaySel = (_time.GetDate(kFALSE) - 20000000); // note: GetDate returns an unsigned integer

  bool found = false;
  for (UInt_t i = 0; i < _iDayFirst.size(); i++) {
    if (_iDayFirst[i] <= iDaySel && iDaySel <= _iDayLast[i]) {
      for (int iPar=0; iPar<_nPar; iPar++) {
        _pSel[iPar] = _p[iPar][i];
      }
      
//      cout << "FOUND: day " << iDaySel << endl;
//      cout << _iDayFirst[i] << " " << _iDayLast[i] << endl;
//      cout << _pSel[0] << " " << _pSel[1] << " " << _pSel[2] << " " << _pSel[3] << endl;

      found = true;
      break;
    }
  }

  if (!found)
    return 0;

  double chi2max;
  
  if (_nPar==4) {
  
    chi2max=_pSel[0] + _pSel[1] * pow(etaMod, _pSel[2]) * (1. + pow(_pSel[3] * etaMod, 2));

//    cout << "chi2max " << chi2max << endl;

  }
  else if (_nPar==5) {
    chi2max=_pSel[0] + _pSel[1] * pow(etaMod, _pSel[2]) * (1. + pow(_pSel[3] * etaMod, _pSel[4]));
  }
  else {
    chi2max=-1;
  }

//  cout << "chi2 " << chi2 << endl;
  
  if (chi2 > chi2max) {
  
//    cout << "KO" << endl;
    
    return 0;

  }

//  cout << "OK" << endl;

  return CUTOK;

}
1	/*
2	* TrkChi2DeflTimeCut.cpp
3	*
4	* Created on: 6-oct-2009
5	* Author: S. Ricciarini
6	*/
7
8	#include "TrkChi2DeflTimeCut.h"
9
10	double TrkChi2DeflTimeCut::_GetChi2(PamLevel2 *event) {
11	return event->GetTrack(0)->GetTrkTrack()->chi2;
12	}
13
14	int TrkChi2DeflTimeCut::Check(PamLevel2 *event) {
15
16	double chi2 = _GetChi2(event);
17	if (chi2 < 0.)
18	return 0;
19	Double_t etaMod = fabs(event->GetTrack(0)->GetTrkTrack()->GetDeflection());
20
21	_time.Set(event->GetOrbitalInfo()->absTime, kFALSE, 0, kFALSE);
22	// TTimestamp::GetDate() returns date in format YYYYMMDD so to compare it
23	// with 6-digits dates (YYMMDD) we must subtract 20000000
24
25	Int_t iDaySel = (_time.GetDate(kFALSE) - 20000000); // note: GetDate returns an unsigned integer
26
27	bool found = false;
28	for (UInt_t i = 0; i < _iDayFirst.size(); i++) {
29	if (_iDayFirst[i] <= iDaySel && iDaySel <= _iDayLast[i]) {
30	for (int iPar=0; iPar<_nPar; iPar++) {
31	_pSel[iPar] = _p[iPar][i];
32	}
33
34	// cout << "FOUND: day " << iDaySel << endl;
35	// cout << _iDayFirst[i] << " " << _iDayLast[i] << endl;
36	// cout << _pSel[0] << " " << _pSel[1] << " " << _pSel[2] << " " << _pSel[3] << endl;
37
38	found = true;
39	break;
40	}
41	}
42
43	if (!found)
44	return 0;
45
46	double chi2max;
47
48	if (_nPar==4) {
49
50	chi2max=_pSel[0] + _pSel[1] * pow(etaMod, _pSel[2]) * (1. + pow(_pSel[3] * etaMod, 2));
51
52	// cout << "chi2max " << chi2max << endl;
53
54	}
55	else if (_nPar==5) {
56	chi2max=_pSel[0] + _pSel[1] * pow(etaMod, _pSel[2]) * (1. + pow(_pSel[3] * etaMod, _pSel[4]));
57	}
58	else {
59	chi2max=-1;
60	}
61
62	// cout << "chi2 " << chi2 << endl;
63
64	if (chi2 > chi2max) {
65
66	// cout << "KO" << endl;
67
68	return 0;
69
70	}
71
72	// cout << "OK" << endl;
73
74	return CUTOK;
75
76	}