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	pam-fi	1.2	/*
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	pam-fi	1.3	// TrkTrack *track = event->GetTrack(0)->GetTrkTrack();
15			if(event->GetNTracks(_trkAlg)==0)return 0;
16			ExtTrack *track = event->GetTrack(0,_trkAlg)->GetExtTrack();
17	pam-fi	1.2	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			}