MiscCuts/OBTPktNumCut/OBTPktNumCut.cpp

/*
 * OBTPktNumCut.cpp
 *
 *  Created on: 18/mar/2010
 *      Author: Nicola Mori
 */

#include "OBTPktNumCut.h"

OBTPktNumCut::OBTPktNumCut(const char *cutName, vector<pair<unsigned int, unsigned int> > &list) :
  PamCut(cutName), _list(list) {
}

OBTPktNumCut::OBTPktNumCut(const char *cutName, TString listFileName) :
  PamCut(cutName) {
  ifstream listFile(listFileName);

  pair<unsigned int, unsigned int> buffer;

  /* File layout:
   *          OBT      PktNum */
  while (listFile >> buffer.first >> buffer.second) {
    _list.push_back(buffer);
  }

  listFile.close();
}

int OBTPktNumCut::Check(PamLevel2 *event) {

  if (_list.size() == 0)
    return 0; // No more couples left in the list

  unsigned int OBT = event->GetOrbitalInfo()->OBT;
  unsigned int pktNum = event->GetOrbitalInfo()->pkt_num;

  vector<pair<unsigned int, unsigned int> >::iterator currPair = _list.begin();
  while (currPair->first != OBT && currPair->second != pktNum){
    currPair++;
    if (currPair > _list.end())
      return 0; // Event not found in the list
  }

  _list.erase(currPair); // Don't search anymore for the current pair (already found)

  return CUTOK;
}
1	/*
2	* OBTPktNumCut.cpp
3	*
4	* Created on: 18/mar/2010
5	* Author: Nicola Mori
6	*/
7
8	#include "OBTPktNumCut.h"
9
10	OBTPktNumCut::OBTPktNumCut(const char *cutName, vector<pair<unsigned int, unsigned int> > &list) :
11	PamCut(cutName), _list(list) {
12	}
13
14	OBTPktNumCut::OBTPktNumCut(const char *cutName, TString listFileName) :
15	PamCut(cutName) {
16	ifstream listFile(listFileName);
17
18	pair<unsigned int, unsigned int> buffer;
19
20	/* File layout:
21	* OBT PktNum */
22	while (listFile >> buffer.first >> buffer.second) {
23	_list.push_back(buffer);
24	}
25
26	listFile.close();
27	}
28
29	int OBTPktNumCut::Check(PamLevel2 *event) {
30
31	if (_list.size() == 0)
32	return 0; // No more couples left in the list
33
34	unsigned int OBT = event->GetOrbitalInfo()->OBT;
35	unsigned int pktNum = event->GetOrbitalInfo()->pkt_num;
36
37	vector<pair<unsigned int, unsigned int> >::iterator currPair = _list.begin();
38	while (currPair->first != OBT && currPair->second != pktNum){
39	currPair++;
40	if (currPair > _list.end())
41	return 0; // Event not found in the list
42	}
43
44	_list.erase(currPair); // Don't search anymore for the current pair (already found)
45
46	return CUTOK;
47	}