TrkCuts/TrkDedxNucleiCut/TrkDedxNucleiCut.cpp

/*
 * TrkDedxNucleiCut.cpp
 *
 *  Created on: 8/jul/2010
 *      Author: Nicola Mori
 */

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

#include "TrkDedxNucleiCut.h"

int TrkDedxNucleiCut::Check(PamLevel2 *event) {

//  TrkTrack *track = event->GetTrack(0)->GetTrkTrack();
    if(event->GetNTracks(_trkAlg)==0)return 0;    
    ExtTrack *track = event->GetTrack(0,_trkAlg)->GetExtTrack();
  float R = track->GetRigidity();
  float dEdx = track->GetDEDX();

  if (dEdx < 3.7 + 4.8 / (R * R))
    return 0;

  return CUTOK;
}
1	/*
2	* TrkDedxNucleiCut.cpp
3	*
4	* Created on: 8/jul/2010
5	* Author: Nicola Mori
6	*/
7
8	/! @file TrkDedxNucleiCut.cpp The TrkDedxNucleiCut class implementation file. /
9
10	#include "TrkDedxNucleiCut.h"
11
12	int TrkDedxNucleiCut::Check(PamLevel2 *event) {
13
14	// TrkTrack *track = event->GetTrack(0)->GetTrkTrack();
15	if(event->GetNTracks(_trkAlg)==0)return 0;
16	ExtTrack *track = event->GetTrack(0,_trkAlg)->GetExtTrack();
17	float R = track->GetRigidity();
18	float dEdx = track->GetDEDX();
19
20	if (dEdx < 3.7 + 4.8 / (R * R))
21	return 0;
22
23	return CUTOK;
24	}