TofCuts/TofTopS1Cut/TofTopS1Cut.cpp

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

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

#include "TofTopS1Cut.h"

int TofTopS1Cut::Check(PamLevel2 *event) {

  // Array initialization
  for (int ipl = 0; ipl < 2; ipl++) {
    _nHitPmtPl[ipl] = 0;
  }

  Int_t pmtId = -1;
  Int_t layerId = 10; // valid values are 0 to 5
  
  // Compute the number of hit PMT on each plane
  for (int ih = 0; ih < event->GetToFLevel2()->npmt(); ih++) { // Int_t mandatory
    pmtId = event->GetToFLevel2()->GetToFPMT(ih)->pmt_id;
    layerId = event->GetToFLevel2()->GetPlaneIndex(pmtId);
    if (0 <= layerId && layerId < 2)
      _nHitPmtPl[layerId]++;
  }

  // CUT on hit PMTs
  for (int ipl = 0; ipl < 2; ipl++) { // check done only for S1 layers
    if (!(_nHitPmtPl[ipl] < 3))
      return 0;
  }

  return CUTOK;
}
1	/*
2	* TofTopS1Cut.cpp
3	*
4	* Created on: 13-mar-2009
5	* Author: Sergio Ricciarini, Nicola Mori
6	*/
7
8	/! @file TofTopS1Cut.cpp The TofTopS1Cut class implementation file /
9
10	#include "TofTopS1Cut.h"
11
12	int TofTopS1Cut::Check(PamLevel2 *event) {
13
14	// Array initialization
15	for (int ipl = 0; ipl < 2; ipl++) {
16	_nHitPmtPl[ipl] = 0;
17	}
18
19	Int_t pmtId = -1;
20	Int_t layerId = 10; // valid values are 0 to 5
21
22	// Compute the number of hit PMT on each plane
23	for (int ih = 0; ih < event->GetToFLevel2()->npmt(); ih++) { // Int_t mandatory
24	pmtId = event->GetToFLevel2()->GetToFPMT(ih)->pmt_id;
25	layerId = event->GetToFLevel2()->GetPlaneIndex(pmtId);
26	if (0 <= layerId && layerId < 2)
27	_nHitPmtPl[layerId]++;
28	}
29
30	// CUT on hit PMTs
31	for (int ipl = 0; ipl < 2; ipl++) { // check done only for S1 layers
32	if (!(_nHitPmtPl[ipl] < 3))
33	return 0;
34	}
35
36	return CUTOK;
37	}