TofCuts/TofDedxS1Cut/TofDedxS1Cut.cpp

/*
 * TofDedxS1Cut.cpp
 *
 *  Created on: 13-mar-2009
 *      Author: Sergio Ricciarini, Nicola Mori
 */

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

#include "TofDedxS1Cut.h"

int TofDedxS1Cut::Check(PamLevel2 *event) {

  int iTrackTof=0; // index initialized to 0 (TOF stand-alone track)
  if (! (_trackTof==0) ) {
    iTrackTof=*_trackTof;
  }

  // for each layer of S1: check that there is one and only one hit paddle and identify it; if this condition is satisfied, do the check on DeDx for the paddle (otherwise no check is done)

  for (int ipl=0; ipl<2; ipl++) {
  
    int nhitpad = 0; // number of hit paddles
    int ihitpad = 100; // identifier of hit paddle; maximum valid value is 24 (taking into account each possible TOF layer)
    float dedxpad;

    for (int ipad=0; ipad<TOFNPADLAYER[ipl]; ipad++) {
      if (event->GetToFLevel2()->HitPaddle(ipl,ipad)) { // HitPaddle is true if and only if there are both TDC signal for that pad AND they are both 'good'
        ihitpad=event->GetToFLevel2()->GetPaddleid(ipl,ipad);
        nhitpad++;
      }
    }

    if (nhitpad == 1) { // there is one and only one hit pad in the layer
      int satwarnpad;
      event->GetToFLevel2()->GetdEdxPaddle(iTrackTof,ihitpad,100,dedxpad,satwarnpad);
      // NOTE: if dedxpad=0 then there is no dE/dL signal (i.e. no ADC signals to evaluate dE/dL): this means that ADC signals are under threshold and therefore event is good
      if ( ! (dedxpad < _maxDedxS1) ) {
//        cout << "return on check " << ipl << " hitpad " << ihitpad << " track id " << *_trackTof << " dedxpad " << dedxpad << endl;
        return ipl;
      }
    }
          
  }

  return CUTOK;
  
}
1	pam-fi	1.1	/*
2			* TofDedxS1Cut.cpp
3			*
4			* Created on: 13-mar-2009
5			* Author: Sergio Ricciarini, Nicola Mori
6			*/
7
8			/! @file TofDedxS1Cut.cpp The TofDedxS1Cut class implementation file /
9
10			#include "TofDedxS1Cut.h"
11
12			int TofDedxS1Cut::Check(PamLevel2 *event) {
13
14			int iTrackTof=0; // index initialized to 0 (TOF stand-alone track)
15			if (! (_trackTof==0) ) {
16			iTrackTof=*_trackTof;
17			}
18
19			// for each layer of S1: check that there is one and only one hit paddle and identify it; if this condition is satisfied, do the check on DeDx for the paddle (otherwise no check is done)
20
21			for (int ipl=0; ipl<2; ipl++) {
22
23			int nhitpad = 0; // number of hit paddles
24			int ihitpad = 100; // identifier of hit paddle; maximum valid value is 24 (taking into account each possible TOF layer)
25			float dedxpad;
26
27			for (int ipad=0; ipad<TOFNPADLAYER[ipl]; ipad++) {
28			if (event->GetToFLevel2()->HitPaddle(ipl,ipad)) { // HitPaddle is true if and only if there are both TDC signal for that pad AND they are both 'good'
29			ihitpad=event->GetToFLevel2()->GetPaddleid(ipl,ipad);
30			nhitpad++;
31			}
32			}
33
34			if (nhitpad == 1) { // there is one and only one hit pad in the layer
35			int satwarnpad;
36			event->GetToFLevel2()->GetdEdxPaddle(iTrackTof,ihitpad,100,dedxpad,satwarnpad);
37			// NOTE: if dedxpad=0 then there is no dE/dL signal (i.e. no ADC signals to evaluate dE/dL): this means that ADC signals are under threshold and therefore event is good
38			if ( ! (dedxpad < _maxDedxS1) ) {
39			// cout << "return on check " << ipl << " hitpad " << ihitpad << " track id " << *_trackTof << " dedxpad " << dedxpad << endl;
40			return ipl;
41			}
42			}
43
44			}
45
46			return CUTOK;
47
48			}