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) {

  int trkSeqNo;
  if(_standAlone)
    trkSeqNo = 0;
  else
    trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;

  float beta = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12];

  //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.06)) - 2.)
    return 0;

  if (dEdx > 1 / (0.25 * (beta - 0.12)))
    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	int trkSeqNo;
15	if(_standAlone)
16	trkSeqNo = 0;
17	else
18	trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;
19
20	float beta = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12];
21
22	//float beta = event->GetTrack(0)->GetToFTrack()->beta[12];
23	//int trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;
24	float dEdx = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100)
25	+ event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
26	dEdx /= 4;
27
28	if (dEdx < 1 / (0.2 * (beta - 0.06)) - 2.)
29	return 0;
30
31	if (dEdx > 1 / (0.25 * (beta - 0.12)))
32	return 1;
33
34	if (dEdx < 11. * beta - 7.7)
35	return 2;
36
37	return CUTOK;
38	}