TrkCuts/TrkChi2QualCut/TrkChi2QualCut.cpp

/*
 * TrkChi2QualCut.cpp
 *
 *  Created on: 13-mar-2009
 *      Author: Sergio Ricciarini, Nicola Mori
 */

/*! @file TrkChi2QualCut.cpp The TrkChi2QualCut class implementation file */

#include "TrkChi2QualCut.h"

int TrkChi2QualCut::Check(PamLevel2 *event) {

  TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();

  float p0 = 1.111588e+00;
  float p1 = 1.707656e+00;
  float p2 = 1.489693e-01;
  float def_0 = 0.07;
  float chi2m025_0 = p0 + fabs(def_0) * p1 + def_0 * def_0 * p2;

  float chi2m025 = p0 + fabs(trkTrack->GetDeflection()) * p1 + pow(trkTrack->GetDeflection(), 2) * p2;

  float chi2m = pow(chi2m025 - chi2m025_0 + pow(6.4, 0.25), 4.); // cut at ~90% of chi^2 distribution

  if (!(trkTrack->chi2 < chi2m))
    return 0;
  return CUTOK;

}
1	pam-fi	1.1	/*
2			* TrkChi2QualCut.cpp
3			*
4			* Created on: 13-mar-2009
5			* Author: Sergio Ricciarini, Nicola Mori
6			*/
7
8			/! @file TrkChi2QualCut.cpp The TrkChi2QualCut class implementation file /
9
10			#include "TrkChi2QualCut.h"
11
12			int TrkChi2QualCut::Check(PamLevel2 *event) {
13
14			TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();
15
16			float p0 = 1.111588e+00;
17			float p1 = 1.707656e+00;
18			float p2 = 1.489693e-01;
19			float def_0 = 0.07;
20			float chi2m025_0 = p0 + fabs(def_0) * p1 + def_0 * def_0 * p2;
21
22			float chi2m025 = p0 + fabs(trkTrack->GetDeflection()) * p1 + pow(trkTrack->GetDeflection(), 2) * p2;
23
24			float chi2m = pow(chi2m025 - chi2m025_0 + pow(6.4, 0.25), 4.); // cut at ~90% of chi^2 distribution
25
26			if (!(trkTrack->chi2 < chi2m))
27			return 0;
28			return CUTOK;
29
30			}