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	/*
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	}