TofCuts/TofDedxHeCut/TofDedxHeCut.cpp

/*
 * TofDedxHeCut.cpp
 *
 *  Created on: 03/ago/2009
 *      Author: Nicola Mori
 */

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

#include "TofDedxHeCut.h"

int TofDedxHeCut::Check(PamLevel2 *event) {

  float beta = event->GetTrack(0)->GetToFTrack()->beta[12];
  int trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;
  float dEdx = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100)
      + event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
  dEdx /= 4;

  if (dEdx < 1 / (0.2 * (beta - 0.05)) - 2.)
    return 0;

  if (dEdx > 1 / (0.25 * (beta - 0.1)))
    return 1;

  if (dEdx < 11. * beta - 7.7)
    return 2;

  return CUTOK;
}
1	pam-fi	1.1	/*
2			* TofDedxHeCut.cpp
3			*
4			* Created on: 03/ago/2009
5			* Author: Nicola Mori
6			*/
7
8			/! @file TofDedxHeCut.cpp The TofDedxHeCut class implementation file. /
9
10			#include "TofDedxHeCut.h"
11
12			int TofDedxHeCut::Check(PamLevel2 *event) {
13
14			float beta = event->GetTrack(0)->GetToFTrack()->beta[12];
15			int trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;
16			float dEdx = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100)
17			+ event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
18			dEdx /= 4;
19
20			if (dEdx < 1 / (0.2 * (beta - 0.05)) - 2.)
21			return 0;
22
23			if (dEdx > 1 / (0.25 * (beta - 0.1)))
24			return 1;
25
26			if (dEdx < 11. * beta - 7.7)
27			return 2;
28
29			return CUTOK;
30			}