/[PAMELA software]/PamCut/TofCuts/TofDedxHeCut/TofDedxHeCut.cpp
ViewVC logotype

Diff of /PamCut/TofCuts/TofDedxHeCut/TofDedxHeCut.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pam-fi, Wed Aug 5 14:03:21 2009 UTC revision 1.2 by pam-fi, Fri Aug 14 10:30:49 2009 UTC
# Line 11  Line 11 
11    
12  int TofDedxHeCut::Check(PamLevel2 *event) {  int TofDedxHeCut::Check(PamLevel2 *event) {
13    
14    float beta = event->GetTrack(0)->GetToFTrack()->beta[12];    int trkSeqNo;
15    int trkSeqNo = event->GetTrack(0)->GetToFTrack()->trkseqno;    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)    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);        + event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100) + event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
26    dEdx /= 4;    dEdx /= 4;
27    
28    if (dEdx < 1 / (0.2 * (beta - 0.05)) - 2.)    if (dEdx < 1 / (0.2 * (beta - 0.06)) - 2.)
29      return 0;      return 0;
30    
31    if (dEdx > 1 / (0.25 * (beta - 0.1)))    if (dEdx > 1 / (0.25 * (beta - 0.12)))
32      return 1;      return 1;
33    
34    if (dEdx < 11. * beta - 7.7)    if (dEdx < 11. * beta - 7.7)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23