TrkCuts/TrkChi2DeflCut/TrkChi2DeflCut.cpp

/*
 * TrkChi2DeflCut.cpp
 *
 *  Created on: 7-apr-2009
 *      Author: Nicola Mori
 */

#include "TrkChi2DeflCut.h"

int TrkChi2DeflCut::Check(PamLevel2 *event) {

  float chi2 = event->GetTrack(0)->GetTrkTrack()->chi2;
  if (chi2 < 0.)
    return 0;
  float eta2 = event->GetTrack(0)->GetTrkTrack()->GetDeflection(); //Sign doesn't matter
  eta2 *= eta2;  // Now eta2 is deflection^2

#ifdef DEBUGPAMCUT

  h_trk_chi2_defl[0]->Fill(event->GetTrack(0)->GetTrkTrack()->GetDeflection(), pow(chi2,0.25));

#endif

  if (chi2 > _p0 + _p1 * eta2 + _p2 * eta2 * eta2)
    return 0;

#ifdef DEBUGPAMCUT

  h_trk_chi2_defl[1]->Fill(event->GetTrack(0)->GetTrkTrack()->GetDeflection(), pow(chi2,0.25));

#endif

  return CUTOK;

}
1	pam-fi	1.1	/*
2			* TrkChi2DeflCut.cpp
3			*
4			* Created on: 7-apr-2009
5			* Author: Nicola Mori
6			*/
7
8			#include "TrkChi2DeflCut.h"
9
10			int TrkChi2DeflCut::Check(PamLevel2 *event) {
11
12			float chi2 = event->GetTrack(0)->GetTrkTrack()->chi2;
13			if (chi2 < 0.)
14			return 0;
15			float eta2 = event->GetTrack(0)->GetTrkTrack()->GetDeflection(); //Sign doesn't matter
16			eta2 *= eta2; // Now eta2 is deflection^2
17
18			#ifdef DEBUGPAMCUT
19
20			h_trk_chi2_defl[0]->Fill(event->GetTrack(0)->GetTrkTrack()->GetDeflection(), pow(chi2,0.25));
21
22			#endif
23
24			if (chi2 > _p0 + _p1 * eta2 + _p2 * eta2 * eta2)
25			return 0;
26
27			#ifdef DEBUGPAMCUT
28
29			h_trk_chi2_defl[1]->Fill(event->GetTrack(0)->GetTrkTrack()->GetDeflection(), pow(chi2,0.25));
30
31			#endif
32
33			return CUTOK;
34
35			}